BZOJ3680 吊打xxx

爬山。

爬山算法即是模拟爬山的过程，随机选择一个位置爬山，每次朝着更高
的方向移动，直到到达山顶，即每次都在临近的空间中选择最优解作为
当前解，直到局部最优解。这样算法会陷入局部最优解，能否得到全局
最优解取决于初始点的位置。初始点若选择在全局最优解附近，则就可
能得到全局最优解。

http://www.matrix67.com/blog/archives/422

By：大奕哥

#include<bits/stdc++.h>
using namespace std;
int n;
double ansx,ansy;
struct node{
    double x,y;
    int w;
}p[10005];
double sqr(double x){return x*x;}
double dis(double x,double y,node b){
    return sqrt(sqr(x-b.x)+sqr(y-b.y));
}
void hillclimb()
{
    double t=1000,x,y;
    for(int i=1;i<=n;++i)
    ansx+=p[i].x*p[i].w,ansy+=p[i].y*p[i].w;
    ansx/=n;ansy/=n;
    while(t>1e-9)
    {
        x=y=0;
        for(int i=1;i<=n;++i)
        {
            x+=(p[i].x-ansx)*p[i].w/dis(ansx,ansy,p[i]);
            y+=(p[i].y-ansy)*p[i].w/dis(ansx,ansy,p[i]);
        }
        ansx+=x*t;ansy+=y*t;
        if(t>0.5)t*=0.5;
        else t*=0.97;
    }
    printf("%.3lf %.3lf",ansx,ansy);
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;++i)
        scanf("%lf%lf%d",&p[i].x,&p[i].y,&p[i].w);
    hillclimb();
    return 0;
}