第七次作业 Millier-Rabin素性测试算法

1. 算法原理
   见教材或PPT
2. 程序输入：
   正奇数n,满足 $n - 1 = 2^st, s >= 1, t为奇数$，检测轮数固定为20；
   （t、s是根据n来求出的，虽然比较明显，这里还是说一下=。=）
3. 程序输出：
   若n是素数，输出YES；反之输出NO;

PS:

1. 实验报告中要有完整的算法流程描述和样例输入输出；
2. 读入输出方式与前几次相同；
3. 算法中用到的求模幂的部分，要求要利用前面试验中的快速幂算法；
4. 每一轮都需要重新选取随机数 $a \in (0, n)$;
5. 请将功能整合到函数，便于之后实验中使用。