Harder Gcd Problem

题目链接

题意:求1-n里有多少对数满足GCD(x,y)>1

首先我们来讨论一下什么样的数在1-n中不符合题目要求，首先1是不符合的，如果n足够大的话那么对于偶数来说毋庸置疑是满足的，然后是奇数，小于n/2的奇数也是可满足的，大于n/2的奇数如果存在某奇数*一奇数等于此奇数，那么也是可满足的，否则是不可满足的。其实上面都是废话
我们不难发现如果某个数不符合要求，则该数是大于n/2的数，且除了1和它本身没有其他的因数，即素数，我们知道任何一个合数都可以由多个素数相乘得到，也就是说，任何一个合数都可以找到一个不与它互质的素数，两者的因数即为找到的素数，所以我们要找的数对满足的格式是：素数，该素数的倍数。回归上文，不符合题目要求的数有1和大于n/2的素数，所以我们先要找到满足小于或等于n/2的最大素数。

本题多组样例，所以不能每一次都求解满足的素数，所以可用先预处理求得这题范围下的所有素数，可以采用欧拉筛完成，时间复杂度为O(n)。

欧拉筛：

bool is_prime[N];
int prime[N];
int cnt;

void get_prime(int maxm) {
	is_prime[1] = is_prime[0] = true;
	for (int i = 2; i <= maxm; i++) {
		if (!is_prime[i]) {
			prime[cnt++] = i;
		}
		for (int j = 0; j < cnt && i * prime[j] < maxm; j++) {
			is_prime[i * prime[j]] = true;
			if (i % prime[j]==0) break;
		}
	}
}

然后我们又遇到一个问题，我们应该先从较小的素数开始寻找数对还是从较大的素数开始呢？其实带入一下实况就很容易确定了。
比如说n=16
首先，小于或等于n/2的素数有2，3，5，7
从小到大:
(2,4),(6,8),(10,12),(14,16),(3,9)
从大到小
(7,14),(5,15),(3,6),(9,12),(2,4),(8,10);

从小到大开始最终有5对，而从大到小开始有6对。

同时还要讨论一下某素数的倍数有奇数个或偶数个小于或等于n的情况，当有偶数个符合显然随便匹配都可以，但当个数为奇数时却存在较优的选择，如果个数为1，则一个数对也构不成，当个数为大于1的奇数时，因为除2外的素数都是奇数，所以他们的倍数的奇偶规律为:当倍数为奇数倍是为奇数，当倍数为偶数倍时是偶数，而偶数是素数2的倍数，所以当我们留下一个偶数的价值大于留下奇数的价值，这也侧面说明了为什么要从较大的素数开始寻找数对。
完(๑′ᴗ‵๑)

AC代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>

using namespace std;

const int N = 2e5 + 50;

struct node {
	int x, y;
}date[N];

vector<node> vec;

bool is_prime[N];
int prime[N];
int vis[N];
int arr[N];
int cnt;

void get_prime(int maxm) {
	is_prime[1] = is_prime[0] = true;
	for (int i = 2; i <= maxm; i++) {
		if (!is_prime[i]) {
			prime[cnt++] = i;
		}
		for (int j = 0; j < cnt && i * prime[j] < maxm; j++) {
			is_prime[i * prime[j]] = true;
			if (i % prime[j]==0) break;
		}
	}
}

int main() {
	int t;
	scanf("%d", &t);
	get_prime(N);
	while (t--) {
		int n;
		vec.clear();
		memset(vis, 0, sizeof(vis));
		scanf("%d", &n);
		int start=0;
		for (int i = 0; i < cnt; i++) {
			if (prime[i] > n / 2) {
				start = i - 1;
				break;
			}
		}
		for (int i = start; i >= 0; i--) {
			int pos = 0;
			for (int j = 1; j; j++) {
				int temp = j * prime[i];
				if (temp <= n && !vis[temp]) {
					arr[++pos] = temp;
				}
				else if(temp>n) {
					break;
				}
			}
			if (pos % 2 == 0) {
				for (int k = 1; k < pos; k += 2) {
					node temp;
					temp.x = arr[k];
					temp.y = arr[k + 1];
					vis[temp.x] = 1;
					vis[temp.y] = 1;
					vec.push_back(temp);
				}
			}
			else {
				if (pos > 1) {//留下偶数
					node temp;
					temp.x = arr[1];
					temp.y = arr[pos];
					vis[temp.x] = 1;
					vis[temp.y] = 1;
					vec.push_back(temp);
				}
				for (int k = 3; k < pos; k += 2) {
					node temp;
					temp.x = arr[k];
					temp.y = arr[k + 1];
					vis[temp.x] = 1;
					vis[temp.y] = 1;
					vec.push_back(temp);
				}
			}
		}
		printf("%d\n", vec.size());
		for (int i = 0; i < vec.size(); i++) {
			printf("%d %d\n", vec[i].x, vec[i].y);
		}
	}
}