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: 1011
    User: HansBug
    Language: Pascal
    Result: Accepted
    Time:2296 ms
    Memory:1984 kb
****************************************************************/
 
const st=50;
var
    i,j,k,l,m,n:longint;  
    p:extended;
    b,c:array[0..100010] of longint;
    a:array[0..100010] of extended;
begin
    readln(n,p);
    for i:=1 to n do read(b[i]);
    readln;
    for i:=1 to n do c[i]:=trunc(i*p);
    for i:=1 to n do
        begin
            if a[i]<>0 then a[i]:=0;
            if c[i]<>0 then
                begin
                    if i<=st then
                        begin
                            for j:=1 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]/2))/(i-(c[i-st]/2)));
                            for j:=c[i-st]+1 to c[i] do a[i]:=a[i]+b[j]/(i-j);
                            a[i]:=a[i]*b[i];
                        end;
                end;
        end;
    for i:=1 to n do writeln(a[i]:0:6);
    readln;
end.