This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A100361 #20 Mar 02 2022 10:44:16 %S A100361 0,1,2,4,6,16,18,54,58,100,120,504,1302,3234,14748,16102,22782,34656, %T A100361 64764,70866,194940,274074,313344,331416,354640 %N A100361 Numbers k such that 2^k - k + 1 is prime. %C A100361 a(21) > 150000. - _Giovanni Resta_, Mar 18 2014 %C A100361 a(26) > 5*10^5. - _Robert Price_, Oct 13 2014 %p A100361 A100361:=n->`if`(isprime(2^n-n+1), n, NULL): seq(A100361(n), n=0..10^3); # _Wesley Ivan Hurt_, Oct 13 2014 %t A100361 {ta={{0}}, tb={{0}}};Do[g=n;s=2^n-n+1; If[PrimeQ[s], Print[n];ta=Append[ta, n]; tb=Append[tb, s]], {n, 1, 10000}];{ta, tb, g} %o A100361 (PARI) is(n)=ispseudoprime(2^n-n+1) \\ _Charles R Greathouse IV_, Feb 20 2017 %Y A100361 Cf. A001580, A052007, A048744, A100357-A100359, A061421. %K A100361 nonn,hard,more %O A100361 1,3 %A A100361 _Labos Elemer_, Nov 19 2004 %E A100361 a(15)-a(20) from _Giovanni Resta_, Mar 18 2014 %E A100361 a(21)-a(25) from _Robert Price_, Oct 13 2014