Lưu ý: Các code mẫu dưới đây chỉ mang tính tham khảo và có thể không AC được bài tập này

Code mẫu của ladpro98

#include <bits/stdc++.h>
#define long long long

using namespace std;
map<long, int> F;
int m, n, k;
long res;

void Factor(long x) {
    int lim = int(sqrt(x)); int cnt;
    for(int i = 1; i <= lim; i++)
    if (x % i == 0) {
        if (F[i] < k) {
            cnt = ++F[i];
            if (cnt >= k) res = max(res, (long)i);
        }
        if ((long)i * i != x) {
            if (F[x / i] < k) {
                cnt = ++F[x / i];
                if (cnt >= k) res = max(res, x / i);
            }
        }
    }
}

int main()
{
    long ai;
    scanf("%d %d", &n, &k);
    for(int i = 1; i <= n; i++) {
        scanf("%lld", &ai);
        Factor(ai);
    }
    printf("%lld", res);
    return 0;
}

Code mẫu của skyvn97

#include<bits/stdc++.h>
#define MAX   100100
#define FOR(i,a,b) for (int i=(a);i<=(b);i=i+1)
#define REP(i,n) for (int i=0;i<(n);i=i+1)
#define FORE(i,v) for (__typeof((v).begin()) i=(v).begin();i!=(v).end();i++)
#define fi   first
#define se   second
using namespace std;
typedef long long ll;
typedef pair<ll,int> ii;
bool notprime[MAX];
vector<int> prime;
int n,k;
ll res,a[MAX];
vector<ii> ana[MAX];
set<ll> check;
void eratosthene(void) {
    REP(i,2) notprime[i]=true;
    FOR(i,2,MAX-1) if (!notprime[i]) {
        prime.push_back(i);
        for (int j=2*i;j<MAX;j=j+i) notprime[j]=true;
    }
    //cerr<<prime.size()<<endl;
}
void init(void) {
    scanf("%d%d",&n,&k);
    FOR(i,1,n) cin>>a[i];
    res=1;
}
void analysis(ll x,vector<ii> &ana) {
    FORE(it,prime) {
        if (*it>x) break;
        int t=0;
        while (x%*it==0) {
            t++;
            x/=*it;
        }
        if (t>0) ana.push_back(ii(*it,t));
    }
    if (x>1) ana.push_back(ii(x,1));
    //FORE(it,ana) printf("%d %d\n",it->fi,it->se);
}
void update(ll x) {
    if (x<=res) return;
    if (check.find(x)!=check.end()) return;
    check.insert(x);
    int t=0;
    FOR(i,1,n) if (a[i]%x==0) {
        t++;
        if (t>=k) {
            res=x;
            return;
        }
    }
}
void backtrack(const vector<ii> &ana,int id,ll cur) {
    REP(i,ana[id].se+1) {
        if (i>0) cur*=ana[id].fi;
        if (id==ana.size()-1) update(cur);
        else backtrack(ana,id+1,cur);
    }
}
void process(void) {
    FOR(i,1,n) analysis(a[i],ana[i]);
    FOR(i,1,n) backtrack(ana[i],0,1);
    cout<<res;
}
int main(void) {
    eratosthene();
    init();
    process();
    return 0;
}