Wagstaff primes: Difference between revisions

Content added Content deleted
Line 88: Line 88:
def wagstaff(N):
def wagstaff(N):
""" find first N Wagstaff primes """
""" find first N Wagstaff primes """
pri, wcount = 2, 0
pri, wcount = 1, 0
while wcount < N:
while wcount < N:
pri +=1
pri += 2
if isprime(pri):
if isprime(pri):
wag = (2**pri + 1) // 3
wag = (2**pri + 1) // 3