login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A123477 Expansion of (1 - b(q)) / 3 in powers of q where b(q) is a cubic AGM theta function. 5
1, 0, -2, 1, 0, 0, 2, 0, -2, 0, 0, -2, 2, 0, 0, 1, 0, 0, 2, 0, -4, 0, 0, 0, 1, 0, -2, 2, 0, 0, 2, 0, 0, 0, 0, -2, 2, 0, -4, 0, 0, 0, 2, 0, 0, 0, 0, -2, 3, 0, 0, 2, 0, 0, 0, 0, -4, 0, 0, 0, 2, 0, -4, 1, 0, 0, 2, 0, 0, 0, 0, 0, 2, 0, -2, 2, 0, 0, 2, 0, -2, 0, 0, -4, 0, 0, 0, 0, 0, 0, 4, 0, -4, 0, 0, 0, 2, 0, 0, 1, 0, 0, 2, 0, 0 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,3
COMMENTS
Cubic AGM theta functions: a(q) (see A004016), b(q) (A005928), c(q) (A005882).
Denoted by lambda(n) on page 4 (1.7) in Kassel and Reutenauer arXiv:1610.07793. - Michael Somos, Dec 10 2017
LINKS
Christian Kassel and Christophe Reutenauer, Complete determination of the zeta function of the Hilbert scheme of n points on a two-dimensional torus, arXiv preprint arXiv:1610.07793 [math.NT], 2016.
FORMULA
Moebius transform is period 9 sequence [1, -1, -3, 1, -1, 3, 1, -1, 0, ...].
a(n) is multiplicative and a(p^e) = -2 if p = 3 and e>0, a(p^e) = e+1 if p == 1 (mod 6), a(p^e) = (1+(-1)^e)/2 if p == 2, 5 (mod 6).
a(3*n + 2) = 0. a(3*n + 1) = A033687(n), a(3*n) = -2*A002324(n).
-3*a(n) = A005928(n) unless n=0. |a(n)| = A113063(n).
EXAMPLE
G.f. = q - 2*q^3 + q^4 + 2*q^7 - 2*q^9 - 2*q^12 + 2*q^13 + q^16 + 2*q^19 + ...
MAPLE
A123477 := proc(n)
local a, pe, p, e;
a := 1;
for pe in ifactors(n)[2] do
p := op(1, pe) ;
e := op(2, pe) ;
if modp(p, 6) = 1 then
a := a*(e+1) ;
elif modp(p, 6) in {2, 5} then
a := a*(1+(-1)^e)/2 ;
elif e > 0 then
a := -2*a ;
end if;
end do:
a ;
end proc:
seq(A123477(n), n=1..100) ; # R. J. Mathar, Feb 22 2021
MATHEMATICA
a[ n_] := If[ n < 1, 0, DivisorSum[ n, {1, -1, -3, 1, -1, 3, 1, -1, 0} [[Mod[#, 9, 1]]] &]]; (* Michael Somos, Dec 10 2017 *)
PROG
(PARI) {a(n) = if( n<1, 0, sumdiv(n, d, [0, 1, -1, -3, 1, -1, 3, 1, -1] [d%9+1]))};
(PARI) {a(n) = my(A, p, e); if( n<1, 0, A = factor(n); prod(k=1, matsize(A)[1], [p, e] = A[k, ]; if( p==3, -2, p%6==1, e+1, !(e%2))))};
CROSSREFS
Sequence in context: A239393 A256637 A113063 * A035225 A298931 A035219
KEYWORD
sign,mult
AUTHOR
Michael Somos, Sep 27 2006
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 22:04 EDT 2024. Contains 371254 sequences. (Running on oeis4.)