P1337 [JSOI2004]平衡点 吊打XXX
模拟退火
初始温度\(T_0\) 终止温度\(T_k\) 温度变化率\(d\)
\(T_k\)略大于0，\(d\)略小于1
当前状态\(x,y\) 当前解\(E\) 当前最优解\(minE\) 当前温度\(T\)
新状态\(nx,ny\) 新解\(nE\) 新解与当前解差值\(\Delta E = nE-E\)
\(if\)新解比当前解更优\(nE
当前状态\(x,y\)移动到 \(nx,ny\)
当前解\(E\)移动到\(nE\)
\(else\) 我们有 \(e^{\frac{\Delta E}{kT}}\) 的概率使当前状态移动到\(nx,ny\)(其中\(k\)为随机数)
当前温度\(T\)降至\(T*d\)
最后\(T\)降到等于低于\(T_k\)时，算法结束。
Tips：
1.\(nx,ny\)可通过将\(x,y\)加上一个\([-T,T]内的随机数(T*(Rand()*2-1))\)得到。
其中Rand()函数如下：
inline double Rand(){
return double(rand())/double(RAND_MAX);
}2.判断当前状态是否能移动到新状态
(本题求的是最小值，因此当\(\Delta E = nE-E\)<0时更优)
inline bool Accept(double delta,double T){
return delta<0||Rand()3.这是个看脸的算法，这是个玄学的算法。（参见代码中的八位大质数和八遍退火）
4.别学我拿两个退火对拍。。。
AC CODE：(Warning ： 本题不保证每次都能AC)
#include
#include
#include
#include
#include
#include
using namespace std;
const int N = 1000 + 10;
int n;
struct node{
double x,y,w;
}a[N];
node Ans;
double minE;
inline double dis(double x,double y,double _x,double _y){
return sqrt((x-_x)*(x-_x)+(y-_y)*(y-_y));
}
inline double cal(double x,double y){
double ans=0;
for(int i=1;i<=n;i++){
double _x=a[i].x,_y=a[i].y;
ans+=a[i].w*dis(x,y,_x,_y);
}
if(ansTk){
double nx=x+T*(Rand()*2-1),ny=y+T*(Rand()*2-1);
//新坐标 T*(Rand()*2-1) 生成[-T,T]内的实数
double nE=cal(nx,ny);//新解
if(Accept(nE-E,T)){//转移
x=nx,y=ny;
E=nE;
}
T*=d;//get low
}
for(int i=1;i<=1000;i++){
x=Ans.x+T*(Rand()*2-1),y=Ans.y+T*(Rand()*2-1);
cal(x,y);
}
}
int main(){
// freopen("data.in","r",stdin);
// freopen("sol.out","w",stdout);
srand(19260817);srand(rand());srand(rand());
scanf("%d",&n);
node _;_.x=0,_.y=0;
for(int i=1;i<=n;i++){
double x,y,w;scanf("%lf%lf%lf",&x,&y,&w);
a[i]=(node){x,y,w};
_.x+=x,_.y+=y;
}
_.x/=n,_.y/=n;
double T0=2333,Tk=1e-3,d=1-7e-2;
SA(_,T0,Tk,d);
SA(Ans,T0,Tk,d);
SA(Ans,T0,Tk,d);
SA(Ans,T0,Tk,d);
SA(Ans,T0,Tk,d);
SA(Ans,T0,Tk,d);
SA(Ans,T0,Tk,d);
SA(Ans,T0,Tk,d);
printf("%.3lf %.3lf\n",Ans.x,Ans.y);
return 0;
}