BZOJ 1011 [HNOI2008]遥远的行星

题意：

其实没啥题意，算星球受的重力，允许有5%的误差。

解题思路：

其实也没什么解题思路，一开始以为是递推，发现好像暴力推太慢了，后来看题解发现2000个以内可以暴力算，超过两千个就a[i]*a[j]/j-i中j取平均值。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<cmath>
#include<cctype>
#include<ctime>
using namespace std;
#define INF (2139062143)
#define F (100000007)
#define MAXN (100000+10)
#define MAXAJ (35000)
#define MAXMi (10000000)
#define eps (1e-6) 

int n,i,j;
double a[MAXN],s[MAXN],f[MAXN]={0.0};
double A;

int main(){
	
	scanf("%d%lf",&n,&A);
	a[0]=s[0]=0;
	for(i=1;i<=n;i++) {scanf("%lf",&a[i]); s[i]=a[i]+s[i-1];}
	
	int Q=min(2000,n);
	
	for (i=1;i<=Q;i++) {
		double maxj=A*(double)i+eps;
		for(j=1;j<=maxj;j++) f[i]+=a[j]/(double)(i-j);
		f[i]*=a[i];
	}
	for (i=Q+1;i<=n;i++){
		double maxj=A*double(i)+eps;
		int Ai=floor(maxj);
		f[i]=a[i]*s[Ai]/double(i-(1+Ai)/2);
	}
	for(i=1;i<=n;i++) printf("%.6lf\n",f[i]);
	return 0;
}