login
This site is supported by donations to The OEIS Foundation.

 

Logo

Annual appeal: Please make a donation to keep the OEIS running! Over 6000 articles have referenced us, often saying "we discovered this result with the help of the OEIS".
Other ways to donate

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A047817 Denominators of Hurwitz numbers H_n (coefficients in expansion of Weierstrass P-function). 4
10, 10, 130, 170, 10, 130, 290, 170, 4810, 410, 10, 2210, 530, 290, 7930, 170, 10, 351130, 10, 6970, 3770, 890, 10, 214370, 1010, 530, 524290, 557090, 10, 325130, 10, 170, 130, 1370, 290, 5969210, 1490, 10, 1081730, 6970, 10, 3770 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

Hurwitz showed (see Katz, eqn. 9) that a(n) = product of the prime p = 2 and the primes p of the form 4*k + 1 such that p - 1 divides 4*n. It follows that a(n) is a divisibility sequence, that is, if n | m then a(n) | a(m). - Peter Bala, Jan 08 2014

REFERENCES

F. Lemmermeyer, Reciprocity Laws, Springer-Verlag, 2000; see p. 276.

LINKS

T. D. Noe, Table of n, a(n) for n = 1..1000

L. Carlitz, The coefficients of the lemniscate function, Math. Comp., 16 (1962), 475-478.

A. Hurwitz, Über die Entwicklungskoeffizienten der lemniskatischen Funktionen, Math. Ann., 51 (1899), 196-226; Mathematische Werke. Vols. 1 and 2, Birkhaeuser, Basel, 1962-1963, see Vol. 2, No. LXVII.

A. Hurwitz, Über die Entwicklungskoeffizienten der lemniskatischen Funktionen, Math. Ann., 51 (1899), 196-226; Mathematische Werke. Vols. 1 and 2, Birkhaeuser, Basel, 1962-1963, see Vol. 2, No. LXVII. [Annotated scanned copy]

N. M. Katz, The Congruences of Clausen - von Staudt and Kummer for Bernoulli-Hurwitz Numbers, Mathematische Annalen 216, 1-4 (1975)

Index to divisibility sequences

FORMULA

Let P be the Weierstrass P-function satisfying P'^2 = 4*P^3 - 4*P. Then P(z) = 1/z^2 + Sum_{n=1..infinity} 2^(4n)*H_n*z^(4n-2)/(4n*(4n-2)!).

Sum_{ (r, s) != (0, 0) } 1/(r+si)^(4n) = (2w)^(4n)*H_n/(4n)! where w = 2 * Integral_{0..1} dx/(sqrt(1-x^4)).

See PARI line for recurrence.

EXAMPLE

Hurwitz numbers H_1, H_2, ... = 1/10, 3/10, 567/130, 43659/170, 392931/10, ...

MAPLE

H := proc(n) local k; option remember; if n = 1 then 1/10 else 3*add((4*k - 1)*(4*n - 4*k - 1)*binomial(4*n, 4*k)*H(k)*H(n - k), k = 1 .. n - 1)/( (2*n - 3)*(16*n^2 - 1)) fi; end;

a := n -> denom(H(n));

# Implementation based on Hurwitz's extension of Clausen's theorem:

GenClausen := proc(n) local k, S; map(k->k+1, numtheory[divisors](n));

    S := select(p-> isprime(p) and p mod 4 = 1, %);

    if S <> {} then 2*mul(k, k=S) else NULL fi end:

A047817_list := proc(n) local i; seq(GenClausen(i), i=1..4*n) end;

A047817_list(42); # Peter Luschny, Oct 03 2011

# Implementation based on Weierstrass's P-function:

c := n -> (n*(4*n-2)!/(2^(4*n-2)))*coeff(series(WeierstrassP(z, 4, 0), z, 4*n+2), z, 4*n-2); a := n -> denom(c(n)); seq(a(n), n=1..42); # Peter Luschny, Aug 18 2014

MATHEMATICA

a[1] = 1/10; a[n_] := a[n] = (3/(2*n - 3)/(16*n^2 - 1))* Sum[(4*k - 1)*(4*n - 4*k - 1)*Binomial[4*n, 4*k]*a[k]* a[n - k], {k, 1, n - 1}]; Denominator[ Table[a[n], {n, 1, 42}]] (* Jean-François Alcover, Oct 18 2011, after PARI *)

a[ n_] := If[ n < 1, 0, Denominator[ 2^(-4 n) (4 n)! SeriesCoefficient[ 1 - x WeierstrassZeta[ x, {4, 0}], {x, 0, 4 n}]]]; (* Michael Somos, Mar 05 2015 *)

PROG

(PARI) do(lim)=v=vector(lim); v[1]=1/10; for(n=2, lim, v[n]=3/(2*n-3)/(16*n^2-1)*sum(k=1, n-1, (4*k-1)*(4*n-4*k-1)*binomial(4*n, 4*k)*v[k]*v[n-k])) \\ Henri Cohen, Mar 18 2002

CROSSREFS

For numerators see A002306.

Cf. A160014.

Sequence in context: A285608 A285612 A287600 * A065243 A105750 A220449

Adjacent sequences:  A047814 A047815 A047816 * A047818 A047819 A047820

KEYWORD

nonn,easy,nice,frac

AUTHOR

N. J. A. Sloane

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy .

Last modified November 21 06:31 EST 2017. Contains 294989 sequences.