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.

A256890 Triangle T(n,k) = t(n-k, k); t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), where f(x) = x + 2.

Original entry on oeis.org

1, 2, 2, 4, 12, 4, 8, 52, 52, 8, 16, 196, 416, 196, 16, 32, 684, 2644, 2644, 684, 32, 64, 2276, 14680, 26440, 14680, 2276, 64, 128, 7340, 74652, 220280, 220280, 74652, 7340, 128, 256, 23172, 357328, 1623964, 2643360, 1623964, 357328, 23172, 256, 512, 72076, 1637860, 10978444, 27227908, 27227908, 10978444, 1637860, 72076, 512
Offset: 0

Views

Author

Dale Gerdemann, Apr 12 2015

Keywords

Comments

Related triangles may be found by varying the function f(x). If f(x) is a linear function, it can be parameterized as f(x) = a*x + b. With different values for a and b, the following triangles are obtained:
a\b 1.......2.......3.......4.......5.......6
The row sums of these, and similarly constructed number triangles, are shown in the following table:
a\b 1.......2.......3.......4.......5.......6.......7.......8.......9
The formula can be further generalized to: t(n,m) = f(m+s)*t(n-1,m) + f(n-s)*t(n,m-1), where f(x) = a*x + b. The following table specifies triangles with nonzero values for s (given after the slash).
a\b 0 1 2 3
-2 A130595/1
-1
0
With the absolute value, f(x) = |x|, one obtains A038221/3, A038234/4, A038247/5, A038260/6, A038273/7, A038286/8, A038299/9 (with value for s after the slash).
If f(x) = A000045(x) (Fibonacci) and s = 1, the result is A010048 (Fibonomial).
In the notation of Carlitz and Scoville, this is the triangle of generalized Eulerian numbers A(r, s | alpha, beta) with alpha = beta = 2. Also the array A(2,1,4) in the notation of Hwang et al. (see page 31). - Peter Bala, Dec 27 2019

Examples

			Array, t(n, k), begins as:
   1,    2,      4,        8,        16,         32,          64, ...;
   2,   12,     52,      196,       684,       2276,        7340, ...;
   4,   52,    416,     2644,     14680,      74652,      357328, ...;
   8,  196,   2644,    26440,    220280,    1623964,    10978444, ...;
  16,  684,  14680,   220280,   2643360,   27227908,   251195000, ...;
  32, 2276,  74652,  1623964,  27227908,  381190712,  4677894984, ...;
  64, 7340, 357328, 10978444, 251195000, 4677894984, 74846319744, ...;
Triangle, T(n, k), begins as:
    1;
    2,     2;
    4,    12,      4;
    8,    52,     52,       8;
   16,   196,    416,     196,      16;
   32,   684,   2644,    2644,     684,      32;
   64,  2276,  14680,   26440,   14680,    2276,     64;
  128,  7340,  74652,  220280,  220280,   74652,   7340,   128;
  256, 23172, 357328, 1623964, 2643360, 1623964, 357328, 23172,   256;
		

Crossrefs

Programs

  • Magma
    A256890:= func< n,k | (&+[(-1)^(k-j)*Binomial(j+3,j)*Binomial(n+4,k-j)*(j+2)^n: j in [0..k]]) >;
    [A256890(n,k): k in [0..n], n in [0..10]]; // G. C. Greubel, Oct 18 2022
    
  • Mathematica
    Table[Sum[(-1)^(k-j)*Binomial[j+3, j] Binomial[n+4, k-j] (j+2)^n, {j,0,k}], {n,0, 9}, {k,0,n}]//Flatten (* Michael De Vlieger, Dec 27 2019 *)
  • PARI
    t(n,m) = if ((n<0) || (m<0), 0, if ((n==0) && (m==0), 1, (m+2)*t(n-1, m) + (n+2)*t(n, m-1)));
    tabl(nn) = {for (n=0, nn, for (k=0, n, print1(t(n-k, k), ", ");); print(););} \\ Michel Marcus, Apr 14 2015
    
  • SageMath
    def A256890(n,k): return sum((-1)^(k-j)*Binomial(j+3,j)*Binomial(n+4,k-j)*(j+2)^n for j in range(k+1))
    flatten([[A256890(n,k) for k in range(n+1)] for n in range(11)]) # G. C. Greubel, Oct 18 2022

Formula

