[nowcoder5669H]Harder Gcd Problem

题目相当于问1-n中最多能选出多少对不互素无交集的二元组，并要求方案

构造：将所有数放入其最小质因子对应的集合，若素数p所对应的集合元素个数为奇数且$p e 2$且$2ple n$，那么就将$2p$从2对应的集合移到p对应的集合，最终每一个集合中选择$frac{|S|}{2}$（下取整）对即可

关于这一构造的正确性证明：

1.合法性：素数p所对应的集合中的每一个元素都有公素因子p，因此两两不互素；同时每一个数最多只存在于一个集合，因此不会重复

2.最优性：首先对于$2p>n$的素数p和1都无法参与答案，那么剩余元素的最小质因子一定在$[2,frac{n}{2}]$中，而通过构造中的调整一定可以使得$(2,frac{n}{2}]$中的集合元素个数都为偶数，因此最多剩下1个数，不存在更多的匹配

#include<bits/stdc++.h>
using namespace std;
#define N 200005
vector<int>v,ans[N];
int E,t,n,p[N],mn[N],vis[N],bl[N];
int main(){
    for(int i=2;i<N-4;i++){
        if (!vis[i]){
            p[++p[0]]=i;
            mn[i]=i;
        }
        for(int j=1;(j<=p[0])&&(i*p[j]<N-4);j++){
            vis[i*p[j]]=1;
            mn[i*p[j]]=p[j];
            if (i%p[j]==0)break;
        }
    }
    scanf("%d",&t);
    while (t--){
        scanf("%d",&n);
        for(int i=1;i<=n;i++)vis[i]=0;
        for(int i=2;i<=n;i++){
            vis[mn[i]]++;
            bl[i]=mn[i];
        }
        for(int i=1;(i<=p[0])&&(p[i]<=n);i++)
            if ((vis[p[i]]&1)&&(p[i]<=n/2))bl[2*p[i]]=p[i];
        for(int i=1;i<=n;i++)ans[i].clear();
        for(int i=2;i<=n;i++)ans[bl[i]].push_back(i);
        int sum=0;
        for(int i=2;i<=n;i++)sum+=ans[i].size()/2;
        printf("%d
",sum);
        for(int i=1;(i<=p[0])&&(p[i]<=n);i++)
            for(int j=0;j+1<ans[p[i]].size();j+=2)printf("%d %d
",ans[p[i]][j],ans[p[i]][j+1]);
    }
}