Description
把1~n分成若干组,每组两个不互素的数,问最多可以分多少组
Input
一个整数n(1<=n<=1e5)
Output
输出最多可以分的组数和每组的两个数
Sample Input
6
Sample Output
2
6 3
2 4
Solution
1和大于n/2的素数肯定无法和任何一个数组成一组,那么只要有一种构造方案使得其他的数全部配对(奇数就扔一个)则一定是最优的
从小到大枚举素因子d,把所有还没被用的且是d的倍数的数全部拿出来,如果是偶数那么这些数两两配对即可,如果是奇数那么把2d拿出来,这样最后就剩下若干2d,这些数任意两个的最大公因子至少是2所以也可以随便组,故剩下的数也可以两两配对,进而除1和大于n/2的素数外其他所有数都配对了,即为最优方案
Code
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
#define INF 0x3f3f3f3f
#define maxn 111111
int mark[maxn],vis[maxn],a[maxn],ans[maxn][2];
void get_prime(int n=100000)
{
memset(mark,0,sizeof(mark));
for(int i=2;i<=n;i++)
if(!mark[i])
for(int j=2*i;j<=n;j+=i)
mark[j]=1;
}
int main()
{
get_prime();
int n;
while(~scanf("%d",&n))
{
memset(vis,0,sizeof(vis));
vis[1]=1;
for(int i=n/2+1;i<=n;i++)
if(!mark[i])vis[i]=1;
int res=0;
for(int i=3;i<=n/2;i++)
if(!mark[i])
{
int temp=0;
for(int j=i;j<=n;j+=i)
if(!vis[j])a[temp++]=j,vis[j]=1;
swap(a[0],a[1]);
if(temp%2)vis[2*i]=0;
for(int j=temp%2;j<temp;j+=2)ans[res][0]=a[j],ans[res++][1]=a[j+1];
}
int temp=0;
for(int i=1;i<=n;i++)
if(!vis[i])a[temp++]=i;
for(int i=0;i+1<temp;i+=2)ans[res][0]=a[i],ans[res++][1]=a[i+1];
printf("%d\n",res);
for(int i=0;i<res;i++)printf("%d %d\n",ans[i][0],ans[i][1]);
}
return 0;
}