cp's OEIS Frontend

A098845 Numbers k such that 4^k - 2^k - 1 is prime.

%I A098845 #64 Apr 03 2023 10:36:10
%S A098845 2,4,5,9,10,18,38,45,50,57,108,161,208,224,225,240,354,597,634,1008,
%T A098845 1080,1468,1525,1560,3298,3329,3846,4129,5430,8616,11834,12988,14610,
%U A098845 43401,45306,53776,54449,67497,74025,122449,136845,142896,164541,171157,187668,274054,316944,349296
%N A098845 Numbers k such that 4^k - 2^k - 1 is prime.
%C A098845 All primes certified using PFGW from primeform group. - _Pierre CAMI_, Mar 07 2005
%C A098845 No terms 2, 3, 7, 12, 13 or 15 (mod 20) except 2. - _Robert Israel_, Dec 08 2015, updated by _Fabrice Lavier_, Jan 10 2019
%C A098845 Using such "Goldilocks" primes (a term coined by Mike Hamburg) as modulus facilitates use of Karatsuba multiplication in elliptic-curve cryptography. - _Francois R. Grieu_, Mar 25 2021
%H A098845 Chris Caldwell, <a href="https://t5k.org/primes/page.php?id=80257">The largest known primes</a>
%H A098845 Mike Hamburg, <a href="https://eprint.iacr.org/2015/625">Ed448-Goldilocks, a new elliptic curve</a>, Cryptology ePrint Archive, Report 2015/625.
%p A098845 select(t -> isprime(4^t-2^t-1), [$1..1000]); # _Robert Israel_, Dec 08 2015
%t A098845 Select[Range[15000], PrimeQ[4^# - 2^# - 1] &] (* _Vincenzo Librandi_, Dec 08 2015 *)
%o A098845 (Magma) [n: n in [0..1000] | IsPrime(2^n*(2^n-1)-1)]; // _Vincenzo Librandi_, Dec 08 2015
%o A098845 (PARI) for(n=1, 1e3, if(ispseudoprime(4^n-2^n-1), print1(n, ", "))) \\ _Altug Alkan_, Dec 08 2015
%o A098845 (Python)
%o A098845 from sympy import isprime
%o A098845 for n in range(1,1000):
%o A098845     if isprime(4**n-2**n-1):
%o A098845         print(n, end=', ') # _Stefano Spezia_, Jan 11 2019
%Y A098845 Cf. similar sequences listed in A265481.
%K A098845 nonn
%O A098845 1,1
%A A098845 _Pierre CAMI_, Oct 10 2004; extended several times: Jun 01 2005, Jun 19 2006, May 03 2007
%E A098845 Extended to a(44) = 349296 (2^698592 - 2^349296 - 1 is a 210298-digit certified prime) by _Pierre CAMI_, Jan 11 2009
%E A098845 Definition simplified by _Pierre CAMI_, May 10 2012
%E A098845 a(30) corrected by _Robert Israel_, Dec 14 2015
%E A098845 4 missing terms between a(41) = 136845 and what is now a(46) = 274054 added by _Fabrice Lavier_, Jan 10 2019