【bzoj3680】【吊打XXX】【模拟退火】

Description

gty又虐了一场比赛，被虐的蒟蒻们决定吊打gty。gty见大势不好机智的分出了n个分身，但还是被人多势众的蒟蒻抓住了。蒟蒻们将
n个gty吊在n根绳子上，每根绳子穿过天台的一个洞。这n根绳子有一个公共的绳结x。吊好gty后蒟蒻们发现由于每个gty重力不同，绳
结x在移动。蒟蒻wangxz脑洞大开的决定计算出x最后停留处的坐标，由于他太弱了决定向你求助。
不计摩擦，不计能量损失，由于gty足够矮所以不会掉到地上。

Input

输入第一行为一个正整数n(1<=n<=10000)，表示gty的数目。
接下来n行,每行三个整数xi，yi，wi，表示第i个gty的横坐标，纵坐标和重力。
对于20%的数据，gty排列成一条直线。
对于50%的数据，1<=n<=1000。
对于100%的数据，1<=n<=10000,-100000<=xi,yi<=100000

Output

输出1行两个浮点数（保留到小数点后3位），表示最终x的横、纵坐标。

Sample Input

3
0 0 1
0 2 1
1 1 1

Sample Output

0.577 1.000

题解：模拟退火裸题。首先选中点，然后随机的往四周跳。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<ctime>
#include<cstdlib> 
using namespace std;
struct use{double x,y,g;}p[10100];
use a;int n;
double minn=999999999999999ll;
double dis(use x,use y){return sqrt((x.x-y.x)*(x.x-y.x)+(x.y-y.y)*(x.y-y.y));}
double cal(use x)
{
    double temp=0;
    for(int i=1;i<=n;i++) temp+=p[i].g*dis(x,p[i]);
    if(temp<minn) {a=x;minn=temp;}
    return temp;
}
double ra(){return rand()%1000/1000.0;}
void work(double t)
{
   use now=a; 
   while (t>0.01)
   {
      use temp;
      temp.x=now.x+t*(ra()*2-1);
      temp.y=now.y+t*(ra()*2-1);
      double d=cal(now)-cal(temp);
      if (d>0||exp(d/t)>ra()) now=temp;
      t*=0.984;
   }
    for(int i=1;i<=1000;i++)
    {
        use temp;
        temp.x=a.x+t*(ra()*2-1);
        temp.y=a.y+t*(ra()*2-1);
        cal(temp);
    }   
}
int main()
{
  scanf("%d",&n);
  srand(23333);
  for (int i=1;i<=n;i++)
   {
     scanf("%lf%lf%lf",&p[i].x,&p[i].y,&p[i].g);
     a.x+=p[i].x;a.y+=p[i].y;
   }
   a.x/=n;a.y/=n;
   work(500000);
   printf("%.3lf %.3lf",a.x,a.y);
}