Little Pony and Alohomora Part I
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 416 Accepted Submission(s): 220
Problem Description
Trixie is a female unicorn pony and traveling magician, having been rumored to be "the most magical unicorn in all of Equestria". One day, Trixie arrives in Ponyville.
There is n boxes on the stage. There is a random key in each box and each box has a key match it. As a unicorn, Trixie can use "Alohomora"(the unlocking spell) to open any box.
Trixie would like to open all the boxes, also Trixie would like to minimize the spells she used. Help Trixie to calculate the expected number of spells she need to used in order to open all the boxes necessary.
![](https://i-blog.csdnimg.cn/blog_migrate/fffaf0ebad9d8e371fcd55995aa490e3.png)
There is n boxes on the stage. There is a random key in each box and each box has a key match it. As a unicorn, Trixie can use "Alohomora"(the unlocking spell) to open any box.
Trixie would like to open all the boxes, also Trixie would like to minimize the spells she used. Help Trixie to calculate the expected number of spells she need to used in order to open all the boxes necessary.
Input
Input contains multiple test cases (less than 100).
For each test case, only one line contains one integer n (1<=n<=1000000000).
For each test case, only one line contains one integer n (1<=n<=1000000000).
Output
For each case, output the corresponding result, rounded to 4 digits after the decimal point.
Sample Input
1 2
Sample Output
1.0000 1.5000HintThe second sample: There are two different permutations when n equal to 2: (1, 2) and (2, 1). Trixie need to use 2 spells in (1, 2) and 1 spell in (2, 1).
思路:写出前几个样例可以发现
f(1)=1.0000
f(2)=1.5000
f(3)=1.8333
f(4)=2.0833
发现了什么 ? 没错f(n)=f(n-1)+1/n
那么就是f(n)=1/1+1/2+1/3+...1/n O(n)递推过去肯定TLE
所以对这个数列求极限,可以得到f(n)=ln(n)+C(C是欧拉常数,等于0.577215664901,当然这道题要求精确到小数点后四位,所以取0.5772即可)
注意前面的项不能用这个极限,这个是n很大才成立的~
代码:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define N 100010
#define C 0.577215664901
double f[N];
int main()
{
long long n;
f[0]=0;
for(int i=1;i<=100000;i++)
{
double b=i;
f[i]=f[i-1]+1/b;
}
while(~scanf("%lld",&n))
{
if(n>100000)
printf("%.4lf\n",log((double)n)+C);
else printf("%.4lf\n",f[n]);
}
return 0;
}