约数（又称因数）

约数

唯一分解定理：任何一个大于1的数都可以被分解成有限个质数乘积的形式 $N=\prod_{i=1}^{m}p_{i}{c_i} $

$N$的正约数个数为：$(c_1+1)\times (c_2+1)\times ...\times (c_m+1)={\prod }_{i=1}^m(c_i+1)$
$N^M$的正约数个数为：$(M\times c_1+1)\times (M\times c_2+1)\times ...\times (M\times c_m+1)={\prod }_{i=1}^m(M\times c_i+1)$
$N$的所有正约数之和为：$\sigma (n)=(1+p_1+p_1^2+...+p_1^{c_1})\times ...\times (1+p_m+p_m^2+...+p_m^{c_m})$ （每一块乘积是一个等比数列）

随机数据下，约数个数的期望是$O({\ln }_{}n)$。一个数的约数个数上界是$2\sqrt{N}$。

例题：

LightOJ-1336 Sigma Function

题意：求$1\sim n$中，$\sigma (n)$是偶数的数的个数。
思路：（主要要自己去推导分析）根据公式：$\sigma (n)=(1+p_1+p_1^2+...+p_1^{c_1})\times ...\times (1+p_m+p_m^2+...+p_m^{c_m})$，当乘积的每一项都为奇数时，$\sigma (n)$是奇数。当$p=2$时，$c$为任意值；当$p!=2$时，$c$为偶数 时，$\sigma (n)$为奇数。当2及其它质数的指数都取偶数时（即平方数），有$\sqrt{n}$个数；当2的指数取奇数，其他质数的指数取偶数时，有$\sqrt{n/2}$个。

#include <bits/stdc++.h>
using namespace std;

#define int long long
int T,n,cs;

signed main(){
	scanf("%lld",&T);
	while(T--){
		scanf("%lld",&n);
		int ans=n;
		ans-=(int)sqrt(n);
		ans-=(int)sqrt(n/2);
		printf("Case %lld: %lld\n",++cs,ans);
	}
}

求一个数的正约数

复杂度：$O(\sqrt{N})$.

for(int i=1;i*i<=n;++i){
	if(n%i==0){
		printf("%d ",i);
		if(n/i!=i) printf("%d ",n/i);
	}
}

如果能知道一个大数（1e18）的质因子分解情况，可以直接dfs枚举幂次，搞出这个大数的所有因子，效率更快！

例题：E2. Divisible Numbers (hard version)

多次求解一个很大的数的正约数个数

比一个一个$O(\sqrt{N})$去求解要更快。

//先欧拉筛筛出 1~根号n 的质因子 
void prime(int n){
	for(int i=2;i<=n;++i){
		if(!vis[i]){
			vis[i]=1;
			pri[++cnt]=i;
		}
		for(int j=1;pri[j]<=n/i;++j){
			vis[pri[j]*i]=1;
			if(i%pri[j]==0) break;
		}
	}
}

//然后利用定理结论 计算正约数的个数 
//根据筛出来的质数对数n进行质因数分解
int cal(int n){
	int ans=1;
	for(int i=1;i<=cnt;++i){
		if(pri[i]>n) break;
		if(n%pri[i]==0){
			int num=0;
			while(n%pri[i]==0) num++,n/=pri[i];
			ans*=(num+1);
		}
	}
	if(n>1) ans*=2;
	return ans;
} 

例题：
LightOJ-1341 Aladdin and the Flying Carpet
题意：多组输入，每次给定$a,b$，求解$a$的因子数，要求因子大于$b$，且小于等于$\sqrt{(}a)$，若$\sqrt{(}a)\times \sqrt{(}a)==a$，$\sqrt{(}a)$不取。

#include <bits/stdc++.h>
using namespace std;

#define int long long
const int N=1e6+10;
int cnt=0,pri[N];
bool vis[N];

void prime(int n){
	for(int i=2;i<=n;++i){
		if(!vis[i]){
			vis[i]=1;
			pri[++cnt]=i;
		}
		for(int j=1;pri[j]<=n/i;++j){
			vis[pri[j]*i]=1;
			if(i%pri[j]==0) break;
		}
	}
}

int cal(int n){
	int ans=1;
	for(int i=1;i<=cnt;++i){
		if(pri[i]>n) break;
		if(n%pri[i]==0){
			int num=0;
			while(n%pri[i]==0) num++,n/=pri[i];
			ans*=(num+1);
		}
	}
	if(n>1) ans*=2;
	return ans/2;
} 

