Harder Gcd Problem
解题思路
直接给个大佬的链接吧。
大佬的博客
代码
/*
从大到小安排质数
这个构造真的妙
*/
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mx=2e5+5;
int prime[mx];
bool noprime[mx],f[mx];
int a[mx];
struct node{
int x,y;
}b[mx];
int cnt=0;
void db(){
for(int i=2;i<mx;i++){
if(!noprime[i]) prime[cnt++]=i;
for(int j=0;j<cnt;j++){
if(i*prime[j]>=mx) break;
noprime[i*prime[j]]=1;
if(i%prime[j]==0){
break;
}
}
}
prime[cnt++]=1000000;
return ;
}
int main(){
ios::sync_with_stdio(0);
db();
int t,n,ans;cin>>t;
while(t--){
cin>>n;
ans=0;
int op=upper_bound(prime,prime+cnt,n/2)-prime-1;
for(int i=op;i>=0;i--){
int sum=0;
for(int j=prime[i];j<=n;j+=prime[i]){
if(j==prime[i]*2) continue;
if(j<=n&&f[j]==0){
a[sum++]=j;
f[j]=1;
}
}
if(sum&1) a[sum++]=prime[i]*2,f[prime[i]*2]=1;
for(int j=0;j<sum;j+=2){
b[ans].x=a[j];b[ans++].y=a[j+1];
}
}
cout<<ans<<"\n";
for(int i=0;i<ans;i++){
cout<<b[i].x<<" "<<b[i].y<<"\n";
}
for(int i=1;i<=n;i++) f[i]=0;
}
return 0;
}