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 17 results. Next

A257993 Least gap in the partition having Heinz number n; index of the least prime not dividing n.

Original entry on oeis.org

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

Views

Author

Emeric Deutsch, May 18 2015

Keywords

Comments

The "least gap" of a partition is the least positive integer that is not a part of the partition. For example, the least gap of the partition [7,4,2,2,1] is 3.
We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product(p_j-th prime, j=1...r) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] we get 2*2*3*7*29 = 2436.
In the Maple program the subprogram B yields the partition with Heinz number n.
Sum of least gaps of all partitions of m = A022567(m).
From Antti Karttunen, Aug 22 2016: (Start)
Index of the least prime not dividing n. (After a formula given by Heinz.)
Least k such that A002110(k) does not divide n.
One more than the number of trailing zeros in primorial base representation of n, A049345.
(End)
The least gap is also called the mex (minimal excludant) of the partition. - Gus Wiseman, Apr 20 2021

Examples

			a(18) = 3 because the partition having Heinz number 18 = 2*3*3 is [1,2,2], having least gap equal to 3.
		

References

  • G. E. Andrews and K. Eriksson, Integer Partitions, Cambridge Univ. Press, 2004, Cambridge.
  • Miklós Bóna, A Walk Through Combinatorics, World Scientific Publishing Co., 2002.

Crossrefs

Positions of 1's are A005408.
Positions of 2's are A047235.
The number of gaps is A079067.
The version for crank is A257989.
The triangle counting partitions by this statistic is A264401.
One more than A276084.
The version for greatest difference is A286469 or A286470.
A maximal instead of minimal version is A339662.
Positions of even terms are A342050.
Positions of odd terms are A342051.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A056239 adds up prime indices, row sums of A112798.
A073491 lists numbers with gap-free prime indices.
A238709 counts partitions by sum and least difference.
A333214 lists positions of adjacent unequal prime gaps.
A339737 counts partitions by sum and greatest gap.

Programs

  • Maple
    with(numtheory): a := proc (n) local B, q: B := proc (n) local nn, j, m: nn := op(2, ifactors(n)): for j to nops(nn) do m[j] := op(j, nn) end do: [seq(seq(pi(op(1, m[i])), q = 1 .. op(2, m[i])), i = 1 .. nops(nn))] end proc: for q while member(q, B(n)) = true do  end do: q end proc: seq(a(n), n = 1 .. 150);
    # second Maple program:
    a:= n-> `if`(n=1, 1, (s-> min({$1..(max(s)+1)} minus s))(
            {map(x-> numtheory[pi](x[1]), ifactors(n)[2])[]})):
    seq(a(n), n=1..100);  # Alois P. Heinz, May 09 2016
    # faster:
    A257993 := proc(n) local p, c; c := 1; p := 2;
    while n mod p = 0 do p := nextprime(p); c := c + 1 od: c end:
    seq(A257993(n), n=1..100); # Peter Luschny, Jun 04 2017
  • Mathematica
    A053669[n_] := For[p = 2, True, p = NextPrime[p], If[CoprimeQ[p, n], Return[p]]]; a[n_] := PrimePi[A053669[n]]; Array[a, 100] (* Jean-François Alcover, Nov 28 2016 *)
    Table[k = 1; While[! CoprimeQ[Prime@ k, n], k++]; k, {n, 100}] (* Michael De Vlieger, Jun 22 2017 *)
  • PARI
    a(n) = forprime(p=2,, if (n % p, return(primepi(p)))); \\ Michel Marcus, Jun 22 2017
  • Python
    from sympy import nextprime, primepi
    def a053669(n):
        p = 2
        while True:
            if n%p!=0: return p
            else: p=nextprime(p)
    def a(n): return primepi(a053669(n)) # Indranil Ghosh, May 12 2017
    
  • Scheme
    (define (A257993 n) (let loop ((n n) (i 1)) (let* ((p (A000040 i)) (d (modulo n p))) (if (not (zero? d)) i (loop (/ (- n d) p) (+ 1 i))))))
    ;; Antti Karttunen, Aug 22 2016
    

