One Person Game（ZOJ3329）

题目链接：http://icpc.moe/onlinejudge/showProblem.do?problemCode=3329

题目意思是，给三个骰子，骰子分别有K1，K2，K3个面，另外设定a,b,c。如果摇到的点数正好是a,b,c那么现有的分数清零。如果不是那么计入当前的分数，分数到达n了之后游戏结束，问结束的次数期望

说真的按照先前的思路。必然定义出dp[i]，表示摇到这个点结束的期望。pk是摇到这个点的概率，p0是归零的概率

然后dp[i]=dp[i]=∑(pk*dp[i+k])+dp[0]*p0+1;

然后就不会了，因为dp[i+k]里边也有dp[0]，而且我们求得就是dp[0],因此化成这样就没办法操作了。

然后看了大神的题解。顿时觉得好牛逼。

http://www.cnblogs.com/kuangbin/archive/2012/10/03/2710648.html

用的是代换的思想

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
double p[20];//摇到点数之和为j的概率
double A[550],B[550];//一开始写成505wa了好几次
int t,n,k1,k2,k3,a,b,c;
double p0;
int main()
{
    cin>>t;
    while(t--)
    {
        memset(A,0,sizeof(A));
        memset(B,0,sizeof(B));
        memset(p,0,sizeof(p));

        cin>>n>>k1>>k2>>k3>>a>>b>>c;

        p0=1.0/k1/k2/k3;//摇到的概率

        for(int i=1; i<=k1; i++)
            for(int j=1; j<=k2; j++)
                for(int k=1; k<=k3; k++)
                    if(i!=a||j!=b||k!=c)
                        p[i+k+j]+=p0;
        //每次加上这个概率最终就是摇到这个点的概率

        for(int i=n; i>=0; i--)
        {
            A[i]=p0;
            B[i]=1;
            for(int j=1; j<=k1+k2+k3; j++)
            {
                A[i]+=p[j]*A[i+j];
                B[i]+=p[j]*B[i+j];
            }
        }
        printf("%.16lf\n",B[0]/(1-A[0]));
    }
    return 0;
}