cp's OEIS Frontend

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.

A215775 Numbers k such that 2*k!! + 1 is a prime.

This page as a plain text file.
%I A215775 #24 Aug 14 2024 01:51:45
%S A215775 0,1,2,3,4,5,6,7,8,10,13,17,18,24,28,30,39,51,149,197,355,444,623,908,
%T A215775 1084,1125,1577,1748,2584,2628,2702,3973,5046,6048,6058,7123,9706,
%U A215775 14223,17775
%N A215775 Numbers k such that 2*k!! + 1 is a prime.
%C A215775 The corresponding primes are in A215777.
%e A215775 17 is in the sequence because 2*17!! + 1 = 2*34459425 + 1 = 68918851 is prime.
%t A215775 lst={};Do[If[PrimeQ[2*n!!+1],AppendTo[lst,n]],{n,0,1500}];lst
%o A215775 (Python)
%o A215775 from sympy import isprime
%o A215775 def afind(startat=0, limit=10**9):
%o A215775   df, dfm1 = 1, 1
%o A215775   for k in range(startat):
%o A215775     df, dfm1 = (k+1)*dfm1, df
%o A215775   for k in range(startat, limit+1):
%o A215775     if isprime(2*df+1): print(k, end=", ")
%o A215775     df, dfm1 = (k+1)*dfm1, df
%o A215775 afind(limit=625) # _Michael S. Branicky_, May 16 2021
%Y A215775 Cf. A006882, A215777.
%K A215775 nonn,hard,more
%O A215775 1,3
%A A215775 _Michel Lagneau_, Aug 23 2012
%E A215775 Offset changed to 1 by and a(29)-a(31) from _Jinyuan Wang_, May 15 2021
%E A215775 a(32)-a(33) from _Michael S. Branicky_, May 16 2021
%E A215775 a(34)-a(39) from _Michael S. Branicky_, Aug 13 2024