BZOJ-1011 遥远的行星 乱搞+估

1011: [HNOI2008]遥远的行星
Time Limit: 10 Sec Memory Limit: 162 MBSec Special Judge
Submit: 3277 Solved: 1208
[Submit][Status][Discuss]

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

搞不明白这题要闹哪样...

前几个计算,后几个估算,大体上用平均算..
最近物理刚学万有引力,数学刚学总体密度曲线??

code:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
int read()
{
    int x=0,f=1; char ch=getchar();
    while (ch<'0' || ch>'9') {if (ch=='-0') f=-1; ch=getchar();}
    while (ch>='0' && ch<='9') {x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
#define maxn 100010 
#define eps 1e-8
#define T 2000
int n,m[maxn];
double F[maxn],A;
int main()
{
    n=read(); scanf("%lf",&A);
    for (int i=1; i<=n; i++) m[i]=read();
    for (int i=1; i<=T; i++)
        {
            int r=eps+floor(i*A);
            for (int j=1; j<=r; j++) 
                F[i]+=(double)m[i]*m[j]/(i-j);
        }
    for (int i=T+1; i<=n; i++)
        {
            int l=eps+floor(A*(i-T)),r=eps+floor(A*i);
            for(int j=l+1; j<=r;j ++)
                F[i]+=(double)m[i]*m[j]/(i-j);
            F[i]+=(double)m[i]*(F[i-T]/m[i-T])*(i-T-(double)l/2)/double(i-(double)l/2);
        }
    for (int i=1; i<=n; i++)
        printf("%.6lf\n",F[i]);
    return 0;
}