poj 2096

 1 /**
 2 程序员调bug
思路： 一共有四种情况，
1、 1个原有的bug 在原有的分类中
2、 1个原有的bug 在新的分类中
3、 1个新的bug 在原有的分类中
4、 1个新bug  在新的分类中


 3 **/
 4 #include <iostream>
 5 #include <cstdio>
 6 using namespace std;
 7 double dp[1010][1010];
 8 int main()
 9 {
10     int n,s;
11     cin>>n>>s;
12     for(int i=n;i>=0;i--){
13         for(int j=s;j>=0;j--){
14             if(i==n&&j==s)
15                 continue;
16             double p1 = (double(s-j)*i)/n/(s);
17             double p2 = (double(n-i)*j)/n/(s);
18             double p3 = (double(s-j)*(n-i))/n/(s);
19             double p0 = 1.0-(double(j*i))/n/(s);
20             dp[i][j] = p1*dp[i][j+1]+p2*dp[i+1][j]+p3*dp[i+1][j+1]+1;
21             dp[i][j] /= p0;
22         }
23     }
24     printf("%.4lf\n",dp[0][0]);
25     return 0;
26 }