Formula

a(n) = A000720(A053669(n)). - Alois P. Heinz, May 18 2015
From Antti Karttunen, Aug 22-30 2016: (Start)
a(n) = 1 + A276084(n).
a(n) = A055396(A276086(n)).
A276152(n) = A002110(a(n)).
(End)
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1 + Sum_{k>=1} 1/A002110(k) = 1.705230... (1 + A064648). - Amiram Eldar, Jul 23 2022
a(n) << log n/log log n. - Charles R Greathouse IV, Dec 03 2022

Extensions

A simpler description added to the name by Antti Karttunen, Aug 22 2016

A325405 Heinz numbers of integer partitions y such that the k-th differences of y are distinct for all k >= 0 and are disjoint from the i-th differences for i != k.

Original entry on oeis.org

1, 2, 3, 5, 7, 10, 11, 13, 14, 15, 17, 19, 22, 23, 26, 29, 31, 33, 34, 35, 37, 38, 39, 41, 43, 46, 47, 51, 53, 55, 57, 58, 59, 61, 62, 67, 69, 71, 73, 74, 77, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 101, 103, 106, 107, 109, 111, 113, 115, 118, 119, 122
Offset: 1

Views

Author

Gus Wiseman, May 02 2019

Keywords

Comments

First differs from A325388 in lacking 130.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
The differences of a sequence are defined as if the sequence were increasing, so for example the differences of (6,3,1) are (-3,-2).
The zeroth differences of a sequence are the sequence itself, while the k-th differences for k > 0 are the differences of the (k-1)-th differences.
The enumeration of these partitions by sum is given by A325404.

