A122004 - OEIS

A122004

Primes p that divide A123373(p-1).

2, 3, 5, 23, 29

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,1

LINKS

EXAMPLE

A123373 begins {4, 48, 3598, 924780, 287358579128, ...}.

a(1) = 2 because 2 divides A123373(1) = 4.

a(2) = 3 because 3 divides A123373(2) = 48.

a(3) = 5 because 5 divides A123373(4) = 924780.

PROG

(PARI) is(p) = isprime(p) && sum(j=1, p-1, sum(k=1, p-1, Mod(prime(k), p)^prime(j))) == 0; \\ Jinyuan Wang, Jan 16 2021

(Python)

from sympy import nextprime, prime

p, A122004_list = 2, []

while p < 10**6:

if 0 == sum(pow(prime(i), prime(j), p) for i in range(1, p) for j in range(1, p)) % p:

A122004_list.append(p)

p = nextprime(p) # Chai Wah Wu, Feb 19 2021

CROSSREFS

KEYWORD

nonn,more,hard

AUTHOR

STATUS

approved