实验5 素数

1

小明是一个聪明的孩子，对数论有着很浓烈的兴趣。他发现求1到正整数10n 之间有多少个素数是一个很难的问题，该问题的难以决定于n 值的大小。现在的问题是，告诉你n的值，让你帮助小明计算小于10n的素数的个数值共有多少位？输入数据有若干组，每组数据包含1个整数n（1 < n < 1000000000），若遇到EOF则处理结束。对应每组数据，将小于10^n 的素数的个数值的位数在一行内输出，格式见样本输出。同组数据的输出，其每个尾数之间空一格，行末没有空格。

#include <bits/stdc++.h>
using namespace std;
 
int main(){
  int n = 0;
  while (cin >> n) {
    cout << int(n - log10(n) - log10(log(10)) + 1) << endl;
  }
  return 0;
}

2

1742年，哥德巴赫在教学中发现，每个不小于6的偶数都是两个素数（只能被1和它本身整除的数）之和。如6＝3＋3，14＝3＋11等等。公元1742年6月7日哥德巴赫写信给当时的大数学家欧拉，欧拉在6月30日给他的回信中说，他相信这个猜想是正确的，但他不能证明。叙述如此简单的问题，连欧拉这样首屈一指的数学家都不能证明，这个猜想便引起了许多数学家的注意。从哥德巴赫提出这个猜想至今，许多数学家都不断努力想攻克它，但都没有成功。我们不需要你去证明哥德巴赫猜想。 如果哥德巴赫猜想是正确的，一个（不小于6的）偶数，都是两个素数之和。那么这个偶数能被至少一个素数对表示，如14，即可以表示为14＝3＋11，也可以表示为14=7+7。不同的偶数对应的素数对的数目是不一样的，如偶数6，就只能表示为6=3+3。对于每个给定的偶数，我们希望知道有多少素数对的和等于该偶数。输入： 有多组测试数据。每组测试数据占一行，包含唯一的一个正偶数n.(6 <= n <= 1e8)

输出：对于每个输入的偶数，输出一行包含唯一的一个整数：表示有多少个素数对的和是输入的偶数。

例子输入：

6

14

例子输出：

1

2

这题对时间复杂度和空间复杂度都有要求，不能直接使用线性筛

#include <bits/stdc++.h>

using namespace std;

#define N 10000



int tag[N];

int prime[N];

int len = 0;

//线性筛函数
void findprime() {

    memset(tag,0,sizeof(tag));

    for(int i = 2; i <= N; ++i) {

      if(!tag[i]) prime[++len] = i;

      for(int j = 1; j <= len && i * prime[j] <= N; ++j) {

        tag[i * prime[j]] = 1;

        if(i % prime[j] == 0) break;

      }

    }

}

//判断是否素数
bool isprime(int x) {

  for(int i = 1; x > prime[i] && i <= len; ++i){

    if(x % prime[i] == 0) return false;

  }

  return true;

}



int main(){

  int n = 0, ans = 0;

  findprime();

  while(cin >> n) {

    ans = 0;

    for(int i = 2; i <= n / 2; ++i)

    if(isprime(i) && isprime(n - i))  ++ans;

    cout << ans << endl;

  }

}

3

分拆素数和

把一个偶数拆成两个不同素数的和，有几种拆法呢？

输入：输入包含一个偶数，其值不会超过10000

输出：对应这个偶数，输出其拆成不同素数的个数，每个结果占一行

例子输入：

30

例子输出：

3

这题和第二题大体上类似，但是不能取相同的素数。n的范围不超过10000，可以直接使用线性筛法

#include <bits/stdc++.h>
using namespace std;

bool isPrime(int x) {
  if(x < 2) return false;
  for(int i = 2; i <= sqrt(x); ++i) {
    if(!(x % i))  return false;
  }
  return true;
}

int main(){
  int n = 0, ans = 0;
  while (cin >> n) {
    ans = 0;
    for(int i = 2; i <= n/2; ++i) {
      if(i != n - i && isPrime(i) && isPrime(n - i))  ++ans;
    }
    cout << ans << endl;
  }
  return 0;
}

4

特殊质数：

如果一个数是质数，而且他在质数表中的编号也是质数，称为特殊的质数。比如质数表中第7个质数是17，因为7和17都是质数，所以我们称17是特殊的质数；比如第10个质数是29，因为10不是质数，所以29不是特殊的质数。现给定非负整数n，求<=n的最大的特殊的质数。输入数据有多组：多组输入数据n（0<=n<=10^6）输出小于等于n的最大的特殊的质数，对于给定的n如果不存在这样的质数，输出-1。

例子输入：

0

20

50

60

例子输出：

-1

17

41

59

<br/>

NEFU OJ的1151题

#include <bits/stdc++.h>
using namespace std;
#define N 1000000

int tag[N];
int prime[N];
int len = 0;
//线性筛函数
void findprime() {
    memset(tag,0,sizeof(tag));
    for(int i = 2; i <= N; ++i) {
      if(!tag[i]) prime[++len] = i;
      for(int j = 1; j <= len && i * prime[j] <= N; ++j) {
        tag[i * prime[j]] = 1;
        if(i % prime[j] == 0) break;
      }
    }
}
bool isprime(int x) {
  if(x < 2) return false;
  for(int i = 2; i <= sqrt(x); ++i)
    if(x % i == 0)  return false;
  return true;
}

int main(){
  int n = 0, ans = 0, k = 0;
  findprime();
  while(cin >> n) {
  ans = -1, k = 0;
  for(int i = len; i >= 1; --i) {
    if(prime[i] <= n && isprime(i)) {
      ans = prime[i]; break;
    }
  }
  cout << ans << endl;
  }
}