例题一:分拆素数和(hdu 2098)
把一个偶数拆成两个不同素数之和,有多少种拆法
输入:包含一个正整数,其值不会超过10 000,输入0结束
输出:输出拆成的不同素数对的个数
#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define N 20000
bool isprime[N];///判定i是否为素数
int n;
void init()
{
memset(isprime,1,sizeof(isprime));
isprime[1]=0;
for(int i=2;i<N;i++)
{
if(isprime[i])
{
for(int j=i*i;j<N;j+=i)
isprime[j]=0;
}
}
}
int main()
{
init();
while(scanf("%d",&n)!=EOF)
{
if(n==0)
break;
int num=0;
for(int i=2;i<=n/2+1;i++)
{
if(isprime[i]&&isprime[n-i])
{
if(i!=n-i)///要求两个素数不相同
num++;
}
}
cout<<num<<endl;
}
return 0;
}
例题二 云之遥—素数(nefu 109)
判定一个数是否为素数
输入:数据有多组,每组输入一个N,且N是大于0且小于2 000 000 000 的整数
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
#define N 50000///根号下20亿
bool isprime[N];
int prime[N],nprime;
void doprime()
{
long long i,j;///防止数据相乘后溢出,超出int 型的范围
nprime=0;
memset(isprime,1,sizeof(isprime));
isprime[1]=0;
for(i=2;i<N;i++)
{
if(isprime[i])
{
nprime++;
prime[nprime]=i;
for(j=i*i;j<N;j+=i)
isprime[j]=0;
}
}
}
bool isp(long long n)///基本法求素数
{
int i;
int k=(int)sqrt(double(n));///sqrt函数的原型为double sqrt(double),参数为double类型
for(i=1;prime[i]<=k;i++)///这里是小于等于
if(n%prime[i]==0)
return 0;
return 1;
}
int main()
{
doprime();
long long n;
while(scanf("%lld",&n)!=EOF)
{
if(n==1)
cout<<"NO"<<endl;
else
{
if(isp(n))
cout<<"YES"<<endl;
else
cout<<"NO"<<endl;
}
}
return 0;
}
例题三 素数距离(pku 2689)
题意:求给定区间内的质数距离最小的一对和质数距离最大的一对。
输入:两个正整数L和U,L<U,L和U相差不超过1 000 000
输出:输出两对相邻素数
分析:使用筛选法筛掉[L,U]内的所有非素数,需要知道[L,U]区间的所有非素数的素数因子,21亿内的数或者是素数,或者能被根号21亿内的素数整除,也就是说[L,U]区间内的所有非素数的素因子都在根号21亿内
预先将根号21亿内的所有素数找出来,然后用这些素数去筛掉指定区间的所有非素数
根据素数分布定律,100 000个数中最多有80 000(n/ln(n))个素数
#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define N 50000
#define maxint 0x3f3f3f
long long l,u;
long long prime2[1000005];///存储区间[l,u]内的素数
int prime1[N];///存储根号21亿内的素数
int nprime1,nprime2;///记录所存素数的个数
bool isprime[N];///判定是否为素数
void doprime()
{
long long i,j;
nprime1=0;
memset(isprime,1,sizeof(isprime));
for(i=2;i<N;i++)
{
if(isprime[i])
{
nprime1++;
prime1[nprime1]=i;
for(j=i*i;j<N;j+=i)
isprime[j]=0;
}
}
}
void judge()
{
long long i,j,b;
memset(isprime,1,sizeof(isprime));
for(i=1;i<=nprime1;i++)
{
b=l/prime1[i];
while(b*prime1[i]<l||b<=1)///从第一个大于等于l的数开始筛选
b++;
for(j=b*prime1[i];j<=u;j+=prime1[i])
{
if(j>=l)
isprime[j-l]=0;
}
}
if(l==1)///考录到l==1的情况
isprime[0]=0;
}
void solve()
{
int i;
long long min=maxint;
long long max=-maxint;
long long minl,minr,maxl,maxr;
judge();
nprime2=0;
for(i=0;i<=u-l;i++)///记录区间[l,u]中素数的个数,并将素数顺序存储起来
{
if(isprime[i])
{
nprime2++;
prime2[nprime2]=i+l;
}
}
if(nprime2<=1)
printf("There are no adjacent primes.\n");
else
{
for(i=2;i<=nprime2;i++)
{
int temp=prime2[i]-prime2[i-1];///记录相邻素数的距离
if(temp<min)
{
min=temp;
minl=prime2[i-1];
minr=prime2[i];
}
if(temp>max)
{
max=temp;
maxl=prime2[i-1];
maxr=prime2[i];
}
}
printf("%lld,%lld are closest,%lld,%lld are most distant.\n",minl,minr,maxl,maxr);
}
}
int main()
{
doprime();
while(scanf("%lld%lld",&l,&u)!=EOF)
{
solve();
}
return 0;
}