2020牛客多校第4场 H

给你一个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;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值