Lucas-Lehmer 級数を p まで取得するために、以下のコードを書きました。p はメルセンヌ数の指数です。確認したところ、11、23、29 などの一部の素数 p では機能しないことがわかりました。
コードは次のとおりです。
def ll_series (p):
ll_list=[4]
print 4
for i in range(1, p+1):
ll_list.append((ll_list[i-1]**2 - 2) % (2**p-1))
print(ll_list[i])
return ll_list