A376993 a(n) = A376992(n) - 10^(n-1).

4, 3, 13, 13, 513, 801, 6781, 30721, 40513, 1057513, 515313, 16728501, 78402181, 13617661, 472012281, 64846161, 5481873013, 2459693601, 116852093013, 62611784481, 1234170737761, 1565435686113, 17492477581161, 2254878102513, 16836143444113, 229959946206301

Offset

1,1

Links

Table of n, a(n) for n=1..26.

Mathematica

a[n_]:=(Module[{k=1}, While[!PrimeQ[m=2k^2+2k+1]||IntegerLength[m]<n, k++]; m]-10^(n-1)); Array[a, 15]

Prog

(Python)
from math import isqrt
from itertools import count
from sympy import isprime
def A376993(n):
    for k in count(isqrt(((a:=10**(n-1))<<1)-1>>2)):
        m = 2*k*(k+1)+1
        if m >= a and isprime(m):
            return m-a # Chai Wah Wu, Oct 13 2024

Crossrefs

Cf. A011557, A376933, A376992.

Sequence in context: A120340 A082018 A056477 * A105200 A088933 A019136

Adjacent sequences: A376990 A376991 A376992 * A376994 A376995 A376996

Keyword

nonn,base

Author

Stefano Spezia, Oct 11 2024

Extensions

a(21)-a(26) from Chai Wah Wu, Oct 13 2024

Status

approved