参看资料:
https://www.cnblogs.com/acgoto/p/9693207.html
http://www.cnblogs.com/fzl194/p/9027045.html
1》朴素算法求因子数目以及因子和
#include<bits/stdc++.h>
using namespace std;
int get_num(int n){
int num=0;
for(int i=1;i<=n;++i)
if(n%i==0)num++;
return num;
}
int get_sum(int n){
int tot=0;
for(int i=1;i<=n;++i)
if(n%i==0)tot+=i;
return tot;
}
int main(){
int n;
while(cin>>n){
cout<<get_num(n)<<endl;//求n的因子个数
cout<<get_sum(n)<<endl;//求n的因子数之和
}
return 0;
}
2》约数个数定理:O(√n)。
#include<bits/stdc++.h>
using namespace std;
int get_num(int n){
int tot=1;
for(int i=2;i*i<=n;++i){
if(n%i==0){
int x=0;
while(n%i==0){
n/=i;
x++;
}
tot*=(x+1);
}
}
if(n>1)tot*=2;
return tot;
}
int get_sum(int n){
int tot=1;
for(int i=2;i*i<=n;i++){
if(n%i==0){
int mul=1;
while(n%i==0){
n/=i;
mul*=i;
}
tot*=(mul*i-1)/(i-1);
}
}
if(n>1)tot*=(n+1);
return tot;
}
int main(){
int n;
while(cin>>n){
cout<<get_num(n)<<endl;//求n的因子个数
cout<<get_sum(n)<<endl;//求n的因子数之和
}
return 0;
}
3》线性筛法打表因子数目
int divisor_num[maxn];
void init(int n)
{
for(int i = 1; i <= n; i++)
{
divisor_num[i]++;
for(int j = i * 2; j <= n; j += i)
divisor_num[j]++;
}
}
或
int divisor_num[maxn];
void init(int n)
{
for(int i = 1; i * i <= n; i++)
{
for(int j = i; i * j <= n; j++)
{
divisor_num[i * j] += 2;
if(i == j)divisor_num[i * j]--;
}
}
}