给你一个1~n的序列 让你求出两个完全不同的集合a,b 使得gcd(ai,bi) > 1 且 尽可能的大。
倒序枚举是:因为大一些的质数(x)只能与比它大的匹配(z),而比它大的质数(z)可能会被比它小的质数(y) 匹配,导致x无匹配项。
欧拉筛筛质数。
按上述方式求解就可
代码:
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define CLOCK printf("Time used:%.2lf\n",(double)clock()/CLOCKS_PER_SEC);
#include<bits/stdc++.h>
#define int long long
#define sc scanf
#define pf printf
using namespace std;
typedef pair<int,int> pii;
typedef long long ll;
const int INF = 0x3f3f3f3f;
const double eps = 1e-5;
const int mod = 998244353;
const int N = 200010;
bool book[N];
int is_prime[N],primes[N],cnt;
vector<int> ve,s;
vector<pii> ans;
void ols(int n){
for(int i=2;i<=n;i++){
if(!is_prime[i]) primes[cnt++] = i;
for(int j=0;j<cnt&&primes[j]*i<=n;j++){
is_prime[primes[j]*i] = 1;
if(i % primes[j] == 0) break;
}
}
}
signed main(){
IOS;
#ifdef ddgo
freopen("C:/Users/asus/Desktop/ddgoin.txt","r",stdin);
#endif
ols(200000);
int tt; cin>>tt;
while(tt --){
memset(book,false,sizeof(book));
s.clear(); ans.clear();
int n; cin>>n;
int pos = upper_bound(primes,primes+cnt,n/2)-primes-1;//找到第一个小于n/2的质数
for(int i=pos;i>=0;i--){
ve.clear();
int m = primes[i];
for(int j=1;j*m<=n;j++){
if(!book[j*m]){
ve.push_back(j*m);
book[j*m] = true;
}
}
if((int)ve.size() & 1){
s.push_back(ve[1]);
ans.push_back({ve[0],ve[2]});
for(int j=4;j<(int)ve.size();j+=2) ans.push_back({ve[j-1],ve[j]});
}else{
for(int j=1;j<(int)ve.size();j+=2) ans.push_back({ve[j-1],ve[j]});
}
}
for(int j=1;j<(int)s.size();j+=2) ans.push_back({s[j-1],s[j]});
cout<<(int)ans.size()<<endl;
for(auto i : ans) cout<<i.first<<" "<<i.second<<endl;
}
return 0;
}