Miller-Rabin随机判断素数

#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 b,long long mod)
{
	long long t = a % mod,r = 1;
	while(b)
	{
		if(b & 1) r = multi(r,t,mod);
		t = multi(t,t,mod); b = b>>1;
	}
	return r;
}
bool Miller_Rabin(long long n)
{
	long long time = 10;
	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))
	{
		if(Miller_Rabin(n)) printf("Yes\n");
		else printf("No\n");
	}
	return 0;
}