OFFSET
0,2
LINKS
Karl-Heinz Hofmann, Visualization of n = 0 to 4.
PROG
(PARI) select( {is(n, p=4)=ispower(vecsum([d^p|d<-digits(n)])+n, p)}, [0..10^5])
(Python)
aupto = 7869300
A362954 = []
A000583 = set(fp**4 for fp in range(0, int(aupto**(1/4)+3)))
for n in range(0, aupto+1):
print(A362954) # Karl-Heinz Hofmann, Jun 02 2023
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Will Gosnell and M. F. Hasler, May 09 2023
STATUS
approved