「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
0 0 1
0 2 1
1 1 1
Sample Output
0.577 1.000
题解
此题应该和JSOI平衡点是一样的
但是这种题放oj上啥心态。。。
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 |
#include<iostream> #include<cstdio> #include<cmath> using namespace std; int n; double ansx,ansy; struct data{double x,y;int w;}p[10005]; inline double sqr(double x){return x*x;} inline double dis(double x,double y,data p) { return sqrt(sqr(x-p.x)+sqr(y-p.y)); } void hillclimb() { double t=1000,x,y; for(int i=1;i<=n;i++) ansx+=p[i].x*p[i].w,ansy+=p[i].y*p[i].w; ansx/=n;ansy/=n; while(t>0.00000001) { x=y=0; for(int i=1;i<=n;i++) { x+=(p[i].x-ansx)*p[i].w/dis(ansx,ansy,p[i]); y+=(p[i].y-ansy)*p[i].w/dis(ansx,ansy,p[i]); } ansx+=x*t;ansy+=y*t; if(t>0.5)t*=0.5; else t*=0.97; } printf("%.3lf %.3lf",ansx,ansy); } int main() { scanf("%d",&n); for(int i=1;i<=n;i++) scanf("%lf%lf%d",&p[i].x,&p[i].y,&p[i].w); hillclimb(); return 0; } |
计数器 9856665
OI 课件=>
OI 课件=>
问一下这个题22~27行是一种向最优值的逼近么？
抑或是一种伪随机？
没太看懂，请教一下。