最开始我的想法是把 C 数组的所有质数因子都乘起来在算是否为 K 的倍数,但是后来想了一下直接乘起来可能会数字特别大,而且也是错误的想法,所有质数因子乘起来不一定是 K 的倍数,例如 C 数组的所有质数因子2*3*5,K=4*6*5,最后的结果是 K 为 C 的倍数,但是这样并不能确定 x mod k 的值。
这个代码几乎是看的大神的,但是为了以后可以查看所以按照自己的理解记录下来。统计 C 数组的每个质数因子的最多的个数,如果比K的个数多则可以,少则不行
#include <iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int M=1e6+6;
int c[M],prime[1000],s[M];
bool isprime(int a)
{
for(int i=2;i*i<=a;i++)
{
if(a%i==0) return 0;
}
return 1;
}
int main()
{
int n,k,cnt=0;
for(int i=2;i<=1000;i++)
{
if(isprime(i))
prime[cnt++]=i;
}
scanf("%d%d",&n,&k);
for(int i=0;i<n;i++) scanf("%d",&c[i]);
for(int i=0;i<n;i++)
{
if(c[i]%k==0)
{
puts("Yes");
return 0;
}
for(int j=0;j<cnt&&c[i]>1&&prime[j]*prime[j]<=c[i];j++)
{
int r=0;
while(c[i]%prime[j]==0)
{
r++;
c[i]/=prime[j];
}
s[prime[j]]=max(s[prime[j]],r);
}
if(c[i]>1) s[c[i]]=max(s[c[i]],1);
}
for(int j=0;j<cnt&&k>1&&prime[j]*prime[j]<=k;j++)
{
int r=0;
while(k%prime[j]==0)
{
r++;
k/=prime[j];
}
if(r>s[prime[j]])
{
puts("No");
return 0;
}
}
if(k>1)
if(s[k]>=1) puts("Yes");
else puts("No");
else puts("Yes");
return 0;
}
另外一个
zeppelin_58
#大神的更短的代码
#include <bits/stdc++.h>
using namespace std;
inline int gcd(int a, int b){
while(a > 0 && b > 0){
if(a > b){
a %= b;
}
else{
b %= a;
}
}
return a + b;
}
inline int lcm(int a, int b){
return a / gcd(a, b) * b;
}
int main(){
int n, k, l;
scanf("%d %d", &n, &k);
l = 1;
for(int i = 0; i < n; ++i){
int x;
scanf("%d", &x);
int g = gcd(x, k);
l = lcm(g, l);
}
if(k == l){
cout<<"Yes"<<endl;
}
else{
cout<<"No"<<endl;
}
return 0;
}