基本上代码和黄学长的一毛一样
模拟退火
主要是生成新状态和评价函数搞不清怎么弄。。
自己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 0x3f3f3f3f
using 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);
}
}