cp's OEIS Frontend

A113063 Associated with theta series of hexagonal net with respect to a node.

%I A113063 #34 Jul 25 2025 15:37:47
%S A113063 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,
%T A113063 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,
%U A113063 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
%N A113063 Associated with theta series of hexagonal net with respect to a node.
%C A113063 Denoted by |lambda(n)| on page 4 (1.7) in Kassel and Reutenauer arXiv:1610.07793. - _Michael Somos_, Jun 04 2015
%H A113063 Antti Karttunen, <a href="/A113063/b113063.txt">Table of n, a(n) for n = 1..16384</a>
%H A113063 Christian Kassel and Christophe Reutenauer, <a href="https://arxiv.org/abs/1505.07229v3">The zeta function of the Hilbert scheme of n points on a two-dimensional torus</a>, arXiv:1505.07229v3 [math.AG], 2015. [Note that a later version of this paper has a different title and different contents, and the number-theoretical part of the paper was moved to the publication which is next in this list.]
%H A113063 Christian Kassel and Christophe Reutenauer, <a href="https://arxiv.org/abs/1610.07793">Complete determination of the zeta function of the Hilbert scheme of n points on a two-dimensional torus</a>, arXiv:1610.07793 [math.NT], 2016.
%H A113063 Christian Kassel and Christophe Reutenauer, <a href="https://arxiv.org/abs/2507.15780">Pairs of intertwined integer sequences</a>, arXiv:2507.15780 [math.NT], 2025. See p. 5.
%F A113063 Moebius transform is period 9 sequence [ 1, -1, 1, 1, -1, -1, 1, -1, 0, ...].
%F A113063 a(n) is multiplicative with 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).
%F A113063 a(3*n + 2) = 0. a(3*n + 1) = A033687(n), a(3*n) = 2 * A002324(n).
%F A113063 3 * a(n) = A113062(n) unless n=0.
%F A113063 G.f.: Sum_{k>0} f(x^k) + f(x^(3*k)) where f(x) := x / (1 + x + x^2). - _Michael Somos_, Jun 04 2015
%F A113063 a(n) = |A123477(n)|. - _Michael Somos_, Dec 10 2017
%F A113063 Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 4*Pi/(9*sqrt(3)) = 0.806133... (A121839 - 1). - _Amiram Eldar_, Dec 28 2023
%e A113063 G.f. = x + 2*x^3 + x^4 + 2*x^7 + 2*x^9 + 2*x^12 + 2*x^13 + x^16 + 2*x^19 + ...
%t A113063 a[ n_] := If[ n < 1, 0, DivisorSum[ n, {1, -1, 1, 1, -1, -1, 1, -1, 0} [[Mod[#, 9, 1]]] &]]; (* _Michael Somos_, Jun 04 2015 *)
%t A113063 f[p_, e_] := If[Mod[p, 6] == 1, e+1, (1+(-1)^e)/2]; f[3, e_] := 2; a[n_] := Times @@ f @@@ FactorInteger[n]; a[1] = 1; Array[a, 100] (* _Amiram Eldar_, Sep 05 2023 *)
%o A113063 (PARI) {a(n) = if( n<1, 0, sumdiv(n, d, [0, 1, -1, 1, 1, -1, -1, 1, -1][d%9 + 1]))};
%o A113063 (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))))};
%Y A113063 Cf. A002324, A033687, A113062, A121839, A123477.
%K A113063 nonn,easy,mult
%O A113063 1,3
%A A113063 _Michael Somos_, Oct 13 2005