Harmonic Number LightOJ - 1234

题目：Harmonic Number   LightOJ - 1234

题意： 求   Hn

思路：题意十分的简单，但是看到这道题第一想法就是暴力，后来发现有T次查询，如果暴力指定超时，所以就想到了打表，打表的时间复杂度是O(n),查询时O(1)，但是问题来了数据范围最大到1e8,内存会爆掉的，所以就需要优化一下节约内存，所以可以采用分段打表，100数打一个表，空间是1e6那么内存就不会。

代码：

#include<bits/stdc++.h>
using namespace std;
const int maxn  = 1e6+5;
const int INF = 1e8+5;
double a[maxn];
void init(){//打表
    double tmp = 0;
    for(int i = 1; i <= INF; i++){
        tmp += double(1.0/i);
        if(i % 100  == 0){
            a[i/100] = tmp;
        }
    }
}
int main(){
    int T;
    init();
    cin>>T;
    int ca = 0;
    while(T--){
        int n;
        scanf("%d", &n);
        int t = n / 100;//算出n在那一块
        double ans = a[t];
        for(int i = t * 100+1; i <= n; i++){//剩下的部分
            ans += 1.0/(i*1.0);
        }
        printf("Case %d: %.10lf\n",++ca, ans);
    }
return 0;
}

还有一种方法就是利用公式的方法， 参考博客：Harmonic Number(调和级数+欧拉常数)

公式：f(n)=ln(n)+C+1/(2*n)， 公式中C≈0.57721566490153286060651209;但是使用的时候n不能过小， 10000以上的数就可以用这个公式，剩下的就是打表，10000以内的数打表应该非常容易。

代码：

#include<bits/stdc++.h>
using namespace std;
const double  C = 0.57721566490153286060651209;
const int maxn = 100000;
double a[maxn];
void Init(){
     a[1] = 1;
     for(int i = 2; i <= maxn; i++){
        a[i] = a[i-1]+double(1.0/i);
     }
}
int main(){
    int T;
    Init();
    cin>>T;
    int ca = 0;
    while(T--){
        int n;
        scanf("%d", &n);
        if(n <= maxn)printf("Case %d: %.10lf\n",++ca, a[n]);
        else {
            printf("Case %d: %.10lf\n",++ca, log(n)+C+double(1.0/(2*n)));
        }
    }

return 0;
}