SPOJ FAVDICE_概率论DP

Each dog has its day.

题意：

抛一个n面的骰子，它的n个面都至少出现一次抛骰子的期望次数是多少。

Input

The first line of the input contains an integer t, the number of test cases. t test cases follow.

Each test case consists of a single line containing a single integer N (1 <= N <= 1000) - the number of sides on BuggyD's die.

Output

For each test case, print one line containing the expected number of times BuggyD needs to throw his N-sided die so that each number appears at least once. The expected number must be accurate to 2 decimal digits.

概率问题我实在是太弱了，碰到这种题总是纠结于计算概率或者期望的公式，而实际上这些问题大都可以用动态规划轻易解决，当然用概率论应该也是可以的，不过概率论都还没学完，而且DP的思想我也更熟悉。

动态规划将问题简化成了状态转移、初始状态、最终状态几个难点，同时建立一个好的模型会使这些问题也简单化。

本题定义dp[i]为已经出现了i个不同的面时n个面都出现的期望。那么此时再抛一次骰子，无疑两种结果，仍然出现n个面，出现了i+1个面，这两个事件出现的概率分别是 i / n 和 ( n - i ) / n 同时抛这次骰子也会增加一次抛骰子。

因此 dp[ i ] = i / n * dp[ i ] + (  n - i ) / n * dp[ i+1 ]

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cmath>
using namespace std;
#define mxn 1010
int n;
double dp[mxn];
int main(){
	int cs;
	scanf("%d",&cs);
	while(cs--){
		scanf("%d",&n);
		dp[n]=0.0;
		for(int i=n-1;i>=0;--i)
			dp[i]=dp[i+1]+(double)n/(n-i);
		printf("%.2lf\n",dp[0]);
	}
	return 0;
}