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-10 of 14 results. Next

A257506 Inverse permutation to A257505.

Original entry on oeis.org

1, 3, 6, 2, 10, 15, 21, 28, 36, 45, 55, 5, 66, 78, 91, 9, 105, 4, 120, 136, 153, 14, 171, 190, 210, 231, 253, 276, 300, 325, 351, 378, 406, 435, 465, 496, 528, 561, 595, 630, 666, 703, 741, 780, 820, 861, 903, 20, 946, 990, 1035, 27, 1081, 1128, 1176, 1225, 1275, 1326, 1378, 35, 1431, 1485, 1540, 44, 1596, 54, 1653, 1711, 1770, 65, 1830, 8
Offset: 1

Views

Author

Antti Karttunen, May 06 2015

Keywords

Crossrefs

Inverse: A257505.

Programs

  • Scheme
    (define (A257506 n) (let ((row (A257681 n)) (col (A257679 n))) (* (/ 1 2) (- (expt (+ row col) 2) row col col col -2))))

Formula

a(n) = (1/2) * ((c+r)^2 - r - 3*c + 2), where r = A257681(n), and c = A257679(n).

A255411 Shift factorial base representation of n one digit left (with 0 added to right), increment all nonzero digits by one, then convert back to decimal; Numbers with no digit 1 in their factorial base representation.

Original entry on oeis.org

0, 4, 12, 16, 18, 22, 48, 52, 60, 64, 66, 70, 72, 76, 84, 88, 90, 94, 96, 100, 108, 112, 114, 118, 240, 244, 252, 256, 258, 262, 288, 292, 300, 304, 306, 310, 312, 316, 324, 328, 330, 334, 336, 340, 348, 352, 354, 358, 360, 364, 372, 376, 378, 382, 408, 412, 420, 424, 426, 430, 432, 436, 444
Offset: 0

Views

Author

Antti Karttunen, Apr 16 2015

Keywords

Comments

Nonnegative integers such that the number of ones (A257511) in their factorial base representation (A007623) is zero.
Nonnegative integers such that the least missing nonzero digit (A257079) in their factorial base representation is one.
a(n) can be also directly computed from n by "shifting left" its factorial base representation (that is, by appending one zero to the right, see A153880) and then incrementing all nonzero digits by one, and then converting the resulting (still valid) factorial base number back to decimal. See the examples.
The sequences A227130 and A227132 are closed under a(n), in other words, permutation listed as the a(n)-th entry in tables A060117 & A060118 has the same parity as the n-th entry in those same tables.

Examples

			Factorial base representation (A007623) of 1 is "1", shifting it left yields "10", and when we increment all nonzero digits by one, we get "20", which is the factorial base representation of 4 (as 4 = 2*2! + 0*1!), thus a(1) = 4.
F.b.r. of 2 is "10", shifting it left yields "100", and "200" is f.b.r. of 12, thus a(2) = 12.
F.b.r. of 43 is "1301", shifting it left and incrementing all nonzeros by one yields "24020", which is f.b.r of 340, thus a(43) = 340.
		

Crossrefs

Complement: A256450.
Positions of ones in A257079, fixed points of A257080, positions of zeros in A257511, A257081 and A257261.
Cf. also A227130/A227132, A060117/A060118 and also arrays A257503 & A257505.

