Я прочитал немало описаний в Интернете.
У меня есть идеальный рабочий тест мощности и я думаю, что моя функция get_r() верна.
Я не уверен, как выполнить эту часть алгоритма:
Код: Выделить всё
For a = 1 to square-root(totient(r) * log(n)):
if (X+a)^n != X^n+a (mod X^r − 1,n), output composite
Ниже приведены ссылки на программу, которую я написал для реализации мельника. -тест Рабина и мой (незаконченный) код aks.
Если кто-нибудь сможет объяснить математические расчеты или дать мне немного псевдокода, со мной все будет в порядке.
спасибо
aks.py
miller.py
Подробнее здесь: https://stackoverflow.com/questions/140 ... -in-python