cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

User: Andrew V. Sutherland

Andrew V. Sutherland's wiki page.

Andrew V. Sutherland has authored 27 sequences. Here are the ten most recent ones:

A158494 Boundary area of the T-square fractal.

Original entry on oeis.org

4, 24, 80, 248, 768, 2360, 7200, 21848, 66048, 199160, 599520, 1802648, 5416128, 16264760, 48827040, 146546648, 439771008, 1319575160, 3959249760, 11878797848, 35638490688, 106919666360, 320767387680, 962318940248, 2886990375168, 8661038234360
Offset: 1

Author

Andrew V. Sutherland, Mar 20 2009

Keywords

Comments

Consider the n-th iteration of the T-square fractal (as defined in the links) drawn on an integer lattice scaled so that the shortest edge on the boundary of the fractal has unit length a(n)gives the number of lattice squares in the unshaded region that are adjacent to a square in the shaded region. For n=1 there is a single shaded square and a(1) counts the 4 adjacent unshaded squares. Proposed by Simone Severini.

Crossrefs

Cf. A000392.

Programs

  • Mathematica
    CoefficientList[Series[4*(1 - 5*x^2 + 2*x^3 + 4*x^4)/((1 - x)*(1 - 2*x)*(1 - 3*x)), {x, 0, 30}], x] (* Wesley Ivan Hurt, Jan 20 2017 *)
  • PARI
    a(n)=4*((n==1)+(n==2)*6+(n>=3)*(1-2^(n-1)+23*3^(n-3))) \\ Jaume Oliver Lafont, Mar 22 2009
    
  • PARI
    Vec(4*x*(1-5*x^2+2*x^3+4*x^4) / ((1-x)*(1-2*x)*(1-3*x)) + O(x^30)) \\ Colin Barker, May 22 2017

Formula

a(1)=4, a(2)=24, a(3)=80; for n>3, a(n) = 3*a(n-1) + 2^n - 8.
G.f.: 4*x*(1 - 5*x^2 + 2*x^3 + 4*x^4) / ((1 - x)*(1 - 2*x)*(1 - 3*x)). - Jaume Oliver Lafont, Mar 21 2009
From Colin Barker, May 22 2017: (Start)
a(n) = 4 - 2^(n+1) + 92*3^(n-3) for n>2.
a(n) = 6*a(n-1) - 11*a(n-2) + 6*a(n-3) for n>5. (End)

Extensions

Edited by Charles R Greathouse IV, Oct 28 2009

A146879 Minimal degree of X_1(n).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 2, 2, 2, 4, 2, 5, 3, 4, 4, 7, 4, 5, 6, 6, 6, 11, 6, 12, 8, 10, 10, 12, 8, 18, 12, 14, 12
Offset: 1

Author

Andrew V. Sutherland, Nov 03 2008

Keywords

Comments

a(n) is the least d>0 for which there exists a plane curve f(x,y)=0 of degree d in x or y which is birationally equivalent to the modular curve X_1(n). There exist infinitely many non-isomorphic elliptic curves defined over number fields of degree a(n) which contain a point of order n. a(n)=1 if and only if X_1(n) has genus 0 and these values of n represent the possible finite orders of a point on an elliptic curve over Q.
By Mazur's theorem, these are 1,2,3,4,5,6,7,8,9,10 and 12. a(n)=2 if and only if X_1(n) is elliptic or hyperelliptic, which occurs only for n=11,13,14,15,16 and 18 [Mestre 1981]. The lower bound a(17)>3 follows from [Parent 1999] and the upper bound a(17)<=4 appears (for example) in [Reichert 1986]. a(20)=3 since it cannot be 1 or 2 and an explicit example of degree 3 is known (see below).
From [Jeon-Kim-Schweizer 2006] it follows that this is the only case when a(n)=3. The results a(21)=4 and a(22)=4 then follow from explicit examples [Sutherland 2008]. a(24) is either 4 or 5 and a(n) is not 4 for any n other than 17, 21, 22, or 24 by the results of [Jeon-Kim-Park 2006]. a(23) must be 5, 6, or 7. See [Sutherland 2008] for these and other upper bounds for n <= 50.
For n = 23 to 40, a(n) has been computed by M. Derickx and M. van Hoeij. For n = 41 to 100, upper bounds for a(n) have been computed by M. van Hoeij (see link). - Mark van Hoeij, Apr 17 2012

Examples

			a(20)<=3 because y^3+(x^2+3)y^2+(x^3+4)y+2=0 is an explicit plane model for X_1(20) and a(20)=3 because it is not 1 or 2 (these are all known).
		

Crossrefs

Cf. A029937.

A138552 Returning walks of length 2n on the upper half of the square lattice, distinct under reflections about the y-axis.

Original entry on oeis.org

1, 2, 11, 90, 889, 9723, 113322, 1380522, 17382365, 224573349, 2962117366, 39741658047, 540862505806, 7450655906450, 103713126384420, 1456845308244810, 20627719676855685, 294136002612344145
Offset: 0

Author

Andrew V. Sutherland, Mar 24 2008

Keywords

Comments

