Poj1379(区域内距最近点最远的位置)模拟退火算法实现

题目链接:http://poj.org/problem?id=1379

 

 

题意:给出指定的区域,及平面内的点集,求出一个该区域内到点集中所有点的最短距离最大的一个位置,点.

 

解题:三角剖分问题,可是对于没接触过三角剖分的我,当看到792行的网上找的代码时,当时那个感觉啊```,还好从Disscuss 上看到了可以用随          机化算法模拟退火来解决,下载了个PPT,想通了其中的环节~然后,决定今天把它敲出来,AC之,可敲出来之后以G++提交,RE了,然后悲剧就开        始了,一直RE,期间改成C++提交,CE,痛苦ing~中午午休没睡着,上课前又从床上跳下来,把第一遍提交的代码拉出来,改用C++提交,在以CE        两次为代价下,把那些语法问题改过了,结果,,,AC了~真搞不懂G++,C++,呃,还有VC,看来得换个编译器么?

 

不报怨啦,上代码~

dooder_daodao

1379

Accepted

196K

79MS

C++

1597B

2011-04-20 13:56:45

 

#include<stdio.h> #include<string.h> #include<stdlib.h> #include<time.h> #include<math.h> #define M 1008 //first time: N = 8 ,L = 108 :172ms; //second : N = 4 ,L = 60 :79ms; //初始点集个数和每次迭代中点随机移动次数 #define L 60 #define N 4 const double eps = 1e-3; const double inf = 10e8; const double Pi = acos((double)-1.0);//CE处,用double 显式转化 typedef struct{ double x,y; void get(){ scanf("%lf%lf",&x,&y); } void Randir(double delta,int l,int c){ double d=Pi*(rand())/16384; x+=delta*sin(d); y+=delta*cos(d); } void Randpo(int l,int c) { x=rand()%l+1; y=rand()%c+1; } }Point; double dis[N]; // Point randp[N],p[M]; int n; double delta;//初始温度,用来控制降温过程 /****************/ double Dis(Point p,Point q) { return (p.x-q.x)*(p.x-q.x)+(p.y-q.y)*(p.y-q.y); } int main() { int t,i,j,k,x,y; double tmp,d; Point cur; scanf("%d",&t); srand(time(NULL)); while(t--){ scanf("%d%d%d",&x,&y,&n); for(i=0;i<n;i++) p[i].get(); for(i=0;i<N;i++){ randp[i].Randpo(x,y); dis[i]=inf; for(j=0;j<n;j++){ tmp=Dis(randp[i],p[j]); if(tmp<dis[i]) dis[i]=tmp; } } //能包含矩形区域的最小圆半径 delta=sqrt(1.0*(x*x+y*y))*0.5;//CE处,sqrt 用double 作参数 while(delta>=eps){//冷却至需求点跳出 for(i=0;i<N;i++){ for(j=0;j<L;j++){ cur=randp[i];//对解集中点进行随机移动 cur.Randir(delta,x,y); if(cur.x<0||cur.y<0||cur.x>x||cur.y>y) continue; tmp=inf; for(k=0;k<n;k++){ d=Dis(cur,p[k]); if(d<tmp) tmp=d; } if(tmp>dis[i]){//如果得到更优解 dis[i]=tmp; randp[i]=cur; } } } delta*=0.8;//冷却 } for(i=k=0;i<N;i++){ if(dis[i]>dis[k]) k=i; } printf("The safest point is (%.1f, %.1f)./n",randp[k].x,randp[k].y); } return 0; }