SGU117

题目大意：

给定n个数，求n个数中有多少个数它的m次方能被k整除

0<n,m,k<=10000

看到这种题就应该想到分解质因数……

把k分解质因数，然后把那n个数都分解质因数

乘方就变成指数的乘法

如果一个数每一项的指数>=k对应的项的指数，那么这个数的m次方可被k整除

快速幂应该也是可以过的……

//Lib #include<cstdio> #include<cstring> #include<cstdlib> #include<cmath> #include<ctime> #include<iostream> #include<algorithm> #include<vector> #include<string> #include<queue> using namespace std; //Macro #define rep(i,a,b) for(int i=a,tt=b;i<=tt;++i) #define rrep(i,a,b) for(int i=a,tt=b;i>=tt;--i) #define erep(i,e,x) for(int i=x;i;i=e[i].next) #define irep(i,x) for(__typedef(x.begin()) i=x.begin();i!=x.end();i++) #define read() (strtol(ipos,&ipos,10)) #define sqr(x) ((x)*(x)) #define pb push_back #define PS system("pause"); typedef long long ll; typedef pair<int,int> pii; const int oo=~0U>>1; const double inf=1e20; const double eps=1e-6; string name="",in=".in",out=".out"; //Var int n,m,k; int ori[10008],d[10008],prime[10008],f,g,x,cnt,ans; int vis[10008]; void getprime(int n) { rep(i,2,n) { if(!vis[i])prime[++cnt]=i,vis[i]=i; for(int j=1;j<=cnt&&prime[j]*i<=n;j++) { vis[i*prime[j]]=prime[j]; if(i%prime[j]==0)break; } } } void divide(int x,int p[]) { rep(i,1,cnt) { while(x%prime[i]==0)x/=prime[i],p[i]++; if(x==1)break; } } void Work() { scanf("%d%d%d",&n,&m,&k); divide(k,ori); rep(i,1,n) { scanf("%d",&x); memset(d,0,sizeof d); divide(x,d);bool flag=true; rep(i,1,cnt) { d[i]*=m; if(d[i]<ori[i]){flag=false;break;} } if(flag)ans++; } cout<<ans<<endl; } int main() { // freopen((name+in).c_str(),"r",stdin); // freopen((name+out).c_str(),"w",stdout); getprime(10000); Work(); return 0; }