2018.09.28【BZOJ4318】OSU!（期望DP）

传送门

---

解析：

不错的期望DP入门题。

思路：

首先，要是没有明白期望的线性性。这道题是无法理解的。

一个常识，要维护高次，我们可以将高次展开，利用低次来维护。

所以我们同时维护长度的期望，长度平方的期望，长度立方的期望。

注意这里可能就有人不懂了，长度平方的期望不等于长度的期望的平方。
同样，长度立方的期望不等于长度的期望的立方。

我们用$len{1}_i,len{2}_i,len{3}_i$表示上述三个参数。

则我们将它展开得到$$len{2}_i=len{2}_{i-1}+2\ast len{1}_{i-1}+1$$$$len{3}_i=len{3}_{i-1}+3\ast len{2}_{i-1}+3\ast len{1}_{i-1}+1$$

套上一个概率就行了。
然后我们的答案就是所有$len3$的和，然而$len3$不会用于修改其他的项，我们可以直接在$len3$上维护前缀和。

---

代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define re register
#define gc getchar
#define pc putchar
#define cs const

inline
ll getint(){
	re ll num=0;
	re char c;
	while(!isdigit(c=gc()));num=c^48;
	while(isdigit(c=gc()))num=(num<<1)+(num<<3)+(c^48);
	return num;
}

inline
double getdb(){
	re double x=0,y=1.0;
	re char c;
	while(!isdigit(c=gc()));
	x=c^48;
	while(isdigit(c=gc()))x=(x*10)+(c^48);
	if(c!='.')return x;
	while(isdigit(c=gc()))x+=(y/=10)*(c^48);
	return x;
}

cs int N=100002;
double len1[N],len2[N],len3[N];
int n;
signed main(){
	n=getint();
	for(int re i=1;i<=n;++i){
		double p=getdb();
		len1[i]=p*(len1[i-1]+1);
		len2[i]=p*(len2[i-1]+2*len1[i-1]+1);
		len3[i]=len3[i-1]+p*(3*len2[i-1]+3*len1[i-1]+1);
	}
	printf("%.1f",len3[n]);
	return 0;
}