神奇模拟淬火

#include<bits/stdc++.h>
using namespace std;
int n;
int x[1009];
int y[1009];
int w[1009];

double ansx, ansy, answ;

double engery(double _x, double _y) {
    double r = 0, dx, dy;
    for(int i = 1; i <= n; i++) {
        dx = _x - x[i];
        dy = _y - y[i];
        r += sqrt(dx * dx + dy * dy) * w[i];
    }
    return r;
}

void sa() {
    double t = 1e-15;
    while(t <= 3000) {
        double ex = ansx + (rand() * 2.0 - RAND_MAX)*t;
        double ey = ansy + (rand() * 2.0 - RAND_MAX)*t;
        double ew = engery(ex, ey);
        double de = ew - answ;
        if(de < 0) {
            ansx = ex;
            ansy = ey;
            answ = ew;
        }else if(exp(-de/t) * RAND_MAX > rand()) {
            ansx = ex;
            ansy = ey;
        }
        t*=1.003;
         
    }
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    srand(time(0));
    
    cin >> n;
    for(int i = 1; i <= n; i++) {
        cin >> x[i] >> y[i] >> w[i];
        ansx += x[i];
        ansy += y[i];
    }    
    ansx/=n;
    ansy/=n;
    answ = engery(ansx, ansy);
    while((double)clock()/CLOCKS_PER_SEC<0.95)sa();
    printf("%.3lf %.3lf",ansx,ansy);
    return 0;
}