标签：f1 P1365 期望 f2 len P1654 长度 CF235B include

期望长度

定义

这里期望长度表示一段序列连续长度的期望。具体来说，对于一段序列，每个点都有一个概率连续和断开。求所有连续序列和的期望。

当然，对于以上期望长度的定义，我们只需要求出每个点存在的期望的和即可。但是题目永远不会这么简单。

Osu!

Osu!是一个音乐游戏，玩家需要对音符在恰当时候进行敲击来通关。一次到位的敲击为o，不到位的为x。一段连续到位的敲击，即combo次数为这段序列的长度。

我们接下来讨论的三个题都和这个游戏有关。

---

• level1

一段Osu!序列为一串字符，包括'o','x','?'。其中'o','x'的定义如上，'?'表示此位置有一半的几率为'o'。游戏得分为所有combo次数平方的和。求得分的期望。

也就是我们要求所有序列长度平方的期望和。

期望

期望具有线性性，但不具有积性。这意味着我们无法对求得的期望长度直接平方来得到答案。

并且请注意一点，若一个值的期望为0，并不意味着它的平方的期望为0。这可以帮助我们理解期望的线性性。

期望的平方在大多数情况下并没有什么实际意义。

但是，期望具有线性性。

考虑我们的答案，实际上就是长度平方的期望。考虑往后的转移。（设\(f1\)表示当前期望长度，\(f2\)表示答案，即长度平方期望的和）

根据公式\((len+1)^2=len^2+2*len+1\)若后一位\(i\)为'o',则后一位\(i\)的期望值分别为

\(f1_i=f1_{i-1}+1\)

\(f2_i=f2_{i-1}+2*f1_{i-1}+1\)

即此位\(f2\)的值其实是可以从前一位线性转移来的。所谓线性，就是其幂为1。

同样，考虑第\(i\)位为'x'的情况，\(f1=0\),\(f2\)直接继承前面的答案。

然后我们就可以得到'?'的情况：上述两种情况和除以2.

\(f1_i=\frac{f1_{i-1}+1}2\)

\(f2_i=\frac{2*f2_{i-1}+2*f1_{i-1}+1}{2}\)

于是我们就能完成P1365

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cctype>
using namespace std;
inline int read(){
	int x=0,w=0;char c=getchar();
	while(!isdigit(c))w|=c=='-',c=getchar();
	while(isdigit(c))x=(x<<3)+(x<<1)+(c^48),c=getchar();
	return w?-x:x;
}
namespace star
{
	int n;
	double f,len,p;
	inline void work(){
		n=read();
		char c=getchar();
		while(c!='o' and c!='?' and c!='x')c=getchar(); 
		while(n--){
			if(c=='o')f=(f+2*len+1),len=len+1;
            else if(c=='?')f=(2*(f+len)+1)/2,len=(len+1)/2;
            else len=0;
		}
		printf("%.4lf",f);
	}
}
signed main(){
	star::work();
	return 0;
}

---

• level2

我们发现，其实对于概率任意的情况也可以推出来。上题三种字符其实就是对应概率为1,0.5,0的三种情况。

设\(p\)为该点为'o'的概率，则有:

\(f2_i=f2_{i-1}+p*(2*f1_{i-1}+1)\)

\(f1_i=p*(f1_{i-1}+1)\)

所以上题的代码的核心部分等同于：

	while(n--){
		p=c=='o'?1.0:c=='?'?0.5:0.0;
		f=f+p*(2*len+1);
		len=p*(len+1);
		c=getchar();
	}

于是我们可以完成CF235B

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cctype>
using namespace std;
inline int read(){
	int x=0,w=0;char c=getchar();
	while(!isdigit(c))w|=c=='-',c=getchar();
	while(isdigit(c))x=(x<<3)+(x<<1)+(c^48),c=getchar();
	return w?-x:x;
}
namespace star
{
	int n;
	double f,len,p;
	inline void work(){
		n=read();
		while(n--){
			scanf("%lf",&p);
			f=(f-len*len+(len+1)*(len+1))*p+f*(1-p);
			len=p*(len+1);
		}
		printf("%.8lf",f);
	}
}
signed main(){
	star::work();
	return 0;
}

---

• level3

我们已经完成了对于长度平方的期望和的问题。那么我们就可以解决新的问题：对于答案为所有combo长度立方的和的期望我们怎么求解呢？

根据期望的线性性，我们再维护一个平方的期望即可。

根据公式\((len+1)^3=len^3+3len^2+3len+1\)，我们可以得到以下转移：

\[f3_i=f3_{i-1}+p*(3*(f2_{i-1}+f1_{i-1})+1)\\ f2_i=p*(f2_{i-1}+2*f1_{i-1}+1)\\ f1_i=p*(f1_{i-1}+1); \]

注意！

我承认我的变量名的定义有亿点点毒瘤，因为读者可以清楚地发现在上一题中\(f2\)的转移为\(f2_i=f2_{i-1}+p*(2*f1_{i-1}+1)\)而非当前转移。实际上在定义\(f2\)时我的定义为答案而非二次项的期望，根据期望的线性性，答案是可以继承上一次的答案进行转移的，也就是对于'x'的情况继承\(f2\)而非\(0\)的原因。

在此level中我对\(f2\)重新定义为长度二次幂的期望。希望不要因为我的毒瘤误导大家。

相同的，我对\(f3\)的定义为答案，因此需要继承之前的答案。

于是我们就完成了P1654

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cctype>
using namespace std;
inline int read(){
	int x=0,w=0;char c=getchar();
	while(!isdigit(c))w|=c=='-',c=getchar();
	while(isdigit(c))x=(x<<3)+(x<<1)+(c^48),c=getchar();
	return w?-x:x;
}
namespace star
{
	int n;
	double f1,f2,f3;
	inline void work(){
		n=read();
		double p;
		while(n--){
			scanf("%lf",&p);
			f3=f3+p*(3*(f2+f1)+1);
			f2=p*(f2+2*f1+1);
			f1=p*(f1+1);
		}
		printf("%.1lf",f3);
	}
}
signed main(){
	star::work();
	return 0;
}

总结

期望好神奇。

标签：f1,P1365,期望,f2,len,P1654,长度,CF235B,include
来源： https://www.cnblogs.com/BrotherHood/p/13739292.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。