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.

A081297 Array T(k,n), read by antidiagonals: T(k,n) = ((k+1)^(n+1)-(-k)^(n+1))/(2k+1).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 7, 5, 1, 1, 1, 13, 13, 11, 1, 1, 1, 21, 25, 55, 21, 1, 1, 1, 31, 41, 181, 133, 43, 1, 1, 1, 43, 61, 461, 481, 463, 85, 1, 1, 1, 57, 85, 991, 1281, 2653, 1261, 171, 1, 1, 1, 73, 113, 1891, 2821, 10501, 8425, 4039, 341, 1, 1, 1, 91, 145, 3305
Offset: 0

Views

Author

Paul Barry, Mar 17 2003

Keywords

Comments

Square array of solutions of a family of recurrences.
Rows of the array give solutions to the recurrences a(n)=a(n-1)+k(k-1)a(n-2), a(0)=a(1)=1.
Subarray of array in A072024. - Philippe Deléham, Nov 24 2013

Examples

			Rows begin
  1, 1,  1,  1,   1,    1, ...
  1, 1,  3,  5,  11,   21, ...
  1, 1,  7, 13,  55,  133, ...
  1, 1, 13, 25, 181,  481, ...
  1, 1, 21, 41, 461, 1281, ...
		

Crossrefs

Columns include A002061, A001844, A072025.
Diagonals include A081298, A081299, A081300, A081301, A081302.

Programs

  • Mathematica
    T[n_, k_]:=((n + 1)^(k + 1) - (-n)^(k + 1)) / (2n + 1); Flatten[Table[T[n - k, k], {n, 0, 10}, {k, 0, n}]] (* Indranil Ghosh, Mar 27 2017 *)
  • PARI
    for(k=0, 10, for(n=0, 9, print1(((k+1)^(n+1)-(-k)^(n+1))/(2*k+1), ", "); ); print(); ) \\ Andrew Howroyd, Mar 26 2017
    
  • Python
    def T(n, k): return ((n + 1)**(k + 1) - (-n)**(k + 1)) // (2*n + 1)
    for n in range(11):
        print([T(n - k, k) for k in range(n + 1)]) # Indranil Ghosh, Mar 27 2017

Formula

T(k, n) = ((k+1)^(n+1)-(-k)^(n+1))/(2k+1).
Rows of the array have g.f. 1/((1+kx)(1-(k+1)x)).

Extensions

Name clarified by Andrew Howroyd, Mar 27 2017

A072024 Table by antidiagonals of T(n,k) = ((n+1)^k - (-n)^k)/(2*n+1).

Original entry on oeis.org

0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 3, 1, 1, 0, 1, 5, 7, 1, 1, 0, 1, 11, 13, 13, 1, 1, 0, 1, 21, 55, 25, 21, 1, 1, 0, 1, 43, 133, 181, 41, 31, 1, 1, 0, 1, 85, 463, 481, 461, 61, 43, 1, 1, 0, 1, 171, 1261, 2653, 1281, 991, 85, 57, 1, 1, 0, 1, 341, 4039, 8425, 10501, 2821, 1891, 113, 73, 1, 1, 0
Offset: 0

Views

Author

Henry Bottomley, Jun 06 2002

Keywords

Comments

Rows of the array have g.f. x/((1+k*x)*(1-(k+1)*x)). - Philippe Deléham, Nov 24 2013

Examples

			Rows start:
0 1 1  1   1    1     1      1       1        1 ...
0 1 1  3   5   11    21     43      85      171 ...
0 1 1  7  13   55   133    463    1261     4039 ...
0 1 1 13  25  181   481   2653    8425    40261 ...
0 1 1 21  41  461  1281  10501   36121   246141 ...
0 1 1 31  61  991  2821  32551  117181  1093711 ...
0 1 1 43  85 1891  5461  84883  314245  3879331 ...
0 1 1 57 113 3305  9633 194713  734161 11638089 ...
...
		

Crossrefs

Rows include A057427, A001045, A015441, A053404, A053428, A053430, A065874, etc. Columns include A000004, A000012, A000012, A002061, A001844, A072025, etc.
Cf. A081297.

Programs

  • Magma
    [((k+1)^(n-k) - (-k)^(n-k))/(2*k+1): k in [0..n], n in [0..12]]; // G. C. Greubel, Jan 27 2020
    
  • Maple
    seq(seq( ((k+1)^(n-k) - (-k)^(n-k))/(2*k+1), k=0..n), n=0..12); # G. C. Greubel, Jan 27 2020
  • Mathematica
    T[n_, k_]:= ((n + 1)^k - (-n)^k)/(2n + 1); Flatten[Join[{0}, Table[T[k, n- k], {n, 1, 15}, {k, 0, n}]]] (* Indranil Ghosh, Mar 27 2017 *)
  • PARI
    for(n=0, 10, for(k=0, 9, print1(((n+1)^k-(-n)^k)/(2*n+1), ", "); ); print(); ) \\ Andrew Howroyd, Mar 26 2017
    
  • Sage
    def T(n, k): return ((n+1)^k - (-n)^k)/(2*n+1)
    [[T(k,n-k) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Jan 27 2020

Formula

T(n, k) = T(n, k-1) + n*(n+1)*T(n, k-2) = A060959(A002378(n), k).
T(k, 2n) = (2n+1)*A047969(n, k+1).

A065874 a(n) = (7^(n+1) - (-6)^(n+1))/13.

Original entry on oeis.org

1, 1, 43, 85, 1891, 5461, 84883, 314245, 3879331, 17077621, 180009523, 897269605, 8457669571, 46142992981, 401365114963, 2339370820165, 19196705648611, 117450280095541, 923711917337203, 5856623681349925, 44652524209512451, 290630718826209301, 2166036735625732243
Offset: 0

Views

Author

Len Smiley, Dec 07 2001

Keywords

Comments

A second-order recurrence of promic type (integer roots).
If the number j = A002378(m) is promic (= i(i+1)), then a(n) = a(n-1) + j*a(n-2), a(0) = a(1) = 1 has a closed-form solution involving only powers of integers. The binomial coefficient sum solves the recurrence regardless of promicity (cf. GKP reference).
Hankel transform is := 1,42,0,0,0,0,0,0,0,0,0,0,... - Philippe Deléham, Nov 02 2008

References

  • R. L. Graham, D. E. Knuth, O. Patashnik, "Concrete Mathematics", Addison-Wesley, 1994, p. 204.

Crossrefs

Cf. A001045 (j=2), A015441 (j=6), A053404 (j=12), A053428 (j=20), A053430 (j=30).

Programs

  • Maple
    n->sum(binomial(n-k, k)*(42)^k, k=0..n)
  • Mathematica
    LinearRecurrence[{1,42},{1,1},30] (* Harvey P. Dale, Apr 30 2017 *)
  • PARI
    a(n) = { (7^(n+1) - (-6)^(n+1))/13 } \\ Harry J. Smith, Nov 02 2009

Formula

a(n) = a(n-1) + 42a(n-2); a(0) = a(1) = 1.
G.f.: -1/((6*x+1)*(7*x-1)). - R. J. Mathar, Nov 16 2007
Showing 1-3 of 3 results.