约数相关习题

轻拍那头牛

直接求每个数的约数难求，因为涉及质因数分解，然后利用性质，非常麻烦，但是求每个数的倍数好求，所以我们枚举i求出有多少数字是i的倍数即可。

#include <bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define ll long long
#define x first
#define y second
typedef pair<int,int> PII;
const int N=1e6+10;
int n,m,k,ecnt,ans;
int a[N],cnt[N],s[N];


int main() {
	
	cin>>n;
	
	for(int i=1;i<=n;i++){
		cin>>a[i];
		cnt[a[i]]++; 
	} 
	
	for(int i=1;i<=N;i++){
		for(int j=i;j<=N;j+=i){
			s[j]+=cnt[i];//j是i的倍数 
		} 
	} 
	
	for(int i=1;i<=n;i++)cout<<s[a[i]]-1<<endl; //减去自身 
	
	
	
	

	return 0;
}

约数个数定理

N=p1 ^ c1+p2 ^ c2+…+pk ^ ck

约数个数=(c1+1) * (c2+1) * … * (ck+1)
这个很好理解，就是排列组合把他们打开以后有一个1，这个1就是每个数都有的约数1，其余约数是不同情况的组合

樱花

求x，y的不同组合等价于求(n!)^2约数个数

根据阶乘分解质因数

我们得出n!=p1 ^ c1+ p2 ^ c2 +… pk ^ ck 推出
(n!) ^ 2=p1 ^ 2c1 + p2 ^ 2c2+…+ pk ^ 2ck

最后根据约数定理的出ans=(2c1+1) * (2c2+1) …(2ck+1)

#include <bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define ll long long
#define x first
#define y second
typedef pair<int,int> PII;
const int N=1e6+10,mod=1e9+7;
int n,m,k,ecnt,cnt;
int vis[N],prim[N],st[N];

void get_prime(int n) {//线性筛
    memset(vis,0,sizeof vis);
	for(int i=2;i<=n;i++){
		if(!vis[i])prim[++cnt]=i;
		for(int j=1;1ll*i*prim[j]<=n;j++){//越界中断 
			vis[i*prim[j]]=1; 
			if(i%prim[j]==0)break; 
		} 
	} 
}

int main() {
	cin>>n; 
	get_prime(n);
	ll ans=1; 
	for(int i=1;i<=cnt;i++){
		int p=prim[i];
		int s=0;
		for(int j=n;j;j/=p) s+=j/p;//累加和
		ans=(ll)ans*(2*s+1)%mod; 
	} 
	cout<<ans; 

	return 0;
}