奇奇怪怪突然出戏的奇葩题
原题:
直线上N颗行星,X=i处有行星i,行星J受到行星I的作用力,当且仅当i<=AJ.此时J受到作用力的大小为 Fi->j=
Mi*Mj/(j-i) 其中A为很小的常量,故直观上说每颗行星都只受到距离遥远的行星的作用。请计算每颗行星的受力
,只要结果的相对误差不超过5%即可.
1<=N<=10^5.0.01< a < =0.35,0<=Mi<=10^7
引用chty的一句话:“如何提高智商???这是个永恒的话题”
乱搞题,似乎没法求出正确答案……
对于每个行星都计算一个k=a*i
如果k比较小,直接暴力计算每一个行星对行星i作用的力
如果k比较大,因为a很小所以i很大,所以如果把所有的i换成k/2答案应该在误差范围内
这个很玄……
我看别人的代码是k<=67的时候暴力
看看代码搞把,这题太奇怪了……
代码:
1 #include<iostream> 2 #include<cstdio> 3 #include<algorithm> 4 #include<cstring> 5 #include<cmath> 6 using namespace std; 7 int n; 8 double a,m[110000]; 9 double sum[110000]; 10 int main(){ 11 //懵逼.jpg 12 //freopen("ddd.in","r",stdin); 13 //freopen("bzoj_1011_planet.in","r",stdin); 14 //freopen("bzoj_1011_planet.out","w",stdout); 15 cin>>n>>a; 16 sum[0]=0; 17 for(int i=1;i<=n;i++){ 18 scanf("%lf",&m[i]); 19 sum[i]=sum[i-1]+m[i]; 20 } 21 for(int i=1;i<=n;i++){ 22 int k=int(a*double(i)+1e-8); 23 double ans=0; 24 if(k<=67) 25 for(int j=1;j<=k;j++) 26 ans+=(m[i]*m[j])/double(i-j); 27 else 28 ans=(sum[k]*m[i])/double(i-((1+k)>>1)); 29 printf("%.6lf ",ans); 30 } 31 return 0; 32 }