Description
在比赛的时候,1Y(1 次AC)是很值得高兴的事情。但很多大牛总会因为很弱智的错误先WA 一次,再AC。
而很多时候,这点罚时的差距使得他们与金牌无缘。弱智错误系列中最显著的就是忘记加龙龙。
龙龙外国人叫它作long long,表示64位整数,输入与输出64位整数则可以使用例如
scanf("%lld", &a)与printf("%lld", a)的形式完成。很多图论和动态规划的题目中,
虽然题目说最后输出的答案是32 位的整数,但中间计算的过程有时会超过int,这时我们就要使用龙龙了。
可惜的是,很多同学刚开始学写程序都是用VC的,在VC上是无法使用long long的,我们要用__int64
代替,输入与输出64位整数则可以使用例如scanf("%I64d", &a)与printf("%I64d", a)的形式完成。
但是提交上OJ 的时候,如果使用GCC或G++,都只支持long long,我们在提交时又得按照上边的改回来(的确挺麻烦,窘)。
为了让知道龙龙的同学们记得使用龙龙,不知道的学习使用龙龙,下边有个很简单的函数,希望大家
求出它的返回值:
long long H(int n){
long long res = 0;
int i;
for(i = 1; i <= n; i=i+1 ){
res = (res + n/i);
}
return res;
}
不过直接使用这个函数是会超时的,必须改造这个函数,当然这一定难不到未来的编程高手–你
输入格式
第一行是数字T(T<=1021)表示case数,接下来T 行,每行一个整数n,n是一个32 位整数(保证可以由int 表示)。
输出格式
函数返回值。
输入样例
2
5
10
输出样例
10
27
解题思路:
这道题很容易让人进入一个误区,认为这是数学问题,可以推出一条公式。
明显这里所有的变量都是整型的,就是说10/3=3而不是3.3333,所以答案应该比公式算出来要小。
那么我们就发现了,不能整除的数,得到的商是跟最近一个能整除的商一样的,
例如:100/100=1,100/99=1,100/98=1……100/51=1,100/50=2,那么一共有(100-50)个1相加。
……100/34=2,100/33=3,接着就是(50-33)个2相加。
那么怎么快速得到100,50,33……呢?显然就是100/1=100,100/2=50,100/3=33……
设n/i=a,n/(i+1)=b,就有(b-a)个i相加,当然别忘了加上a,a也是结果中的一部分,
那么当i>a的时候就处理完了,但是如果最后i=a,也就是一个完全平方数,如果按照原来的方式计算,
加上a,再加上(b-a)个i,就重复计算了,所以,完全平方数的最后一步需要单独处理。
最后时间复杂度由原来的O(n)优化成了O(√n)。
当n为25时,计算数据如下:
①.a=25,i=1 25/1 = 25
②.a=12,i=2 25/2 = 12
③.a=8,i=3 25/3 = 8*
④.a=6,i=4 25/4 = 6
⑤.a=5,i=5 25/5 = 5
①.25/6 = 4*(7-6) 1个4
②.25/7 25/8 = 3*(9-7) 2个3
③.25/9 25/10 25/11 25/12 = 2*(13-9) 4个2
④.25/13 25/14 25/15 25/16 25/17 25/18 25/19 25/20
25/21 25/22 25/23 25/24 25/25 = 1*(25-13) 12个1
代码:
#include <stdio.h>
long long H_2(int n)
{
long long ret = 0;
int a , b, i;
for (i = 1, a = n; a > i; i++)
{
ret+= a;//当前商
b = n/(i+1);//下一个商
ret += (a-b) * i;//累加和
a = b;//商移动
}
if (a == i)//加上中间一项
ret += a;
return ret;
}
int main()
{
int T, i, n;
scanf("%d", &T);
for (i = 0; i < T; i++)
{
scanf("%d", &n);
printf("%lld\n", H_2(n));
}
return 0;
}