当年讲概率和期望的时候我没在啊,,都是自学乱讲的,,不知道今天有多少人被搞晕了qwq
概率这个大概都会,期望就是所有概率*贡献的和。反正就是象征性地加起来。
主要是做题。然后有个小结论就是一般来说概率正推期望反推。
因为你在起点的概率是确定的1,但到终点就不知道了。
同样,你从终点到终点期望是已知的0,但从1过去就不知道了。
实际上反例也多,比如高斯消元来解什么的,,实际上还是靠做题积累经验。
放东西问题,因为每块向左向右倒的概率是一样的,可以方便我们来从中间状态入手。
假设我要放第k块东西且左右都放好了。
那期望分三种:向左倒,向右倒,不倒。同时还要加上左右都放好的期望。
不倒的概率是1-l-r,期望值放置数量就是1/(1-l-r)。
所以总期望是
如果我们考虑前i块都放好,枚举中间倒下的位置k,就是
这样是n方,感性理解一下发现:当i逐渐增大,如果j不变,则E(i-j-1)逐渐增大,而E(j-1)却不变,所以答案会越来越坏。所以我们可以感性得到j这个决策点是单调的,可以直接弄个单调的数据结构维护一下(实际上是个下凸壳)
#include<bits/stdc++.h>
using namespace std;
#define in read()
int in{
int cnt=0,f=1;char ch=0;
while(!isdigit(ch)){
ch=getchar();if(ch=='-')f=-1;
}
while(isdigit(ch)){
cnt=cnt*10+ch-48;
ch=getchar();
}return cnt*f;
}
double l,r;
double f[1003];int n;
double F(int a,int b){
return f[b]+f[a-b-1]+(1+f[b]*l+f[a-b-1]*r)/(1-l-r);
}
double solve(){
for(int i=1,j=0;i<=n;i++){
while(j<i-1&&F(i,j)>F(i,j+1))++j;
f[i]=F(i,j);
} return f[n];
}
signed main(){
while(n=in){
scanf("%lf %lf",&l,&r);
printf("%.2lf\n",solve());
}
return 0;
}
反正题上说残图最后能装走就行。我们就不管那么多直接三维动态规划。
打到第i场,胜利了j场,反正一个残图只有-1的背包消耗,开个上界即可。然后整体挪动200位。
最后把合法的加起来。
#include<bits/stdc++.h>
using namespace std;
double f[203][203][403];
int n,l,k;
double p[203];
int a[203];
signed main(){
cin>>n>>l>>k;
for(int i=1;i<=n;i++){
scanf("%lf",&p[i]);p[i]/=100.0;
}
for(int i=1;i<=n;i++)cin>>a[i];
if(k>n)k=n;
f[0][0][k+200]=1;
for(int i=0;i<n;i++){
for(int j=0;j<=i;j++){
for(int kk=200-i;kk<=200+n;kk++){
f[i+1][j][kk]+=f[i][j][kk]*(1-p[i+1]);
f[i+1][j+1][min(kk+a[i+1],200+n)]+=f[i][j][kk]*p[i+1];
}
}
}double ans=0;
for(int j=l;j<=n;j++){
for(int kk=200;kk<=200+n;kk++){
ans+=f[n][j][kk];
}
}printf("%.6lf",ans);
return 0;
}
同样枚举现在已经获得了i种宝石j个洞口,然后倒着推期望。注意要乘上自己对答案有贡献的期望。不能单纯地乘概率。
注意double运算。
#include<bits/stdc++.h>
using namespace std;
#define in read()
double n,s;
double f[1003][1003];
signed main(){
scanf("%lf %lf",&n,&s);
for(int i=n;i>=0;i--){
for(int j=s;j>=0;j--){
if(j==s&&i==n)continue;
if(j+1<=s)f[i][j]+=f[i][j+1]*((i/n)*(1-j/s))/(1-(i/n)*(j/s));
if(i+1<=n)f[i][j]+=f[i+1][j]*((1-i/n)*(j/s))/(1-(i/n)*(j/s));
if(i+1<=n&&j+1<=s)f[i][j]+=f[i+1][j+1]*((1-i/n)*(1-j/s))/(1-(i/n)*(j/s));
f[i][j]+=1/(1-(i/n)*(j/s));
}
}printf("%.4lf",f[0][0]);
return 0;
}