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 }