【jzoj4756】【幻象】

题目大意

每秒他的思绪中都有一定的概率浮现出奇异的幻象,持续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)^2*x*y=>x*a^2+x*y*b^2+2*x*a*y*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;
}