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, 不搞了不搞了。