P1654 OSU!

最新推荐文章于 2024-08-19 11:12:00 发布

Snow_raw

最新推荐文章于 2024-08-19 11:12:00 发布

阅读量325

点赞数

分类专栏：概率与期望

本文链接：https://blog.csdn.net/qq_53504307/article/details/127432430

版权

概率模型字符串处理期望分数转移方程动态规划

关键词由CSDN通过智能技术生成

概率与期望专栏收录该内容

4 篇文章 0 订阅

订阅专栏

P1654 OSU!

题意： 给定一串字符串，每一串字符串中的每一位仅存在 $0$ 和 $1$ 两种情况，且出现 $1$ 的概率为 $p [i]$ ，在这个串中了连续的 $X$ 的 $1$ 可以贡献 $X^3$ 的分数。这 $X$ 个 $1$ 不能被其他连续的 $1$ 所包含，求期望分数。

思路： 首先我们不知道怎么求 $X^3$ 的期望分数。我们来推导一下式子， $X+1)^3$ $=$ $X^3$ $+$ $3*X^2$ $+$ $3 * X$ $+ 1$ ，那么我们可以发现多的新的贡献是 $3*X^2$ $+$ $3 * X$ $+ 1$ ，且此贡献仅在前一位为 $1$ 的情况下，当前为概率为 $p [i]$ 的情况下可以转移获得。若当前位为 $0$ ，既 $(1 - p [i])$ $*$ $an s [i - 1]$ ， $an s [i - 1]$ 表示上一层 $X^3$ 的总分数期望。

因此转移方程：
$X^1$ $=$ $X [i]$ $=$ $(X [i - 1] + 1) * p [i]$
$X^2$ $=$ $Y [i]$ $=$ $(Y [i - 1] + 2 * X [i - 1] + 1) * p [i]$
$X^3$ $=$ $(an s [i - 1] + 3 * Y [i - 1] + 3 * X [i - 1] + 1) * p [i]$

三者表示的都是前一位为 $1$ 情况下所能获得的分数期望。
$an s [i] = (an s [i - 1] + 3 * Y [i - 1] + 3 * X [i - 1] + 1) * p [i] + an s [i - 1] * (1 - p [i]) = an s [i - 1] + (3 * Y [i - 1] + 3 * X [i - 1] + 1) * p [i]$

Code:

#include<bits/stdc++.h>
using namespace std;
#define int long long 
const int N = 1e5+10;
int n;
double dp[N];
double x[N];
double y[N];
signed main(){
	cin>>n;
	for(int i=1;i<=n;i++){
		double xx;
		cin>>xx;
		x[i]=(x[i-1]+1)*xx;
		y[i]=(y[i-1]+2*x[i-1]+1)*xx;
		dp[i]=dp[i-1]+(3*y[i-1]+3*x[i-1]+1)*xx;
	}
	printf("%.1lf",dp[n]); 
	return 0;
}