poj 2096 马克科夫链

/*重点是理解dp数组的含义。Poj 的collect bugs.题意就是给n种bug和s个子系统.发现一个bug需要一天的时间.求期望的天数.主要是dp[i][j]的构造.它表示当前已经在j个子系统内发现n个bug,里最后结果还相差的期望天数.一次推出四个转移方程.如dp[i+1][j]表示新发现了一种bug.则对应的概率为(n-i)*j/(n*s*1.0).表示在还剩n-i的bug里除以总数n.子系统还是没有更新.即为j/s.以此类推/*/#include <stdio.h>double dp[1010][1010];int main(){    int n,s;    while(scanf("%d%d",&n,&s)==2)    {        dp[n][s]=0;        for(int i=n;i>=0;i--)        {            for(int j=s;j>=0;j--)            {                double p1=1.0*i*j/(n*s*1.0);                double p2=1.0*(n-i)*j/(n*s*1.0);                double p3=1.0*i*(s-j)/(n*s*1.0);                double p4=1.0*(n-i)*(s-j)/(n*s*1.0);                if(i==n&&j==s) continue;//上面一个式子必须要转化,不能把待求量当做已知条件.                //dp[i][j]=(dp[i][j]+1)*p1+(dp[i+1][j]+1)*p2+(1+dp[i][j+1])*p3+(1+dp[i+1][j+1])*p4;                dp[i][j]=(dp[i+1][j]*p2+dp[i][j+1]*p3+dp[i+1][j+1]*p4+1)/(1-p1)*1.0;            }        }        printf("%.4lf\n",dp[0][0]);    }    return 0;}