poj 2096 Collecting Bugs（期望DP）

用dp[i][j]表示在i中类型的软件里发现了j种bug，第j种bug有四种可能：

1、一种已经发现的bug在已经发现bug的子系统中找到；

2、一种尚未发现的bug在已经发现bug的子系统中找到；

3、一种已经发现的bug在尚未发现bug的子系统中找到；

4、一种尚未发现的bug在尚未发现bug的子系统中找到；

可以得到状态转移方程：dp[i][j]=dp[i][j]*i/n*j/s+dp[i-1][j]*(n-i+1)/n*j/s+dp[i][j-1]*i/n*(n-j+1)/s+dp[i-1][j-1]*(n-i+1)/n*(s-j+1)/s;

化简这个方程： (1-ij/ns)*dp[i][j]=dp[i-1][j]*j(n-i+1)/ns+dp[i][j-1]*i(n-j+1)/ns+dp[i-1][j-1]*(n-i+1)(s-j+1)/ns;
dp[i][j]=j*(n-i+1)/(ns-ij)*dp[i-1][j]+i*(n-j+1)/(ns-ij)+(n-i+1)*(s-j+1)/(ns-ij);

（以上方程除法优先级高于乘法，ij表示i*j）；

将这个方程逆推出来：dp[i][j]=1.0*j*(n-i)/(n*s-i*j)*dp[i+1][j]+1.0*i*(s-j)/(n*s-i*j)*dp[i][j+1]+1.0*(n-i)*(s-j)/(n*s-i*j)*dp[i+1][j+1]+1.0*n*s/(n*s-i*j);

#include<stdio.h>
#include<string.h>
#define N 1005
double dp[N][N];
int main()
{
    int n,s;
    while(scanf("%d%d",&n,&s)!=EOF)
    {
        int i,j;
        for(i=0;i<=n+1;i++)
        {
            for(j=0;j<=s+1;j++)
                dp[i][j]=0;
        }
        for(i=n;i>=0;i--)
        {
            for(j=s;j>=0;j--)
            {
                if(i==n&&j==s)
                    continue;
                dp[i][j]=1.0*j*(n-i)/(n*s-i*j)*dp[i+1][j]+1.0*i*(s-j)/(n*s-i*j)*dp[i][j+1]+1.0*(n-i)*(s-j)/(n*s-i*j)*dp[i+1][j+1]+1.0*n*s/(n*s-i*j);
            }
        }
        printf("%.4f\n",dp[0][0]);
    }
    return 0;
}