POJ 2096 Collecting Bugs 概率dp

题意：

有n种bug存在于s个系统中，每次可以找到一个bug问在s个系统中找到n种bug且每个系统最少找到一种的概率是多少

分析：

因为bug的数量是无限所以每次都有n*s种可能

可设dp【i，j】表示当前已经在j种系统中发现i种bug了

那么对于dp【i，j】有下面五种后继dp【i，j】 dp【i+1，j】 dp【i，j+1】dp【i+1,j+1】

坑点是输出用 .4f

ACcode:

#pragma warning(disable:4786)//使命名长度不受限制
#pragma comment(linker, "/STACK:102400000,102400000")//手工开栈
#include <map>
#include <set>
#include <queue>
#include <cmath>
#include <stack>
#include <cctype>
#include <cstdio>
#include <cstring>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
#define rd(x) scanf("%d",&x)
#define rd2(x,y) scanf("%d%d",&x,&y)
#define rds(x) scanf("%s",x)
#define rdc(x) scanf("%c",&x)
#define ll long long int
#define ull unsigned long long
#define maxn 1010
#define mod 1000000007
#define INF 0x3f3f3f3f //int 最大值
#define FOR(i,f_start,f_end) for(int i=f_start;i<=f_end;++i)
#define MT(x,i) memset(x,i,sizeof(x))
#define PI  acos(-1.0)
#define E  exp(1)
#define eps 1e-8
ll gcd(ll a,ll b){return b==0?a:gcd(b,a%b);}
ll mul(ll a,ll b,ll p){ll sum=0;for(;b;a=(a+a)%p,b>>=1)if(b&1)sum=(sum+a)%p;return sum;}
inline void Scan(int &x) {
      char c;while((c=getchar())<'0' || c>'9');x=c-'0';
      while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
}
using namespace std;
double dp[maxn][maxn];
int main(){
    int n,loop,cnt=1,s;
    while(scanf("%d%d",&n,&s)!=EOF){
        dp[n][s]=0;
        double t=1.0/(n*s);
        for(int i=n;i>=0;--i)for(int j=s;j>=0;j--)if(i!=n||j!=s)
            dp[i][j]=(dp[i+1][j]*(n-i)*j*t+dp[i][j+1]*i*(s-j)*t+dp[i+1][j+1]*(n-i)*(s-j)*t+1)/(1.0-i*j*t);
        printf("%.4f\n",dp[0][0]);
    }
    return 0;
}
/*
1 2
1 4
29 88
2 5
2 23
1 1
*/