[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=100;
const double threshold=1e-8,delta=0.98,initial_temperature=100;
const double d[4][2]={{0,1},{1,0},{0,-1},{-1,0}};
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=0;
    for(int i=0;i<n;i++) {
        ret+=dist(x,p[i]);
    }
    return ret;
}
int main() {
    scanf("%d
",&n);
    for(int i=0;i<n;i++) {
        scanf("%lf%lf",&p[i].x,&p[i].y);
    }
    Point nowp=p[0];
    double ans=getsum(nowp),t=initial_temperature;
    while(t>threshold) {
        bool finished=false;
        while(!finished) {
            finished=true;
            for(int i=0;i<4;i++) {
                Point nextp;
                nextp.x=nowp.x+d[i][0]*t;
                nextp.y=nowp.y+d[i][1]*t;
                double tmp=getsum(nextp);
                if(tmp<ans) {
                    ans=tmp;
                    nowp=nextp;
                    finished=false;
                }
            }
        }
        t*=delta;
    }
    printf("%.0f
",ans);
    return 0;
}