A365445 a(n) is the index of the least prime that is also a centered n-gonal number, or -1 if none exists.

8, 3, 11, 4, 14, -1, -1, 5, 19, 6, 22, 14, 36, 7, 27, 8, 43, 18, 31, 9, 34, 21, 36, 22, 38, 10, 795, 11, 64, 25, 46, 27, 47, 12, 48, 50, 183, 13, 394, 14, 83, 121, 58, 15, 61, 169, 94, 36, 63, 16, 489, 38, 67, 68, 105, 17, 623, 18, 73, 74, 75, 44, 347

Offset

3,1

Links

Table of n, a(n) for n=3..65.

Ilya Gutkovskiy, Scatterplot of a(n) up to n=10000

Eric Weisstein's World of Mathematics, Centered Polygonal Number

Formula

a(n) = A000720(A365815(n)). - Michel Marcus, Sep 27 2023

Mathematica

Table[PrimePi[SelectFirst[n # (# + 1)/2 + 1 & /@ Range[100], PrimeQ]], {n, 3, 65}] /. (PrimePi[Missing["NotFound"]] -> -1)

Prog

(PARI) a(n) = if ((n==8) || (n==9), return(-1)); my(k=0, p); while (!isprime(p=1+n*k*(k-1)/2), k++); primepi(p); \\ Michel Marcus, Sep 27 2023

Crossrefs

Centered k-gonal primes listed in A276261.

Cf. A000040, A000720, A101321, A365815.

Sequence in context: A069610 A257811 A122159 * A069200 A069218 A248295

Adjacent sequences: A365442 A365443 A365444 * A365446 A365447 A365448

Keyword

sign

Author

Ilya Gutkovskiy, Sep 25 2023

Status

approved