login
A383476
Numbers k such that binomial(2k,k) + k is prime.
0
1, 3, 5, 23, 55, 61, 191, 1933, 3049, 8011, 10589, 58687, 100469
OFFSET
1,2
COMMENTS
Corresponding primes start 3, 23, 257, 8233430727623, 98527218530093856775578873054487.
MATHEMATICA
Select[Range[3500], PrimeQ[Binomial[2*#, #] + #] &] (* Amiram Eldar, Apr 28 2025 *)
PROG
(Magma) [k: k in [1..1000] | IsPrime(Binomial(2*k, k)+k)];
(Python)
from math import comb
from gmpy2 import is_prime
def ok(n): return is_prime(comb(2*n, n) + n)
print([k for k in range(3050) if ok(k)]) # Michael S. Branicky, Apr 28 2025
CROSSREFS
Sequence in context: A148554 A120937 A075307 * A100302 A290384 A023247
KEYWORD
nonn,hard,more
AUTHOR
EXTENSIONS
a(7) inserted and a(10)-a(13) from Michael S. Branicky, Apr 29 2025
STATUS
approved