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.

Previous Showing 41-49 of 49 results.

A109447 Binomial coefficients C(n,k) with n-k odd, read by rows.

Original entry on oeis.org

1, 2, 1, 3, 4, 4, 1, 10, 5, 6, 20, 6, 1, 21, 35, 7, 8, 56, 56, 8, 1, 36, 126, 84, 9, 10, 120, 252, 120, 10, 1, 55, 330, 462, 165, 11, 12, 220, 792, 792, 220, 12, 1, 78, 715, 1716, 1287, 286, 13, 14, 364, 2002, 3432, 2002, 364, 14, 1, 105, 1365, 5005, 6435, 3003, 455, 15
Offset: 1

Views

Author

Philippe Deléham, Aug 27 2005

Keywords

Comments

The same as A119900 without 0's. A reflected version of A034867 or A202064. - Alois P. Heinz, Feb 07 2014
From Vladimir Shevelev, Feb 07 2014: (Start)
Also table of coefficients of polynomials P_1(x)=1, P_2(x)=2, for n>=2, P_(n+1)(x) = 2*P_n(x)+(x-1)* P_(n-1)(x). The polynomials P_n(x)/2^(n-1) are connected with sequences A000045 (x=5), A001045 (x=9), A006130 (x=13), A006131 (x=17), A015440 (x=21), A015441 (x=25), A015442 (x=29), A015443 (x=33), A015445 (x=37), A015446 (x=41), A015447 (x=45), A053404 (x=49); also the polynomials P_n(x) are connected with sequences A000129, A002605, A015518, A063727, A085449, A002532, A083099, A015519, A003683, A002534, A083102, A015520. (End)

Examples

			Starred terms in Pascal's triangle (A007318), read by rows:
1;
1*, 1;
1, 2*, 1;
1*, 3, 3*, 1;
1, 4*, 6, 4*, 1;
1*, 5, 10*, 10, 5*, 1;
1, 6*, 15, 20*, 15, 6*, 1;
1*, 7, 21*, 35, 35*, 21, 7*, 1;
1, 8*, 28, 56*, 70, 56*, 28, 8*, 1;
1*, 9, 36*, 84, 126*, 126, 84*, 36, 9*, 1;
Triangle T(n,k) begins:
1;
2;
1,    3;
4,    4;
1,   10,  5;
6,   20,  6;
1,   21,  35,   7;
8,   56,  56,   8;
1,   36, 126,  84,  9;
10, 120, 252, 120, 10;
		

Crossrefs

Cf. A109446.

Programs

  • Maple
    T:= (n, k)-> binomial(n, 2*k+1-irem(n, 2)):
    seq(seq(T(n, k), k=0..ceil((n-2)/2)), n=1..20);  # Alois P. Heinz, Feb 07 2014
  • Mathematica
    Flatten[ Table[ If[ OddQ[n - k], Binomial[n, k], {}], {n, 0, 15}, {k, 0, n}]] (* Robert G. Wilson v *)

Extensions

More terms from Robert G. Wilson v, Aug 30 2005
Corrected offset by Alois P. Heinz, Feb 07 2014

A160444 Expansion of g.f.: x^2*(1 + x - x^2)/(1 - 2*x^2 - 2*x^4).

Original entry on oeis.org

0, 1, 1, 1, 2, 4, 6, 10, 16, 28, 44, 76, 120, 208, 328, 568, 896, 1552, 2448, 4240, 6688, 11584, 18272, 31648, 49920, 86464, 136384, 236224, 372608, 645376, 1017984, 1763200, 2781184, 4817152, 7598336, 13160704, 20759040, 35955712, 56714752
Offset: 1

Views

Author

Willibald Limbrunner (w.limbrunner(AT)gmx.de), May 14 2009

Keywords

Comments

This sequence is the case k=3 of a family of sequences with recurrences a(2*n+1) = a(2*n) + a(2*n-1), a(2*n+2) = k*a(2*n-1) + a(2*n), a(1)=0, a(2)=1. Values of k, for k >= 0, are given by A057979 (k=0), A158780 (k=1), A002965 (k=2), this sequence (k=3). See "Family of sequences for k" link for other connected sequences.
It seems that the ratio of two successive numbers with even, or two successive numbers with odd, indices approaches sqrt(k) for these sequences as n-> infinity.
This algorithm can be found in a historical figure named "Villardsche Figur" of the 13th century. There you can see a geometrical interpretation.

Crossrefs

