【BZOJ】1011 遥远的行星

题意

直线上$n$颗行星，$x=i$处有行星$i$，行星$j$受到行星$i$的作用力，当且仅当$i<=Aj$，此时$j$受到作用力的大小为 $F_{i\rightarrow j}={M_i\times M_j\over j-i}$，其中$A$为很小的常量，故直观上说每颗行星都只受到距离遥远的行星的作用。
请计算每颗行星的受力，只要结果的相对误差不超过5%即可.

$1\leq n\leq {10}^5$
$0.01\leq a\leq 0.35$

分析

设$g_i$表示第$i$颗行星受$1$到$g_i$行星的作用力，则第$i$颗行星受力

$$f_i=\sum_{j=1}^{g_i}{M_i\times M_j\over i-j}=M_i\sum_{j=1}^{g_i}{M_j\over i-j}$$ 。

直接计算$f_i$的时间复杂度为$O(n)$，总的时间复杂度为$O(n^2)$，很明显会超时。

题目中提及“只要结果的相对误差不超过5%即可”，说明我们求出的值不用特别的精确，那么就考虑用近似算法进行估算。

考虑用递推的方法，把所求尽可能凑出已知量。
假设我们已经求出

$$f_i=M_i\sum_{j=1}^{g_i}{M_j\over i-j}$$
则 $$f_{i+T}=M_{i+T}(\sum_{j=1}^{g_i}{M_j\over i-j}\times {i-j\over i+T-j}+\sum_{j=g_i+1}^{g_{i+T}}{M_j\over i+T-j})$$

对于括号内的第二部分，注意$g$值具有的单调性，所以总共$O(n)$可以解决。
关键处理好括号内的第一部分即可。

第一部分可以这样估算：

$$\sum_{j=1}^{g_i}{M_j\over i-j}\times {i-j\over i+T-j}\approx {f_i\over M_i}\times {i-{g_i\over 2}\over i+T-{g_i\over 2}}$$

所以

$$f_{i+T}\approx M_{i+T}\times ({f_i\over M_i}\times {i-{g_i\over 2}\over i+T-{g_i\over 2}}+\sum_{j=g_i+1}^{g_{i+T}}{M_j\over i+T-j})$$

为了使答案更加精确，我们需求更多的准确值。
对于$i\leq T$，我们直接用暴力算法求解。
对于$i>T$，我们用$f_{i-T}$来估算$f_i$。
$T$的取值应尽可能大，而且暴力算法不会超时。
我的$T$取了$1000$。

还有就是要注意取整不能直接用$int(A*j)$，可以使用$int(A*j+EPS)$，或者使用函数$floor(A*j)$。
浮点误差会导致答案错误。

小结

• 给定很少输入，求$1$到$N$的某个函数值的问题，通常有两种方法：直接求、递推。

• 浮点数取整直接用$(int)a$会因为浮点误差带来结果的偏差，应使用$floor(a)$更精确。

• 关于这类估算的问题
  ①特征：“结果的相对误差不超过5%”之类的话，再给个SPJ；
  ②原则：能精确尽可能精确；
  ③对于不超时的精确部分，就让它保留，对于会超时的精确部分，才考虑估算；
  ④在用递推的方法时，对于小数据可以先暴力求，对于大数据再估算；
  ⑤估算的一种常见方法：取平均数、中位数、众数之类的。

代码

#include <cstdio>
#include <cmath>

const int N=100010;
const int T=1000;
const double EPS=1e-7;

int n;
double a;

double m[N];
int g[N];

double w[N];

int main(void)
{
    scanf("%d%lf",&n,&a);

    for (int i=1;i<=n;i++) scanf("%lf",&m[i]);
    for (int i=1;i<=n;i++) g[i]=floor(a*i);

    for (int i=1;i<=T;i++)
    {
        for (int j=1;j<=g[i];j++)
            w[i]+=m[j]/(i-j);
        w[i]*=m[i];
    }

    double tmp;
    for (int i=1;i<=n-T;i++)
    {
        tmp=0;
        tmp=tmp+w[i]*(i-g[i]/2)/m[i]/(i+T-g[i]/2);
        for (int j=g[i]+1;j<=g[i+T];j++) tmp=tmp+m[j]/(i+T-j);
        w[i+T]=tmp*m[i+T];
    }

    for (int i=1;i<=n;i++)
        printf("%0.6lf\n",w[i]);

    return 0;
}