传送门
3680: 吊打XXX
Time Limit: 10 Sec Memory Limit: 128 MBSec Special Judge
Submit: 6513 Solved: 2221
[Submit][Status][Discuss]
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
HINT
Source
By wangxz
题意:给n个质点的位置(xi,yi)以及质量wi,求所有质点的重心(也就是使 ∑ i = 1 n ( ( x i − x ) ∗ ( x i − x ) + ( y i − y ) ∗ ( y i − y ) ∗ w i ) \sum_{i=1}^n(\sqrt{(xi-x)*(xi-x)+(yi-y)*(yi-y)} * wi) ∑i=1n((xi−x)∗(xi−x)+(yi−y)∗(yi−y)∗wi))最小的(x,y)
题解:模拟退火。
关于调参:srand()里的常数以及降温系数d还有初始温度T0和终止温度Tk都可能影响结果的正确性…而且复杂度比较玄,一般来说多烧几次虽然会提高准确性但是可能会造成TLE…还是随缘调参看脸吧(
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define debug(x) cout<<#x<<" is "<<x<<endl;
const int maxn=1e4+5;
const double inf=1e7;
const double eps=1e-3;
int n;
struct Pot{
double x;
double y;
double ac;
}p[maxn],ans;
double Rand(){
return (double)rand()/RAND_MAX;
}
double cal(struct Pot a){
double ak=0;
for(int i=1;i<=n;i++){
ak+=sqrt((a.x-p[i].x)*(a.x-p[i].x)+(a.y-p[i].y)*(a.y-p[i].y))*p[i].ac;
}
return ak;
}
void SA(){
double t=200;
struct Pot now=ans;
struct Pot newpot;
while(t>eps){
newpot.x=now.x+(2*Rand()-1)*t;
newpot.y=now.y+(2*Rand()-1)*t;
newpot.ac=cal(newpot);
if(newpot.ac<now.ac||exp(-(newpot.ac-now.ac)/t)>Rand())now=newpot;
if(now.ac<ans.ac)ans=now;
t*=0.97;
}
for(int i=1;i<=1000;i++){
newpot.x=ans.x+(2*Rand()-1)*t;
newpot.y=ans.y+(2*Rand()-1)*t;
newpot.ac=cal(newpot);
if(newpot.ac<ans.ac)ans=newpot;
}
}
int main(){
srand(19990716);
scanf("%d",&n);
ans.x=0;
ans.y=0;
for(int i=1;i<=n;i++){
scanf("%lf%lf%lf",&p[i].x,&p[i].y,&p[i].ac);
ans.x+=p[i].x;
ans.y+=p[i].y;
}
ans.x/=n;ans.y/=n;ans.ac=cal(ans);
SA();
printf("%.3f %.3f\n",ans.x,ans.y);
return 0;
}