int T,a,b,cs;

signed main(){
	prime(1000000);
	scanf("%lld",&T);
	while(T--){
		scanf("%lld%lld",&a,&b);
		if(b*b>=a) printf("Case %lld: %lld\n",++cs,0ll);
		else{
			int ans=cal(a);
			for(int i=1;i<b;++i) if(a%i==0) ans--;
			printf("Case %lld: %lld\n",++cs,ans);	
		}
	}
}

求$1\sim n$中每个数的正约数集合

时间复杂度：$O(n{\log }_{}n)$
倍数法的思路是，对于每个数d,它的倍数$2d,3d,4d\dots \dots$等等的约数里就有$d$。所以时间复杂度为 $O(n+\frac{n}{2}+\frac{n}{3}+..+\frac{n}{n})\approx O(n{\log }_{}n)$。
常用与一些于因子有关的计算，如计算 $\sum_{i=1}^{n} \sum_{d|n} d $ 或者 ${\sum }_{i=1}^n{\sum }_{d|n}f(d)$.

#include <bits/stdc++.h>
using namespace std;

int n;
vector<int> factor[100000];

int main(){
	cin>>n;
	for(int i=1;i<=n;++i) {
		for(int j=1;j*i<=n;++j) { //1倍的i，2倍的i，3倍的i，...，都有i
			factor[i*j].push_back(i); //因子从小到大被求出 //写成 factor[i*j].push_back(j)的话，因子从大到小被求出 
		}
	}
    /*另一种写法
    for(int i=1;i<=n;i++){ 
        for(int j=i;j<=n;j+=i){ //1倍的i，2倍的i，3倍的i，...，都有i
            factor[j].push_back(i);
        }
    }
    */
	for(int i=1;i<=n;++i) {
		cout<<i<<": ";
		for(int j=0;j<factor[i].size();++j) cout<<factor[i][j]<<' ';
		cout<<endl;
	}
}

求$1\sim n$中每个数的质因子集合

欧拉筛中求得

vector<int> v[N];
int cnt=0,pri[N];
bool vis[N];

void prime(int n){
	for(int i=2;i<=n;++i){
		if(!vis[i]){
			pri[++cnt]=i;
			v[i].pb(i);
			vis[i]=true;
		}
		for(int j=1;pri[j]<=n/i;++j){
			vis[pri[j]*i]=true;
			v[pri[j]*i]=v[i];
			if(i%pri[j]==0) break;
			v[pri[j]*i].pb(pri[j]);
		}
	}
}

相关例题：
Acwing-97.约数之和
求解$A^B$的所有约数之和，$0<=A,B<=5e7$.
做法：质因数分解$A$，然后运用 约数之和公式+等比数列公式 求解。
有一个细节说明：
在求乘法逆元时，当p[i]-1是MOD的倍数时，即$p[i]$不存在MOD的逆元，即$p[i]-1\equiv 0(\mathrm{mod}MOD)$，即$p[i]\equiv 1(\mathrm{mod}MOD)$，所以 $1+p_i+p_i^2+...+p_i^{c_i}\equiv 1+1+1+...+1\equiv c_i\ast b+1(\mathrm{mod}MOD)$.

#include <bits/stdc++.h>
using namespace std;

#define int long long

const int N=1e2;
const int MOD=9901;
int a,b;
int p[N],c[N],cnt=0;

void divide(int x){
	for(int i=2;i<=sqrt(x);++i){
		if(x%i==0) p[++cnt]=i;
		while(x%i==0){
			x/=i;
			c[cnt]++;
		}
	}
	if(x>1) p[++cnt]=x,c[cnt]=1;
}

int qpow(int a,int b){
	int ans=1;
	while(b){
		if(b&1) ans=ans*a%MOD;
		a=a*a%MOD;
		b>>=1;
	}
	return ans;
}

signed main(){
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);
	cin>>a>>b;
	int ans=1;
	divide(a);
	for(int i=1;i<=cnt;++i){
		if((p[i]-1)%MOD==0){
			ans=(c[i]*b+1)*ans%MOD;
		}
		else{
			ans=(qpow(p[i],c[i]*b+1)-1+MOD)%MOD*qpow(p[i]-1,MOD-2)%MOD*ans%MOD; //指数不能取mod 
		}
	}
	if(a) cout<<ans<<'\n'; //卡0点 
	else cout<<0<<'\n';
}