Miller_Rabin

Posted 2020-10-09

#include <bits/stdc++.h>
using namespace std;
int power(int a,int k,int mod){
int ans=1;
	while(k){
		if(k&1)ans*=a;
		a*=a;
		a%=mod;
		ans%=mod;
		k>>=1;
	}
	return ans;
}
bool isp(int n){
int a;
bool flag=false;
        if(n==2)return true;
	if(n==3215031751)return false;
	for(int i=1;i<=10;i++){
		a=rand()%(n-2)+2;
		if(power(a,n-1,n)!=1)
			return false;
	}
	return true;
}
int main(){
int n,a,b=0;
	srand(time(0));
	cin>>n;
	if(isp(n))
		puts("Yes");
	else
		puts("No");
	return 0;
}