Kuangbin数论专题，Harmonic Number (II)（ LightOJ - 1245，找规律）

kuangbin
题意：题目很简单让你求n/1+n/2+……+n/n;思路：其实这道题和数论没有太多关系，就是一个找规律的题目，一开始我观察这个公式是n(1+1/2+1/3+……+1/n)以为会有一个求和公式，结果它是一个调和级数，不收敛，虽然有一个近似公式约等于欧拉常数+ln（n）但是只能在n趋近于无穷的时候无限逼近，所以只能找其中的规律先求出前sqrt(n)项和：即n/1+n/2+…+n/sqrt(n)再求出后面所以项之和.后面每一项的值小于或者等于sqrt(n),计算值为1到sqrt(n)的项的个数，乘以其项值即可快速得到答案例如:10/1+10/2+10/3+…+10/10sqrt(10) = 3先求出其前三项的和为10/1+10/2+10/3在求出值为1的项的个数为(10/1-10/2)个，分别是(10/10,10/9,10/8,10/7,10/6),值为2个项的个数（10/2-10/3）分别是（10/5,10/4）,在求出值为3即sqrt(10)的项的个数.显然，值为sqrt(10)的项计算了2次,减去一次即可得到答案。当n/(int)sqrt(n) == (int)sqrt(n)时，值为sqrt(n)的值会被计算2次。但是对于n=12，就不需要减,每一个值都被计算了一次
12/1+12/2+12/3+12/4+12/5+12/6+12/7+12/8+12/9+12/10+12/11+12/12
3 2 1
12/3-12/4个 12/2-12/3个 12/1-12/2个
代码：

#include< iostream>
#include< cstdio>
#include< cstdlib>
#include< cstring>
#include< string>
#include< vector>
#include< algorithm>
#include< cmath>
using namespace std;
int main()
{
    int t;
    scanf("%d",&t);
    int C=1;
    while(t--){
        int n;
        scanf("%d",&n);
        long long ans=0;
        int i;
        for(i=1;i<=(int)sqrt((double)n);i++){
            ans+=(n/i);
            if(n/i>n/(i+1)){
                ans+=(((n/i)-(n/(i+1)))*i);
            }
        }
        i--;
        if((n/i)==i){
            ans-=i;
        }
        printf("Case %d: %lld\n",C++,ans);
    }
    return 0;
}