Harmonic Number (II)

博客探讨了如何解决计算Harmonic Number的问题,当数值范围过大时,暴力计算不可行。文章介绍了两种优化方法:1) 利用反比例函数性质,将和分为两部分并减去重叠部分;2) 分别计算前半段和后半段的和,避免超时。给出了两种方法的AC代码实现,并提供了样例输入输出。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

Harmonic Number (II)

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
11
1
2
3
4
5
6
7
8
9
10
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

题意: 计算n/1+n/2+…+n/n的值

思路: 题所给数据范围太大,故用题中所给方法暴力肯定是不行的,所以要找规律,我看其他大佬的博客学会了两种

1.大佬博客

y=n/x,是一个反比例函数,关于y=x对称,所以我们求n/x的和,相当于求图形的面积,可以分成两份一样的面积(在y=x处分,即一份x从1到sqrt(n),一份从sqrt(n)到n),此时x=sqrt(n),但两份面积有重叠,所以要减去重叠部分面积sqrt(n)*sqrt(n)(可能光看这些文字想不明白,所以可去看上述博客里的图片,一目了然)
2.大佬博客
这第二种方法就是找sqrt(n)到n之间n/x的和,前半段可暴力,不会超时,后半段按规律相加即可,上述博客叙述的很详细,我就不再多赘述了

AC代码:

第一种:

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <algorithm>
using namespace std;
int main()
{
	int t;
	scanf("%d",&t);
	int kk=1;
	while(t--)
	{
		int x;
		scanf("%d",&x);
		long long sum=0;
		int m=sqrt(x);
		for(int i=1; i<=m; i++)//前半段面积
			sum+=x/i;
		sum=sum*2-m*m;//*2-重叠面积
		printf("Case %d: %lld\n",kk++,sum);
	}

第二种:

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
int main()
{
	int t,n,i,j=1,k;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d",&n);
		k=sqrt(n);
		long long sum=0;
		for(i=1;i<=k;i++)//前半段和
			sum+=n/i;
		for(i=1;i<=k;i++)//后半段和
		sum+=(n/i-n/(i+1))*i;
		if(n/k==k)//减算重的
		sum-=k;
		printf("Case %d: %lld\n",j++,sum);
	}
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值