G - Harmonic Number (II)

最新推荐文章于 2020-06-25 12:32:18 发布

titi2018815

最新推荐文章于 2020-06-25 12:32:18 发布

阅读量546

点赞数 1

本文链接：https://blog.csdn.net/titi2018815/article/details/82223561

版权

题目：

I was trying to solve problem '1234 - Harmonic Number', I wrote the following code

long long H( int n ) {
  long long res = 0;
  for( int i = 1; i <= n; i++ )
res = res + n / i;
  return res;
}

Yes, my error was that I was using the integer divisions only. However, you are given n, you have to find H(n) as in my code.

Input

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

Each case starts with a line containing an integer n (1 ≤ n < 231).

Output

For each case, print the case number and H(n) calculated by the code.

Sample Input

2147483647

Sample Output

Case 1: 1

Case 2: 3

Case 3: 5

Case 4: 8

Case 5: 10

Case 6: 14

Case 7: 16

Case 8: 20

Case 9: 23

Case 10: 27

Case 11: 46475828386

题意：

让我们求f（n）=n/1+n/2+n/3+......+n/n；同时注意n/i取整；

思路：

首先我们先看数据的范围，n (1 ≤ n < 2 31)，数据范围太大，如果我们按照题目中的代码直接暴力肯定超时，那么，我们就要优化代码；

f=n/x这个函数关于y = x 对称对称点刚好是sqrt(n)，于是就简单了直接求sum+n/i (i*i<n && i >=1)，然后乘以2，再减去i*i即可。

当n==10时，m=sqrt(10)==13，f(10)==（10+5+3）*2-3*3=27；

关于m对称，左右相差m*m，左边10+5+3=18，右边2+2+1+1+1+1+1=9，和左边相差3*3=9，

故，这道题我们就可以简化了；

公式：

m=sqrt（n）；

f（n）=（n/1+n/2+n/3+......+n/k）*2-m*m；（n/k<=m）;

代码如下：

#include<stdio.h>
#include<math.h>

int main()
{
    int t,kk=1;
    scanf("%d",&t);
    while(t--)
    {
        long long n,sum=0;
        scanf("%lld",&n);
        int m=(int)sqrt(n);
        for(int i=1; i<=m; i++)
            sum=sum+n/i;
        sum*=2;
        sum-=m*m;
        printf("Case %d: %lld\n",kk++,sum);
    }
    return 0;
}