题解——约数个数和

题目

首先，我们考虑最容易想到的做法，按照题意模拟
代码如下：

#include <iostream>
#include<cmath>
using namespace std;
int f(int a)//求约数个数
{
	int sum = 0;
	for (int i = 1; i <= a; i++)
	{
		if (a % i == 0) { sum++; }
	}
	return sum;
}
int main()
{
	int n;
	cin >> n;
	int sum = 0;
	for (int i = 1; i <= n; i++)
	{
		sum += f(i);
	}
	cout << sum;
	return 0;
}

不难发现，这样子做虽然能到达目的，但是时间复杂度太高，所以我们先对函数f进行第一次优化

int f(int a)//求约数个数
{
	int sum = 0, i;
	for (i = 1; i * i < a; i++)
	{
		if (a % i == 0) 
		{ 
			sum+=2;
		}
	}
	if (i * i == a) { sum++; }
	return sum;
}

如此，就是把类似6%2=0和6%3=0这两种情况合成为一种情况（2*3=6），降低了一定的时间复杂度，但我们还可以考虑更优的解决方法

引入如下定理

约数个数定理

我们可以考虑利用此定理再次对函数f进行优化，代码如下：

int f(int n)
{
    int res = 1;//约数个数
    if (n == 1)
    {
        return 1;
    }
    for (int i = 2; i * i <= n; i++)
    {
        int k = 0;
        while (n % i == 0)
        {
            n = n / i;
            k++;
        }
        if (k)
        {
            res *= (k + 1);
        }
    }
    if (n != 1) //剩一个质因数
    { 
        res = res * 2;//相当于res*=(k+1),k=1
    }
    if (res == 1)//质数情况
    {
        return 2;
    }
    return res;
}

到这里，利用函数算得每个数的约数个数，然后进行求和的方法已经很难进行实质性的优化了
接下来，我们考虑另一种思路
我们不一定要对每个数的约数个数进行考虑，我们可以对约数本身进行考虑，比如，假设n为6，那么1肯定是1，2，3，4，5，6的因数，即1作为因数会在求和过程中出现6次，2是2，4，6的因数，即2作为因数会出现3次，依次类推，我们不难发现，对每个因数，它出现的次数为n/i,i为对应因数，我们可以考虑计算因数出现次数来求约数个数和，代码如下:

#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
int main()
{
	int n, i, s;
	cin >> n;
	s = n;
	for (i = 2; i <= n; i++)
	{
		s = s + n / i;
	}
	cout << s;
	return 0;
}

当i>n/2时，n/i恒为1，可以利用这点继续进行优化，至此，解决问题的代码的时间复杂度已经降低到非常可观的程度