T(n,k) = t(n-k, k); t(0,0) = 1, t(n,m) = 0 if n < 0 or m < 0 else t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), where f(x) = x + 2.
Sum_{k=0..n} T(n, k) = A001715(n).
T(n,k) = Sum_{j = 0..k} (-1)^(k-j)*binomial(j+3,j)*binomial(n+4,k-j)*(j+2)^n. - Peter Bala, Dec 27 2019
Modified rule of Pascal: T(0,0) = 1, T(n,k) = 0 if k < 0 or k > n else T(n,k) = f(n-k) * T(n-1,k-1) + f(k) * T(n-1,k), where f(x) = x + 2. - Georg Fischer, Nov 11 2021
From G. C. Greubel, Oct 18 2022: (Start)
T(n, n-k) = T(n, k).
T(n, 0) = A000079(n). (End)

A257607 Triangle read by rows: T(n,k) = t(n-k, k); t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), where f(x) = x + 5.

Original entry on oeis.org

1, 5, 5, 25, 60, 25, 125, 535, 535, 125, 625, 4210, 7490, 4210, 625, 3125, 30885, 86110, 86110, 30885, 3125, 15625, 216560, 880735, 1377760, 880735, 216560, 15625, 78125, 1471235, 8330745, 18948695, 18948695, 8330745, 1471235, 78125, 390625, 9764910, 74498800, 234897010, 341076510, 234897010, 74498800, 9764910, 390625
Offset: 0

Views

Author

Dale Gerdemann, May 03 2015

Keywords

Examples

			Triangle begins as:
      1;
      5,       5;
     25,      60,      25;
    125,     535,     535,      125;
    625,    4210,    7490,     4210,      625;
   3125,   30885,   86110,    86110,    30885,    3125;
  15625,  216560,  880735,  1377760,   880735,  216560,   15625;
  78125, 1471235, 8330745, 18948695, 18948695, 8330745, 1471235, 78125;
		

Crossrefs

Similar sequences listed in A256890.

Programs

  • Mathematica
    T[n_, k_, a_, b_]:= T[n, k, a, b]= If[k<0 || k>n, 0, If[n==0, 1, (a*(n-k)+b)*T[n-1, k-1, a, b] + (a*k+b)*T[n-1, k, a, b]]];
    Table[T[n,k,1,5], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Mar 24 2022 *)
  • Sage
    def T(n,k,a,b): # A257607
        if (k<0 or k>n): return 0
        elif (n==0): return 1
        else: return  (a*k+b)*T(n-1,k,a,b) + (a*(n-k)+b)*T(n-1,k-1,a,b)
    flatten([[T(n,k,1,5) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Mar 24 2022

Formula

T(n,k) = t(n-k, k); t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), where f(x) = x + 5.
Sum_{k=0..n} T(n, k) = A049198(n).
From G. C. Greubel, Mar 24 2022: (Start)
T(n, k) = (a*k + b)*T(n-1, k) + (a*(n-k) + b)*T(n-1, k-1), with T(n, 0) = 1, a = 1, and b = 5.
T(n, n-k) = T(n, k).
T(n, 0) = A000351(n).
T(n, 1) = 10*6^n - 5^n*(10 + n).
T(n, 2) = 55*7^n - 10*6^n*(n+10) + 5^n*binomial(n+10, 2). (End)

A049199 Nonsquarefree numbers whose Euler totient function is also nonsquarefree.

Original entry on oeis.org

8, 12, 16, 20, 24, 25, 27, 28, 32, 36, 40, 44, 45, 48, 50, 52, 54, 56, 60, 63, 64, 68, 72, 75, 76, 80, 81, 84, 88, 90, 92, 96, 99, 100, 104, 108, 112, 116, 117, 120, 124, 125, 126, 128, 132, 135, 136, 140, 144, 147, 148, 150, 152, 153, 156, 160, 162, 164, 168, 169
Offset: 1

Views

Author

Keywords

Comments

Prime powers with sufficiently large exponents are in this sequence, including 8, 16, 32, 64, ..., 27, 81, ..., 25, 125.

Examples

			124 = 4*31 is in the sequence because phi(124)=60 and neither 124 nor 60 is squarefree.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[170], !SquareFreeQ[#] && !SquareFreeQ[EulerPhi[#]] &] (* Amiram Eldar, Jul 13 2019 *)
  • PARI
    isok(n) = !issquarefree(n) && ! issquarefree(eulerphi(n)); \\ Michel Marcus, Jul 13 2019

Formula

k such that abs(mu(k)) = 0 and abs(mu(phi(k))) = 0.
Showing 1-3 of 3 results.