BZOJ1011

1011: [HNOI2008]遥远的行星

思路:乱搞题
不要想着算出正确答案……

对于每个行星i,我们计算K=A*i
如果K比较小,我们直接暴力计算每一个行星对行星i的答案贡献
那K比较大呢?

不虚,我们观察计算式中每一个分母
如果K很大,而A又很小,所以i很大
所以我们如果把所有的(i-j)换成(i-k/2),答案差不了多少
然后就可以用前缀和加速计算啦!

代码:

#include 
#include 
#include 
#include 
#include 
#include 
#define eps 1e-8
using namespace std;
const int MAXN=100010;
int n;
double a;
double sum[MAXN];
double m[MAXN];
void init()
{
    cin>>n>>a;
    for (int i=1;i<=n;i++) {
        scanf("%lf",&m[i]);
        sum[i]=sum[i-1]+m[i]; }
    for (int i=1;i<=n;i++) {
        int k=(int)(a*(double)(i)+eps);
        double ans=0;
        if (k<=67) {
            for (int j=1;j<=k;j++)
                ans+=(m[i]*m[j])/(double)(i-j); }
        else ans=(sum[k]*m[i])/(double)(i-((1+k)>>1));
        printf("%.6lf\n",ans); }
    return ;
}
int main()
{

    init();

    return 0;
}

你可能感兴趣的:(BZOJ,乱搞)