用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;
}