问题 B: 要塞任务(数论）

#include<bits/stdc++.h> 
#include <math.h>
using namespace std;
typedef long long ll; 
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();}
    while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();}
    return x*f;
}
const int maxn=1e5+100;
const int inf=0x3fffffff;
const int mod=1e9+7; 
int biaoji[maxn],p[maxn],ol[maxn],cnt,a[maxn];
void P(){
    for(int i=2;i<=maxn;i++){
        if(!biaoji[i]){
            p[++cnt]=i;
            ol[i]=i-1;
        }
        for(int j=1;j<=cnt&&p[j]*i<=maxn;j++){
            biaoji[p[j]*i]=1;
            if(i%p[j]==0){
                break;
            }
        } 
    }
    ol[1]=1;
    int t;
    for(int i=2;i<=maxn;i++){
        if(biaoji[i]){
            ol[i]=i;
            t=i;
            for(int j=1;j<=cnt&&t>1;j++)
            {
                if(i%p[j]==0)
                {
                    ol[i]=ol[i]*(p[j]-1)/p[j];
                    while(t%p[j]==0){
                        t/=p[j];
                    }
                }
                if(p[j]*2>i){
                    break;
                }
            }
        }
    }
}
int t;
void inint(){
    cin>>t;
}
int main(){
    P();
    int ans=-1,n,m;
    int ma=-1e9,tot=0;
    inint();
    while(t--){
        ans=-1;
        ma=-1e9;
        memset(a,0,sizeof(a));
        cin>>n>>m;
        for(int i=1;i<=n;i++){
            int p;
            p=read();
            a[p]++;
            ma=max(ma,p);
        }
        for(int i=1;i<=ma;i++){
            tot=0;
            for(int j=i;j<=ma;j+=i){
                tot+=a[j];
            }
            if(tot>=m){
                ans=max(ans,ol[i]);
            }
        }
        printf("%d\n",ans);
    }
}