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;
}