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 }

poj 2096