[POJ2420]A Star not a Tree?

来源：

Waterloo Local 2002.01.26

题目大意：

找出$n$个点的费马点。

思路：

模拟退火。
首先任取其中一个点（或随机一个坐标）作为基准点，每次向四周找距离为$t$的点，如果找到的点总距离更小，就把该点作为新的基准点。
每次找完后降低温度$t$，重复上述过程，直到温度$t$低于阈值$threshold$。
另外注意在POJ上double类型输出不能用%lf，只能用%f，一开始因为这个WA了好几次。

 #include<cmath>
 #include<cstdio>
 const double inf=1e9;
 const int N=;
 const double threshold=1e-,delta=0.98,initial_temperature=;
 const double d[][]={{,},{,},{,-},{-,}};
 int n;
 struct Point {
     double x,y;
 };
 Point p[N];
 inline double sqr(const double x) {
     return x*x;
 }
 inline double dist(const Point a,const Point b) {
     return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y));
 }
 inline double getsum(const Point x) {
     double ret=;
     for(int i=;i<n;i++) {
         ret+=dist(x,p[i]);
     }
     return ret;
 }
 int main() {
     scanf("%d\n",&n);
     for(int i=;i<n;i++) {
         scanf("%lf%lf",&p[i].x,&p[i].y);
     }
     Point nowp=p[];
     double ans=getsum(nowp),t=initial_temperature;
     while(t>threshold) {
         bool finished=false;
         while(!finished) {
             finished=true;
             for(int i=;i<;i++) {
                 Point nextp;
                 nextp.x=nowp.x+d[i][]*t;
                 nextp.y=nowp.y+d[i][]*t;
                 double tmp=getsum(nextp);
                 if(tmp<ans) {
                     ans=tmp;
                     nowp=nextp;
                     finished=false;
                 }
             }
         }
         t*=delta;
     }
     printf("%.0f\n",ans);
     return ;
 }