[bzoj1011][HNOI2008]遥远的行星【乱搞】

【题目描述】

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

Sample Output

0.000000
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;
}