BZOJ1011

最新推荐文章于 2018-03-26 22:43:30 发布

Tgop_Knight

最新推荐文章于 2018-03-26 22:43:30 发布

阅读量660

点赞数

分类专栏： BZOJ 文章标签：乱搞

本文链接：https://blog.csdn.net/Tgop_Knight/article/details/47425173

版权

BZOJ 专栏收录该内容

11 篇文章 0 订阅

订阅专栏

1011: [HNOI2008]遥远的行星

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

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

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

代码：

#include <cstring>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#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;
}