调和级数求和(打表/欧拉常数) - Harmonic Number LightOJ - 1234

调和级数求和(打表/欧拉常数) - Harmonic Number LightOJ - 1234

题意：

$计算调和级数前n项的和:$

$$H_n=1+\frac{1}{2}+\frac{1}{3}+...+\frac{1}{n}=\sum _{i=1}^n\frac{1}{i}$$

Input

Input starts with an integer T (≤ 10000), denoting the number of test cases.

Each case starts with a line containing an integer n (1 ≤ n ≤ 10⁸).

Output

For each case, print the case number and the nth harmonic number $H_n$. Errors less than 10^-8 will be ignored.

Sample Input

12
1
2
3
4
5
6
7
8
9
90000000
99999999
100000000

Sample Output

Case 1: 1
Case 2: 1.5
Case 3: 1.8333333333
Case 4: 2.0833333333
Case 5: 2.2833333333
Case 6: 2.450
Case 7: 2.5928571429
Case 8: 2.7178571429
Case 9: 2.8289682540
Case 10: 18.8925358988
Case 11: 18.9978964039
Case 12: 18.9978964139

$Timelimit：3000ms，Memorylimit：32768kB$

---

分析：

方法一：打表

$考虑到空间限制，10^8组数据不能够全部打出来，$

$但可以每隔50组数据记录一次，存储2\times 10^6组数据，对于每一个n，再做一个50次以内的递推。$

方法二：欧拉常数

$欧拉常数：是一个主要应用于数论的数学常数。它的定义是调和级数与自然对数的差值的极限。$

$\gamma \approx 0.57721566490153286060651209。$

渐进展开： $$\gamma =H_n-lnn-\frac{1}{2n}+\frac{1}{12n^2}-\frac{1}{120n^4}+...$$

$故当n较大时，我们可近似计算H_n=lnn+\frac{1}{2n}+\gamma$

$由于本题误差在10^8数量级，故当n<10000时，我们递推计算；否则我们直接套公式计算。$

法一代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>

using namespace std;

const int N=2e6+10;

double H[N];
int idx;

void Init()
{
    double tmp=0;
    for(int i=1;i<=1e8;i++)
    {
        tmp+=1.0/i;
        if(i%50==0) H[++idx]=tmp;
    }
}

int main()     
{
    Init();
    
    int T; cin>>T;
    for(int t=1;t<=T;t++)
    {
        int n; scanf("%d",&n);
        int s=n/50;
        double d=0;
        for(int i=s*50+1;i<=n;i++) d+=1.0/i;
        printf("Case %d: %.8lf\n",t,H[s]+d);
    }

    return 0;
}

法二代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>

#define ll long long

using namespace std;

const int N=1e4+10;
const double C=0.57721566490153286060651209;

double H[N];

void Init()
{
    for(int i=1;i<=1e4;i++) H[i]=H[i-1]+1.0/i;
}

int main()     
{
    Init();
    
    int T; cin>>T;
    for(int t=1;t<=T;t++)
    {
        int n; scanf("%d",&n);
        double ans;
        if(n<10000) ans=H[n];
        else ans=log(n)+1.0/(2*n)+C;
        
        printf("Case %d: %.8lf\n",t,ans);
    }

    return 0;
}