Programs

  • Mathematica
    factBaseIntDs[n_] := Module[{m, i, len, dList, currDigit}, i = 1; While[n > i!, i++]; m = n; len = i; dList = Table[0, {len}]; Do[currDigit = 0; While[m >= j!, m = m - j!; currDigit++]; dList[[len - j + 1]] = currDigit, {j, i, 1, -1}]; If[dList[[1]] == 0, dList = Drop[dList, 1]]; dList]; s = Table[FromDigits[factBaseIntDs[n]], {n, 500}]; {0}~Join~Flatten@ Position[s, x_ /; DigitCount[x][[1]] == 0](* Michael De Vlieger, Apr 27 2015, after Alonso del Arte at A007623 *)
    Select[Range[0, 444], ! MemberQ[IntegerDigits[#, MixedRadix[Reverse@ Range@ 12]], 1] &] (* Michael De Vlieger, May 30 2016, Version 10.2 *)
    r = MixedRadix[Reverse@Range[2, 12]]; Table[FromDigits[Map[If[# == 0, 0, # + 1] &, IntegerDigits[n, r]]~Join~{0}, r], {n, 0, 60}] (* Michael De Vlieger, Aug 14 2016, Version 10.2 *)
  • Python
    from sympy import factorial as f
    def a007623(n, p=2): return n if n

    0 else '0' for i in x)[::-1] return 0 if n==0 else sum(int(y[i])*f(i + 1) for i in range(len(y))) print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 20 2017

A256450 Numbers that have at least one 1-digit in their factorial base representation (A007623).

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15, 17, 19, 20, 21, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 53, 54, 55, 56, 57, 58, 59, 61, 62, 63, 65, 67, 68, 69, 71, 73, 74, 75, 77, 78, 79, 80, 81, 82, 83, 85, 86, 87, 89, 91, 92, 93, 95, 97, 98, 99, 101
Offset: 0

Views

Author

Antti Karttunen, Apr 27 2015

Keywords

Comments

Numbers n for which A257679(n) = 1, i.e., numbers n such that the least nonzero digit in their factorial base representation (A007623) is 1.
Involution A225901 maps each term of this sequence to a unique term of A273670, and vice versa.
Starting offset is zero (with a(0) = 1) because it is the most natural offset for the given fast recurrence.

Crossrefs

Complement of A255411.
Cf. A257680 (characteristic function), A273662 (left inverse).
First row of A257503, first column of A257505.
Subsequences: A059590 (apart from its zero-term), A255341, A255342, A255343, A257262, A257263, A258198, A258199.
Cf. also A227187 (numbers with at least one nonleading zero) and A273670, A225901.

Programs

  • Mathematica
    Select[Range@ 101, MemberQ[IntegerDigits[#, MixedRadix[Reverse@ Range@ 12]], 1] &] (* Michael De Vlieger, May 30 2016, Version 10.2 *)
    r = MixedRadix[Reverse@ Range[2, 12]]; Select[Range@ 101, Min[IntegerDigits[#, r] /. 0 -> Nothing] == 1 &]  (* Michael De Vlieger, Aug 14 2016, Version 10.2 *)
  • Python
    def A(n, p=2): return n if n

    =1]) # Indranil Ghosh, Jun 19 2017

Formula

a(0) = 1, and for n >= 1, if A257511(1+a(n-1)) > 0, then a(n) = a(n-1) + 1, otherwise a(n-1) + 2. [In particular, if the previous term is 2k, then the next term is 2k+1, because all odd numbers are members.]
Other identities:
For all n >= 0, A273662(a(n)) = n. [A273662 works as the left inverse for this sequence.]
From Antti Karttunen, May 26 2015: (Start)
Alternative recurrence for the same sequence:
Set k = A258198(n), d = n - A258199(n) and f = A000142(k+1) = (k+1)! If d < f then b(n) = f+d, otherwise b(n) = ((2+floor((d-f)/A258199(n))) * f) + b((d-f) mod A258199(n)). For offset=1 sequence, define a(n) = b(n-1).
(End)

Extensions

Starting offset changed from 1 to 0 by Antti Karttunen, May 30 2016

A062119 a(n) = n! * (n-1).

Original entry on oeis.org

0, 2, 12, 72, 480, 3600, 30240, 282240, 2903040, 32659200, 399168000, 5269017600, 74724249600, 1133317785600, 18307441152000, 313841848320000, 5690998849536000, 108840352997376000, 2189611807358976000, 46225138155356160000, 1021818843434188800000
Offset: 1

Views

Author

Olivier Gérard, Jun 13 2001

Keywords

Comments

For n > 0, a(n) = number of permutations of length n+1 that have 2 predetermined elements nonadjacent; e.g., for n=2, the permutations with, say, 1 and 2 nonadjacent are 132 and 231, therefore a(2)=2. - Jon Perry, Jun 08 2003
Number of multiplications performed when computing the determinant of an n X n matrix by definition. - Mats Granvik, Sep 12 2008
Sum of the length of all cycles (excluding fixed points) in all permutations of [n]. - Olivier Gérard, Oct 23 2012
Number of permutations of n distinct objects (ABC...) 1 (one) times >>("-", A, AB, ABC, ABCD, ABCDE, ..., ABCDEFGHIJK, infinity) and one after the other to resemble motif: A (1) AB (1-1), AAB (2-1), AAAB (3-1), AAAAB (4-1), AAAAAB (5-1), AAAAAAB (6-1), AAAAAAAB (7-1), AAAAAAAAB (8-1) etc.,>> "1(one) fixed point". Example:motif: AAAB (or BBBA) 12 * one (1) fixed point etc. Let: AAAB ................ 'A'BCD 1. 'A'BDC 2. 'A'CBD 3. ACDB 'A'DBC 4. 'A'DCB B'A'CD 5. B'A'DC 6. BCAD 7. BCDA BD'A'C 8. BDCA C'A'BD 9. C'A'DB CB'A'D 10. CBDA CDAB CDBA D'A'BC 11. DACB DB'A'C 12. DBCA DCAB DCBA. - Zerinvary Lajos, Nov 27 2009 (does anybody understand what this is supposed to say? - Joerg Arndt, Jan 10 2015)
a(n) is the number of ways to arrange n books on two bookshelves so that each shelf receives at least one book. - Geoffrey Critzer, Feb 21 2010
a(n) = number whose factorial base representation (A007623) begins with digit {n-1} and is followed by n-1 zeros. Viewed in that base, this sequence looks like this: 0, 10, 200, 3000, 40000, 500000, 6000000, 70000000, 800000000, 9000000000, A0000000000, B00000000000, ... (where "digits" A and B stand for placeholder values 10 and 11 respectively). - Antti Karttunen, May 07 2015

Crossrefs

Column 2 of A257503 (apart from initial zero. Equally, row 2 of A257505).
Cf. A001286 (same sequence divided by 2).
Cf. A001563. - Zerinvary Lajos, Aug 27 2008
Cf. sequences with formula (n + k)*n! listed in A282466.

Programs

Formula

a(n) = n! * (n-1).
E.g.f.: x^2/(1-x)^2. - Geoffrey Critzer, Feb 21 2010
a(n) = 2 * A001286(n).
a(n) = A001563(n) - A000142(n). - Antti Karttunen, May 07 2015, hinted by crossref left by Lajos.
From Amiram Eldar, Jul 11 2020: (Start)
Sum_{n>=2} 1/a(n) = Ei(1) + 2 - e - gamma = A091725 + 2 - A001113 - A001620.
Sum_{n>=2} (-1)^n/a(n) = gamma - Ei(-1) - 1/e = A001620 + A099285 - A068985. (End)

Extensions

Last term a(19) corrected by Harry J. Smith, Aug 02 2009

A276955 Square array A(row,col): A(row,1) = A273670(row-1), and for col > 1, A(row,col) = A153880(A(row,col-1)); Dispersion of factorial base left shift A153880.

Original entry on oeis.org

1, 2, 3, 6, 8, 4, 24, 30, 12, 5, 120, 144, 48, 14, 7, 720, 840, 240, 54, 26, 9, 5040, 5760, 1440, 264, 126, 32, 10, 40320, 45360, 10080, 1560, 744, 150, 36, 11, 362880, 403200, 80640, 10800, 5160, 864, 168, 38, 13, 3628800, 3991680, 725760, 85680, 41040, 5880, 960, 174, 50, 15, 39916800, 43545600, 7257600, 766080, 367920, 46080, 6480, 984, 246, 56, 16
Offset: 1

Views

Author

Antti Karttunen, Sep 22 2016

Keywords

Comments

The square array A(row,col) is read by descending antidiagonals: A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), etc.
When viewed in factorial base (A007623) the terms on each row start all with the same prefix, but with an increasing number of zeros appended to the end. For example, for row 8 (A001344 from a(1)=11 onward), the terms written in factorial base look as: 121, 1210, 12100, 121000, ...

Examples

			The top left {1..9} x {1..18} corner of the array:
   1,  2,   6,   24,   120,    720,    5040,    40320,    362880
   3,  8,  30,  144,   840,   5760,   45360,   403200,   3991680
   4, 12,  48,  240,  1440,  10080,   80640,   725760,   7257600
   5, 14,  54,  264,  1560,  10800,   85680,   766080,   7620480
   7, 26, 126,  744,  5160,  41040,  367920,  3669120,  40279680
   9, 32, 150,  864,  5880,  46080,  408240,  4032000,  43908480
  10, 36, 168,  960,  6480,  50400,  443520,  4354560,  47174400
  11, 38, 174,  984,  6600,  51120,  448560,  4394880,  47537280
  13, 50, 246, 1464, 10200,  81360,  730800,  7297920,  80196480
  15, 56, 270, 1584, 10920,  86400,  771120,  7660800,  83825280
  16, 60, 288, 1680, 11520,  90720,  806400,  7983360,  87091200
  17, 62, 294, 1704, 11640,  91440,  811440,  8023680,  87454080
  18, 72, 360, 2160, 15120, 120960, 1088640, 10886400, 119750400
  19, 74, 366, 2184, 15240, 121680, 1093680, 10926720, 120113280
  20, 78, 384, 2280, 15840, 126000, 1128960, 11249280, 123379200
  21, 80, 390, 2304, 15960, 126720, 1134000, 11289600, 123742080
  22, 84, 408, 2400, 16560, 131040, 1169280, 11612160, 127008000
  23, 86, 414, 2424, 16680, 131760, 1174320, 11652480, 127370880
		

Crossrefs

Inverse permutation: A276956.
Transpose: A276953.
Cf. A276949 (index of column where n appears), A276951 (index of row).
Cf. A153880.
Columns 1-3: A273670, A276932, A276933.
The following lists some of the rows that have their own entries. Pattern present in the factorial base expansion of the terms on that row is given in double quotes:
Row 1: A000142 (from a(1)=1, "1" onward),
Row 2: A001048 (from a(2)=3, "11" onward),
Row 3: A052849 (from a(2)=4, "20" onward).
Row 4: A052649 (from a(1)=5, "21" onward).
Row 5: A108217 (from a(3)=7, "101" onward).
Row 6: A054119 (from a(3)=9, "111" onward).
Row 7: A052572 (from a(2)=10, "120" onward).
Row 8: A001344 (from a(1)=11, "121" onward).
Row 13: A052560 (from a(3)=18, "300" onward).
Row 16: A225658 (from a(1)=21, "311" onward).
Row 20: A276940 (from a(3) = 27, "1011" onward).
Related or similar permutations: A257505, A275848, A273666.
Cf. also arrays A276617, A276588 & A276945.

Programs

Formula

A(row,1) = A273670(row-1), and for col > 1, A(row,col) = A153880(A(row,col-1))
As a composition of other permutations:
a(n) = A275848(A257505(n)).

A257503 Square array A(row,col) read by antidiagonals: A(1,col) = A256450(col-1), and for row > 1, A(row,col) = A255411(A(row-1,col)); Dispersion of factorial base shift A255411 (array transposed).

Original entry on oeis.org

1, 2, 4, 3, 12, 18, 5, 16, 72, 96, 6, 22, 90, 480, 600, 7, 48, 114, 576, 3600, 4320, 8, 52, 360, 696, 4200, 30240, 35280, 9, 60, 378, 2880, 4920, 34560, 282240, 322560, 10, 64, 432, 2976, 25200, 39600, 317520, 2903040, 3265920, 11, 66, 450, 3360, 25800, 241920, 357840, 3225600, 32659200, 36288000, 13, 70, 456, 3456, 28800, 246240, 2540160, 3588480, 35925120, 399168000, 439084800
Offset: 1

Views

Author

Antti Karttunen, Apr 27 2015

Keywords

Comments

The array is read by antidiagonals: A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), etc.
The first row (A256450) contains all the numbers which have at least one 1-digit in their factorial base representation (see A007623), after which the successive rows are obtained from the terms on the row immediately above by shifting their factorial representation one left and then incrementing the nonzero digits in that representation with a factorial base shift-operation A255411.

Examples

			The top left corner of the array:
     1,     2,     3,     5,      6,      7,      8,      9,     10,     11,     13
     4,    12,    16,    22,     48,     52,     60,     64,     66,     70,     76
    18,    72,    90,   114,    360,    378,    432,    450,    456,    474,    498
    96,   480,   576,   696,   2880,   2976,   3360,   3456,   3480,   3576,   3696
   600,  3600,  4200,  4920,  25200,  25800,  28800,  29400,  29520,  30120,  30840
  4320, 30240, 34560, 39600, 241920, 246240, 272160, 276480, 277200, 281520, 286560
  ...
		

Crossrefs

Transpose: A257505.
Inverse permutation: A257504.
Row index: A257679, Column index: A257681.
Row 1: A256450, Row 2: A257692, Row 3: A257693.
Columns 1-3: A001563, A062119, A130744 (without their initial zero-terms).
Column 4: A213167 (without the initial one).
Column 5: A052571 (without initial zeros).
Cf. also permutations A255565 and A255566.
Thematically similar arrays: A083412, A135764, A246278.

Programs

Formula

A(1,col) = A256450(col-1), and for row > 1, A(row,col) = A255411(A(row-1,col)).

Extensions

Formula changed because of the changed starting offset of A256450 - Antti Karttunen, May 30 2016

A257679 The smallest nonzero digit present in the factorial base representation (A007623) of n, 0 if no nonzero digits present.

Original entry on oeis.org

0, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 3, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 1, 2, 1, 3, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 3, 1, 1, 1, 2, 1, 4, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 3, 1, 1, 1, 2, 1, 1
Offset: 0

Views

Author

Antti Karttunen, May 04 2015

Keywords

Comments

a(0) = 0 by convention, because "0" has no nonzero digits present.
a(n) gives the row index of n in array A257503 (equally, the column index for array A257505).

Examples

			Factorial base representation (A007623) of 4 is "20", the smallest digit which is not zero is "2", thus a(4) = 2.
		

Crossrefs

Positions of records: A001563.
Cf. A256450, A257692, A257693 (positions of 1's, 2's and 3's in this sequence).
Cf. also A257079, A246359 and arrays A257503, A257505.

Programs

  • Mathematica
    a[n_] := Module[{k = n, m = 2, rmin = n, r}, While[{k, r} = QuotientRemainder[k, m]; k != 0 || r != 0, If[0 < r < rmin, rmin = r]; m++]; rmin]; Array[a, 100, 0] (* Amiram Eldar, Jan 23 2024 *)
  • Python
    def A(n, p=2):
        return n if n
  • Scheme
    (define (A257679 n) (let loop ((n n) (i 2) (mind 0)) (if (zero? n) mind (let ((d (modulo n i))) (loop (/ (- n d) i) (+ 1 i) (cond ((zero? mind) d) ((zero? d) mind) (else (min d mind))))))))
    ;; Alternative implementations based on given recurrences, using memoizing definec-macro:
    (definec (A257679 n) (if (zero? (A257687 n)) (A099563 n) (min (A099563 n) (A257679 (A257687 n)))))
    (definec (A257679 n) (cond ((zero? n) n) ((= 1 (A257680 n)) 1) (else (+ 1 (A257679 (A257684 n))))))
    

Formula

If A257687(n) = 0, then a(n) = A099563(n), otherwise a(n) = min(A099563(n), a(A257687(n))).
In other words, if n is either zero or one of the terms of A051683, then a(n) = A099563(n) [the most significant digit of its f.b.r.], otherwise take the minimum of the most significant digit and a(A257687(n)) [value computed by recursing with a smaller value obtained by discarding that most significant digit].
a(0) = 0, and for n >= 1: if A257680(n) = 1, then a(n) = 1, otherwise 1 + a(A257684(n)).
Other identities:
For all n >= 0, a(A001563(n)) = n. [n * n! gives the first position where n appears. Note also that the "digits" (placeholders) in factorial base representation may get arbitrarily large values.]
For all n >= 0, a(2n+1) = 1 [because all odd numbers end with digit 1 in factorial base].

A265609 Array read by ascending antidiagonals: A(n,k) the rising factorial, also known as Pochhammer symbol, for n >= 0 and k >= 0.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 2, 2, 0, 1, 3, 6, 6, 0, 1, 4, 12, 24, 24, 0, 1, 5, 20, 60, 120, 120, 0, 1, 6, 30, 120, 360, 720, 720, 0, 1, 7, 42, 210, 840, 2520, 5040, 5040, 0, 1, 8, 56, 336, 1680, 6720, 20160, 40320, 40320, 0
Offset: 0

Views

Author

Peter Luschny, Dec 19 2015

Keywords

Comments

The Pochhammer function is defined P(x,n) = x*(x+1)*...*(x+n-1). By convention P(0,0) = 1.
From Antti Karttunen, Dec 19 2015: (Start)
Apart from the initial row of zeros, if we discard the leftmost column and divide the rest of terms A(n,k) with (n+k) [where k is now the once-decremented column index of the new, shifted position] we get the same array back. See the given recursive formula.
When the numbers in array are viewed in factorial base (A007623), certain repeating patterns can be discerned, at least in a few of the topmost rows. See comment in A001710 and arrays A265890, A265892. (End)
A(n,k) is the k-th moment (about 0) of a gamma (Erlang) distribution with shape parameter n and rate parameter 1. - Geoffrey Critzer, Dec 24 2018

Examples

			Square array A(n,k) [where n=row, k=column] is read by ascending antidiagonals as:
A(0,0), A(1,0), A(0,1), A(2,0), A(1,1), A(0,2), A(3,0), A(2,1), A(1,2), A(0,3), ...
Array starts:
n\k [0  1   2    3     4      5        6         7          8]
--------------------------------------------------------------
[0] [1, 0,  0,   0,    0,     0,       0,        0,         0]
[1] [1, 1,  2,   6,   24,   120,     720,     5040,     40320]
[2] [1, 2,  6,  24,  120,   720,    5040,    40320,    362880]
[3] [1, 3, 12,  60,  360,  2520,   20160,   181440,   1814400]
[4] [1, 4, 20, 120,  840,  6720,   60480,   604800,   6652800]
[5] [1, 5, 30, 210, 1680, 15120,  151200,  1663200,  19958400]
[6] [1, 6, 42, 336, 3024, 30240,  332640,  3991680,  51891840]
[7] [1, 7, 56, 504, 5040, 55440,  665280,  8648640, 121080960]
[8] [1, 8, 72, 720, 7920, 95040, 1235520, 17297280, 259459200]
.
Seen as a triangle, T(n, k) = Pochhammer(n - k, k), the first few rows are:
   [0] 1;
   [1] 1, 0;
   [2] 1, 1,  0;
   [3] 1, 2,  2,   0;
   [4] 1, 3,  6,   6,    0;
   [5] 1, 4, 12,  24,   24,    0;
   [6] 1, 5, 20,  60,  120,  120,     0;
   [7] 1, 6, 30, 120,  360,  720,   720,     0;
   [8] 1, 7, 42, 210,  840, 2520,  5040,  5040,     0;
   [9] 1, 8, 56, 336, 1680, 6720, 20160, 40320, 40320, 0.
		

References

  • Ronald L. Graham, Donald E. Knuth and Oren Patashnik, Concrete Mathematics, Addison-Wesley, 1994.
  • H. S. Wall, Analytic Theory of Continued Fractions, Chelsea 1973, p. 355.

Crossrefs

Triangle giving terms only up to column k=n: A124320.
Row 0: A000007, row 1: A000142, row 3: A001710 (from k=1 onward, shifted two terms left).
Column 0: A000012, column 1: A001477, column 2: A002378, columns 3-7: A007531, A052762, A052787, A053625, A159083 (shifted 2 .. 6 terms left respectively, i.e. without the extra initial zeros), column 8: A239035.
Row sums of the triangle: A000522.
A(n, n) = A000407(n-1) for n>0.
2^n*A(1/2,n) = A001147(n).
Cf. also A007623, A008279 (falling factorial), A173333, A257505, A265890, A265892.

Programs

  • Maple
    for n from 0 to 8 do seq(pochhammer(n,k), k=0..8) od;
  • Mathematica
    Table[Pochhammer[n, k], {n, 0, 8}, {k, 0, 8}]
  • Sage
    for n in (0..8): print([rising_factorial(n,k) for k in (0..8)])
    
  • Scheme
    (define (A265609 n) (A265609bi (A025581 n) (A002262 n)))
    (define (A265609bi row col) (if (zero? col) 1 (* (+ row col -1) (A265609bi row (- col 1)))))
    ;; Antti Karttunen, Dec 19 2015

Formula

A(n,k) = Gamma(n+k)/Gamma(n) for n > 0 and n^k for n=0.
A(n,k) = Sum_{j=0..k} n^j*S1(k,j), S1(n,k) the Stirling cycle numbers A132393(n,k).
A(n,k) = (k-1)!/(Sum_{j=0..k-1} (-1)^j*binomial(k-1, j)/(j+n)) for n >= 1, k >= 1.
A(n,k) = (n+k-1)*A(n,k-1) for k >= 1, A(n,0) = 1. - Antti Karttunen, Dec 19 2015
E.g.f. for row k: 1/(1-x)^k. - Geoffrey Critzer, Dec 24 2018
A(n, k) = FallingFactorial(n + k - 1, k). - Peter Luschny, Mar 22 2022
G.f. for row n as a continued fraction of Stieltjes type: 1/(1 - n*x/(1 - x/(1 - (n+1)*x/(1 - 2*x/(1 - (n+2)*x/(1 - 3*x/(1 - ... ))))))). See Wall, Chapter XVIII, equation 92.5. Cf. A226513. - Peter Bala, Aug 27 2023

A257681 Column index of A257503: If A257685(n) = 0, then a(n) = A257682(n), otherwise a(n) = a(A257685(n)).

Original entry on oeis.org

0, 1, 2, 3, 1, 4, 5, 6, 7, 8, 9, 10, 2, 11, 12, 13, 3, 14, 1, 15, 16, 17, 4, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 5, 43, 44, 45, 6, 46, 47, 48, 49, 50, 51, 52, 7, 53, 54, 55, 8, 56, 9, 57, 58, 59, 10, 60, 2
Offset: 0

Views

Author

Antti Karttunen, May 04 2015

Keywords

Comments

a(0) = 0 by convention, because 0 is not present in arrays A257503 and A257505.

Crossrefs

Column index of A257503, row index of A257505.
Cf. also A257679 (gives the other index).

Formula

If A257685(n) = 0, then a(n) = A257682(n), otherwise a(n) = a(A257685(n)).

A052571 E.g.f. x^3/(1-x)^2.

Original entry on oeis.org

0, 0, 0, 6, 48, 360, 2880, 25200, 241920, 2540160, 29030400, 359251200, 4790016000, 68497228800, 1046139494400, 16999766784000, 292919058432000, 5335311421440000, 102437979291648000, 2067966706950144000
Offset: 0

Views

Author

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Keywords

Comments

For n >= 3, a(n) = number whose factorial base representation (A007623) begins with digit {n-2} followed by n-1 zeros. Viewed in that base, this sequence looks like this: 0, 0, 0, 100, 2000, 30000, 400000, 5000000, 60000000, 700000000, 8000000000, 90000000000, A00000000000, B000000000000, ... (where "digits" A and B stand for placeholder values 10 and 11 respectively). - Antti Karttunen, May 07 2015

Crossrefs

Column 5 of A257503 (apart from zero terms. Equally, row 5 of A257505).
Cf. sequences with formula (n + k)*n! listed in A282466.

Programs

  • Magma
    [0,0] cat [n*(n+1)*(n+2)*Factorial(n): n in [0..20]]; // Vincenzo Librandi, Oct 11 2011
    
  • Maple
    spec := [S,{S=Prod(Z,Z,Z,Sequence(Z),Sequence(Z))},labeled]: seq(combstruct[count](spec,size=n), n=0..20);
    [seq (n*(n+1)*(n+2)*n!,n=0..17)]; # Zerinvary Lajos, Nov 25 2006
    a:=n->add((n!),j=1..n-2):seq(a(n), n=0..21); # Zerinvary Lajos, Aug 27 2008
    G(x):=x^3/(1-x)^2: f[0]:=G(x): for n from 1 to 21 do f[n]:=diff(f[n-1],x) od: x:=0: seq(f[n],n=0..19); # Zerinvary Lajos, Apr 01 2009
  • Mathematica
    Table[Sum[n!, {i, 3, n}], {n, 0, 19}] (* Zerinvary Lajos, Jul 12 2009 *)
    With[{nn=20},CoefficientList[Series[x^3/(1-x)^2,{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Feb 27 2025 *)
  • Scheme
    (define (A052571 n) (if (< n 2) 0 (* (- n 2) (A000142 n)))) ;; Antti Karttunen, May 07 2015

Formula

E.g.f.: x^3/(-1+x)^2.
Recurrence: {a(1)=0, a(0)=0, a(2)=0, a(3)=6, (1-n^2)*a(n)+(-2+n)*a(n+1)=0}.
For n >= 2, a(n) = (n-2)*n!.
a(n+2) = n*(n+1)*(n+2)*n!. - Zerinvary Lajos, Nov 25 2006
a(n) = 3*A090672(n-2) = 6*A005990(n-2). - Zerinvary Lajos, May 11 2007
From Amiram Eldar, Jan 14 2021: (Start)
Sum_{n>=3} 1/a(n) = 9/4 - e - gamma/2 + Ei(1)/2 = 9/4 - A001113 - (1/2)*A001620 + (1/2)*A091725.
Sum_{n>=3} (-1)^(n+1)/a(n) = -1/4 + gamma/2 - Ei(-1)/2 = -1/4 + (1/2)*A001620 + (1/2)*A099285. (End)
Showing 1-10 of 14 results. Next