种花小游戏 随机化搜索
植物大战僵尸这款游戏中，还有个特别有意思的赚钱方式——种花（能长金币的花）。
种出来的金币需要玩家点击才能得到，或者，玩家可以购买一只蜗牛来帮助捡金币。然而，蜗牛爬得慢是众所周知的。所以，场上有若干金币时，蜗牛总是喜欢以最少的行程来捡走所有的金币。
现在告诉你场上n个金币所在位置的坐标，以及蜗牛所在位置，让你求出蜗牛捡走所有金币的最小行程。
比较神奇，数据量是16，dfs可以过；随机化搜索，很稳定（所有的点都相同的时间）；
1 #include<iostream> 2 #include<cstdio> 3 #include<ctime> 4 #include<cstdlib> 5 #include<cmath> 6 using namespace std; 7 #define up(i,p,n) for(int i=p;i<=n;i++) 8 #define down(i,n,p) for(int i=n;i>=p;i--) 9 namespace OI{ 10 const int maxn=20,inf=100000000; 11 int n,a[maxn][2],q[maxn]; 12 double d[maxn][maxn]; 13 void swa(int i,int j){ 14 for(;i<=j;i++,j--){ 15 swap(q[i],q[j]); 16 } 17 } 18 double slove(){ 19 up(i,1,100){ 20 int x=rand()%n+1; 21 int y=rand()%n+1; 22 if(x==n)x=n-1;if(y==n)y=n-1; 23 swap(q[x],q[y]); 24 } 25 up(i,1,n)up(j,i+1,n-1) 26 if(d[q[i]][q[i+1]]+d[q[j-1]][q[j]]>d[q[i]][q[j-1]]+d[q[i+1]][q[j]])swa(i+1,j-1); 27 double t=0; 28 up(i,1,n-1)t+=d[q[i]][q[i+1]]; 29 return t; 30 } 31 int squ(int x){return x*x;} 32 void init(){ 33 srand(time(0)); 34 scanf("%d",&n);n++; 35 up(i,1,n)scanf("%d%d",&a[i][0],&a[i][1]); 36 up(i,1,n)up(j,1,n) 37 d[i][j]=sqrt(1.0*squ(a[i][0]-a[j][0])+squ(a[j][1]-a[i][1])); 38 //up(i,1,n)up(j,1,n)printf("%d %d %.2lf\n",i,j,d[i][j]); 39 up(i,1,n)q[i]=i; 40 double ans=(double)inf; 41 up(i,1,50000){ 42 double t=slove(); 43 if(t<ans)ans=t; 44 } 45 printf("%.2lf\n",ans); 46 } 47 } 48 int main(){ 49 //freopen("1.in","r",stdin); 50 //freopen("1.out","w",stdout); 51 OI::init(); 52 return 0; 53 }
