【poj2096】 Collecting Bugs

http://poj.org/problem?id=2096 (题目链接)

题意

　　有一个程序，其中有s个子结构，每个子结构出bug的概率相等。bug总共分成n类，每种bug出现的概率相等。每天找出一个bug，求所有子结构都出现bug并且每种bug都出现过所需要的期望天数。

Solution

　　终于领会了期望dp的一点点。。。

细节

　　%.4lf用G++交会Wa。。

代码

// poj2096
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<queue>
#define LL long long
#define inf 1<<30
#define Pi acos(-1.0)
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
using namespace std;
 
const int maxn=1010;
double f[maxn][maxn];
int n,m;
 
int main() {
	scanf("%d%d",&n,&m);
	f[n][m]=0;
	for (int i=n;i>=0;i--)
		for (int j=m;j>=0;j--) {
			if (i+1<=n) f[i][j]+=f[i+1][j]*(1-(double)i/n)*(double)j/m;
			if (j+1<=m) f[i][j]+=f[i][j+1]*(1-(double)j/m)*(double)i/n;
			if (i+1<=n && j+1<=m) f[i][j]+=f[i+1][j+1]*(1-(double)i/n)*(1-(double)j/m);
			if (i!=n || j!=m) f[i][j]=(f[i][j]+1)/(1-(double)i/n*(double)j/m);
		}
	printf("%.4lf",f[0][0]);
	return 0;
}