poj3842-效率啊

这是一道算法选择不好会TLE的题目,同时这道题目涉及到的内容也非常丰富。至少我学到了一些知识。

首先题目的内容很简单,就是最多给你7个数字,问你用着7个数字能组成多少个素数。

 

由于原来搞过线性筛选素数,于是就将9999999以内的素数保存下来,实际上这个过程大概需要半秒钟吧。然后枚举枚举每个素数,看看这个素数用到的数字是不是在给的7个里面,如果是就+1.代码很快敲完了,发现交上去TLE。

后来想想也是,最多200组数据,如果每次都是9999999,那么每次都要枚举30万个素数,肯定超时,想不好算法就下手很容易吃力不讨好。。。。

 

如果这样不行,那么只能对这几个数字进行全排列,并且判断每一个是不是素数,走这个思路,好像不会超时,同时也避免了9999999这种重复数字,但是组合的数字很少的情况。

算法:

想找出3400以内的素数,这个很快。为什么要找3400以内的呢?因为如果要判断9999999以内的数是不是素数,只需要让这个数对sqrt(9999999)以内的素数取余数,看看是不是0就行了。

然后求这几个数字的组合,因为题目并没有说所有数字都要用,可以1个2个。。。。

然后对每个组合求全排列

然后对每个全排列求数,判断是不是素数,同时这里要注意,如果重复了,就不要判了,还要判重。(比如9999999,如果组合3个数,有很多个999,这些就没必要全排列了,因为已经算过一次了,直接pass)

最后得出结果,输出有多少个。

 

用到知识点

求素数

求组合

求全排列

这些知识点放到后面介绍吧。。。。。这里直接上代码了,如果想了解如何求快速素数,求组合,求全排列,请看下一篇。。。

代码:

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <set>
using namespace std;
char ch[10];
int num[10];
int zuhe[10];
bool f[3400];
int prim[1500];
int pNum = 0;

void Prime()//get primer
{
	for (int i = 2; i < 3400; ++ i)
	{
		if (! f[i])
		{
			prim[pNum ++] = i;
		}
		for (int j = 0; (j < pNum) && (prim[j] * i < 3400); ++j)
		{
			f[prim[j] * i] = true;
			if (i % prim[j] == 0)
			{
				break;
			}
		}
	}
}
bool IsPrime(int tmp)//judge it is a primer or not
{
	if (tmp < 2)
	{
		return false;
	}
	for (int i = 0;prim[i] * prim[i] <= tmp; ++ i)
	{
		if (tmp % prim[i] == 0)
		{
			return false;
		}
	}
	return true;
}
int getNum(int len)//get the num for the digits
{
	int sum = 0;
	for (int i = 0;i < len; ++ i)
	{
		sum = sum * 10 + zuhe[i];
	}
	return sum;
}
void Char2Int()// char to digits
{
	for (int i = 0; i < strlen(ch); ++i)
	{
		num[i] = ch[i] - '0';
	}
}
int main()
{
	int n;
	Prime();//init 
	while (scanf("%d", &n) != EOF)
	{
		while (n --)
		{
			int times = 0;
			scanf("%s", ch);
			Char2Int();
			int len = strlen(ch);
			sort(num, num + len);	//sort	
			set <int> SS;//avoid repeating 
			for (int i = 1; i < (1 << len); ++ i)// get zuhe
			{
				
				int len1 = 0;
				for (int j = 0; j < len; ++ j)
				{
					if (i & (1 << j))
					{
						zuhe[len1 ++] = num[j];
					}
				}
				do 
				{
					int tmp = getNum(len1);
					if (! SS.count(tmp) && IsPrime(tmp))
					{
						SS.insert(tmp);
						times ++;
					}
				} while (next_permutation(zuhe, zuhe + len1));//qu pai lie 

			}
			printf("%d\n", times);
		}
	}
}


 

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值