2020年牛客暑期多校训练（第四场）harder gcd problem

 题目链接：https://ac.nowcoder.com/acm/contest/5669/H 

输入

2
4
10

输出

1
2 4
4
3 9
5 10
8 2
4 6

题解： 首先先排除1，和大于n/2的质数，因为它们是不能配对的，然后我们考虑将剩下的质数倒序记录下它们未被记录的倍数（因为正序的话，有些数记录不到它们的倍数），然后如果一个数记录的数的数目为偶数（这里包括它们本身），这些数都可以两两匹配，如果为奇数，取出一个偶数即可，剩余的依旧可以配对，最后把取出的偶数两两配对即可。

 

#include<bits/stdc++.h>
using namespace std;
#define ll long long
typedef pair<int,int> P;
const int N = 2e5+5;
int prime[N],tot;
bool st[N];
vector<int> v[N];
void get_primes(int n) {
    for(int i=2;i<=n;i++) {
        if(!st[i]) prime[tot++] = i;
        for(int j=0;prime[j] <= n/i;j++) {
            st[i * prime[j]] = true;
            if(i % prime[j] == 0) break;
        }
    }
}
int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    get_primes(N);
    int t;
    cin >> t;
    while(t--) {
        memset(st,false,sizeof(st));
        vector<int> p;
        vector<P> res;
        int n,cnt = 0;
        cin >> n;
        for(int i=0;i<=n;i++) v[i].clear();
        int m = upper_bound(prime,prime+tot,n/2)-prime-1;
        for(int i=m;i>=0;i--) {
            for(int j=1;j <= n/prime[i];j++) {
                 if(!st[j * prime[i]]) {
                     v[prime[i]].push_back(j * prime[i]);
                     st[j * prime[i]] = true;
                 }
            }
            int len = v[prime[i]].size();
            if(len & 1) {
                 res.push_back({v[prime[i]][0],v[prime[i]][2]});
                 p.push_back(v[prime[i]][1]);
                 for(int j=4;j<len;j+=2) {
                      res.push_back({v[prime[i]][j],v[prime[i]][j-1]});
                 }
            }else {
                for(int j=1;j<len;j+=2) {
                    res.push_back({v[prime[i]][j],v[prime[i]][j-1]});
                }
            }
        }
        for(int i=1;i<p.size();i+=2) {
            res.push_back({p[i],p[i-1]});
        }
        cout << res.size() << endl;
        for(int i=0;i<res.size();i++)
        cout << res[i].first << ' ' << res[i].second << endl;
    }
    return 0;
}