【Luogu P1337】平衡点 / 吊打XXX

题目链接

题目描述

如图：有n个重物，每个重物系在一条足够长的绳子上。每条绳子自上而下穿过桌面上的洞，然后系在一起。图中X处就是公共的绳结。假设绳子是完全弹性的（不会造成能量损失），桌子足够高（因而重物不会垂到地上），且忽略所有的摩擦。
问绳结X最终平衡于何处。
注意：桌面上的洞都比绳结X小得多，所以即使某个重物特别重，绳结X也不可能穿过桌面上的洞掉下来，最多是卡在某个洞口处。

题解

退火板子题。
系统最后处于势能最低状态。(最稳定)
可以直接理解为应让所有点到答案点的距离乘上重量的和最小。

但这题我硬是调了一晚上的参，AC率还不高。
脸太黑没办法TAT 。(是不是不该用时间种子来生成随机数QAQ)

代码如下:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<cstring>
#include<queue>
using namespace std;
typedef double db;
const int N=1e3+10;
struct point {
    db x,y;
}P[N];
db G[N];
const db INF=1e99;
point ans;
db Ans;int n;
const int T=100000;
const db delta=0.99;
const db eps=1e-15;
inline db dist(point a,point b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
inline db Get(point Q)
{
    register db res=0;
    for(register int i=1;i<=n;++i) res+=dist(Q,P[i])*G[i];
    return res;
}
inline db Rand(){
    return ((rand()<<1)-RAND_MAX);
}
inline void SA(db XX,db YY)
{
    ans.x=XX;ans.y=YY;
    Ans=Get(ans);register db t=T;
    point u;u.x=XX;u.y=YY;
    db pans=Ans;
    while(t>eps){
        bool flag=1;
        while(flag){
            flag=0;
            point v;
            v.x=u.x+Rand()*t;v.y=u.y+Rand()*t;
            register db res=Get(v);
            if(res<pans){
                flag=1;u=v;pans=res;
                if(pans<Ans) ans=u,Ans=pans;
            }
            else if(rand()>exp((res-pans)/t)*RAND_MAX){
                flag=1;u=v;pans=res;
            }
        }
        t*=delta;
    }
    for(register int i=1;i<=5000;++i){
        point v;
        v.x=ans.x+Rand()*t;v.y=ans.y+Rand()*t;
        register db res=Get(v);
        if(res<Ans) Ans=res,ans=v;
    }
}
int main()
{
    Ans=INF;scanf("%d",&n);
    srand(time(NULL));
    point K;K.x=0;K.y=0;
    for(register int i=1;i<=n;++i) {
        scanf("%lf %lf %lf",&P[i].x,&P[i].y,&G[i]);
        K.x+=P[i].x;K.y+=P[i].y;
    }
    K.x/=n;K.y/=n;
    SA(K.x,K.y);
    printf("%.3lf %.3lf\n",ans.x,ans.y);
}