Collecint Bugs

题目：
http://poj.org/problem?id=2096

一个软件有$s$个子系统，会产生$n$种$bug$。某人一天能够发现一个$bug$，这个$bug$属于一个子系统的概率为$\frac{1}{s}$，属于某一类$bug$的概率为$\frac{1}{n}$。问发现$n$种$bug$，并且每个子系统也都发现$bug$的天数的期望。

思路：
$dp[i][j]$表示已经找到$i$种$bug$，并且这些$bug$属于$j$个系统，达到目标状态的天数的期望。有
$$\begin{array}{rl}dp[i][j]& =1+\frac{i}{n}\ast \frac{j}{s}dp[i][j]+\frac{i}{n}\ast \frac{s-j}{s}dp[i][j+1]\\ & +\frac{n-i}{n}\ast \frac{j}{s}dp[i+1][j]+\frac{n-i}{n}\ast \frac{s-j}{s}dp[i+1][j+1]\\ (1-\frac{i\ast j}{s\ast n})\ast dp[i][j]& =1+\frac{i}{n}\ast \frac{s-j}{s}dp[i][j+1]+\frac{n-i}{n}\ast \frac{j}{s}dp[i+1][j]\\ & +\frac{n-i}{n}\ast \frac{s-j}{s}dp[i+1][j+1]\end{array}$$
$dp[n][s]=0$

#include"cstdlib"
#include"cstdio"
#include"cstring"
#include"cmath"
#include"queue"
#include"algorithm"
#include"iostream"
#define eps 1e-5
using namespace std;
double dp[1020][1020];
int main()
{
    
    
    int n,s;
    while(cin>>n>>s)
    {
    
    
        memset(dp,0,sizeof(dp));
        for(int i=n; i>=0; i--)
        {
    
    
            for(int j=s; j>=0; j--)
            {
    
    
                if(i==n&&j==s) continue;
                double tep=1.0;
                if(i<n) tep+=(n*1.0-i)*j/(n*s)*dp[i+1][j];
                if(j<s) tep+=(s*1.0-j)*i/(n*s)*dp[i][j+1];
                if(i<n&&j<s) tep+=(n*1.0-i)*(s*1.0-j)/(n*s)*dp[i+1][j+1];
                dp[i][j]=tep/(1-(i*j*1.0)/(s*n));
            }
        }
        printf("%.4f\n",dp[0][0]);
    }
    return 0;
}