Miller Rabin算法：大质数判断

问题概述：判断一个数n是不是质数(n<=10^18)

输入样例：                              对应输出：

7                                              Yes

费马小定理：如果p是质数，且a,p互质，那么a^(p-1)%p==1

Miller-Rabin算法的理论基础：如果p是一个大于2的质数，先将p-1表示成2^s*r的形式(r是奇数)，令a是和n互素的任

意整数，那么a^r%p==1或对某个j(0<=j<=s-1,j∈Z)，等式a^(2^j*r)%p==±1成立

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<algorithm>
#include<iostream>
#include<math.h>
using namespace std;
#define LL long long
LL Multi(LL a, LL b, LL mod)
{
	LL ans = 0;
	a %= mod;
	while(b)
	{
		if(b%2==1)  ans = (ans+a)%mod, b--;
		else  a = (a+a)%mod, b /= 2;
	}
	return ans;
}

LL Pow(LL a, LL b, LL mod)
{
	LL ans = 1;
	a %= mod;
	while(b)
	{
		if(b&1)  ans = Multi(ans, a, mod), b--;
		else  a = Multi(a, a, mod), b /= 2;
	}
	return ans;
}

int Miller_Rabin(LL n)
{
	int i, j, k;
	LL a, x, y, mod;
	if(n==2)  return 1;
	if(n<2 || n%2==0)  return 0;
	k = 0, mod = n-1;
	while(mod%2==0)
	{
		k++;
		mod /= 2;
	}
	for(i=1;i<=10;i++)
	{
		a = rand()%(n-1)+1;			//随机生成一个小于n的正整数
		x = Pow(a, mod, n);
		y = 0;
		for(j=1;j<=k;j++)
		{
			y = Multi(x, x, n);
			if(y==1 && x!=1 && x!=n-1)
				return 0;
			x = y;
		}
		if(y!=1)
			return 0;
	}
	return 1;
}

int main(void)
{
	LL n;
	scanf("%I64d",&n);
	if(Miller_Rabin(n))
		printf("Yes\n");
	else
		printf("No\n");
	return 0;
}