ZCMU-1769-约数问题

1769: 约数问题

Time Limit: 1 Sec   Memory Limit: 128 MB
Submit: 27   Solved: 15
[ Submit][ Status][ Web Board]

Description

Input

多行输入，每行一个整数 N（0 < N < 1000000）。

Output

每行输出一个数，为整数M，即f(1)到f(N)的累加和。

Sample Input

3

Sample Output

5

【解析】
这道题在网上看了有两种方法，其实有共同之处，一种是这么想的因为他让你算的是f1到fn的约数的总和，比如说
1-10，约数为1的其实就是10/1,约数为2的其实就是10/2,约数为3的其实就是10/3，就这样类推，我们只用枚举1到
n/2就可以了，后面的n/2到n肯定就都是只有一个约数的。代码附上。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
long long facs(int m)
{
   long long Sum=0;
   int i;
   for(i=1;i<=m/2;i++)
   {
       Sum+=m/i;
   }
   Sum+=m-i+1;
   return Sum;
}
int main()
{
    int n;
    long long sum;
    while(~scanf("%d",&n))
    {
        sum=facs(n);
        printf("%lld\n",sum);
    }
    return 0;
}

还有一种就是大佬的高效算法，其实是这样子的，我们要算n的因数，我们画一个函数的双曲线xy=n，双曲线上所有
的整数点的横坐标都是n的因数，双曲线又是关于y=x对称我们只用考虑一半的双曲线的横纵坐标就是它的因子总数了
对于小于n的整数对应的双曲线在xy=n的下面，所以我们只需要求出xy=n的一半面积覆盖的整数点乘2就可以了，我们
从1到根号n进行枚举，每一个横坐标对应的点的个数其实就是(n/i)-i+1,这里减去i是减去在y=x下方的点，比如我们
可以这么理解，拿6举个例子，6/2-2+1为2，表示算上的点是2,3和2,2本来是还应该有1,2，不过因为在算2的时候，y
是小于x的了，这里我们只算一半的面积最后乘2的时候，是因为横纵坐标都要算上去。我们在算对称轴上的点的时候
我们乘2其实已经算了两次，所以我们要减去一次，就是最后的(Sum<<1)-i+1。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
long long facs(int m)
{
   long long Sum=0;
   int i;
   for(i=1;i*i<=m;i++)
   {
       Sum+=m/i-i+1;
   }
   return ((Sum<<1)-i+1);
}
int main()
{
    int n;
    long long sum;
    while(~scanf("%d",&n))
    {
        sum=facs(n);
        printf("%lld\n",sum);
    }
    return 0;
}