求约数之和

先看题

这个题可以欧拉筛再分解

AC代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
ll t,n,sumn[5000008],cnt;
ll prime[5000008],flag[5000008];
void olas(ll n){//欧拉筛
	flag[1]=1;
	for(int i=2;i<=n;i++){
		if(flag[i]==0) prime[++cnt]=i;
		for(int j=1;j<=cnt&&prime[j]*i<=n;j++){
			flag[prime[j]*i]=1;
			if(i%prime[j]==0) break; 
		}
	}
}
void fun(ll n){//分解过程
	for(int i=1;i<=cnt&&prime[i]<=n;i++){
		while(n%prime[i]==0){
			sumn[i]++;
			n/=prime[i]; 
		}
	}
	if(n!=1){
		sumn[++cnt]++;
		prime[cnt]=n;
	}
}
int main(){
	olas(5000008);
	cin>>n;
	for(int i=1;i<=n;i++){
		cin>>t;
		fun(t);
	}
	ll ans=1;
	for(int i=1;i<=cnt;i++){
		if(sumn!=0){
			ll sum=1,fac=1;
			for(int j=1;j<=sumn[i];j++){
				fac*=prime[i]%1000000007;
				sum+=fac%1000000007;
			}
			ans*=sum%1000000007;
		}
	}
	cout<<ans;
	return 0;
}