NYOJ-468 Fibonacci数列（三）【大素数判断】

题目链接：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;
}