[spoj Favorite Dice ][期望dp]

（1）https://vjudge.net/problem/SPOJ-FAVDICE

题意：有一个n面的骰子，每一面朝上的概率相同，求所有面都朝上过至少一次的总次数期望。

题解：令dp[i]表示 i 面满足条件的期望次数，则有 dp[i]=$sum_{j=1}^{i-1}$(Pj*(本次操作对最终期望的贡献+dp[i]))+$sum_{j=1}^{n-(i-1)}$(Qj*(本次操作对最终期望的贡献+dp[i-1]))，其中Pj表示出现已经出现过数字的概率，由于是等概率事件，所以这里所有的Pj都是1/n，同理Qj都是1/n，这道题目求的是次数，所以每次操作对最终期望的贡献都是1，所以这道题写下来就是dp[i]=$sum_{j=1}^{i-1}$((1/n+dp[i])+$sum_{j=1}^{n-(i-1)}$((1/n))，化简整理可得转移方程dp[i]=dp[i-1]+n/(n-(i-1))

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<vector>
#include<queue>
#include<map>
using namespace std;
//#define io_test
#define debug(x) cout<<x<<"$$$$"<<endl;
typedef long long ll;
double dp[1005];
int main()
{
#ifdef io_test
    freopen("in.txt","r",stdin);
    freopen("out.txt","w",stdout);
#endif // io_test
    int t;
    scanf("%d",&t);
    while(t--){
        int n;
        scanf("%d",&n);
        dp[1]=1;
        for(int i=2;i<=n;i++){
            dp[i]=dp[i-1]+1.0*n/(n-(i-1));
           // debug(dp[i]);
        }
        printf("%.2lf
",dp[n]);
    }
    return 0;
}

（2）http://2050.acmclub.cn/contests/contest_showproblem.php?pid=1001&cid=2

题意：有n+m条路，其中n条是正确的路，走每条路消耗的时间为ai，正确的路通往终点，剩下m条为错误的路，走每条路消耗的时间为bi，求到终点用时的期望是否大于y

题解：令dp[i]表示第i次找到正确的路的期望时间，显然dp[0]=0，而最后要求的时间期望就是dp[1]，则效仿上题有dp[1]=$sum_{i=1}^{m}$(Pi*(bi+dp[1]))+$sum_{i=1}^{n}$(Qi*(ai+dp[0]))，由于是等概率事件，所以Pi和Qi都是1/(n+m)，所以dp[1]=$sum_{i=1}^{m}$((1/(n+m))*(bi+dp[1]))+$sum_{i=1}^{n}$((1/(n+m)*ai))，化简得到dp[1]=($sum_{i=1}^{n}$ai+$sum_{i=1}^{m}$bi)/n

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<vector>
#include<queue>
#include<map>
using namespace std;
//#define io_test
#define debug(x) cout<<x<<"$$$$"<<endl;
typedef long long ll;
int a[15],b[15];
int main()
{
#ifdef io_test
    freopen("in.txt","r",stdin);
    freopen("out.txt","w",stdout);
#endif // io_test
    int t;
    scanf("%d",&t);
    while(t--){
        int n,m,y;
        ll sum1=0;
        ll sum2=0;
        scanf("%d%d%d",&n,&m,&y);
        for(int i=1;i<=n;i++){scanf("%d",&a[i]);sum1+=a[i];}
        for(int i=1;i<=m;i++){scanf("%d",&b[i]);sum2+=b[i];}
        if(sum1+sum2>n*y){
            printf("Wait
");
        }
        else{
            printf("Go
");
        }
    }
    return 0;
}