LightOJ1234 Harmonic Number 调和级数求和

Posted czsharecode

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了LightOJ1234 Harmonic Number 调和级数求和相关的知识,希望对你有一定的参考价值。

【题目】

技术分享图片

 

【预备知识】

技术分享图片,其中r是欧拉常数,const double r= 0.57721566490153286060651209;

这个等式在n很大 的时候 比较精确。

 

【解法】可以在 n较小的时候,比如n<1e6时,直接用预处理的打表O(1)求值,在n比较 大的时候,运用以上公式,此时要减去 1/(2*n)加以修正。

#include<iostream>
#include<cmath>
using namespace std;
const double euler= 0.57721566490153286060651209;
const int maxn = 1e6;
double a[maxn];
int cas = 1;
int main(){
    long long n;
    a[1] = 1;
    for(int i=2; i<maxn; i++){
        a[i] = a[i-1] + 1.0 / i;
    }
    int t;
    cin>>t;
    while(t--){
        cin>>n;
        if(n < maxn){
            printf("Case %d: %.10lf
",cas++,a[n]);
            continue;
        }
        double ans = log(1+n) + euler - 1.0/(2*n);
        printf("Case %d: %.10lf
",cas++,ans);
    }
    return 0;
} 

 

 

【分块打表】

虽然1e8的表打不出来,但1e6的表很好打,所以每隔100个数记录一次前缀和。到时用的时候,O(1)取出最接近n的前缀和,余下不足100个数暴力 求和即可。

 

#include<iostream>
#include<cmath>
using namespace std;
const double euler= 0.57721566490153286060651209;
const int maxn = 1e8+100;

double a[maxn/100];
int count = 1;

int cas = 1;
int main(){
    long long n;
    a[0] = 0;
    double s = 0;
    for(int i=1; i<maxn; i++){
        s += 1.0/i;
        if( i % 100 == 0){
            a[count++] = s;
        }
    }
    int t;
    cin>>t;
    while(t--){
        double ans = 0;
        cin>>n;
        int num = n / 100;//对应a[num]
        ans += a[num];
         for(long long i=num * 100 + 1; i<=n; i++){
             ans += 1.0/i;
         } 
         printf("Case %d: %.10lf
", cas++, ans);
        
    }
    return 0;
} 

 

以上是关于LightOJ1234 Harmonic Number 调和级数求和的主要内容,如果未能解决你的问题,请参考以下文章

LightOJ 1234Harmonic Number

Harmonic Number LightOJ - 1234 ??????????????????

LightOJ1234 Harmonic Number 调和级数求和

I - Harmonic Number LightOJ - 1234 (分段打表+暴力)

Light oj 1234 - Harmonic Number

1245 - Harmonic Number (II)---LightOJ1245