Harder Gcd Problem(构造)

Harder Gcd Problem

解题思路

直接给个大佬的链接吧。
大佬的博客

代码

/*
从大到小安排质数
这个构造真的妙 
*/
#include 
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;
}