bzoj1011: [HNOI2008]遥远的行星

这题好神，完全接受不能啊；

在看题解之前，我千方百计想着怎么精确地算出结果，但苦思之后决定看题解，没想到题解是估计近似；

个人题解：题目中的要求是结果误差不能超过5%，这个条件很关键；

可以设置一个范围，在这个范围之内精确计算，超过这个范围就近似；

虽然我的感觉是很可能被卡（我想的是有可能通过1 1 1 1000000 1 1 1这样的方式卡掉），

但实际上AC了，所以在考试时要大胆的假设，这种估算法在某些题比较关键；

代码其实都类似：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cstdlib>
#include<ctime>
#include<vector>
#include<algorithm>
#include<queue>
#include<map>
#include<vector>
#include<cmath>
using namespace std;
#define LL long long
#define up(i,j,n) for(int i=j;i<=n;i++)
#define down(i,n,j) for(int i=n;i>=j;i--)
const double eps=1e-8;
const int maxn=101000;
int n;
double m[maxn],A,f[maxn],pre[maxn];
void init(){
    freopen("1.in","r",stdin);
    freopen("1.out","w",stdout);
    scanf("%d%lf",&n,&A);
    up(i,1,n)scanf("%lf",&m[i]);
}
void work(){
    int T=20;
    for(int i=1;i<=T;i++){
        int t=(int)(eps+floor(A*i));
        for(int j=1;j<=t;j++)f[i]+=m[j]/(i-j)*m[i];
    }
    for(int i=T+1;i<=n;i++){
        int r=(int)(eps+floor(A*i));
        int l=(int)(eps+floor(A*(i-T)));
        for(int j=l+1;j<=r;j++)f[i]+=m[j]*m[i]/(i-j);
        f[i]+=m[i]*(f[i-T]/m[i-T])*(i-T-1.0*l/2)/(i-1.0*l/2);
    }
    up(i,1,n)printf("%.6lf\n",f[i]);
}
int main(){
    init();
    work();
}