做题感悟:这题是在学习了区间筛法后才做的,学习了区间筛素数后,这题真 so easy !
解题思路:可以用两次筛法,也可以直接用区间平移。
代码:
#include<stdio.h>
#include<iostream>
#include<map>
#include<stack>
#include<string>
#include<string.h>
#include<stdlib.h>
#include<math.h>
#include<vector>
#include<queue>
#include<algorithm>
using namespace std ;
#define LEN sizeof(struct node)
#define lld long long int
const double PI = 3.1415926535898 ;
const double INF = 99999999 ;
const long long mod= 1000 ;
const int MX = 1000005 ;
bool is_prime[MX] ;
int prime[MX],num ;
void search_prime(int L,int U)
{
num=0 ;
int n=sqrt(1.0*U),j ;
int d=U-L+1 ; // 区间长度
memset(is_prime,false,sizeof(is_prime)) ;// 先默认全部为素数
for(int i=( L%2 != 0 ) ;i<d ;i+=2) // 把区间中的偶数筛掉
is_prime[i]=true ;
for(int i=3 ;i<=n ;i+=2)
{
if(i>L&&is_prime[i-L])// 如果 i 大于 L且is_prime[i-l] 为真则必为合数
continue ;
j=(L/i)*i ;// j 为最接近L 的和数且为 i的倍数
if(j<L)
j+=i ;
if(j==i) j+=i ;
j-=L ;
for( ;j<d ;j+=i) // j每次加 i 均为合数 :j=(L/i)*i+i+i……
is_prime[j]=true ;
}
if(L<=1) is_prime[1-L]=true ;// 0 1 均不是素数
if(L<=2) is_prime[2-L]=false ;
for(int i=0 ;i<d ;i++)
if(!is_prime[i])
prime[num++]=i+L ;
}
void find()
{
int min=9999999,max=0 ;
int a1=-1,a2=-1,a3=-1,a4=-1 ;
for(int i=1 ;i<num ;i++)
{
int mx=prime[i]-prime[i-1] ;
if(mx<min)
{
a1=prime[i-1] ;
a2=prime[i] ;
min=mx ;
}
if(mx>max)
{
a3=prime[i-1] ;
a4=prime[i] ;
max=mx ;
}
}
if(a1==-1)
printf("There are no adjacent primes.\n") ;
else printf("%d,%d are closest, %d,%d are most distant.\n",a1,a2,a3,a4) ;
}
int main()
{
int L,U ;
while(~scanf("%d%d",&L,&U))
{
search_prime(L,U) ;
find() ;
}
return 0 ;
}
代码(二次筛法):
#include<stdio.h>
#include<iostream>
#include<map>
#include<stack>
#include<string>
#include<string.h>
#include<stdlib.h>
#include<math.h>
#include<vector>
#include<queue>
#include<algorithm>
using namespace std ;
#define LEN sizeof(struct node)
#define lld long long int
const double PI = 3.1415926535898 ;
const double INF = 99999999 ;
const long long mod= 1000 ;
const int MX = 1000005 ;
int num=0 ;
bool is_prime[MX] ;
__int64 prime[MX] ;
__int64 prime2[MX] ;
void init() // 用筛法筛出第一批素数,用于筛指定区间的素数
{
num=0 ;
memset(is_prime,false,sizeof(is_prime)) ;
for(__int64 i=2 ;i<MX ;i++)
{
if(!is_prime[i])
prime[num++]=i ;
for(int j=0 ;j<num&&i*prime[j]<MX ;j++)
{
is_prime[i*prime[j]]=true ;
if(i%prime[j]==0)
break ;
}
}
}
void search_prime(__int64 L,__int64 U)
{
memset(is_prime,false,sizeof(is_prime)) ;
for(int i=0 ;i<num&&prime[i]*prime[i]<=U ;i++)
{
__int64 p=prime[i] ;
for(__int64 j=(L+p-1)/p ;p<=U/j ;j++)
if(j>1)
is_prime[j*p-L]=true ;
}
if(L<=1) is_prime[1-L]=true ;
if(L<=2) is_prime[2-L]=false ;// 巧妙
int nx=0 ;
for(__int64 i=L ;i<=U ;i++)
if(!is_prime[i-L])
prime2[nx++]=i ; // 存区间的素数
__int64 min=99999999,max=0,a1=-1,a2=-1,a3=-1,a4=-1 ;
for(int i=1 ;i<nx ;i++) // 寻找答案
{
__int64 mx=prime2[i]-prime2[i-1] ;
if(min>mx)
{
a1=prime2[i-1] ;
a2=prime2[i] ;
min=mx ;
}
if(max<mx)
{
a3=prime2[i-1] ;
a4=prime2[i] ;
max=mx ;
}
}
if(a1!=-1)
printf("%I64d,%I64d are closest, %I64d,%I64d are most distant.\n",a1,a2,a3,a4) ;
else
printf("There are no adjacent primes.\n") ;
}
int main()
{
__int64 L,U ;
init() ;
while(~scanf("%I64d%I64d",&L,&U))
{
search_prime(L,U) ;
}
return 0 ;
}