OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..39
EXAMPLE
12563307224 is a term since both 12563307224 and 12563307225 are in A377386: 12563307224/A034968(12563307224) = 369509036, 369509036/A034968(369509036) = 9723922 and 9723922/A034968(9723922) = 373997 are integers, and 12563307225/A034968(12563307225) = 358951635, 358951635/A034968(358951635) = 7976703 and 7976703/A034968(7976703) = 257313 are integers.
PROG
(PARI) fdigsum(n) = {my(k = n, m = 2, r, s = 0); while([k, r] = divrem(k, m); k != 0 || r != 0, s += r; m++); s; }
is1(k) = {my(f = fdigsum(k), f2, m); if(k % f, return(0)); m = k/f; f2 = fdigsum(m); !(m % f2) && !((m/f2) % fdigsum(m/f2)); }
lista(kmax) = {my(q1 = is1(1), q2); for(k = 2, kmax, q2 = is1(k); if(q1 && q2, print1(k-1, ", ")); q1 = q2); }
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Amiram Eldar, Oct 29 2024
STATUS
approved