Under reasonable assumptions, a(n)=E[X^{2n}] where the random variable X is the unitarized Frobenius trace X=a_p/sqrt(p) (as p varies) of a genus 2 curve whose Jacobian is isogenous to the product of two elliptic curves, exactly one of which has complex multiplication.

Examples

			a(2) = 11 because EEWW, EWEW, EWWE, EWNS, ENSW, ENWS, NEWS, NESW, NSEW, NSNS, NNSS are all the walks of length 4 on the upper half of the square lattice that are distinct under reflections about the y-axis.
		

Crossrefs

Programs

  • Mathematica
    CoefficientList[Series[(3 Pi-2 Pi Sqrt[1-4x]-2EllipticE[16 x])/(8Pi x), {x, 0, 20}], x] (* Benedict W. J. Irwin, Jul 13 2016 *)

Formula

a(n) = (A000891(n) + A000108(n))/2.
G.f.: (3*Pi-2*Pi*sqrt(1-4*x)-2*EllipticE(16*x))/(8*Pi*x). - Benedict W. J. Irwin, Jul 13 2016
a(n) ~ 16^n*n^(-2)/Pi. - Ilya Gutkovskiy, Jul 13 2016
Recurrence: n*(n+1)^2*(3*n - 2)*a(n) = 2*n*(2*n - 1)*(15*n^2 - n - 4)*a(n-1) - 8*(2*n - 3)*(2*n - 1)^2*(3*n + 1)*a(n-2). - Vaclav Kotesovec, Jul 14 2016

A138548 Central moment sequence of tr(A^6) in USp(6).

Original entry on oeis.org

1, 0, 5, 1, 63, 46, 1135, 1800, 25431, 66232, 666387, 2397605, 19650565, 87187842, 633498229, 3214996309, 21829972815, 120665223560, 790528831099, 4613644505799, 29715748525937, 179604102525370, 1149406514424945
Offset: 0

Author

Andrew V. Sutherland, Mar 24 2008

Keywords

Comments

If A is a random matrix in the compact group USp(6) (6x6 complex matrices that are unitary and symplectic), then a(n)=E[(tr(A^6)+1)^n] is the n-th central moment of the trace of A^6, since E[tr(A^6)] = -1 (see A138546).

Examples

			a(5) = 46 because E[(tr(A^6)+1)^5] = 46 for a random matrix A in USp(6).
		

Crossrefs

Formula

mgf is A(z)=e^zF(z) where F(z) is the mgf of A138546.

A138545 Central moment sequence of tr(A^4) in USp(6).

Original entry on oeis.org

1, 0, 3, 1, 27, 26, 385, 708, 7231, 20296, 164277, 608565, 4286161, 19021302, 123867107, 617758729, 3862576095, 20774382552, 127548675709, 720773229015, 4401180707397, 25709943020830, 157204921750191, 939751281408962
Offset: 0

Author

Andrew V. Sutherland, Mar 24 2008

Keywords

Comments

If A is a random matrix in the compact group USp(6) (6x6 complex matrices which are unitary and symplectic), then a(n)=E[(tr(A^4+1))^n] is the n-th central moment of the trace of A^4, since E[tr(A^4)] = -1 (see A138544).

Examples

			a(5) = 26 because E[(tr(A^4)+1)^5] = 26 for a random matrix A in USp(6).
		

Crossrefs

Cf. A138544.

Formula

mgf is A(z)=e^zF(z) where F(z) is the mgf of A138544.

A138356 Moment sequence of t^2 coefficient in det(tI-A) for random matrix A in USp(4).

Original entry on oeis.org

1, 1, 2, 4, 10, 27, 82, 268, 940, 3476, 13448, 53968, 223412, 949535, 4128594, 18310972, 82645012, 378851428, 1760998280, 8288679056, 39457907128, 189784872428, 921472827272, 4512940614960, 22279014978544, 110797225212112
Offset: 0

Author

Andrew V. Sutherland, Mar 17 2008

Keywords

Comments

Let the random variable X be the coefficient of t^2 in the characteristic polynomial det(tI-A) of a random matrix in USp(4) (4 X 4 complex matrices that are unitary and symplectic). Then a(n) = E[X^n].
Let L_p(T) be the L-polynomial (numerator of the zeta function) of a genus 2 curve C. Under a generalized Sato-Tate conjecture, for almost all C,
a(n) is the n-th moment of the coefficient of t^2 in L_p(t/sqrt(p)), as p varies.
See A095922 for central moments.

Examples

			a(3) = 4 because E[X^3] = 4 for X the t^2 coeff of det(tI-A) in USp(4).
a(3) = 1*2^3*(1*1-0^2) + 3*2^2*(0*0-1^2) + 3*2^1*(1*2-0^2) + 1*2^0*(0*0-2^2) = 8 - 12 + 12 - 4 = 4.
		

Crossrefs

Formula

a(n) = (1/2)Integral_{x=0..Pi,y=0..Pi}(4cos(x)cos(y)+2)^n(2cos(x)-2cos(y))^2(2/Pi*sin^2(x))(2/Pi*sin^2(y))dxdy.
a(n) = Sum_{i=0..n}binomial(n,i)2^{n-i}*(A126120(i)*A126120(i+2)-A126120(i+1)^2).

