【BZOJ3680】吊打XXX（模拟退火）

题面

BZOJ

题解

模拟退火。。。

就是模拟退火

然后这题有毒

各种调参数之后终于\(AC\)了。。

这种题就是玄学呀。。。

温度要调大

最后跑完还要向四周多\(rand\)几次

保证能够找到最优解。。。

#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>1e-9||exp(dlt/T)>Rand())now=nw;

        T*=0.97;

    }

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

    {

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

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

        Calc(now);

    }

}

int main()

{

    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(1000);

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

    return 0;

}