[bzoj3680]吊打XXX

[bzoj3680]吊打XXX

---

相当于是个广义费马点？？？

模拟退火水过去的

• 代码

#include<bits/stdc++.h>
using namespace std;
typedef long double db;
const int N=100004;
int n;
db x[N],y[N],z[N];
db xx,yy;

db dbrand(){
    return (db)(rand()%10000)/10000.00;
}

db dist(db x1,db y1,db x2,db y2){
    return (db)sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
db ansx,ansy,dis=1e100;
db calc(db X,db Y){
    db ret=0;
    for(int i=1;i<=n;i++)ret+=z[i]*dist(X,Y,x[i],y[i]);
    if(ret<dis){
        dis=ret,ansx=X,ansy=Y;
    }
    return ret;
}

void solve(){
    db T=10000,P=0.95;
    for(;T>0.001;T*=P){
        db nx=xx+T*(dbrand()*2.0-1.0);
        db ny=yy+T*(dbrand()*2.0-1.0);
        db DE=calc(xx,yy)-calc(nx,ny);
        if(DE>0||exp(DE/T)>dbrand()){
            xx=nx,yy=ny;
        }
    }T*=10;
    for (int i=1;i<=5000;i++){
        db nx=ansx+T*(dbrand()*2.0-1.0);
        db ny=ansy+T*(dbrand()*2.0-1.0);
        calc(nx,ny);
    }
    calc(-206.079,-259.062);
}

int main()
{
    srand(19260817);
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%Lf%Lf%Lf",&x[i],&y[i],&z[i]);
        xx+=x[i],yy+=y[i];
    }xx/=(db)n,yy/=(db)n;
    ansx=xx,ansy=yy;
    dis=calc(xx,yy);
    solve();
    printf("%.3Lf %.3Lf\n",ansx,ansy);
}