A138549 Moment sequence of t^2 coefficient in det(tI-A) for random matrix A in USp(6).

Original entry on oeis.org

1, 1, 2, 5, 16, 62, 282, 1459, 8375, 52323, 350676, 2493846, 18659787, 145918295, 1186129168, 9978055080, 86545684565, 771571356565, 7051538798490, 65913863945775, 628919704903746, 6114899366942556, 60492393411513722
Offset: 0

Author

Andrew V. Sutherland, Mar 24 2008

Keywords

Comments

Let the random variable X be the coefficient of t^2 in the characteristic polynomial det(tI-A) of a random matrix in USp(6) (6x6 complex matrices that are unitary and symplectic). Then a(n) = E[X^n].
Let L_p(T) be the L-polynomial (numerator of the zeta function) of a genus 3 curve C. Under a generalized Sato-Tate conjecture, for almost all C, a(n) is the n-th moment of the coefficient of t^2 in L_p(t/sqrt(p)), as p varies.
See A138550 for central moments.

Examples

			a(3) = 5 because E[X^3] = 5 for X the t^2 coeff of det(tI-A) in USp(6).
		

Crossrefs

Formula

See Prop. 12 of first Kedlaya-Sutherland reference.

A138551 Moment sequence of t^3 coefficient in det(tI-A) for random matrix A in USp(6).

Original entry on oeis.org

1, 0, 2, 0, 23, 0, 684, 0, 34760, 0, 2493096, 0, 228253267, 0, 25091028820, 0, 3179942075960, 0, 451649016238160, 0, 70421753109861592, 0, 11869050034269797984, 0, 2136758627313217104448, 0
Offset: 0

Author

Andrew V. Sutherland, Mar 24 2008

Keywords

Comments

Let the random variable X be the coefficient of t^3 in the characteristic polynomial det(tI-A) of a random matrix in USp(6) (6x6 complex matrices that are unitary and symplectic). Then a(n) = E[X^n].
Let L_p(T) be the L-polynomial (numerator of the zeta function) of a genus 3 curve C. Under a generalized Sato-Tate conjecture, for almost all C, a(n) is the n-th moment of the coefficient of t^3 in L_p(t/sqrt(p)), as p varies.

Examples

			a(4) = 23 because E[X^4] = 23 for X the t^3 coeff of det(tI-A) in USp(6).
		

Crossrefs

Formula

See Prop. 12 of Kedlaya-Sutherland.

A138546 Moment sequence of tr(A^5) in USp(6).

Original entry on oeis.org

1, 0, 4, 0, 42, 0, 660, 0, 12810, 0, 281736, 0, 6727644, 0, 170316432, 0, 4504487130, 0, 123255492360, 0, 3465702008340, 0, 99645553785960, 0, 2918768920720380, 0, 86852063374902000, 0, 2619552500788984200, 0, 79939673971478231760, 0
Offset: 0

Author

Andrew V. Sutherland, Mar 24 2008

Keywords

Comments

If A is a random matrix in the compact group USp(6) (6 X 6 complex matrices that are unitary and symplectic), then a(n) = E[(tr(A^5))^n] is the n-th moment of the trace of A^5.

Examples

			a(4) = 42 because E[(tr(A^5))^4] = 42 for a random matrix A in USp(6).
		

Crossrefs

Cf. A138540.

Formula

mgf is A(z) = det[F_{i+j-2}(z)], 1<=i,j<=3, where F_m(z) = Sum_j binomial(m,j)(B_{(2j-m)/5}(z)-B_{(2j-m+2)/5}(z)) and B_v(z)=0 for non-integer v and otherwise B_v(z)=I_v(2z) with I_v(z) the hyperbolic Bessel function (of the first kind) of order v.

A138547 Moment sequence of tr(A^6) in USp(6).

Original entry on oeis.org

1, -1, 6, -15, 90, -310, 1860, -7455, 44730, -195426, 1172556, -5416026, 32496156, -156061620, 936369720, -4628393055, 27770358330, -140348412490, 842090474940, -4331544836190, 25989269017140, -135614951248140, 813689707488840, -4296741195214650, 25780447171287900
Offset: 0

Author

Andrew V. Sutherland, Mar 24 2008

Keywords

Comments

If A is a random matrix in the compact group USp(6) (6 X 6 complex matrices that are unitary and symplectic), then a(n) = E[(tr(A^6))^n] is the n-th moment of the trace of A^6. See A138547 for central moments.

Examples

			a(3) = -15 because E[(tr(A^6))^3] = -15 for a random matrix A in USp(6).
		

Crossrefs

Formula

mgf is A(z) = det[F_{i+j-2}(z)], 1<=i,j<=3, where F_m(z) = Sum_j binomial(m,j)(B_{(2j-m)/6}(z)-B_{(2j-m+2)/6}(z)) and B_v(z)=0 for non-integer v and otherwise B_v(z)=I_v(2z) with I_v(z) is the hyperbolic Bessel function (of the first kind) of order v.