Programs

  • Magma
    I:=[0,1,1,1]; [n le 4 select I[n] else 2*(Self(n-2) +Self(n-4)): n in [1..40]]; // G. C. Greubel, Feb 18 2023
    
  • Mathematica
    LinearRecurrence[{0,2,0,2}, {0,1,1,1}, 40] (* G. C. Greubel, Feb 18 2023 *)
  • SageMath
    @CachedFunction
    def a(n): # a = A160444
        if (n<5): return ((n+1)//3)
        else: return 2*(a(n-2) + a(n-4))
    [a(n) for n in range(1, 41)] # G. C. Greubel, Feb 18 2023

Formula

a(n) = 2*a(n-2) + 2*a(n-4).
a(2*n+1) = A002605(n).
a(2*n) = A026150(n-1).

Extensions

Edited by R. J. Mathar, May 14 2009

A163305 Numerators of fractions in the approximation of the square root of 5 satisfying: a(n)= (a(n-1)+ c)/(a(n-1)+1); with c=5 and a(1)=0. Also product of the powers of two and five times the Fibonacci numbers.

Original entry on oeis.org

0, 5, 10, 40, 120, 400, 1280, 4160, 13440, 43520, 140800, 455680, 1474560, 4771840, 15441920, 49971200, 161710080, 523304960, 1693450240, 5480120320, 17734041600, 57388564480, 185713295360, 600980848640, 1944814878720
Offset: 1

Views

Author

Mark Dols, Jul 24 2009

Keywords

Comments

For denominators see: A084057 (= product of Lucas numbers (excluding first number (2)), and powers of 2).

Crossrefs

Programs

Formula

From Colin Barker, Jun 20 2012: (Start)
a(n) = 2*a(n-1) + 4*a(n-2).
G.f.: 5*x^2/(1-2*x-4*x^2). (End)
a(n) = 5*A063727(n-1). - R. J. Mathar, Mar 08 2021

Extensions

More terms from Franklin T. Adams-Watters, Aug 06 2009

A171648 a(1) = 1, a(n) = 2*a(n-1) if n is even; a(n) = a(n-1)*Fibonacci((n+1)/2)/Fibonacci((n-1)/2) if n is odd.

Original entry on oeis.org

1, 2, 2, 4, 8, 16, 24, 48, 80, 160, 256, 512, 832, 1664, 2688, 5376, 8704, 17408, 28160, 56320, 91136, 182272, 294912, 589824, 954368, 1908736, 3088384, 6176768, 9994240, 19988480, 32342016, 64684032, 104660992, 209321984, 338690048, 677380096, 1096024064
Offset: 1

Views

Author

Gary W. Adamson, Dec 13 2009

Keywords

Comments

a(n)/a(n-1) apparently tends to phi = A001622 if n=odd; e.g. a(21)/a(20) = 91136/56320 = 1.61818...
a(n)/a(n-2) apparently tends to 1+sqrt(5) = 3.236...= A134945; where a(21)/a(19) = 91136/28160 = 3.23636...
a(1)=1, a(2)=2, a(3)=2, for n>3 a(n)=2*a(n-1) if n is even and a(n)=2*(a(n-1)-a(n-2)+a(n-3)) if n is odd. - Vincenzo Librandi, Dec 06 2010

Examples

			a(8) = 48 = 2*a(7) = 2*24. a(9) = 80 = (5/3)*48 since Fibonacci(5) = 5 and Fibonacci(4) = 3.
		

Crossrefs

Cf. A063727 (bisection), A103435 (bisection).

Programs

  • PARI
    Vec(x*(1+2*x)/(1-2*x^2-4*x^4) + O(x^50)) \\ Colin Barker, Aug 02 2016

Formula

a(1) = 1, a(n) = 2*a(n-1) if n is even; a(n) = a(n-1)*A000045((n+1)/2)/A000045((n-1)/2) if n is odd.
From Colin Barker, Aug 02 2016: (Start)
a(n) = 2*a(n-2) + 4*a(n-4) for n>4.
G.f.: x*(1+2*x) / (1-2*x^2-4*x^4).
(End)

Extensions

Defined "F", removed abundant parentheses, added punctuation to examples, added a factor to the definition, corrected a(13) and added more terms - R. J. Mathar, Dec 15 2009

A206800 Riordan array (1/(1-3*x+x^2), x*(1-x)/(1-3*x+x^2)).

Original entry on oeis.org

1, 3, 1, 8, 5, 1, 21, 19, 7, 1, 55, 65, 34, 9, 1, 144, 210, 141, 53, 11, 1, 377, 654, 534, 257, 76, 13, 1, 987, 1985, 1905, 1111, 421, 103, 15, 1, 2584, 5911, 6512, 4447, 2041, 641, 134, 17, 1, 6765, 17345, 21557, 16837, 9038, 3440, 925, 169, 19, 1
Offset: 0

Views

Author

Philippe Deléham, Feb 12 2012

Keywords

Examples

			Triangle begins :
1
3, 1
8, 5, 1
21, 19, 7, 1
55, 65, 34, 9, 1
144, 210, 141, 53, 11, 1
377, 654, 534, 257, 76, 13, 1
987, 1985, 1905, 1111, 421, 103, 15, 1
2584, 5911, 6512, 4447, 2041, 641, 134, 17, 1
6765, 17345, 21557, 16837, 9038, 3440, 925, 169, 19, 1
Triangle (0,3,-1/3,1/3,0,0,0,0,0,...) DELTA (1,0,-1/3,1/3,0,0,0,0,...) begins :
1
0, 1
0, 3, 1
0, 8, 5, 1
0, 21, 19, 7, 1
0, 55, 65, 34, 9, 1...
		

References

  • Subtriangle of the triangle given by (0, 3, -1/3, 1/3, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, -1/3, 1/3, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938.
  • Antidiagonal sums are A072264(n).

Crossrefs

Formula

T(n,k) = 3*T(n-1,k) + T(n-1,k-1) - T(n-2,k) - T(n-2,k-1).
G.f.: 1/(1-(y+3)*x+(y+1)*x^2).
Sum_{k, 0<=k<=n} T(n,k)*x^k = (-1)^n* A015587(n+1), (-1)^n*A190953(n+1), (-1)^n*A015566(n+1), (-1)*A189800(n+1), (-1)^n*A015541(n+1), (-1)^n*A085939(n+1), (-1)^n*A015523(n+1), (-1)^n*A063727(n), (-1)^n*A006130(n), A077957(n), A000045(n+1), A000079(n), A001906(n+1), A007070(n), A116415(n), A084326(n+1), A190974(n+1), A190978(n+1), A190984(n+1), A190990(n+1), A190872(n) for x = -12, -11, -10, -9, -8, -7, -6, -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8 respectively.

A231729 Triangular array read by rows: row n shows the coefficients of the polynomial u(n) = c(0) + c(1)*x + ... + c(2n-1)*x^(2n-1) which is the denominator of the n-th convergent of the continued fraction [k, k, k, ... ], where k = (x^2 + 1)/(x + 1).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 4, 5, 3, 1, 1, 3, 7, 11, 11, 9, 5, 1, 1, 5, 15, 29, 35, 32, 22, 13, 7, 1, 1, 8, 28, 62, 90, 103, 91, 65, 37, 17, 9, 1, 1, 13, 53, 134, 226, 296, 302, 257, 183, 110, 56, 21, 11, 1, 1, 21, 97, 273, 521, 775, 915, 903, 743, 523, 319, 167
Offset: 1

Views

Author

Clark Kimberling, Nov 13 2013

Keywords

Comments

Sum of numbers in row n: 2*A063727(n). Left edge: A000045 (Fibonacci numbers).

Examples

			First 3 rows:
1 . . . 1
1 . . . 1 . . . 1 . . . 1
2 . . . 4 . . . 5 . . . 3 . . . 1 . . . 1 . . . 1
First 3 polynomials:  1 + x, 1 + x + x^2 + x^3, 2 + 4*x + 5*x^2 + 3*x^3 + x^4 + x^5.
		

Crossrefs

Programs

  • Mathematica
    t[n_] := t[n] = Table[(1 + x^2)/(1 + x), {k, 0, n}];
    b = Table[Factor[Convergents[t[n]]], {n, 0, 10}];
    p[x_, n_] := p[x, n] = Last[Expand[Denominator[b]]][[n]];
    u = Table[p[x, n], {n, 1, 10}]
    v = CoefficientList[u, x]; Flatten[v]

A247584 a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + 3*a(n-5) with a(0) = a(1) = a(2) = a(3) = a(4) = 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 13, 43, 113, 253, 509, 969, 1849, 3719, 8009, 18027, 40897, 91257, 198697, 423777, 894081, 1886011, 4007301, 8594411, 18560081, 40181493, 86872293, 187197193, 402060793, 861827743, 1846685729, 3960390059, 8504658049, 18283290609, 39325827729
Offset: 0

Views

Author

Alexander Samokrutov, Sep 20 2014

Keywords

Comments

a(n)/a(n-1) tends to 2.1486... = 1 + 2^(1/5), the real root of the polynomial x^5 - 5*x^4 + 10*x^3 - 10*x^2 + 5*x - 3.
If x^5 = 2 and n >= 0, then there are unique integers a, b, c, d, g such that (1 + x)^n = a + b*x + c*x^2 + d*x^3 + g*x^4. The coefficient a is a(n) (from A052102). - Alexander Samokrutov, Jul 11 2015
If x=a(n), y=a(n+1), z=a(n+2), s=a(n+3), t=a(n+4) then x, y, z, s, t satisfies Diophantine equation (see link). - Alexander Samokrutov, Jul 11 2015

Crossrefs

Cf. A005531.

Programs

  • Magma
    [n le 5 select 1 else 5*Self(n-1) -10*Self(n-2) +10*Self(n-3) -5*Self(n-4) +3*Self(n-5): n in [1..40]]; // Vincenzo Librandi, Jul 11 2015
    
  • Maple
    m:=50; S:=series( (1-x)^4/(1 -5*x +10*x^2 -10*x^3 +5*x^4 -3*x^5), x, m+1):
    seq(coeff(S, x, j), j=0..m); # G. C. Greubel, Apr 15 2021
  • Mathematica
    LinearRecurrence[{5,-10,10,-5,3}, {1,1,1,1,1}, 50] (* Vincenzo Librandi, Jul 11 2015 *)
  • Maxima
    makelist(sum(2^k*binomial(n,5*k), k, 0, floor(n/5)), n, 0, 50); /* Alexander Samokrutov, Jul 11 2015 */
    
  • PARI
    Vec((1-x)^4/(1-5*x+10*x^2-10*x^3+5*x^4-3*x^5) + O(x^100)) \\ Colin Barker, Sep 22 2014
    
  • Sage
    [sum(2^j*binomial(n, 5*j) for j in (0..n//5)) for n in (0..50)] # G. C. Greubel, Apr 15 2021

Formula

a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + 3*a(n-5).
a(n) = Sum_{k=0...floor(n/5)} (2^k*binomial(n,5*k)). - Alexander Samokrutov, Jul 11 2015
G.f.: (1-x)^4/(1 -5*x +10*x^2 -10*x^3 +5*x^4 -3*x^5). - Colin Barker, Sep 22 2014

A099092 Riordan array (1,2+4x).

Original entry on oeis.org

1, 0, 2, 0, 4, 4, 0, 0, 16, 8, 0, 0, 16, 48, 16, 0, 0, 0, 96, 128, 32, 0, 0, 0, 64, 384, 320, 64, 0, 0, 0, 0, 512, 1280, 768, 128, 0, 0, 0, 0, 256, 2560, 3840, 1792, 256, 0, 0, 0, 0, 0, 2560, 10240, 10752, 4096, 512, 0, 0, 0, 0, 0, 1024, 15360, 35840, 28672, 9216, 1024, 0, 0, 0
Offset: 0

Views

Author

Paul Barry, Sep 25 2004

Keywords

Comments

Row sums are A063727. Diagonal sums are A052907.
The Riordan array (1, s+tx) defines T(n,k) = binomial(k,n-k)*s^k*(t/s)^(n-k). The row sums satisfy a(n) = s*a(n-1) + t*a(n-2) and the diagonal sums satisfy a(n) = s*a(n-2) + t*a(n-3).

Examples

			Rows begin
  {1},
  {0,  2},
  {0,  4,  4},
  {0,  0, 16,  8},
  {0,  0, 16, 48, 16}, ...
		

Crossrefs

Formula

Number triangle T(n,k) = binomial(k, n-k)*2^n; columns have g.f. (2x+4x^2)^k.
T(n,k) = A113953(n,k)*2^k = A026729(n,k)*2^n. - Philippe Deléham, Dec 11 2008

A292847 a(n) is the smallest odd prime of the form ((1 + sqrt(2*n))^k - (1 - sqrt(2*n))^k)/(2*sqrt(2*n)).

Original entry on oeis.org

5, 7, 101, 11, 13, 269, 17, 19, 509, 23, 709, 821, 29, 31, 46957, 55399, 37, 168846239, 41, 43, 9177868096974864412935432937651459122761, 47, 485329129, 2789, 53, 3229, 3461, 59, 61, 1563353111, 139237612541, 67, 5021, 71, 73, 484639, 6221, 79, 6869, 83, 7549
Offset: 1

Views

Author

XU Pingya, Sep 24 2017

Keywords

Examples

			For k = {1, 2, 3, 4, 5}, ((1 + sqrt(6))^k - (1 - sqrt(6))^k)/(2*sqrt(6)) = {1, 2, 9, 28, 101}. 101 is odd prime, so a(3) = 101.
		

Crossrefs

Programs

  • Mathematica
    g[n_, k_] := ((1 + Sqrt[n])^k - (1 - Sqrt[n])^k)/(2Sqrt[n]);
    Table[k = 3; While[! PrimeQ[Expand@g[2n, k]], k++]; Expand@g[2n, k], {n, 41}]
  • PARI
    g(n,k) = ([0,1;2*n-1,2]^k*[0;1])[1,1]
    a(n) = for(k=3,oo,if(ispseudoprime(g(n,k)),return(g(n,k)))) \\ Jason Yuen, Apr 12 2025

Formula

When 2*n + 3 = p is prime, a(n) = p.
Previous Showing 41-49 of 49 results.