Light oj 1234 - Harmonic Number

shiyuankongbu

于 2013-05-20 15:56:44 发布

阅读量1.8k

点赞数

分类专栏： .

本文链接：https://blog.csdn.net/shiyuankongbu/article/details/8950862

版权

. 专栏收录该内容

21 篇文章 0 订阅

订阅专栏

本文介绍了Lightoj1234-HarmonicNumber问题的解决思路，通过使用欧拉-马斯刻若尼常数和特定公式来计算调和数。文中还探讨了当输入数值较大时的处理方法。

摘要由CSDN通过智能技术生成

题目：Light oj 1234 - Harmonic Number

思路：。。。。

没思路。。。目测应该有公式。。。forum里面的预处理的方法，目测不大会。。。。

所以公式就公式吧

Euler–Mascheroni constant

定义 gamma = 0.57721566490153286060651209008240243104215933593992

$\gamma = \lim_{n \rightarrow \infty } \left(\sum_{k=1}^n \frac{1}{k} - \ln(n) \right)=\int_1^\infty\left({1\over\lfloor x\rfloor}-{1\over x}\right)\,dx.$

详细介绍：here

开始看公式的时候不仔细，看到n很大的时候取 log(n)+gamma 就可以搞定了，结果WA了，目测应该是刚好第八位不够精确 = =

然后仔细的往下看了看，= = 好多无聊的公式啊

then 我试了试第三个公式，过了....

可是这个题，有木有非公式的做法呢？

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <iostream>
using namespace std;
#define maxn 1000010
#define gamma 0.57721566490153286060651209008240243104215933593992
double sum[maxn];
void init()
{
    sum[0]=0;
    for(int i=1;i<maxn;i++)
        sum[i]=sum[i-1]+1.0/i;
}
int main()
{
    init();
    int t;
    scanf("%d",&t);
    for(int cases=1;cases<=t;cases++)
    {
        int n;
        scanf("%d",&n);
        printf("Case %d: ",cases);
        if(n<maxn)
            printf("%.8lf\n",sum[n]);
        else
            printf("%.8lf\n",(log(n)+log(n+1))/2+gamma);
    }
    return 0;
}