趣题12:质数拆分

己知1000以内的质数集合A,如41=2+3+17+19,这样的质数可以拆分为两个或两个以上的不同的质数之和称为超质数,求A中最大的超质数?


#include <stdio.h>
#include <stdlib.h>

int s[200];
int prime[200];
int n;

struct node
{
	int num[200];
	int n;
};

node bit[200000];

bool isprime(int n)
{
	if(n == 1) return false;

	for(int i = 2; i <n; ++i)
	{
		if(n % i == 0)
			return false;
	}

	return true;
}

void check(int k)
{
	n = 0;

	for(int i = 2; i < 1000; ++i)
	{
		if(isprime(i))
		{
			prime[n] = i;
			++n;
		}
	}

	s[0] = 2;

	for(int i = 1; i < n; ++i)
	{
		s[i] = s[i-1] + prime[i];
	}

	memset(bit, 0, sizeof(node)*200000);

	k = k > n ? n : k;

	for(int i = 0; i < k; ++i)
	{
		if(bit[prime[i]].n == 0)
		{
			bit[prime[i]].n = 1;
			bit[prime[i]].num[0] = prime[i];
		}

		for(int j = s[i-1]; j >= 2; --j)
		{
			if((bit[j].n > 1 && j != prime[i]) || (j != prime[i] && bit[j].n == 1) )
			{
				for(int k = 0; k < bit[j].n; ++k)
				{
					bit[j+prime[i]].num[k] = bit[j].num[k];
				}

				bit[j+prime[i]].n = bit[j].n+1;
				bit[j+prime[i]].num[bit[j].n] = prime[i];
			}

		}
	}

	for(int m = 0; m <= 1000; ++m)
	{
		if(bit[m].n > 1 && isprime(m))
		{
			printf("%d = ", m);

			for(int t = 0; t < bit[m].n-1; ++t)
			{
				printf("%d+", bit[m].num[t]);
			}

			if(bit[m].n - 1 >= 0)
			{
				printf("%d\n", bit[m].num[bit[m].n - 1]);
			}
		}
	}
}

int main()
{
	check(4);

	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值