Codeforces Round #579 (Div. 3)--Common Divisors（gcd，约数个数，质因数分解）

题目链接

题意：
给定个数为n的一串数字序列，找出所有数字的公共因子个数。

题目条件：
1<=n<=4*10 ^5，1<=ai<=10 ^12。

分析：
经过仔细思考，不难想到首先应当计算出所有数字的最大公约数，只需对所有数字算一遍gcd即可，然后就想一想，其实就是求所有数字最大公约数的因子个数，如果直接暴力从1~n判断是否整除会超时，那么只需要进行质因数分解即可，这里质数可以先打素数表（ai最大1e12,只要开1e6内存即可），再质因数分解，也可以直接试除法质因数分解，两种方法复杂度相差不大，打好的素数表质因数分解的更快些，但实测试除法用时更短些，可能是第一种方法打素数表时耗时过多。然后再利用约数个数公式(c1+1)(c2+1)(c3+1)…(cm+1)，即可得到最终答案。

素数打表质因数分解主要代码：

int n;
VI primes(1e6+100),v(1e6+100);
int cnt=0;
void Ai_prime(){
	for(int i=2;i<=1e6;i++){
		if(!v[i])
		primes[cnt++]=i;
		else
		continue;
		for(int j=i;j<1e6/i;j++){
			v[i*j]=1;
		}
	}
}
ll gcd(ll a,ll b){
	return b?gcd(b,a%b):a;
}
void solve(){
    vector<ll> a(n);
    REP1(i,n){
    	cin>>a[i];
    }
    sort(ALL(a));
	ll maxgcd;
	if(n>1){
		maxgcd=gcd(a[0],a[1]);
    for(int i=2;i<n;i++){
    	if(maxgcd==1)
    	break;
    	maxgcd=gcd(a[i],maxgcd);
    }
	}else
    maxgcd=a[0];
    ll sum=1;
    for(int i=0;i<cnt&&primes[i]<=sqrt(maxgcd);i++){
    	int temp=0;
    	while(maxgcd%primes[i]==0){
    		temp++;
    		maxgcd/=primes[i];
    	}
    	if(temp)
    	sum*=temp+1;
    }
    if(maxgcd>1)
    sum*=2;
    cout<<sum<<endl;
}
int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    Ai_prime();
//    freopen("debug\\in.text","r",stdin);
//    freopen("debug\\out1.text","w",stdout);
    while(cin>>n)
    solve();
    return 0;
}

试除法质因数分解主要代码：

int n;
ll divide(ll n,VI& p,VI& c){
	int m=0;
	ll sum=1;
	for(int i=2;i<=sqrt(n);i++){
		if(n%i==0){
			p[++m]=i,c[m]=0;
			while(n%i==0)n/=i,c[m]++;
			sum*=c[m]+1;
		}
	}
	if(n>1)
	p[++m]=n,c[m]=1,sum*=2;
	return sum;
}
ll gcd(ll a,ll b){
	return b?gcd(b,a%b):a;
}
void solve(){
    vector<ll> a(n);
    VI p(1e6+100),c(1e6+100);
    REP1(i,n){
    	cin>>a[i];
    }
    sort(ALL(a));
	ll maxgcd;
	if(n>1){
		maxgcd=gcd(a[0],a[1]);
    for(int i=2;i<n;i++){
    	if(maxgcd==1)
    	break;
    	maxgcd=gcd(a[i],maxgcd);
    }
	}else
    maxgcd=a[0];
    ll sum=divide(maxgcd,p,c);
    cout<<sum<<endl;
}
int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
//    freopen("debug\\in.text","r",stdin);
//    freopen("debug\\out1.text","w",stdout);
    while(cin>>n)
    solve();
    return 0;
}