bzoj3680: 吊打XXX

Posted chenyushuo

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了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
 
题解:
puts("nan nan")得AC
求的东西其实是广义带权费马点,然后裸上模拟退火。
code:
 1 #include<cstdio>
 2 #include<iostream>
 3 #include<cmath>
 4 #include<cstring>
 5 #include<algorithm>
 6 #define sqr(x) ((x)*(x))
 7 using namespace std;
 8 char ch;
 9 bool ok;
10 void read(int &x){
11     for (ok=0,ch=getchar();!isdigit(ch);ch=getchar()) if (ch==-) ok=1;
12     for (x=0;isdigit(ch);x=x*10+ch-0,ch=getchar());
13     if (ok) x=-x;
14 }
15 void read(double &x){
16     for (ok=0,ch=getchar();!isdigit(ch);ch=getchar()) if (ch==-) ok=1;
17     for (x=0;isdigit(ch);x=x*10+ch-0,ch=getchar());
18     if (ok) x=-x;
19 }
20 const int maxn=10005;
21 const double pi=3.14159265358979;
22 int n;
23 double maxv=-100000,tx,ty,g[maxn];
24 struct Point{
25     double x,y;
26     void init(){read(x),read(y),tx+=x,ty+=y,maxv=max(maxv,x),maxv=max(maxv,y);}
27 }point[maxn],nxt,ans;
28 double getrand(){return rand()%100000/100000.0;}
29 double calc(Point t){
30     double res=0;
31     for (int i=1;i<=n;i++){
32         double dis=sqrt(sqr(point[i].x-t.x)+sqr(point[i].y-t.y));
33         res+=dis*g[i];
34     }
35     return res;
36 }
37 void sa(){
38     for (double T=maxv,t1=calc(ans),t2;T>0.0000001;T*=0.99){
39         double deg=getrand()*pi*2,len=pow(T,2.0/3);
40         nxt.x=ans.x+cos(deg)*len,nxt.y=ans.y+sin(deg)*len;
41         t2=calc(nxt);
42         if (exp((t1-t2)/T)>getrand()) ans=nxt,t1=t2;
43     }
44     printf("%.3lf %.3lf\n",ans.x,ans.y);
45 }
46 int main(){
47     srand(19990617);
48     read(n);
49     for (int i=1;i<=n;i++) point[i].init(),read(g[i]);
50     ans=(Point){tx/n,ty/n},sa();
51     return 0;
52 }

 

以上是关于bzoj3680: 吊打XXX的主要内容,如果未能解决你的问题,请参考以下文章

BZOJ3680吊打XXX(模拟退火)

bzoj3680: 吊打XXX

bzoj3680 吊打XXX

BZOJ3680吊打XXX 模拟退火

BZOJ3680吊打xxx [模拟退火]

BZOJ 3680: 吊打XXX模拟退火算法裸题学习,爬山算法学习