标签:
/*数论题 考察唯一分解定理 当然用到一些技巧*/ #include<iostream> #include<cstdio> #include<cstring> #define maxn 30010 using namespace std; int n,m1,m2,prime[maxn],num,mi[maxn],S,s[maxn],ans=0x7fffffff,cnt; bool f[maxn]; void Get_prime(int x) { for(int i=2;i<=x;i++) { if(f[i]==0)prime[++num]=i; for(int j=1;j<=num;j++) { if(i*prime[j]>x)break; f[i*prime[j]]=1; if(i%prime[j]==0)break; } } } void Get_mi(int x,int y) { if(x<2)return;//m1==1 的时候会re for(int i=1;prime[i]<=x;i++) { while(x%prime[i]==0)mi[i]++,x/=prime[i]; if(x==1)break; } for(int i=1;i<=m2;i++) mi[i]*=y; } int cla(int a,int b) { if(a<b)return 0; int t=a/b; if(t*b==a)return t; else return t+1; } int main() { scanf("%d%d%d",&n,&m1,&m2); Get_prime(m1); Get_mi(m1,m2); for(int k=1;k<=n;k++) { scanf("%d",&S);int falg=0,A=0; for(int i=1;prime[i]<=m1;i++)//相当于只分解m1 那这m1的因子来分解Si { cnt=0;if(!prime[i])break; while(S%prime[i]==0)cnt++,S/=prime[i]; if(mi[i]&&!cnt) { falg=1;break; } if(mi[i]==0&&cnt==0)continue; A=max(A,cla(mi[i],cnt)); } if(falg==0)ans=min(ans,A); } if(ans>=0x7fffff)printf("-1\n"); else printf("%d\n",ans); return 0; }
标签:
原文地址:http://www.cnblogs.com/yanlifneg/p/5745037.html