Time Limit:1000MS Memory Limit:32768KB 64bit IO Format:%lld & %llu
Description
Given a dice with n sides, you have to find the expected number of times you have to throw that dice to see all its faces at least once. Assume that the dice is fair, that means when you throw the dice, the probability of occurring any face is equal.
For example, for a fair two sided coin, the result is 3. Because when you first throw the coin, you will definitely see a new face. If you throw the coin again, the chance of getting the opposite side is 0.5, and the chance of getting the same side is 0.5. So, the result is
1 + (1 + 0.5 * (1 + 0.5 * ...))
= 2 + 0.5 + 0.52 + 0.53 + ...
= 2 + 1 = 3
Input
Input starts with an integer T (≤ 100), denoting the number of test cases.
Each case starts with a line containing an integer n (1 ≤ n ≤ 105).
Output
For each case, print the case number and the expected number of times you have to throw the dice to see all its faces at least once. Errors less than 10-6 will be ignored.
Sample Input
5
1
2
3
6
100
Sample Output
Case 1: 1
Case 2: 3
Case 3: 5.5
Case 4: 14.7
Case 5: 518.7377517640
故设dp[i]为在已经扔出了i个不同面的情况下扔出n个不同面的期望次数,dp[n]=0,答案为dp[0]
则\(dp[i]=dp[i]*\frac{i}{n}+dp[i+1]*\frac{n-i}{n}+1\)
移项得：\(dp[i]=dp[i+1]+\frac{n}{n-i}\)
#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
#define N 100010 int main()
{
int T,n,iCase=;
double dp[N];
scanf("%d",&T);
while(T--)
{
scanf("%d",&n);
dp[n]=;
for(int i=n-;i>=;i--) dp[i]=dp[i+]+n*1.0/(n-i);
printf("Case %d: %.7f\n",iCase++,dp[]);
}
return ;
}