bzoj4318: OSU! (概率+dp)

题解

f[i] f [ i ] 为以 i i 结尾的后缀1长度
g[i] 为以 i i 结尾的后缀1长度的平方
h[i] 为以 i i 的答案

f[i]=(f[i1]+1)p[i]
因为 (x+1)2=x2+2x+1 ( x + 1 ) 2 = x 2 + 2 x + 1 ,可知 g[i]=(g[i1]+2f[i1]+1)p[i] g [ i ] = ( g [ i − 1 ] + 2 ∗ f [ i − 1 ] + 1 ) ∗ p [ i ]
因为 (x+1)3=x3+3x2+3x+1 ( x + 1 ) 3 = x 3 + 3 x 2 + 3 x + 1 ,可知 h[i]=(h[i1]+3g[i1]+3f[i1]+1)p[i]+h[i1](1p[i]) h [ i ] = ( h [ i − 1 ] + 3 ∗ g [ i − 1 ] + 3 ∗ f [ i − 1 ] + 1 ) ∗ p [ i ] + h [ i − 1 ] ∗ ( 1 − p [ i ] ) (h为所求,因此要加上以0结尾的答案)

代码

#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
using namespace std; 
#define N 100010
int n,T;
double f[N],g[N],h[N],p[N];
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%lf",&p[i]);
        f[i]=(f[i-1]+1)*p[i];
        g[i]=(g[i-1]+2*f[i-1]+1)*p[i];
        h[i]=(h[i-1]+3*g[i-1]+3*f[i-1]+1)*p[i]+h[i-1]*(1-p[i]);
    }
    printf("%.1lf",h[n]);
    return 0;
} 
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值