Examples

			The sequence of terms together with their prime indices begins:
    1: {}
    2: {1}
    3: {2}
    5: {3}
    7: {4}
   10: {1,3}
   11: {5}
   13: {6}
   14: {1,4}
   15: {2,3}
   17: {7}
   19: {8}
   22: {1,5}
   23: {9}
   26: {1,6}
   29: {10}
   31: {11}
   33: {2,5}
   34: {1,7}
   35: {3,4}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],UnsameQ@@Join@@Table[Differences[primeMS[#],k],{k,0,PrimeOmega[#]}]&]

A295370 Number of permutations of [n] avoiding three consecutive terms in arithmetic progression.

Original entry on oeis.org

1, 1, 2, 4, 18, 80, 482, 3280, 26244, 231148, 2320130, 25238348, 302834694, 3909539452, 54761642704, 816758411516, 13076340876500, 221396129723368, 3985720881222850, 75503196628737920, 1510373288335622576, 31634502738658957588, 696162960370556156224, 15978760340940405262668
Offset: 0

Views

Author

Alois P. Heinz, Nov 20 2017

Keywords

Comments

These are permutations of n whose second-differences are nonzero. - Gus Wiseman, Jun 03 2019

Examples

			a(3) = 4: 132, 213, 231, 312.
a(4) = 18: 1243, 1324, 1342, 1423, 2134, 2143, 2314, 2413, 2431, 3124, 3142, 3241, 3412, 3421, 4132, 4213, 4231, 4312.
		

Crossrefs

Programs

  • Maple
    b:= proc(s, j, k) option remember; `if`(s={}, 1,
          add(`if`(k=0 or 2*j<>i+k, b(s minus {i}, i,
              `if`(2*i-j in s, j, 0)), 0), i=s))
        end:
    a:= n-> b({$1..n}, 0$2):
    seq(a(n), n=0..12);
  • Mathematica
    Table[Length[Select[Permutations[Range[n]],!MemberQ[Differences[#,2],0]&]],{n,0,5}] (* Gus Wiseman, Jun 03 2019 *)
    b[s_, j_, k_] := b[s, j, k] = If[s == {}, 1, Sum[If[k == 0 || 2*j != i + k, b[s~Complement~{i}, i, If[MemberQ[s, 2*i - j ], j, 0]], 0], {i, s}]];
    a[n_] := a[n] = b[Range[n], 0, 0];
    Table[Print[n, " ", a[n]]; a[n], {n, 0, 16}] (* Jean-François Alcover, Nov 20 2023, after Alois P. Heinz *)

Extensions

a(22)-a(23) from Vaclav Kotesovec, Mar 22 2022

A325404 Number of reversed integer partitions y of n such that the k-th differences of y are distinct for all k >= 0 and are disjoint from the i-th differences for i != k.

Original entry on oeis.org

1, 1, 1, 1, 2, 3, 2, 4, 4, 4, 5, 7, 5, 11, 12, 11, 12, 20, 15, 24, 22, 27, 28, 37, 28, 45, 43, 48, 50, 66, 58, 79, 72, 84, 87, 112, 106, 135, 128, 158, 147, 186, 180, 218, 220, 265, 246, 304, 303, 354, 340, 412, 418, 471, 463, 538, 543, 642, 600, 711, 755
Offset: 0

Views

Author

Gus Wiseman, May 02 2019

Keywords

Comments

The differences of a sequence are defined as if the sequence were increasing, so for example the differences of (6,3,1) are (-3,-2).
The zeroth differences of a sequence are the sequence itself, while the k-th differences for k > 0 are the differences of the (k-1)-th differences.
The Heinz numbers of these partitions are given by A325405.

Examples

			The a(1) = 1 through a(12) = 5 reversed partitions (A = 10, B = 11, C = 12):
  (1)  (2)  (3)  (4)   (5)   (6)   (7)   (8)   (9)   (A)   (B)    (C)
                 (13)  (14)  (15)  (16)  (17)  (18)  (19)  (29)   (39)
                       (23)        (25)  (26)  (27)  (28)  (38)   (57)
                                   (34)  (35)  (45)  (37)  (47)   (1B)
                                                     (46)  (56)   (2A)
                                                           (1A)
                                                           (146)
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[Reverse/@IntegerPartitions[n],UnsameQ@@Join@@Table[Differences[#,k],{k,0,Length[#]}]&]],{n,0,30}]

A355526 Maximal difference between adjacent prime indices of n, or k if n is the k-th prime.

Original entry on oeis.org

1, 2, 0, 3, 1, 4, 0, 0, 2, 5, 1, 6, 3, 1, 0, 7, 1, 8, 2, 2, 4, 9, 1, 0, 5, 0, 3, 10, 1, 11, 0, 3, 6, 1, 1, 12, 7, 4, 2, 13, 2, 14, 4, 1, 8, 15, 1, 0, 2, 5, 5, 16, 1, 2, 3, 6, 9, 17, 1, 18, 10, 2, 0, 3, 3, 19, 6, 7, 2, 20, 1, 21, 11, 1, 7, 1, 4, 22, 2, 0, 12
Offset: 2

Views

Author

Gus Wiseman, Jul 10 2022

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Examples

			The prime indices of 9842 are {1,4,8,12}, with differences (3,4,4), so a(9842) = 4.
		

Crossrefs

Crossrefs found in the link are not repeated here.
Positions of first appearances are 4 followed by A000040.
Positions of 0's are A025475, minimal version A013929.
Positions of 1's are 2 followed by A066312, minimal version A355527.
Triangle A238710 counts m such that A056239(m) = n and a(m) = k.
Prepending 0 to the prime indices gives A286469, minimal version A355528.
See also A286470, minimal version A355524.
The minimal version is A355525, triangle A238709.
The augmented version is A355532.
A001522 counts partitions with a fixed point (unproved), ranked by A352827.
A287352, A355533, A355534, A355536 list the differences of prime indices.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[If[PrimeQ[n],PrimePi[n],Max@@Differences[primeMS[n]]],{n,2,100}]

A355525 Minimal difference between adjacent prime indices of n, or k if n is the k-th prime.

Original entry on oeis.org

1, 2, 0, 3, 1, 4, 0, 0, 2, 5, 0, 6, 3, 1, 0, 7, 0, 8, 0, 2, 4, 9, 0, 0, 5, 0, 0, 10, 1, 11, 0, 3, 6, 1, 0, 12, 7, 4, 0, 13, 1, 14, 0, 0, 8, 15, 0, 0, 0, 5, 0, 16, 0, 2, 0, 6, 9, 17, 0, 18, 10, 0, 0, 3, 1, 19, 0, 7, 1, 20, 0, 21, 11, 0, 0, 1, 1, 22, 0, 0, 12
Offset: 2

Views

Author

Gus Wiseman, Jul 10 2022

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Examples

			The prime indices of 9842 are {1,4,8,12}, with differences (3,4,4), so a(9842) = 3.
		

Crossrefs

Crossrefs found in the link are not repeated here.
Positions of first appearances are 4 followed by A000040.
Positions of 0's are A013929, see also A130091.
Triangle A238709 counts m such that A056239(m) = n and a(m) = k.
For maximal instead of minimal difference we have A286470.
Positions of terms > 1 are A325160, also A325161.
See also A355524, A355528.
Positions of 1's are A355527.
A001522 counts partitions with a fixed point (unproved), ranked by A352827.
A238352 counts partitions by fixed points, rank statistic A352822.
A287352, A355533, A355534, A355536 list the differences of prime indices.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[If[PrimeQ[n],PrimePi[n],Min@@Differences[primeMS[n]]],{n,2,100}]

A325406 Triangle read by rows where T(n,k) is the number of reversed integer partitions of n with k distinct differences of any degree.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 2, 0, 0, 1, 2, 2, 0, 0, 1, 1, 3, 2, 0, 0, 1, 4, 2, 3, 1, 0, 0, 1, 1, 5, 5, 2, 1, 0, 0, 1, 3, 5, 6, 3, 3, 1, 0, 0, 1, 3, 4, 8, 7, 1, 4, 2, 0, 0, 1, 3, 6, 11, 7, 5, 2, 4, 2, 1, 0, 1, 1, 6, 13, 8, 9, 9, 0, 4, 3, 1, 0, 1, 6, 7, 11, 12, 9
Offset: 0

Views

Author

Gus Wiseman, May 03 2019

Keywords

Comments

The differences of a sequence are defined as if the sequence were increasing, so for example the differences of (6,3,1) are (-3,-2).
The zeroth differences of a sequence are the sequence itself, while the k-th differences for k > 0 are the differences of the (k-1)-th differences. The distinct differences of any degree are the union of the k-th differences for all k >= 0. For example, the k-th differences of (1,1,2,4) for k = 0...3 are:
(1,1,2,4)
(0,1,2)
(1,1)
(0)
so there are a total of 4 distinct differences of any degree, namely {0,1,2,4}.

Examples

			Triangle begins:
  1
  0  1
  0  1  1
  0  1  2  0
  0  1  2  2  0
  0  1  1  3  2  0
  0  1  4  2  3  1  0
  0  1  1  5  5  2  1  0
  0  1  3  5  6  3  3  1  0
  0  1  3  4  8  7  1  4  2  0
  0  1  3  6 11  7  5  2  4  2  1
  0  1  1  6 13  8  9  9  0  4  3  1
  0  1  6  7 11 12  9 10  8  4  3  2  2
  0  1  1  7 18  9 14 19  5 10  3  5  4  1
  0  1  3  9 17  9 22 20 15  9  7  6  5  4  1
  0  1  4  8 22 11 16 24 22 19 10 11  2  8  7  2
  0  1  4 10 23 15 24 23 27 27 12 14 11  8  8  5  5
Row n = 8 counts the following partitions:
  (8)  (44)        (17)       (116)     (134)   (1133)   (111122)
       (2222)      (26)       (125)     (233)   (11123)
       (11111111)  (35)       (1115)    (1223)  (11222)
                   (224)      (1124)
                   (1111112)  (11114)
                              (111113)
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[Reverse/@IntegerPartitions[n],Length[Union@@Table[Differences[#,i],{i,0,Length[#]}]]==k&]],{n,0,16},{k,0,n}]

A325849 Number of strict compositions of n with no three consecutive parts in arithmetic progression.

Original entry on oeis.org

1, 1, 1, 3, 3, 5, 9, 13, 19, 23, 51, 57, 91, 117, 179, 283, 381, 531, 737, 1017, 1335, 2259, 2745, 3983, 5289, 7367, 9413, 13155, 19461, 25129, 33997, 45633, 61225, 80481, 107091, 137475, 205243, 253997, 345527, 447003, 604919, 768331, 1026167, 1299227
Offset: 0

Views

Author

Gus Wiseman, May 31 2019

Keywords

Comments

A composition of n is a finite sequence of positive integers with sum n. a(n) is the number of strict compositions of n with no two of their adjacent first-differences equal, or with no 0's in their second-differences.

Examples

			The a(1) = 1 through a(8) = 19 compositions:
  (1)  (2)  (3)   (4)   (5)   (6)    (7)    (8)
            (12)  (13)  (14)  (15)   (16)   (17)
            (21)  (31)  (23)  (24)   (25)   (26)
                        (32)  (42)   (34)   (35)
                        (41)  (51)   (43)   (53)
                              (132)  (52)   (62)
                              (213)  (61)   (71)
                              (231)  (124)  (125)
                              (312)  (142)  (134)
                                     (214)  (143)
                                     (241)  (152)
                                     (412)  (215)
                                     (421)  (251)
                                            (314)
                                            (341)
                                            (413)
                                            (431)
                                            (512)
                                            (521)
		

Crossrefs

The non-strict case is A238423.

Programs

  • Mathematica
    Table[Length[Select[Join@@Permutations/@Select[IntegerPartitions[n],UnsameQ@@#&],!MemberQ[Differences[#,2],0]&]],{n,0,30}]

A325852 Number of (strict) integer partitions of n whose differences of all degrees are nonzero.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 3, 5, 6, 6, 9, 11, 11, 15, 19, 19, 26, 31, 31, 41, 49, 53, 62, 75, 81, 97, 112, 124, 145, 171, 175, 215, 244, 274, 307, 344, 388, 446, 497, 561, 599, 700, 779, 881, 981, 1054, 1184, 1340, 1500, 1669, 1767, 2031, 2237, 2486, 2765, 2946, 3300
Offset: 0

Views

Author

Gus Wiseman, May 31 2019

Keywords

Comments

The differences of a sequence are defined as if the sequence were increasing, so for example the differences of (6,3,1) are (-3,-2). The zeroth differences are the sequence itself, while k-th differences for k > 0 are the differences of the (k-1)-th differences. The differences of all degrees of a sequence are the union of its zeroth through m-th differences, where m is the length of the sequence.

Examples

			The a(1) = 1 through a(11) = 11 partitions (A = 10, B = 11):
  (1)  (2)  (3)   (4)   (5)   (6)   (7)    (8)    (9)    (A)    (B)
            (21)  (31)  (32)  (42)  (43)   (53)   (54)   (64)   (65)
                        (41)  (51)  (52)   (62)   (63)   (73)   (74)
                                    (61)   (71)   (72)   (82)   (83)
                                    (421)  (431)  (81)   (91)   (92)
                                           (521)  (621)  (532)  (A1)
                                                         (541)  (542)
                                                         (631)  (632)
                                                         (721)  (641)
                                                                (731)
                                                                (821)
		

Crossrefs

The case for only degrees > 1 is A325874.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],!MemberQ[Union@@Table[Differences[#,i],{i,Length[#]}],0]&]],{n,0,30}]

A339737 Triangle read by rows where T(n,k) is the number of integer partitions of n with greatest gap k.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 2, 0, 1, 0, 2, 1, 1, 1, 0, 3, 1, 1, 1, 1, 0, 4, 1, 2, 2, 1, 1, 0, 5, 1, 3, 2, 2, 1, 1, 0, 6, 2, 3, 4, 3, 2, 1, 1, 0, 8, 2, 4, 5, 4, 3, 2, 1, 1, 0, 10, 2, 5, 7, 6, 5, 3, 2, 1, 1, 0, 12, 3, 6, 8, 9, 6, 5, 3, 2, 1, 1, 0, 15, 3, 8, 11, 11, 10, 7, 5, 3, 2, 1, 1, 0
Offset: 0

Views

Author

Gus Wiseman, Apr 20 2021

Keywords

Comments

We define the greatest gap of a partition to be the greatest nonnegative integer less than the greatest part and not in the partition.

Examples

			Triangle begins:
   1
   1   0
   1   1   0
   2   0   1   0
   2   1   1   1   0
   3   1   1   1   1   0
   4   1   2   2   1   1   0
   5   1   3   2   2   1   1   0
   6   2   3   4   3   2   1   1   0
   8   2   4   5   4   3   2   1   1   0
  10   2   5   7   6   5   3   2   1   1   0
  12   3   6   8   9   6   5   3   2   1   1   0
  15   3   8  11  11  10   7   5   3   2   1   1   0
  18   4   9  13  15  13  10   7   5   3   2   1   1   0
  22   5  10  17  19  18  14  11   7   5   3   2   1   1   0
  27   5  13  20  24  23  20  14  11   7   5   3   2   1   1   0
For example, row n = 9 counts the following partitions:
  (3321)       (432)   (333)      (54)      (522)    (63)    (72)   (81)  (9)
  (22221)      (3222)  (4311)     (441)     (531)    (621)   (711)
  (32211)              (33111)    (4221)    (5211)   (6111)
  (222111)             (3111111)  (42111)   (51111)
  (321111)                        (411111)
  (2211111)
  (21111111)
  (111111111)
		

Crossrefs

Column k = 0 is A000009.
Row sums are A000041.
Central diagonal is A000041.
Column k = 1 is A087897.
The version for least gap is A264401, with Heinz number encoding A257993.
The version for greatest difference is A286469 or A286470.
An encoding (of greatest gap) using Heinz numbers is A339662.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A015723 counts strict partitions with a selected part.
A048004 counts compositions by greatest part.
A056239 adds up prime indices, row sums of A112798.
A064391 is the version for crank.
A064428 counts partitions of nonnegative crank.
A073491 list numbers with gap-free prime indices.
A107428 counts gap-free compositions.
A238709/A238710 counts partitions by least/greatest difference.
A342050/A342051 have prime indices with odd/even least gap.

Programs

  • Mathematica
    maxgap[q_]:=Max@@Complement[Range[0,If[q=={},0,Max[q]]],q];
    Table[Length[Select[IntegerPartitions[n],maxgap[#]==k&]],{n,0,15},{k,0,n}]
  • PARI
    S(n,k)={if(k>n, O(x*x^n), x^k*(S(n-k,k+1) + 1)/(1 - x^k))}
    ColGf(k,n) = {(k==0) + S(n,k+1)/prod(j=1, k-1, 1 - x^j + O(x^max(1,n-k)))}
    A(n,m=n)={Mat(vector(m+1, k, Col(ColGf(k-1,n), -(n+1))))}
    { my(M=A(10)); for(i=1, #M, print(M[i,1..i])) } \\ Andrew Howroyd, Jan 13 2024

Extensions

Offset corrected by Andrew Howroyd, Jan 13 2024
Showing 1-10 of 17 results. Next