OFFSET
1,4
COMMENTS
Previous name was: Numerators of the third row of the Akiyama-Tanigawa algorithm (or transformation) applied to A001008(n+1)/A002805(n+1).
Successive rows:
3/2, 11/6, 25/12, 137/60, 49/20, 363/140, 761/280, 7129/2520, ...;
1/6, 1/5, 1/5, 4/21, 5/28, 1/6, 7/45, 8/55, 3/22, ...;
-1/30, 0, ...;
-1/30.
First column denominators: 2,3,6,30,30,... = A051717(n+1).
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..2500
FORMULA
a(n) = -a(-n) for all n in Z. - Michael Somos, Aug 01 2017
From Amiram Eldar, Nov 17 2022: (Start)
Multiplicative with a(2) = 1, a(2^e) = 2^e for e > 1, a(3^e) = 3^(e-1), and a(p^e) = p^e for p >= 5.
Sum_{k=1..n} a(k) ~ (49/144) * n^2. (End)
Dirichlet g.f.: zeta(s-1)*(1-1/2^s+2/4^s)*(1-2/3^s). - Amiram Eldar, Jan 05 2023
MATHEMATICA
a[1, n_] := HarmonicNumber[n+1]; a[n_, m_] := a[n, m] = m*(a[n-1, m]-a[n-1, m+1]); Table[a[3, m] // Numerator, {m, 1, 72}] (* Jean-François Alcover, Feb 11 2014 *)
a[ n_] := n / {1, 2, 3, 1, 1, 6, 1, 1, 3, 2, 1, 3}[[Mod[n, 12, 1]]]; (* Michael Somos, Aug 01 2017 *)
PROG
(PARI) {a(n) = if( n<0, -a(-n), numerator( 2*n / ((n+2) * (n+3))))}; /* Michael Somos, Aug 01 2017 */
(Magma) [Numerator(2*n/((n+2)*(n+3))): n in [1..50]]; // G. C. Greubel, Aug 07 2018
CROSSREFS
KEYWORD
nonn,frac,mult
AUTHOR
Paul Curtz, Feb 03 2014
EXTENSIONS
New name using Somos's Pari code from Joerg Arndt, May 27 2018
Keyword:mult added by Andrew Howroyd, Jul 31 2018
STATUS
approved