吊打XXX

3680: 吊打XXX

Time Limit: 10 Sec  Memory Limit: 128 MBSec  Special Judge
Submit: 2670  Solved: 984
[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

pre.cjk { font-family: "Droid Sans Fallback", monospace }
p { margin-bottom: 0.25cm; line-height: 120% }
a:link { }

模拟退火算法:
①设定一个初始温度t(在不会TLE的情况下尽量大),设定一个初始答案。
②用随机函数获得新解,计算答案的增量da。
③若da>0,即新解更优,则用新解更新最优解和当前解,否则获得一个随机概率
(-1,1),若这个随机概率小于exp(da/t),则用这个解更新当前解(注意不是更新最优解)。
④t乘以一个设定的常数(0,1),降温。
⑤一直重复做上面的步骤,直到温度已经降为设定的最小温度。
⑥以当前温度在最优解的附近再探测几次,这样可能发现更优的解。
 1 #include<set>
 2 #include<map>
 3 #include<queue>
 4 #include<stack>
 5 #include<ctime>
 6 #include<cmath>
 7 #include<string>
 8 #include<vector>
 9 #include<cstdio>
10 #include<cstdlib>
11 #include<cstring>
12 #include<iostream>
13 #include<algorithm>
14 #define maxn 10010
15 using namespace std;
16 struct data{
17   double x,y,w;
18 }gty[maxn];
19 int n;
20 double ansx=0.0,ansy=0.0,MIN=199999999999999999;
21 double cale(double x,double y){
22   double ret=0;
23   for(int i=1;i<=n;i++)
24     ret+=gty[i].w*sqrt((x-gty[i].x)*(x-gty[i].x)+(y-gty[i].y)*(y-gty[i].y));
25   if(ret<MIN){ansx=x,ansy=y,MIN=ret;}
26   return ret;
27 }
28 double RAND(){return rand()%1000/1000.0;}
29 void SA(double t){
30   double nowx=ansx,nowy=ansy;
31   while(t>0.001){
32     double tmpx,tmpy;
33     tmpx=nowx+t*(RAND()*2-1);
34     tmpy=nowy+t*(RAND()*2-1);
35     double da=cale(nowx,nowy)-cale(tmpx,tmpy);
36     if(da>0 || exp(da/t)>RAND()) nowx=tmpx,nowy=tmpy;
37     t*=0.993;
38   }
39   for(int i=1;i<=1000;i++){
40     double tmpx,tmpy;
41     tmpx=ansx+t*(RAND()*2-1);
42     tmpy=ansy+t*(RAND()*2-1);
43     cale(tmpx,tmpy);
44   }
45   printf("%.3lf %.3lf",ansx,ansy);
46 }
47 int main()
48 {
49   freopen("!.in","r",stdin);
50   freopen("!.out","w",stdout);
51   srand(23333);
52   scanf("%d",&n);
53   for(int i=1;i<=n;i++)
54     scanf("%lf%lf%lf",&gty[i].x,&gty[i].y,&gty[i].w),ansx+=gty[i].x,ansy+=gty[i].y;
55   ansx/=n,ansy/=n;
56   SA(5000000);
57   return 0;
58 }
时间: 2024-10-26 01:52:08

吊打XXX的相关文章

【BZOJ3680】吊打XXX 模拟退火

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

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

3680: 吊打XXX Time Limit: 10 Sec  Memory Limit: 128 MBSec  Special JudgeSubmit: 3192  Solved: 1198[Submit][Status][Discuss] Description gty又虐了一场比赛,被虐的蒟蒻们决定吊打gty.gty见大势不好机智的分出了n个分身,但还是被人多势众的蒟蒻抓住了.蒟蒻们将 n个gty吊在n根绳子上,每根绳子穿过天台的一个洞.这n根绳子有一个公共的绳结x.吊好gty后蒟蒻们发现

【BZOJ3680】吊打XXX(模拟退火)

[BZOJ3680]吊打XXX(模拟退火) 题面 BZOJ 题解 模拟退火... 就是模拟退火 然后这题有毒 各种调参数之后终于\(AC\)了.. 这种题就是玄学呀... 温度要调大 最后跑完还要向四周多\(rand\)几次 保证能够找到最优解... #include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<cmath> #include<

P1337 [JSOI2004]平衡点 吊打XXX - 模拟退火

P1337 [JSOI2004]平衡点 吊打XXX 模拟退火 初始温度\(T_0\) 终止温度\(T_k\) 温度变化率\(d\) \(T_k\)略大于0,\(d\)略小于1 当前状态\(x,y\) 当前解\(E\) 当前最优解\(minE\) 当前温度\(T\) 新状态\(nx,ny\) 新解\(nE\) 新解与当前解差值\(\Delta E = nE-E\) \(if\)新解比当前解更优\(nE<E\) 当前状态\(x,y\)移动到 \(nx,ny\) 当前解\(E\)移动到\(nE\) $

BZOJ 3680 吊打XXX

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

bzoj3680 吊打XXX

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

BZOJ 3680 吊打XXX 模拟退火

首先这题应该改名叫吊打出题人 题目大意:给定n个质点,求重心 这n个质点的重心满足Σ(重心到点i的距离)*g[i]最小 模拟退火的裸题 尼玛交了两篇 死活过不去 各种改参数 最后发现是我的INF不够大 尼玛! 这题INF开0x3f妥妥过不去...起码要max_of _long_long附近才可以 最后写了10188MS,BZOJ倒数第一--这也是种艺术啊0.0 #include<cmath> #include<cstdio> #include<cstring> #inc

【bzoj3680】吊打XXX 随机化

题目描述 gty又虐了一场比赛,被虐的蒟蒻们决定吊打gty.gty见大势不好机智的分出了n个分身,但还是被人多势众的蒟蒻抓住了.蒟蒻们将n个gty吊在n根绳子上,每根绳子穿过天台的一个洞.这n根绳子有一个公共的绳结x.吊好gty后蒟蒻们发现由于每个gty重力不同,绳结x在移动.蒟蒻wangxz脑洞大开的决定计算出x最后停留处的坐标,由于他太弱了决定向你求助.不计摩擦,不计能量损失,由于gty足够矮所以不会掉到地上. 输入 输入第一行为一个正整数n(1<=n<=10000),表示gty的数目.接

B3680 吊打xxx 物理???

看到一道很有意思的题,这个题简直有毒,是一道物理题...好像得用模拟退火...但显然我太弱了不会模拟退火,只能用正交分解暴力... 每次沿着力的方向走一定的距离,假如转头了,则走的步长就减小一点. 不过这里有一个坑,就是假如每次二分,就会错...0.7或0.9就可以.有可能走过了回不来了吧. 题面: Description gty又虐了一场比赛,被虐的蒟蒻们决定吊打gty.gty见大势不好机智的分出了n个分身,但还是被人多势众的蒟蒻抓住了.蒟蒻们将 n个gty吊在n根绳子上,每根绳子穿过天台的一