449 C. Jzzhu and Apples

最新推荐文章于 2018-10-27 08:26:00 发布

pizzaaaaa

最新推荐文章于 2018-10-27 08:26:00 发布

阅读量399

点赞数

CC 4.0 BY-SA版权

分类专栏： codeforces

本文链接：https://blog.youkuaiyun.com/u011277193/article/details/38363691

codeforces 专栏收录该内容

37 篇文章

订阅专栏

本文介绍了一种使用贪心策略的素数匹配算法，通过线性素数筛法预先计算素数，并利用这些素数及其倍数进行匹配。该算法特别针对大素数的倍数，若倍数个数为奇数，则将2倍素数值空开留给2来匹配，以此达到最优匹配。

/*
 贪心算法,对于大素数的倍数,如果倍数的个数是奇数个,可以考虑将2*素数这个值空开,留给2来匹配
 */
#include <cstdio>
#include <cstring>
using namespace std;

//线性素数筛法
const int N = 100010;
int prime[N] = {0},num_prime = 0;
int isNotPrime[N] = {1, 1};
void produce_prime()
{
    for(int i = 2 ; i < N ; i ++)
    {
		if(! isNotPrime[i])
	 		prime[num_prime ++]=i;
		for(int j = 0 ; j < num_prime && i * prime[j] <  N ; j ++)
        {
            
            isNotPrime[i * prime[j]] = 1;
	  		if( !(i % prime[j] ) )
				break;
		}
	}
}

int match[N];
bool vis[N];
int main()
{
    produce_prime();
    int n;
    scanf("%d",&n);
    int save[N];
    int tmp[N];
    int cnt=0;
    for(int i=0;i<num_prime;i++)
    {
        if(prime[cnt]*2>n)  break;
        save[cnt++]=prime[i];
    }
    memset(vis,false,sizeof(vis));
    memset(match,0,sizeof(match));
    int ans=0;
    for(int i=cnt-1;i>=0;i--)
    {
        int c=0;
        for(int j=1;;j++)
        {
            if(prime[i]*j>n)  break;
            if(!vis[prime[i]*j])
            {
                tmp[c++]=prime[i]*j;
            }
        }
        if(c%2==0)
        {
            for(int j=0;j<c;j+=2)
            {
                ans++;
                match[tmp[j]]=tmp[j+1];
                vis[tmp[j]]=true;
                vis[tmp[j+1]]=true;
            }
        }
        else if (c%2==1&&c>1)
        {
            ans++;
            match[tmp[0]]=tmp[2];
            vis[tmp[0]]=true;
            vis[tmp[2]]=true;
            for(int j=3;j<c;j+=2)
            {
                ans++;
                match[tmp[j]]=tmp[j+1];
                vis[tmp[j]]=true;
                vis[tmp[j+1]]=true;
            }
        }
    }
    printf("%d\n",ans);
    for(int i=2;i<=n;i++)
    {
        if(match[i]!=0)
        {
            printf("%d %d\n",i,match[i]);
        }
    }
    return 0;
}