[JSOI2004]平衡点 / 吊打XXX

考虑模拟退火。

题目要我们找到一个点，使得整个系统平衡。

这个要求等价于让我们找到一个点，使得系统总能量最小。

我们退火出一个点，然后计算其能量即可。

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5 #include<queue>
 6 #include<cstdlib>
 7 #include<cmath>
 8 #define N 10005
 9 using namespace std;
10 int read()
11 {
12     int x=0,f=1;char ch=getchar();
13     while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();}
14     while(ch>=‘0‘&&ch<=‘9‘){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
15     return x*f;
16 }
17 struct node
18 {
19     int x,y,w;
20 }a[N];
21 int n;
22 double t,ansx,ansy,sx,sy,ans=1e18;
23 const double delta=0.993;
24 double calc(double x,double y)
25 {
26     double res=0;
27     for(int i=1;i<=n;i++)
28     {
29         double dtx=x-a[i].x,dty=y-a[i].y;
30         res+=sqrt((dtx*dtx)+(dty*dty))*a[i].w;
31     }
32     return res;
33 }
34 void SA()
35 {
36     double x=ansx,y=ansy;
37     t=2137;
38     while(t>1e-14)
39     {
40         double X=x+((rand()<<1)-RAND_MAX)*t;
41         double Y=y+((rand()<<1)-RAND_MAX)*t;
42         double now=calc(X,Y);
43         double D=now-ans;
44         if(D<0)
45         {
46             x=X;y=Y;
47             ansx=x;ansy=y;ans=now;
48         }
49         else if(exp(-D/t)*RAND_MAX>rand())x=X,y=Y;
50         t*=delta;
51     }
52 }
53 void solve()
54 {
55     ansx=(double)sx/n;ansy=(double)sy/n;
56     SA();
57     SA();
58     SA();
59     SA();SA();
60 }
61 int main()
62 {
63     srand(44042137);
64     n=read();
65     for(int i=1;i<=n;i++)
66     {
67         a[i].x=read();a[i].y=read();a[i].w=read();
68         sx+=a[i].x;sy+=a[i].y;
69     }
70     solve();
71     printf("%.3lf %.3lf\n",ansx,ansy);
72     return 0;
73 }

原文地址：https://www.cnblogs.com/szmssf/p/11555752.html