OFFSET
1,1
COMMENTS
The only known prime terms are 2 and 4793; there are no others < 1.5*10^10.
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..53
Carlos Rivera, Puzzle 1101. Primes p that divide the concatenation of sigma(p-1) and sigma(p+1), The Prime Puzzles & Problems Connection, September 2022.
EXAMPLE
2 is a term since sigma(1) = 1, sigma(3) = 4 and 2 divides 14.
MATHEMATICA
Select[Range[2, 10^6], Divisible[FromDigits[Join @@ IntegerDigits @ DivisorSigma[1, # + {-1, 1}]], #] &] (* Amiram Eldar, Sep 09 2022 *)
PROG
(Python)
from sympy import divisor_sigma
def sigma(n): return divisor_sigma(n, 1)
def ok(n): return n > 1 and int(str(sigma(n-1)) + str(sigma(n+1)))%n == 0
print([k for k in range(10**5) if ok(k)])
(PARI) isok(k) = (eval(concat(Str(sigma(k-1)), Str(sigma(k+1)))) % k) == 0; \\ Michel Marcus, Sep 09 2022
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Michael S. Branicky and Giorgos Kalogeropoulos, Sep 08 2022
STATUS
approved
