【Luogu1337】平衡点（模拟退火）

题面

洛谷

题解

和BZOJ3680吊打XXX是一样的。。

但是数据很强呀。。

疯狂调参

各种WA。。。

很无奈呀。。。。

#include<iostream>

#include<cstdio>

#include<cstdlib>

#include<cstring>

#include<cmath>

#include<algorithm>

#include<set>

#include<map>

#include<vector>

#include<queue>

#include<ctime>

using namespace std;

#define ll long long

#define RG register

#define MAX 20000

#define sqr(x) ((x)*(x))

inline int read()

{

    RG int x=0,t=1;RG char ch=getchar();

    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();

    if(ch=='-')t=-1,ch=getchar();

    while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();

    return x*t;

}

struct Node{double x,y,w;}p[MAX],ans;

double mm=1e100;

int n;

double Dis(Node a,Node b){return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y));}

double Calc(Node P)

{

    double ret=0;

    for(int i=1;i<=n;++i)ret+=p[i].w*Dis(P,p[i]);

    if(ret<mm)mm=ret,ans=P;

    return ret;

}

double Rand(){return rand()%100000/100000.00;}

void SA(double T)

{

    Node now=ans;

    while(T>1e-3)

    {

        Node nw;

        nw.x=now.x+T*(2*Rand()-1);

        nw.y=now.y+T*(2*Rand()-1);

        double dlt=Calc(now)-Calc(nw);

        if(dlt>0||exp(dlt/T)>Rand())now=nw;

        T*=0.99;

    }

	for(int i=1;i<=5000;++i)

	{

        now.x=ans.x+T*(2*Rand()-1);

        now.y=ans.y+T*(2*Rand()-1);

        Calc(now);

    }

}

int main()

{

	srand(233);

    n=read();

    for(int i=1;i<=n;++i)ans.x+=p[i].x=read(),ans.y+=p[i].y=read(),p[i].w=read();

    ans.x/=n;ans.y/=n;Calc(ans);

	SA(100000);

    printf("%.3lf %.3lf\n",ans.x,ans.y);

    return 0;

}