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.

Showing 1-3 of 3 results.

A333345 Decimal expansion of (11 + sqrt(85))/2.

Original entry on oeis.org

1, 0, 1, 0, 9, 7, 7, 2, 2, 2, 8, 6, 4, 6, 4, 4, 3, 6, 5, 5, 0, 0, 1, 1, 3, 7, 1, 4, 0, 8, 8, 1, 3, 9, 6, 5, 7, 8, 6, 2, 3, 4, 0, 2, 5, 2, 4, 3, 6, 1, 2, 3, 2, 0, 0, 4, 0, 0, 3, 8, 7, 6, 1, 0, 2, 7, 2, 1, 3, 3, 5, 5, 1, 3, 4, 0, 0, 9, 3, 7, 7, 3, 0, 3, 8, 3, 9, 4, 7, 0, 4, 5, 3, 9, 6, 6, 4, 0, 2, 8, 2, 4, 7, 0, 1, 6, 9, 9
Offset: 2

Views

Author

Kevin Ryde, Mar 15 2020

Keywords

Comments

This constant is Heuberger and Wagner's lambda. They consider the number of maximum matchings a tree of n vertices may have, and show that the largest number of maximum matchings (A333347) grows as O(lambda^(n/7)) (see A333346 for the 7th root). Lambda is the larger eigenvalue of matrix M = [8,3/5,3] which is raised to a power when counting matchings in a chain of "C" parts in the trees (their lemma 6.2).
Apart from the first digit the same as A176522. - R. J. Mathar, Apr 03 2020

Examples

			10.1097722286...
		

Crossrefs

Sequences growing as this power: A147841, A190872, A333344.
Cf. A333346 (seventh root), A176522.

Programs

  • Mathematica
    With[{$MaxExtraPrecision = 1000}, First@ RealDigits[(11 + Sqrt[85])/2, 10, 105]] (* Michael De Vlieger, Mar 15 2020 *)
  • PARI
    (11 + sqrt(85))/2 \\ Michel Marcus, May 21 2020

Formula

Equals continued fraction [10; 9] = 10 + 1/(9 + 1/(9 + 1/(9 + 1/...))). - Peter Luschny, Mar 15 2020

A333347 Largest number of maximum matchings in a tree of n vertices.

Original entry on oeis.org

1, 1, 1, 2, 3, 4, 5, 8, 11, 15, 21, 30, 41, 56, 81, 112, 153, 216, 303, 418, 571, 819, 1133, 1560, 2187, 3063, 4235, 5832, 8280, 11455, 15807, 22140, 30966, 42823, 59049, 83709, 115808, 160083, 224100, 313059, 432992, 597861, 846279, 1170793, 1618650, 2268000, 3164955
Offset: 0

Views

Author

Kevin Ryde, Mar 15 2020

Keywords

Comments

Heuberger and Wagner consider how many different maximum matchings a tree of n vertices may have. They determine the unique tree (free tree) of n vertices with the largest number of maximum matchings, or at n=6 and n=34 the two trees with equal largest number. a(n) is the largest number of maximum matchings. They show that a(n) grows as O(1.391...^n), where the power is ((11 + sqrt(85))/2)^(1/7) = A333346.
They note an algebraic interpretation too, that a(n) is the largest possible absolute value of the product of the nonzero eigenvalues of the adjacency matrix of a tree of n vertices. This is simply that, in the usual way, a term +- m*x^j in the characteristic polynomial of that matrix means there are m matchings which have j vertices unmatched. The smallest j with a nonzero m is the maximum matchings, and that m is also the product of the nonzero roots.
In Heuberger and Wagner's Sage code, optimal_m(n) is a(n) for the general case tree forms. Their general case symbolic calculations are in terms of lambda = (11 + sqrt(85))/2 = A333345 and its quadratic conjugate lambdabar = (11 - sqrt(85))/2 (called alpha and alphabar in the code). The resulting coefficients give constants c_0 through c_6 in their paper for a(n) -> c_{n mod 7} * lambda^(n/7) (theorem 1.2).
The combinations of powers of lambda and lambdabar occurring are linear recurrences. Recurrence coefficients can be found from a symbolic calculation, or from explicit values and an upper bound on recurrence orders from the patterns of branch lengths and powers. Each case n mod 7 is a recurrence of order up to 44. The simplest is G_k = A190872(k) for n == 1 (mod 7) in the formulas below. Other cases are G variants, and possible additional terms growing slower than G.
The full recurrence for all n is order 574 applying at n=31 onwards (after the last initial exception at n=30). See the links for recurrence coefficients and generating function.

Crossrefs

Cf. A190872, A333345, A333346 (growth power), A333348 (matching number).

Formula

For n == 0 (mod 7) and k = n/7 >= 1, a(n) = 8*A190872(k) - 7*A190872(k-1).
For n == 1 (mod 7) and k = (n-1)/7, a(n) = A190872(k+1). [Heuberger and Wagner theorem 3.3 (1) and lemma 6.2 (2)]
For n == 4 (mod 7) and k = (n-4)/7, a(n) = 3*A333344(k). [Heuberger and Wagner theorem 3.3 (4) and lemma 6.2 (2)]

A147841 a(n) = 11*a(n-1) - 9*a(n-2) with a(0)=1, a(1)=9.

Original entry on oeis.org

1, 9, 90, 909, 9189, 92898, 939177, 9494865, 95990922, 970446357, 9810991629, 99186890706, 1002756873105, 10137643587801, 102489267607866, 1036143151396317, 10475171256888693, 105901595463208770, 1070641008783298233, 10823936737447401633, 109427535032871733866, 1106287454724562457829, 11184314186674341431325
Offset: 0

Views

Author

Philippe Deléham, Nov 14 2008

Keywords

Crossrefs

Cf. A147703, A190872, A333344, A333345 (growth power).

Programs

  • Maple
    A147841:= n-> simplify( 3^n*(ChebyshevU(n, 11/6) - (2/3)*ChebyshevU(n-1, 11/6)) ):
    seq(A147841(n), n=0..25); # G. C. Greubel, May 28 2020
  • Mathematica
    Table[3^n*(ChebyshevU[n, 11/6] - (2/3)*ChebyshevU[n-1, 11/6]), {n,0,25}] (* G. C. Greubel, May 28 2020 *)
    LinearRecurrence[{11,-9},{1,9},30] (* Harvey P. Dale, Feb 28 2023 *)
  • PARI
    a(n) = polcoeff(lift(('x-2)*Mod('x,'x^2-11*'x+9)^n), 1); \\ Kevin Ryde, Apr 11 2020

Formula

a(n) = Sum_{k=0..n} A147703(n,k)*8^k.
G.f.: (1-2*x)/(1 -11*x +9*x^2).
a(n) = 9*A333344(n-1) = A190872(n+1) - 2*A190872(n) = A333344(n) - A190872(n). - Kevin Ryde, Apr 11 2020
a(n) = 3^n*(ChebyshevU(n, 11/6) - (2/3)*ChebyshevU(n-1, 11/6)). - G. C. Greubel, May 28 2020
E.g.f.: exp(11*x/2)*(85*cosh(sqrt(85)*x/2) + 7*sqrt(85)*sinh(sqrt(85)*x/2))/85. - Stefano Spezia, Mar 02 2023

Extensions

Entries corrected by Paolo P. Lava, Nov 18 2008
Terms a(18) onward added by G. C. Greubel, May 28 2020
Showing 1-3 of 3 results.