sgu113

最新推荐文章于 2016-04-14 18:59:16 发布

原创最新推荐文章于 2016-04-14 18:59:16 发布 · 840 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#sgu113 #nearly prime numbers #欧拉筛法 #素数

SGU 专栏收录该内容

20 篇文章

订阅专栏

SGU113 Nearly Prime Numbers

题目大意：

求一个数N是否能由两个质数p1*p2得到(p1可以=p2)

输入：

数据组数Q(1<=Q<=10)

Q个数，为题目要求的数字N(1<=n<=10^9)

输出：

对于每个数字，如果找得到，则输出"Yes"，否则输出"No"，用回车隔开。

样例输入：

样例输出：

Yes

看似比较水的一道题目(实际也比较水)。

首先我们不妨设p1<=p2，然后用欧拉筛打出10^6以内的质数表。

枚举p1看看p2是不是质数(如果比10^6大就用性质：因子只有两个判断一下，如果比10^6小，嘿嘿，不要忘了自己手上的表！)

输出答案即可。

注意事项：

1.请仔细阅读我的翻译题目(原本题目让你把Yes和No摆一排输出，你要这样就等着WA吧)。

2.sqrt后的数字转double是因为LINUX的需要(SGU使用WINDOWS，所以不用改，当然改了也不会错)。

3.不要犯大小写这种傻逼错误。。。

下面附上我的代码：

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <stdbool.h>
#include <math.h>
#define SQRT(a) ((int)(sqrt((double)(a))))
int n;
int prime[10001],total;
bool check[100001];
void eular() //欧拉筛
{
  int i,j;
  check[0]=check[1]=false;
  for (i=2;i<=100000;i++)
    {
	if (check[i])
	  prime[++total]=i;
	for (j=1;j<=total;j++)
	  {
	  if (i*prime[j]>100000)
	    break;
	  check[i*prime[j]]=false;
	  if (!i%prime[j])
	    break;
	  }
	}
  return ;
}
bool prim(int a) //判断a是否是质数
{
  int i,maxa=SQRT(a);
  for (i=1;prime[i]<=SQRT(a);i++)
    if (a%prime[i]==0)
      return false;
  return true;
}
bool work(int a)
{
  int p1,p2,i,j;
  int maxa=SQRT(a);
  for (i=1;prime[i]<=maxa;i++)
    {
    p1=prime[i];
    if (a%p1==0)
      {
	  p2=a/p1;
	  if (p2<100000 && check[p2])
	    return true;
	  else if (p2>100000 && prim(p2))
	    return true;
	  return false;
	  }
    }
  return false;
}
void init()
{
  int i;
  memset(check,true,sizeof(check));
  eular();
  scanf("%d",&n);
  for (;n;n--)
    {
    scanf("%d",&i);
    if (work(i))
      printf("Yes\n"); //请千万记住输出回车
    else
      printf("No\n"); //请千万记住输出回车
	}
  return ;
}
int main()
{
  init();
  return 0;
}