1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
#include<stdio.h>
#include<stdlib.h>
//计算a*b(mod n) 
long long PowerMulti(long long a, long long b, long long n){
    long ans = 0;
    while(b){
        if(b&1){
            ans = (ans+a)%n;
        }
        b >>= 1;
        a = (a*2)%n;
    }
    return ans;
}
//快速幂，计算a^b(mod n) 
long long PowerMod(long long a, long long b, long long n){
    long long ans = 1;
    while(b){
        if(b&1){
            ans = PowerMulti(ans, a, n);
        }
        b >>= 1;
        a = PowerMulti(a, a, n);
    }
    return ans;
}
int Miller_Rabin(long long n){
    long long u = n-1, x, y, a, k;
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX