【BZOJ1011】遥远的行星（？？？）

题面

BZOJ
洛谷

题解

大概就是分个块，然后每块取平均数算贡献啥的。
BZOJ上过不去？？？

#include<iostream>
#include<cstdio>
using namespace std;
#define MAX 100100
int n,sqr,blk[MAX];
double a,A,m[MAX],s[MAX];
double calc(int j,int a)
{
	double ret=0;
	if(a<=2000){for(int i=sqr+1;i<=a;++i)ret+=m[j]*m[i]/(j-i);return ret;}
	for(int i=1;i<blk[a];++i)ret+=s[i]*m[j]/(j-i*sqr+sqr/2);
	for(int i=(blk[a]-1)*sqr+1;i<=a;++i)ret+=m[j]*m[i]/(j-i);
	return ret;
}
int main()
{
	ios::sync_with_stdio(false);
	cin>>n>>A;sqr=50;blk[0]=1;
	for(int i=1;i<=n;++i)cin>>m[i];
	for(int i=1;i<=n;++i)s[(i-1)/sqr+1]+=m[i],blk[i]=(i-1)/sqr+1;
	for(int i=1;i<=n;++i)printf("%.10lf
",calc(i,i*A));
	return 0;
}