1011: [HNOI2008]遥远的行星

1011: [HNOI2008]遥远的行星

Time Limit: 10 Sec  Memory Limit: 162 MBSec  Special Judge
Submit: 2241  Solved: 793
[Submit][Status]

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

 题解：呵呵呵呵呵。。。感觉自己真心捉鸡——

从来就只知道只有物理竞赛化学竞赛里面会用近似处理，没想到这里也近似。

大致思路就是先进行分段处理，然后小的部分暴力，大的部分近似处理，有效的利用了题目中所给的误差不超过5%呵呵哒^_^

接下来开始引用（感觉省选现场能A掉此题的都是神OTL，这个近似真心略诡异啊）

设f[i]为第i个行星受到的引力，g[i]表示第i个行星所受引力的最后来源行星

然后小范围（我是取 α * i < 9）暴力搞，大范围就乱搞，忽略除数（i - j）的差距

 /**************************************************************
     Problem:
     User: HansBug
     Language: Pascal
     Result: Accepted
     Time: ms
     Memory: kb
 ****************************************************************/

 const st=;
 var
     i,j,k,l,m,n:longint;
     p:extended;
     b,c:array[..] of longint;
     a:array[..] of extended;
 begin
     readln(n,p);
     for i:= to n do read(b[i]);
     readln;
     for i:= to n do c[i]:=trunc(i*p);
     for i:= to n do
         begin
             if a[i]<> then a[i]:=;
             if c[i]<> then
                 begin
                     if i<=st then
                         begin
                             for j:= to c[i] do a[i]:=a[i]+b[j]/(i-j);
                             a[i]:=a[i]*b[i];
                         end
                     else
                         begin
                             a[i]:=(a[i-st]/b[i-st])*((i-st-(c[i-st]/))/(i-(c[i-st]/)));
                             for j:=c[i-st]+ to c[i] do a[i]:=a[i]+b[j]/(i-j);
                             a[i]:=a[i]*b[i];
                         end;
                 end;
         end;
     for i:= to n do writeln(a[i]::);
     readln;
 end.