HDU 3076 ssworld VS DDD

这是一道很简答的概率DP，但是我刚开始做的时候自己一个人想 还是一头雾水，，orz，看了别人的博客才晓得原来是这么回事，我最开始的思路是错的，不过可以过很多数据，就不说了，讲讲大神的思路，设p1,p2,p表示第一个人赢，第二个人赢，两者平局。当你要计算谁可能赢的时候，你要计算一局里面赢得概率，那么ssworld 想要赢一局他的骰子的点数必须要多，那么当他是2点的时候 DDD只能是1点，当他是3点的时候他只能是1或2，，以此类推。。这样就可以求出 两者赢一局的概率，那么在整个过程之中某人赢一局的概率是怎样的？无疑对N局来说我赢一局， 你们，没赢，那就说明我们剩下的都是平局，哈哈，所以再设片q1，q2表示整个过程之中赢一局的概率那就是对ssworld 我们有：p1+p1*p+p1*p^2+......+p*p1^n(n->无穷)最后就等于p1/(1-p).所以q1=p1/（1-p）讨论一下p就行，最后就是算他们的概率了设dp[i][j]表示ssworld 赢j，次DDD赢了i次， 他的由来有两方面组成：ssworld 从j-1次赢了一次到j，DDD从j-1次赢了一次到j次。那么就有了状态方程。最后想加就好。最后杭电上的血量反了，颠倒过来就好

#include<iostream>#include<cstring>#include<stdio.h>using namespace std;double a[10],b[10];double f[2005][2005];//表示第一个赢j次第二个赢i次int main(){    int n,m;    int i,j,k;    double p1,p2,p;    while(cin>>m>>n)    {        p1=0,p2=0,p=0;        memset(a,0,sizeof(a));        memset(b,0,sizeof(b));        memset(f,0,sizeof(f));        for(i=1;i<=6;i++)            scanf("%lf",&a[i]);        for(i=1;i<=6;i++)            scanf("%lf",&b[i]);        for(i=2;i<=6;i++)            for(j=1;j<i;j++)            p1+=a[i]*b[j];        for(j=2;j<=6;j++)            for(i=1;i<j;i++)            p2+=a[i]*b[j];        p=1-p1-p2;        double q1;        double q2;       //cout<<p<<endl;        if(p==1){q1=0,q2=0;}        else            {                q1=p1/(1-p);                q2=p2/(1-p);            }           // cout<<p1<<" "<<p2<<endl;        f[0][0]=1;        for(j=0;j<m;j++)        for(i=0;i<=n;i++)            {                if(j>0)f[i][j]+=f[i][j-1]*q1;                if(i>0)f[i][j]+=f[i-1][j]*q2;            }            double max=0;        for(i=0;i<n;i++)            max+=f[i][m-1]*q1;            //double M=0.5;            //printf("%.6lf\n",M);            if(max>1)max=1;        printf("%.6lf\n",max);    }    return 0;}