但是我们还可以观察到另一个性质,比如说6,他的因子有1,2,3,6,其中1*6=6,2*3=6。所以我们可以看到一个显而易见的优化,那就是只枚举上述两个式子中较小的数字,也就是1(其实不算),和2。举个例子,4,我们就要枚举到sqrt(4)。所以我们又得到了一个更为高效的程序:bool isPrime(int Number){
for (int i=2;i<Number;++i){
if (Number % i == 0){
return true;
}
}
return false;
}
# include <cmath>
bool isPrime(int Number){
for (int i=2;i<=sqrt(Number);++i){
if (Number % i == 0){
return true;
}
}
return false;
}
# include <cmath>
bool isPrime(int Number){
if ((Number & 1) == 0){
return false;
}
for (int i=3;i<=sqrt(Number); i += 2){
if (Number % i == 0){
return true;
}
}
return false;
}
# include <cstring>
const int MAXN = 10001;
bool Primes[MAXN]; //质数表
int PrimesSize = 0; //质数表的大小
void MakePrimeTables(int MaxSize){
if (MaxSize <= PrimesSize){
//无须继续打表,因为已经打过了
return;
}
else{
memset(Primes,true,sizeof(Primes));
Primes[0] = false;
Primes[1] = false;
int NowPos = 2;
while (NowPos < MaxSize){
if (Primes[NowPos]){
int TimesNumber = NowPos + NowPos;
while (TimesNumber <= MaxSize){
Primes[TimesNumber] = false;
TimesNumber += NowPos;
}
}
}
}
}
int MAXN = 100001;
int MAXP = 10001;
bool Primes[MAXN];
int PrimeNumber[MAXP];
void GetPrime()
{
memset(Primes,true,sizeof(Primes));
memset(PrimeNumber,0,sizeof(PrimeNumber));
for(int i=2;i<=MAX;i++)
{
if(Prime[i]) PrimeNumber[PrimeNumberSize++]=i; // 说明i是一个素数
for(int j=0; j<PrimeNumberSize && i*PrimeNumber[j]<=MAX; j++) // 遍历所有的质因子
{
Primes[i*PrimeNumber[j]]=false;
if(i%PrimeNumber[j]==0) break; // 如果i能够整除pri[j]那么i*pri[j+1]就一定被pri[j]数整除。
}
}
}
这里用到了随机数,所以可能有一点不够稳定,而这个算法的时间复杂度大概是O(log2n),运用了快速幂的优化之后相当不错。//****************************************************************
// Miller_Rabin 算法进行素数测试
//速度快,而且可以判断 <2^63的数
//****************************************************************
const int S = 20; //随机算法判定次数,S越大,判错概率越小
//计算 (a*b)%c. a,b都是long long的数,直接相乘可能溢出的
// a,b,c <2^63
long long mult_mod(long long a, long long b, long long c)
{
a %= c;
b %= c;
long long ret = 0;
while(b)
{
if(b & 1)
{
ret += a;
ret %= c;
}
a <<= 1;
if(a >= c)a %= c;
b >>= 1;
}
return ret;
}
//计算 x^n %c
long long pow_mod(long long x, long long n, long long mod) //x^n%c
{
//Fast-Pow
if(n == 1)return x % mod;
x %= mod;
long long tmp = x;
long long ret = 1;
while(n)
{
if(n & 1) ret = mult_mod(ret, tmp, mod);
tmp = mult_mod(tmp, tmp, mod);
n >>= 1;
}
return ret;
}
//以a为基,n-1=x*2^t a^(n-1)=1(mod n) 验证n是不是合数
//一定是合数返回true,不一定返回false
//利用费尔马小定理:
bool check(long long a, long long n, long long x, long long t)
{
/*a^x%n*/
long long ret = pow_mod(a, x, n);
long long last = ret;
for(int i = 1; i <= t; i++)
{
/*ret*ret%n*/
ret = mult_mod(ret, ret, n);
if(ret == 1 && last != 1 && last != n - 1) return true; //合数
last = ret;
}
if(ret != 1) return true;
return false;
}
// Miller_Rabin()算法素数判定
//是素数返回true.(可能是伪素数,但概率极小)
//合数返回false;
bool Miller_Rabin(long long n)
{
/*防止误判*/
if(n < 2)return false;
if(n == 2)return true;
/*提高效率的特判*/
if((n & 1) == 0) return false; //偶数
long long x = n - 1;
long long t = 0;
while((x & 1) == 0)
{
x >>= 1;
t++;
}
for(int i = 0; i < S; i++)
{
/*
下面这里是以随机数来测试质数
但是,可以直接测试某一些数,就直接100%零错率做到又快又好
测试,但仍然要排除一些强伪素数
*/
long long a = rand() % (n - 1) + 1;
if(check(a, n, x, t))
return false;//合数
}
return true;
}