【题目描述】
Description
直线上N颗行星,X=i处有行星i,行星J受到行星I的作用力,当且仅当i<=AJ.此时J受到作用力的大小为 Fi->j=
Mi*Mj/(j-i) 其中A为很小的常量,故直观上说每颗行星都只受到距离遥远的行星的作用。请计算每颗行星的受力
,只要结果的相对误差不超过5%即可.
Input
第一行两个整数N和A. 1<=N<=10^5.0.01< a < =0.35,接下来N行输入N个行星的质量Mi,保证0<=Mi<=10^7
Output
N行,依次输出各行星的受力情况
Sample Input
5 0.3
3
5
6
2
4
3
5
6
2
4
Sample Output
0.000000
0.000000
0.000000
1.968750
2.976000
0.000000
0.000000
1.968750
2.976000
HINT
精确结果应该为0 0 0 2 3,但样例输出的结果误差不超过5%,也算对
Source
【题解】
近的暴力,原点将距离取平均值算? == 反正A了
但这样一片全是0一头有一个比较大的数的时候会挂,因此要分段算(虽然我没写 QWQ)
/* -------------- user Vanisher problem bzoj-1011 ----------------*/ # include <bits/stdc++.h> # define ll long long # define N 100100 # define T 0 # define eps 1e-8 using namespace std; int read(){ int tmp=0, fh=1; char ch=getchar(); while (ch<'0'||ch>'9'){if (ch=='-') fh=-1; ch=getchar();} while (ch>='0'&&ch<='9'){tmp=tmp*10+ch-'0'; ch=getchar();} return tmp*fh; } int n; double a,w[N],s[N]; int main(){ n=read(); scanf("%lf",&a); for (int i=1; i<=n; i++){ scanf("%lf",&w[i]); s[i]=s[i-1]+w[i]; } for (int i=1; i<=n; i++){ double now=0; if (i<=2000){ for (int j=1; j<=a*i+eps; j++) now=now+w[i]*w[j]/(i-j); } else { int l=1, r=(int)(a*i+eps)-T; if (l<=r){ double mid=(l+r)/2.0; now=now+s[r]*w[i]/(i-mid); } } printf("%.8lf ",now); } return 0; }