数论---------质数距离

给定两个整数L和U,你需要在闭区间[L,U]内找到距离最接近的两个相邻质数C1和C2(即C2-C1是最小的),如果存在相同距离的其他相邻质数对,则输出第一对。
同时,你还需要找到距离最远的两个相邻质数D1和D2(即D1-D2是最大的),如果存在相同距离的其他相邻质数对,则输出第一对。
输入格式
每行输入两个整数L和U,其中L和U的差值不会超过1000000。
输出格式
对于每个L和U ,输出一个结果,结果占一行。
结果包括距离最近的相邻质数对和距离最远的相邻质数对。(具体格式参照样例)
如果L和U之间不存在质数对,则输出“There are no adjacent primes.”。
数据范围
1≤L<U≤231−11≤L<U≤231−1
输入样例:
2 17
14 17

输出样例:
2,3 are closest, 7,11 are most distant.
There are no adjacent primes.

思路:线性筛法只能筛选出1-n的数,但是要筛选出很大的数现在,而根据性质,每一个很大的数n都存在一个小于等于根号n的质因子,所以先筛选出1到根号n的所有质因子,然后拿这些质因子筛选出那些范围内的合数,然后再依次枚举比较

#include <cstring>
#include <cstdio>
#include <algorithm>
#include <iostream>
using namespace std;
const int N = 1000010;
typedef long long LL;
int primes[N], cnt;
bool st[N];
void init(int n){
 memset(st, 0, sizeof st);
 cnt = 0;
 for (int i = 2; i <= n; i ++){
  if (!st[i])   primes[cnt ++] = i;
  for (int j = 0; primes[j] * i <= n; j ++){
   st[i * primes[j]] = true;
   if (i % primes[j] == 0)   break;
  }
 }
}
int main(){
 int l, r;
 while(cin >> l >> r){
  init(50000);
  memset(st, 0, sizeof st);
   for(int i = 0; i < cnt; i ++){
   LL p = primes[i];
   for (LL j = max(2 * p, (l + p - 1) / p * p); j <= r; j += p)
   st[j - l] = true;
  }
   cnt = 0;
    for (int i = 0; i <= r - l; i ++)
   if (!st[i] && i + l >= 2)
   primes[cnt ++] = l + i;
    if (cnt < 2)   puts("There are no adjacent primes.");
   else{
    int minp = 0, maxp = 0;
    for (int i = 0; i + 1 < cnt; i ++){
     int d = primes[i + 1] - primes[i];
     if (d < primes[minp + 1] - primes[minp])   minp = i;
     if (d > primes[maxp + 1] - primes[maxp])   maxp = i;
    }
         printf("%d,%d are closest, %d,%d are most distant.\n",
                primes[minp], primes[minp + 1],
                primes[maxp], primes[maxp + 1]);
   }
 }
 return 0;
}
发布了106 篇原创文章 · 获赞 67 · 访问量 5430

猜你喜欢

转载自blog.csdn.net/qq_45772483/article/details/104940293