50分求助

#include<bits/stdc++.h>
using namespace std;
int gcd(int n,int m){
	if(n%m==0)return m;
	return gcd(m,n%m);
} 
int a[30010],b[30010],my[30010],s[30010];
int main(){
	int n,m1,m2;
	cin>>n>>m1>>m2;
	for(int i=1;i<=n;i++)scanf("%d",&s[i]);
	int M=m1;
	for(int i=2;i<=m2;i++)M*=m1;
	int m=0;
	for(int i=2;i<=m1;i++){
		if(b[i]==0){
			a[++m]=i;
			b[i]=i;
		}
		for(int j=1;j<=m;j++){
			if(a[j]<=b[i]&&i*a[j]<=m1)b[i*a[j]]=a[j];
			else break;
		}
	}
	int sw=0;
	while(m1>1){
		int y=b[m1];
		my[++sw]=y;
		while(m1%y==0)m1/=y;
	} 
	int ans=1e9;
	for(int i=1;i<=n;i++){
		bool flag=true;
		for(int j=1;j<=sw;j++){
			if(s[i]%my[j]!=0){
				flag =false;
				break;
			}
		}
		if(flag){
			long long ys=gcd(M,s[i]);
			int k=0;
			long long s=1;
			while(s%M!=0){
				k++;
				s*=ys;
			}
			ans=min(ans,k);
		}
	}
	if(ans==1e9)cout<<-1<<endl;
	else cout<<ans<<endl;
	return 0;
}