OFFSET
1,1
EXAMPLE
a(1) = 2 and a(2) = 3 are Niven numbers with derivatives 2' = 1 and 3' = 1 which are Niven numbers, and a(1) + a(2) = 2 + 3 = 5 which is Niven number and 5' = 1 is also a Niven number.
PROG
(Magma) f:=func<n |n le 1 select 0 else n*(&+[Factorisation(n)[i][2] / Factorisation(n)[i][1]: i in [1..#Factorisation(n)]])>; a:=[]; niven:=func<n|n mod &+Intseq(n) eq 0 >; ff:=func<n|niven(n) and niven(Floor(f(n)))>; a:=[2]; for n in [2..58] do k:=3; while k in a or not(ff(k) and ff(k+a[n-1])) do k:=k+1; end while; Append(~a, k); end for; a;
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Marius A. Burtea, Nov 20 2021
STATUS
approved