做了几道期望题,头昏脑涨......
首先期望这个东西个人理解就是sigma概率*贡献,它满足线性性:对于两个不相关的随机变量x,y,E(x±y)=E(x)0±E(y),E(xy)=E(x)E(y),E(x/y)=E(x)/E(y).
然后还有一些有关期望奇怪的公式,不会。。。
例题:洛谷 P1297 [国家集训队]单选错位 https://www.luogu.org/problemnew/show/P1297
推一波可发现错位的某个位置对的期望就是它和它错位的位置之中选项更多的那个的选项数分之1,再根据期望线性性加起来即可。
#include<bits/stdc++.h>
#define db double
using namespace std;
const int N=1e7+10;
int a[N],n,A,B,C;
db ans=0;
int main()
{
int nw,nx;
scanf("%d%d%d%d%d",&n,&A,&B,&C,a+1);
for (int i=2;i<=n;i++)
a[i] = ((long long)a[i-1] * A + B) % 100000001;
for (int i=1;i<=n;i++)
a[i] = a[i] % C + 1;
for(int i=1;i<=n;i++)
{
nw=a[i];
if(i!=n)nx=a[i+1];
else nx=a[1];
if(nw<=nx)ans+=1.0/nx;
else ans+=1.0/nw;
}
printf("%.3lf",ans);
}
洛谷P1365 WJMZBMR打osu! / Easy 和 P1654 OSU!
https://www.luogu.org/problemnew/show/P1365
https://www.luogu.org/problemnew/show/P1654
两题基本相同,T1维护两个数组分别表示到该位期望长度以及答案数组,T2维护三个数组表示到该位期望长度,到该位期望(长度平方),答案数组。
转移即可。之前一直搞不清啥是平方的期望、期望的平方,但其实从期望定义上就可以说明,一个是p*(a^2),一个是(p*a)^2,所以注意运算的时候不能混淆,搞不清时就把期望写成概率*贡献的形式。
#include<bits/stdc++.h>
#define db double
using namespace std;
const int N=5e5+10;
int n;
char s[N];
db f[N],g[N];
int main()
{
scanf("%d",&n);
scanf("%s",s+1);
f[0]=0,g[0]=0;
for(int i=1;i<=n;i++)
{
if(s[i]=='x')f[i]=f[i-1],g[i]=0;
if(s[i]=='o')f[i]=f[i-1]+2*g[i-1]+1,g[i]=g[i-1]+1;
if(s[i]=='?')
{
g[i]=0.5*(g[i-1]+1.0);
f[i]=0.5*f[i-1]+0.5*(f[i-1]+1+2*g[i-1]);
}
}
printf("%.4lf",f[n]);
}
#include<bits/stdc++.h>
#define db double
using namespace std;
const int N=5e5+10;
int n;db a[N];
db f[N],g[N],h[N];
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
scanf("%lf",&a[i]);
f[0]=0,g[0]=0,h[0]=0;
for(int i=1;i<=n;i++)
{
g[i]=a[i]*(g[i-1]+1.0);
h[i]=a[i]*(h[i-1]+2*g[i-1]+1);
f[i]=(1.0-a[i])*f[i-1]+a[i]*(f[i-1]+3*h[i-1]+3*g[i-1]+1);
}
printf("%.1lf",f[n]);
}