费马素性检测

题目要求：

关键算法：快速模指数算法，用来求指数很大时，模p的值（p是一个很大的数）
相关代码

import random
def quick_algorithm(a, b, c):  # 求a**b mod c
    a = a % c
    ans = 1
    while b != 0:
        if b & 1:
            ans = (ans * a) % c
        b >>= 1
        a = (a * a) % c
    return ans

def Isprime(p,k):
    for i in range(k):
        a = random.randint(2,p - 2)
        if quick_algorithm(a,p - 1, p) != 1:
            return 0
        return 1

p是要检测的数，k是自己要检测的次数！！
！！本文仅供自己参考！！
（另外quick_algorithm函数借鉴了csdn一位大佬的代码;自己也写了一个，但没这个简单，就不贴上来了）