题意:求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);
}
}
}