题意简述
给出正整数n,你要把1-n之间的正整数分成尽可能多组,使得每一组两个数的最大公约数大于1;输出能分成最多组的个数,并按任意顺序输出每组的两个数.
很妙的一道题。
首先我们考虑去处理每个质数的倍数,那么这个质数一定是小于等于 n/2 的,不然它在 n 的范围内是不会有倍数的。
那么我们先把 $1~n/2$ 范围内的所有质数筛出来。
枚举质数
然后我们要怎么用质数去处理答案呢?
首先我们从大到小枚举这些质数,然后去枚举它的倍数。
然鹅这样复杂度不会炸么?不会。 $O(\sigma_{i=2}^{n/2}n/i)$ 的复杂度吧,而且还远达不到。
信仰一下应该是过得了的,实际上确实也过了。
处理倍数
然后我们再考虑如果这些倍数该怎么处理。
我们首先看这些倍数有没有被使用过,没有的话就入栈。
最优性理解
那么为什么能用就入栈的解法最优呢?
其实很简单,因为当前枚举那个数可以被当前的质数整除,那么该数与当前枚举质数的倍数去匹配肯定不会产生更劣的结果。
那么如果有数字剩下来呢?
如果说当前枚举质数的未使用过的倍数有奇数个,那么我们可以将第 2 个倍数与最后一个换一下,然后弹出栈,标记为未使用。
为什么这样做能达到最优?
考虑第二个倍数一定是 2 的倍数,那么把他弹出的话就可以和 2 的倍数去匹配,而且我们每次遇到奇数个的情况都是弹 2 的倍数,那么这些多余的数就可以凑一块儿了。
1 //by Judge
2 #include<vector>
3 #include<cstdio>
4 #include<iostream>
5 #define P make_pair
6 using namespace std;
7 const int M=1e5+11;
8 char sr[1<<21],z[20];int C=-1,Z;
9 inline void Ot(){fwrite(sr,1,C+1,stdout),C=-1;}
10 inline void print(int x,char chr='\n'){
11 if(C>1<<20)Ot();if(x<0)sr[++C]=45,x=-x;
12 while(z[++Z]=x%10+48,x/=10);
13 while(sr[++C]=z[Z],--Z);sr[++C]=chr;
14 } int n,cnt,ans,top,prim[M],is[M],vis[M],num[M];
15 pair<int,int> tmp[M];
16 inline void prep(){
17 for(int i=2;i<=n>>1;++i){
18 if(!is[i]) prim[++cnt]=i;
19 for(int j=1;j<=cnt&&i*prim[j]<=n/2;++j){
20 is[i*prim[j]]|=1;
21 if(i%prim[j]==0) break;
22 }
23 }
24 }
25 int main(){
26 scanf("%d",&n),prep();
27 for(int i=cnt;i;--i){ top=0;
28 for(int j=prim[i];j<=n;j+=prim[i])
29 if(!vis[j]) num[++top]=j,vis[j]|=1;
30 if(top&1) swap(num[2],num[top]),vis[num[top]]=0,--top;
31 for(int j=1;j<=top;j+=2) tmp[++ans]=P(num[j],num[j+1]);
32 } print(ans);
33 for(int i=1;i<=ans;++i)
34 print(tmp[i].first,' '),
35 print(tmp[i].second);
36 return Ot(),0;
37 }