OFFSET
1,2
COMMENTS
As 1" = 0 and p" = 0 when p is prime, 1 and every prime are terms, hence A008578 is a subsequence. - Bernard Schott, Oct 12 2021
EXAMPLE
8 is in the sequence since 8" = 16 and 8 | 16.
MAPLE
d:= n-> n*add(i[2]/i[1], i=ifactors(n)[2]):
q:= n-> is(irem(d(d(n)), n)=0):
select(q, [$1..250])[]; # Alois P. Heinz, Oct 15 2021
PROG
(PARI) ad(n) = vecsum([n/f[1]*f[2]|f<-factor(n+!n)~]); \\ A003415
isok(k) = !(ad(ad(k)) % k); \\ Michel Marcus, Oct 10 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Oct 09 2021
STATUS
approved