题目
首先,我们考虑最容易想到的做法,按照题意模拟
代码如下:
#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,可以利用这点继续进行优化,至此,解决问题的代码的时间复杂度已经降低到非常可观的程度