OFFSET
1,1
COMMENTS
The auxiliary sequence defined by b(n)=sum_{d|n} A002322(d) starts 1,2,3,4,5,6,7,6,9,10,11,10,13,14,11,10,17,18,19,16,...
The auxiliary sequence is A141258. [Reinhard Zumkeller, Feb 17 2012]
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1000
W. D. Banks and F. Luca, On integers with a special divisibility property, Archivum Mathematicum (BRNO) 42 (2006) pp 31-42.
FORMULA
n such that (sum_{d|n} A002322(d)) | n.
MATHEMATICA
Select[ Range[100000], Divisible[#, s = Total[ CarmichaelLambda /@ Divisors[#]]] && s < # &] (* Jean-François Alcover, Jun 24 2013 *)
PROG
(PARI) lambda(p, alpha)={ if(p>=3 || alpha<=2, return(p^(alpha-1)*(p-1)), return(2^(alpha-2)) ; ) ; } A002322(n)={ local(pf, rmax, resul) ; if(n==1, return(1) ) ; pf=factor(n) ; rmax=matsize(pf)[1] ; resul= lambda(pf[1, 1], pf[1, 2]) ; for(r=2, rmax, resul=lcm(resul, lambda(pf[r, 1], pf[r, 2])) ; ) ; return(resul) ; } b(n)={ sumdiv(n, d, A002322(d)) ; } { for(n=1, 120000, l=b(n) ; if( l != 1 && l != n && n%l==0, print1(n, ", ") ) ; ) ; }
CROSSREFS
KEYWORD
nonn
AUTHOR
R. J. Mathar, Jul 29 2007
STATUS
approved