【jzoj4756】【幻象】

最新推荐文章于 2016-09-09 18:47:56 发布

原创最新推荐文章于 2016-09-09 18:47:56 发布 · 440 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#jzoj4756 #幻象

jzoj 同时被 2 个专栏收录

275 篇文章

订阅专栏

分治

4 篇文章

订阅专栏

本文介绍了一种计算在特定时间内，基于概率出现的幻象及其产生的幻象值的方法。使用了动态规划思想来求解每秒产生幻象值的期望，并提供了一种分治策略实现nlogn的时间复杂度。

题目大意

每秒他的思绪中都有一定的概率浮现出奇异的幻象,持续x秒的幻象将产生 x^2的幻象值。他想知道,在N秒内他期望产生的幻象值是多少。

解题思路

设L[i]为第i秒幻象的持续时间的期望.显然L[i]= (L[i-1] + 1) * a[i]%设f[i]表示前i秒的答案f[i] = f[i-1] + ((L[i-1] + 1)2–L[i-1]2) * a[i]%时间复杂度为O(N)

这里提供一个nlogn的做法，我们考虑分治。设过中点左边大小为a，发生概率为x，右边大小为b，发生概率为y。贡献为(a+b)^2xy=>xa^2+xyb^2+2xay*b所以我们可以枚举左端点，右边的贡献可以预处理。

code

#include<set>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define LF double
#define LL long long
#define max(a,b) ((a>b)?a:b)
#define min(a,b) ((a>b)?b:a)
#define fo(i,j,k) for(LL i=j;i<=k;i++)
#define fd(i,j,k) for(LL i=j;i>=k;i--)
using namespace std;
int const maxn=1000000,maxm=100,inf=2147483647;
int n,m;
double a[maxn+10],ans;
int read(){
    char ch=getchar();int val=0;
    for(;(ch<'0')||(ch>'9');ch=getchar());
    for(;(ch>='0')&&(ch<='9');val=val*10+ch-'0',ch=getchar());
    return val;
}
void dfs(LL l,LL r){
    if(l==r){
        ans+=a[l]*(1-a[l-1])*(1-a[l+1]);
        return;
    }
    LL m=(l+r)/2;double tmp=1,p=0,pr=0,prr=0;
    fo(i,m+1,r){
        tmp*=a[i];
        p+=tmp*(1-a[i+1]);
        pr+=tmp*(1-a[i+1])*(i-m);
        prr+=tmp*(1-a[i+1])*(i-m)*(i-m);
    }tmp=1;
    fd(i,m,l){
        tmp*=a[i];
        ans+=tmp*(1-a[i-1])*(m-i+1)*(m-i+1)*p;
        ans+=tmp*(1-a[i-1])*(m-i+1)*pr*2;
        ans+=tmp*(1-a[i-1])*prr;
    }
    dfs(l,m);dfs(m+1,r);
}
int main(){
    freopen("d.in","r",stdin);
    freopen("d.out","w",stdout);
    n=read();
    fo(i,1,n)
        a[i]=read()/100.0;
    dfs(1,n);
    printf("%.1lf",ans);
    return 0;
}