OFFSET
1,2
COMMENTS
tau() = A000005().
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MAPLE
with(numtheory); isA141551 := proc(n) local dvs, d; dvs := divisors(n) ; for d in dvs do if not tau(tau(d)) in dvs then RETURN(false): fi; od: RETURN(true); end:
t1:=[]; for n from 1 to 60000 do if isA141551(n) then t1:=[op(t1), n]; fi; od:
MATHEMATICA
aQ[n_] := AllTrue[Divisors[n], Divisible[n, DivisorSigma[0, DivisorSigma[0, #]]] &]; Select[Range[1000], aQ] (* Amiram Eldar, Jul 08 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Sep 12 2008
STATUS
approved