洛谷 P1403 [1，n]所有数的约数个数之和

https://www.luogu.org/problemnew/show/P1403

题目描述

科学家们在Samuel星球上的探险得到了丰富的能源储备，这使得空间站中大型计算机“Samuel II”的长时间运算成为了可能。由于在去年一年的辛苦工作取得了不错的成绩，小联被允许用“Samuel II”进行数学研究。

小联最近在研究和约数有关的问题，他统计每个正数N的约数的个数，并以f(N)来表示。例如12的约数有1、2、3、4、6、12。因此f(12)=6。下表给出了一些f(N)的取值：

f(n)表示n的约数个数，现在给出n，要求求出f(1)到f(n)的总和。

输入输出格式

输入格式：

 

输入一行，一个整数n

 

输出格式：

 

输出一个整数，表示总和

 

输入输出样例

输入样例#1： 复制

3

输出样例#1： 复制

5

说明

【数据范围】

20%N<=5000

100%N<=1000000

思路：[1，n]中所有数含有约数i的数的个数有n/i个(向下取整)，因此我们可以得到一个O(n)算法。

#include<iostream>
using namespace std;

int main()
{
    int n,cnt=0;
    cin>>n;
    for(int i=1;i<=n;i++)//[1,n]中的所有数含有约数i的有n/i个(向下取整)
        cnt+=n/i;
    cout<<cnt<<endl;
    return 0;
}

改进：有很多情况下的n/i是一样的，根据这个情况可以再优化，比如n=10的情况：

n=10	i=1	i=2	i=3	i=4	i=5	i=6	i=7	i=8	i=9	i=10
n/i的值	10	5	3	2	2	1	1	1	1	1

我们可以发现：对于所有满足n/i=t的i的个数=n/(n/i)-min(i)+1 不妨带入验证一下，i=2时，结果为1；i=3时，结果为1；i=4时，结果为2；i=6时，结果为5，那么我们可以得到一个更加快速的算法：(时间复杂度大概在O(根号n))

#include<iostream>
using namespace std;

int main() //[1,n]中约数有i的个数为 (int)n/i
{          //很多情况下的 n/i 的值是一样的 这样的个数有n/(n/i)-min(i)+1个
    int n,cnt=0;
    cin>>n;
    int j;
    for(int i=1;i<=n;i=j+1)
    {
        j=n/(n/i);
        cnt+=(n/i)*(j-i+1);
    }
    cout<<cnt<<endl;
    return 0;
}

这道题还有更快的算法， 但是难度太大了orz， 不搞了不搞了。