poj 2420 A Star not a Tree?

萌萌哒传送门

基本上代码和黄学长的一毛一样

模拟退火

主要是生成新状态和评价函数搞不清怎么弄。。

自己YY的一个每次跑出来答案不一样T_T

有人说是以要移动的点为中心，以2t为长的正方形中的随机一个移动位置来求最优值。。

看了代码发现他那个讲道理也不算模拟退火，因为没有概率接受非更优解。。

然而为甚么黄学长在生成新状态没有用随机化。。。

最后得出结论：大神思路太diao吾等蒟蒻只配%%%

代码（有谁看懂了给解释下好么。。）：

#include<set>#include<map>#include<cmath>#include<ctime>#include<queue>#include<stack>#include<cstdio>#include<vector>#include<cstdlib>#include<cstring>#include<iostream>#include<algorithm>#define g getchar()#define ll long long#define inf 0x3f3f3f3fusing namespace std;inline ll read(){ll x=0,f=1;char ch=g;for(;ch<'0'||ch>'9';ch=g)if(ch=='-')f=-1;for(;ch>='0'&&ch<='9';ch=g)x=x*10+ch-'0';return x*f;}inline void out(ll x){int a[25],wei=0;if(x<0)putchar('-'),x=-x;for(;x;x/=10)a[++wei]=x%10;if(wei==0){puts("0");return;}for(int j=wei;j>=1;--j)putchar('0'+a[j]);putchar('\n');}struct re{int x,y;}p[105];int n;inline double sqr(double x){return x*x;}double dis(double x,double y,re p){return sqrt(sqr(x-p.x)+sqr(y-p.y));}double getsum(double x,double y){double tmp=0;for(int i=1;i<=n;++i)tmp+=dis(x,y,p[i]);return tmp;}int main(){srand(time(0));while(scanf("%d",&n)==1){double xx=0,yy=0;double ans=1e20;double t=100000;for(int i=1;i<=n;++i){p[i].x=read();p[i].y=read();xx+=p[i].x;yy+=p[i].y;}xx/=n;yy/=n;ans=getsum(xx,yy);double tmp,x,y;for(;t>0.02;t*=0.9){x=y=0;for(int i=1;i<=n;++i){                       //生成新状态x+=(p[i].x-xx)/dis(xx,yy,p[i]);y+=(p[i].y-yy)/dis(xx,yy,p[i]);}tmp=getsum(xx+x*t,yy+y*t);if(tmp<ans){ans=tmp;xx+=x*t;yy+=y*t;}else if(log(tmp-ans)/t<(rand()%10000)/10000.0){          //随机化判定是否接受移动ans=tmp;xx+=x*t;yy+=y*t;}}printf("%.0lf",ans);}}

萌萌哒传送门2