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-2 of 2 results.

A144018 Triangle T(n,k), n >= 1, 1 <= k <= n, read by rows, where sequence a_k of column k has a_k(0)=0, followed by (k+1)-fold 1 and a_k(n) shifts k places left under Euler transform.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 4, 2, 1, 1, 9, 3, 2, 1, 1, 20, 6, 3, 2, 1, 1, 48, 10, 5, 3, 2, 1, 1, 115, 20, 8, 5, 3, 2, 1, 1, 286, 36, 14, 7, 5, 3, 2, 1, 1, 719, 72, 23, 12, 7, 5, 3, 2, 1, 1, 1842, 137, 40, 18, 11, 7, 5, 3, 2, 1, 1, 4766, 275, 69, 30, 16, 11, 7, 5, 3, 2, 1, 1, 12486, 541, 121, 47, 25, 15, 11, 7, 5, 3, 2, 1, 1
Offset: 1

Views

Author

Alois P. Heinz, Sep 07 2008

Keywords

Examples

			T(5,1) = ([1,2,4]*[1,1,4] + [1]*[1]*4 + [1,2]*[1,1]*2 + [1,3]*[1,2]*1)/4 = 36/4 = 9.
Triangle begins:
    1;
    1,  1;
    2,  1,  1;
    4,  2,  1,  1;
    9,  3,  2,  1, 1;
   20,  6,  3,  2, 1, 1;
   48, 10,  5,  3, 2, 1, 1;
  115, 20,  8,  5, 3, 2, 1, 1;
  286, 36, 14,  7, 5, 3, 2, 1, 1;
  719, 72, 23, 12, 7, 5, 3, 2, 1, 1;
		

Crossrefs

T(2n,n) gives A000041(n).
Cf. A316074.

Programs

  • Maple
    etrk:= proc(p) proc(n, k) option remember; `if`(n=0, 1,
             add(add(d*p(d, k), d=numtheory[divisors](j))*
             procname(n-j, k), j=1..n)/n)
           end end:
    B:= etrk(T):
    T:= (n, k)-> `if`(n<=k, `if`(n=0, 0, 1), B(n-k, k)):
    seq(seq(T(n, k), k=1..n), n=1..14);
  • Mathematica
    etrk[p_] := Module[{f}, f[n_, k_] := f[n, k] = If[n == 0, 1, (Sum[Sum[d*p[d, k], {d, Divisors[j]}]*f[n-j, k], {j, 1, n-1}] + Sum[d*p[d, k], {d, Divisors[n]}])/n]; f]; b = etrk[t]; t[n_, k_] := If[n <= k, If[n == 0, 0, 1], b[n-k, k]]; Table[t[n, k], {n, 1, 13}, {k, 1, n}] // Flatten (* Jean-François Alcover, Aug 01 2013, after Alois P. Heinz *)

A316082 Sequence shifts left ten places under Weigh transform with a(n) = signum(n) for n<10.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 3, 4, 5, 6, 8, 10, 12, 15, 19, 24, 31, 40, 51, 67, 87, 114, 149, 195, 255, 334, 438, 574, 753, 989, 1300, 1711, 2255, 2973, 3927, 5190, 6867, 9092, 12046, 15971, 21188, 28124, 37351, 49631, 65981, 87763, 116794
Offset: 0

Views

Author

Alois P. Heinz, Jun 23 2018

Keywords

Crossrefs

Column k=10 of A316074.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          add(binomial(a(i), j)*b(n-i*j, i-1), j=0..n/i)))
        end:
    a:= n-> (k-> `if`(n
    				

Formula

G.f.: x + x^2 + x^3 + x^4 + x^5 + x^6 + x^7 + x^8 + x^9 + x^10 * Product_{n>=1} (1 + x^n)^a(n). - Ilya Gutkovskiy, May 09 2019
Showing 1-2 of 2 results.