[poj2096] Collecting Bugs【概率dp 数学期望】

传送门:http://poj.org/problem?id=2096

题面很长,大意就是说,有n种bug,s种系统,每一个bug只能属于n中bug中的一种,也只能属于s种系统中的一种。一天能找一个bug,问找到的bug涵盖所有种类的bug与所有种类的系统期望需要几天。

令f(i, j)为找到了i种bug,j种系统期望的天数,那么今天再找一个bug,有4种情况:

①,bug种类为已找到的i种中的一种,系统种类为已找到的j种中的一种,则概率p1 = (i / n) * (j / s)

②,bug种类为未找到的(n - i)种中的一种,系统种类为已找到的j种中的一种,则概率p2 = ((n - i) / n) * (j / s)

③,bug种类为已找到的i种中的一种,系统种类为未找到的(s - j)种中的一种,则概率p3 = (i / n) * ((s - j) / s)

④,bug种类为未找到的(n - i)种中的一种,系统种类为未找到的(s - j)种中的一种,则概率p3 = ((n - i) / n) * ((s - j) / s)

则有f(i, j) = f(i, j) * p1 + f(i + 1, j) * p2 + f(i, j + 1) * p3 + f(i + 1, j + 1) * p4 + 1

即f(i, j) = ( f(i + 1, j) * p2 + f(i, j + 1) * p3 + f(i + 1, j + 1) * p4 + 1 ) / (1 - p1)

#include <cstdio>

const int maxn = 1005, maxs = 1005;

int n, s;
double f[maxn][maxs];

int main(void) {
	scanf("%d%d", &n, &s);
	for (int i = n; ~i; --i) {
		for (int j = s; ~j; --j) {
			if (i == n && j == s) {
				 continue;
			}
			f[i][j] = (f[i + 1][j] * (1 - (double)i / (double)n) * ((double)j / (double)s) +
					  f[i][j + 1] * ((double)i / (double)n) * (1 - (double)j / (double)s) +
					  f[i + 1][j + 1] * (1 - (double)i / (double)n) * (1 - (double)j / (double)s) + 1) / (1 - (double)(i * j) / (double)(n * s));
		}
	}
	printf("%.4f\n", f[0][0]);
	return 0;
}

  

posted @ 2016-11-29 21:02  ciao_sora  阅读(156)  评论(0编辑  收藏  举报