OFFSET
1,132
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..65537
FORMULA
A133943(n) = -a(n), for n >= 2.
MAPLE
A133944 := proc(n) local divs, k, i, a ; divs := convert(numtheory[divisors](n), list) ; a := 0 ; for i from 1 to nops(divs) do k := op(i, divs) ; if k-1 in divs or k+1 in divs then a := a+numtheory[mobius](k) ; fi ; od: RETURN(a) ; end: seq(A133944(n), n=1..120) ; # R. J. Mathar, Oct 21 2007
PROG
(PARI) A133944(n) = sumdiv(n, d, (!if((1==d), (n%2), (n%(d-1))&&(n%(d+1))))*moebius(d)); \\ Antti Karttunen, Sep 02 2018
CROSSREFS
KEYWORD
sign
AUTHOR
Leroy Quet, Sep 30 2007
EXTENSIONS
More terms from R. J. Mathar, Oct 21 2007
STATUS
approved