【BZOJ3680】吊打XXX（模拟退火）
题面
题解
模拟退火。。。
就是模拟退火
然后这题有毒
各种调参数之后终于\(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;
}