BestCoder Round #18 1001 primes problem

Primes Problem

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 659    Accepted Submission(s): 302

Problem Description

Given a number n, please count how many tuple(p1, p2, p3) satisfied that p1<=p2<=p3, p1,p2,p3 are primes and p1 + p2 + p3 = n.

 

Input

Multiple test cases(less than 100), for each test case, the only line indicates the positive integer n(n≤10000) .

 

Output

For each test case, print the number of ways.

 

Sample Input

  
  

   
   3
9
  
  

 

Sample Output

  
  

   
   0
2
题意很简单啦，给出一个n问有多少组合使得p1+p2+p3=n，p1<=p2<=p3 且p1、p2,p3都是素数,显然这是一个关于素数的题目，题意很直接。
很容易想到暴力，但是作为BC的题目自然是要考虑很多因素，我们来分析一下时间复杂度。
对于n的最大值10000，很容易打个表就知道有1229个素数，那么我们考虑极端数据情况下，如果直接先把表打好之后，对每个n，进行p1,p2,p3的三重循环暴力枚举，
这样的话运算次数会是1229的三次方，即10^9运算，而1000MS的时间限制也差不多就是10^9运算，在再加上最多100组测试数据，这样肯定会超时！
所以，不能暴力枚举，应该进行优化，我们注意到其实p1,p2定出来之后，那么需要考虑的p3也确定了为n-p2-p1，那么我们可以将p3的循环改为测试可能的p3是不是素数，
这里用朴素的素数判定法（时间复杂度为n^1/2）而且p3至少为n/3，这样复杂度就变成了1229*1229*（n/3）^1/2 再算上100组数据，和打表的时间应该是勉强能过的！
好了，下面帖代码！

   
   
#include<cstdio>
#include<iostream>
#include<string.h>
int p[2000],tag[10100];int cnt;
void get(int n)  //这个线性筛法好久没用了，刚写的时候还不会写了，O(∩_∩)O哈！，要多练啊！
{
	cnt=0;
	for(int i=2;i<=n;i++)
	{	
		if(!tag[i]) p[cnt++]=i;
		for(int j=0;j<cnt&&p[j]*i<=n;j++)
		{
			tag[i*p[j]]=1;
			if(i%p[j]==0)
			break;
		}
	}
}
bool check(int num)
{	bool flag=true;
	for(int i=2;i*i<=num;i++)
	{
		if(num%i==0)
		{
			flag=false;
			break;
		}
	}
	return flag;
}
int main()
{
	int n;int ans=0;
	get(10000);
	while(scanf("%d",&n)!=EOF)
	{	ans=0;memset(tag,0,sizeof(tag));
	for(int i=0;i<cnt;i++)
	{	
		for(int j=i;j<cnt;j++)
		{
			if(p[i]+p[j]>n)
			break;
			if(check(n-p[i]-p[j])&&(n-p[i]-p[j])>=p[j])
			{
				//printf("%d %d %d\n",p[i],p[j],n-p[i]-p[j]);
				ans++;
			}
		}
	}
	printf("%d\n",ans); 
	}
} 
-------------以下是超时代码！----------------------------
   
   

   
   
#include<cstdio>
#include<iostream>
#include<string.h>
int p[8000],tag[10100];int cnt;
void get(int n)
{
	 cnt=0;
	for(int i=2;i<=n;i++)
	{	
		if(!tag[i]) p[cnt++]=i;
		for(int j=0;j<cnt&&p[j]*i<=n;j++)
		{
			tag[i*p[j]]=1;
			if(i%p[j]==0)
			break;
		}
	}
}
int main()
{
	int n;int ans=0;
	get(10000);
	while(scanf("%d",&n)!=EOF)
	{	ans=0;memset(tag,0,sizeof(tag));
	for(int i=0;i<cnt;i++)
	{	
		for(int j=i;j<cnt;j++)
		{
			if(p[i]+p[j]>n)
			break;
			for(int k=j;k<cnt;k++)
			{
				if(p[i]+p[j]+p[k]>n)
				break;
				if(p[i]+p[j]+p[k]==n)
				{
					//printf("%d %d %d\n",p[i],p[j],p[k]);
				ans++;
				}
				
			}
		}
	}
	printf("%d\n",ans); 
	}
}