题目链接:http://acm.nyist.net/JudgeOnline/problem.php?pid=468
解题思路:
这道题考察的是大素数的判断,但是有一个知识点。可以在黑书的221页找到。
就是斐波那契数列从第5项开始,如果它的项数为素数,那么它就是斐波那契素数。因为斐波那契数列为1,1,2,3,5。其实第二项为1,而1不是素数,所以不是斐波那契素数,第4项为3,虽然3和1,2互素,但是它的项数是4,4不是素数,所以这是一个特殊情况。
综上,第二项和第四项是特殊情况,需要特殊考虑,其他情况均满足这个性质。
斐波那契素数和素数不同:
斐波那契素数指:该项的值与前面所有的斐波那契数互素。
而素数是指只能被1和本身整除。
所以,斐波那契素数不一定是素数。
大素数判断模板,果断水过。。。。
代码如下:
#include<stdio.h>
#include<algorithm>
#include<iostream>
using namespace std;
long long multi(long long a, long long b, long long mod)
{
long long temp = a, sum = 0;
while(b)
{
if(b & 1) sum = (sum + temp) % mod;
temp = (temp + temp) % mod;
b = b >> 1;
}
return sum;
}
long long Modular_exponent(long long a, long long x, long long mod)
{
long long t = a % mod, r = 1;
while(x)
{
if(x & 1) r = multi(r, t, mod);
t = multi(t, t, mod); x = x >> 1;
}
return r;
}
bool Miller_Rabin(long long n)
{
long long time = 20;
if(n < 2) return false;
if(n == 2) return true;
bool flag = false;
for(long long k = 0; k <= time; ++k)
{
flag = false;
long long d = n - 1, r = 0, t, a = rand() % (n - 2) + 2;
while((d & 1) == 0)
{
d = d>>1;
r++;
}
t = Modular_exponent(a, d, n);
if(t == 1 || t == n-1) {flag = true; continue;}
for(long long i = 1; i < r; i++)
{
t = multi(t, t, n);
if(t == 1) {flag = false; return flag;}
if(t == n-1) {flag = true; break;}
}
if(!flag) break;
}
return flag;
}
int main()
{
long long n;
while(~scanf("%lld",&n))
puts((n == 4 || (n != 2 && Miller_Rabin(n))) ? "Yes" : "No");
return 0;
}