【BZOJ】【P3680】【吊打XXX】【题解】【模拟退火】

传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=3680

模拟退火果然是好东西……不过常数调了好久……

Code：

#include<cmath>
#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
const double eps=1e-9;
const int maxn=1e4+10;
int n;
double a[maxn][2];
double w[maxn];
double dis(double x,double y,int i){
	return sqrt((x-a[i][0])*(x-a[i][0])+(y-a[i][1])*(y-a[i][1]));
}
double dist(double x,double y){
	double ans=0;
	for(int i=1;i<=n;i++)
	ans+=dis(x,y,i)*w[i];
	return ans;
}
double rnd(){
	return (double)rand()/RAND_MAX;
}
int main(){
	srand(10086);
	while(scanf("%d",&n)==1){
		double x=0,y=0;
		for(int i=1;i<=n;i++)
		scanf("%lf%lf%lf",&a[i][0],&a[i][1],&w[i]),
		x=max(x,a[i][0]),y=max(y,a[i][1]);
		x++;y++;
		double T=1e3,ans=dist(x,y);
		while(T>eps){
			double dx=0,dy=0;
			for(int i=1;i<=n;i++){
				dx+=(a[i][0]-x)/dis(x,y,i)*w[i];
				dy+=(a[i][1]-y)/dis(x,y,i)*w[i];
			}
			double res=dist(x+dx*T,y+dy*T);
			if(res<ans){
				ans=res;x+=dx*T;y+=dy*T;
			}else{
				double test=pow(2.718281828,(ans-res)/T);
				if(rnd()<test){
					ans=res;x+=dx*T;y+=dy*T;
				}
			}
			T*=0.98;
		}printf("%.3lf %.3lf\n",x,y);
	}
	return 0;
}