A122004 Primes p that divide A123373(p-1).
2, 3, 5, 23, 29
Offset: 1
Examples
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.
Programs
-
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