A076954 -id:A076954 - cp's OEIS frontend

A008278 Reflected triangle of Stirling numbers of 2nd kind, S(n,n-k+1), n >= 1, 1 <= k <= n.

1, 1, 1, 1, 3, 1, 1, 6, 7, 1, 1, 10, 25, 15, 1, 1, 15, 65, 90, 31, 1, 1, 21, 140, 350, 301, 63, 1, 1, 28, 266, 1050, 1701, 966, 127, 1, 1, 36, 462, 2646, 6951, 7770, 3025, 255, 1, 1, 45, 750, 5880, 22827, 42525, 34105, 9330, 511, 1

Offset: 1

N. J. A. Sloane

The n-th row also gives the coefficients of the sigma polynomial of the empty graph \bar K_n. - Eric W. Weisstein, Apr 07 2017
The n-th row also gives the coefficients of the independence polynomial of the (n-1)-triangular honeycomb bishop graph. - Eric W. Weisstein, Apr 03 2018
From Gus Wiseman, Aug 11 2020: (Start)
Conjecture: also the number of divisors of the superprimorial A006939(n - 1) that have 0 <= k <= n distinct prime factors, all appearing with distinct multiplicities. For example, row n = 4 counts the following divisors of 360:
  1  2  12  360
     3  18
     4  20
     5  24
     8  40
     9  45
        72
Equivalently, T(n,k) is the number of length-n vectors 0 <= v_i <= i with k nonzero values, all of which are distinct.
Crossrefs:
A006939 lists superprimorials or Chernoff numbers.
A022915 counts permutations of prime indices of superprimorials.
A076954 can be used instead of A006939.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A336420 is the version counting all prime factors, not just distinct ones.
Cf. A000005, A027423, A095149, A124010, A317829, A327498, A336419, A336421.
(End)
From Leonidas Liponis, Aug 26 2024: (Start)
It appears that this sequence is related to the combinatorial form of Faà di Bruno's formula. Specifically, the number of terms for the n-th derivative of a composite function y = f(g(x)) matches the number of partitions of n.
For example, consider the case where g(x) = e^x, in which all derivatives of g(x) are equal. The first 5 rows of A008278 appear as the factors of derivatives of f(x), highlighted here in brackets:
          dy/dx = [ 1 ] * f'(e^x) * e^x
          d^2y/dx^2 = [ 1 ] * f''(e^x) * e^{2x} + [ 1 ] * f'(e^x) * e^x
          d^3y/dx^3 = [ 1 ] * f'''(e^x) * e^{3x} + [ 3 ] * f''(e^x) * e^{2x} + [ 1 ] * f'(e^x) * e^x
          d^4y/dx^4 = [ 1 ] * f''''(e^x) * e^{4x} + [ 6 ] * f'''(e^x) * e^{3x} + [ 7 ] * f''(e^x) * e^{2x} + [ 1 ] * f'(e^x) * e^x
          d^5y/dx^5 = [ 1 ] * f'''''(e^x) * e^{5x} + [ 10 ] * f''''(e^x) * e^{4x} + [ 25 ] * f'''(e^x) * e^{3x} + [ 15 ] * f''(e^x) * e^{2x} + [ 1 ] * f'(e^x) * e^x
This pattern is observed in Mathematica for the first 10 cases, using the code below.
(End)

			The e.g.f. of [0,0,1,7,25,65,...], the k=3 column of A008278, but with offset n=0, is exp(x)*(1*(x^2)/2! + 4*(x^3)/3! + 3*(x^4)/4!).
Triangle starts:
  1;
  1,  1;
  1,  3,   1;
  1,  6,   7,    1;
  1, 10,  25,   15,    1;
  1, 15,  65,   90,   31,    1;
  1, 21, 140,  350,  301,   63,    1;
  1, 28, 266, 1050, 1701,  966,  127,   1;
  1, 36, 462, 2646, 6951, 7770, 3025, 255, 1;
  ...

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 835.
F. N. David, M. G. Kendall and D. E. Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 223.
R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics, Addison-Wesley, 2nd ed., 1994.

T. D. Noe, Rows n = 0..100 of triangle, flattened
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
Noureddine Chair, Exact two-point resistance, and the simple random walk on the complete graph minus N edges, Ann. Phys. 327, No. 12, 3116-3129 (2012), eq. (27).
Xi Chen, Bishal Deb, Alexander Dyachenko, Tomack Gilmore, and Alan D. Sokal, Coefficientwise total positivity of some matrices defined by linear recurrences, arXiv:2012.03629 [math.CO], 2020.
T. Copeland, Generators, Inversion, and Matrix, Binomial, and Integral Transforms, 2015
U. N. Katugampola, A new Fractional Derivative and its Mellin Transform, arXiv preprint arXiv:1106.0965 [math.CA], 2011.
Eric Weisstein's World of Mathematics, Bell Polynomial
Eric Weisstein's World of Mathematics, Empty Graph
Eric Weisstein's World of Mathematics, Independence Polynomial
Eric Weisstein's World of Mathematics, Sigma Polynomial
Eric Weisstein's World of Mathematics, Stirling Number of the Second Kind

See A008277 and A048993, which are the main entries for this triangle of numbers.

Cf. A094262, A008277, A008276, A003422, A000166, A000110, A000204, A000045, A000108, A213735.

a008278 n k = a008278_tabl !! (n-1) !! (k-1)
a008278_row n = a008278_tabl !! (n-1)
a008278_tabl = iterate st2 [1] where
  st2 row = zipWith (+) ([0] ++ row') (row ++ [0])
            where row' = reverse $ zipWith (*) [1..] $ reverse row
-- Reinhard Zumkeller, Jun 22 2013

rows = 10; Flatten[Table[StirlingS2[n, k], {n, 1, rows}, {k, n, 1, -1}]] (* Jean-François Alcover, Nov 17 2011, *)
Table[CoefficientList[x^n BellB[n, 1/x], x], {n, 10}] // Flatten (* Eric W. Weisstein, Apr 05 2017 *)
n = 5; Grid[Prepend[Transpose[{Range[1, n], Table[D[f[Exp[x]], {x, i}], {i, 1, n}]}], {"Order","Derivative"}], Frame -> All, Spacings -> {2, 1}] (* Leonidas Liponis, Aug 27 2024 *)

for(n=1,10,for(k=1,n,print1(stirling(n,n-k+1,2),", "))) \\ Hugo Pfoertner, Aug 30 2020

T(n, k)=0 if n < k, T(n, 0)=0, T(1, 1)=1, T(n, k) = (n-k+1)*T(n-1, k-1) + T(n-1, k) otherwise.
O.g.f. for the k-th column: 1/(1-x) if k=1 and A(k,x):=((x^k)/(1-x)^(2*k+1))*Sum_{m=0..k-1} A008517(k,m+1)*x^m if k >= 2. A008517 is the second-order Eulerian triangle. Cf. p. 257, eq. (6.43) of the R. L. Graham et al. book. - Wolfdieter Lang, Oct 14 2005
E.g.f. for the k-th column (with offset n=0): E(k,x):=exp(x)*Sum_{m=0..k-1} A112493(k-1,m)*(x^(k-1+m))/(k-1+m)! if k >= 1. - Wolfdieter Lang, Oct 14 2005
a(n) = abs(A213735(n-1)). - Hugo Pfoertner, Sep 07 2020

Name edited by Gus Wiseman, Aug 11 2020

A317829 Number of set partitions of multiset {1, 2, 2, 3, 3, 3, ..., n X n}.

Original entry on oeis.org

1, 1, 4, 52, 2776, 695541, 927908528, 7303437156115, 371421772559819369, 132348505150329265211927, 355539706668772869353964510735, 7698296698535929906799439134946965681, 1428662247641961794158621629098030994429958386, 2405509035205023556420199819453960482395657232596725626

Offset: 0

Antti Karttunen, Aug 10 2018

nonn

Number of factorizations of the superprimorial A006939(n) into factors > 1. - Gus Wiseman, Aug 21 2020

			For n = 2 we have a multiset {1, 2, 2} which can be partitioned as {{1}, {2}, {2}} or {{1, 2}, {2}} or {{1}, {2, 2}} or {{1, 2, 2}}, thus a(2) = 4.

Index entries for sequences related to partitions

Cf. A033312, A317826.
Subsequence of A317828.
A000142 counts submultisets of the same multiset.
A022915 counts permutations of the same multiset.
A337069 is the strict case.
A001055 counts factorizations.
A006939 lists superprimorials or Chernoff numbers.
A076716 counts factorizations of factorials.
A076954 can be used instead of A006939 (cf. A307895, A325337).
A181818 lists products of superprimorials, with complement A336426.
Cf. A000178, A002110, A022559, A027423, A124010, A303279, A318284, A322583, A336417, A336496.

g:= proc(n, k) option remember; uses numtheory; `if`(n>k, 0, 1)+
     `if`(isprime(n), 0, add(`if`(d>k or max(factorset(n/d))>d, 0,
        g(n/d, d)), d=divisors(n) minus {1, n}))
    end:
a:= n-> g(mul(ithprime(i)^i, i=1..n)$2):
seq(a(n), n=0..5);  # Alois P. Heinz, Jul 26 2020

chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Length[facs[chern[n]]],{n,3}] (* Gus Wiseman, Aug 21 2020 *)

\\ See A318284 for count.
a(n) = {if(n==0, 1, count(vector(n,i,i)))} \\ Andrew Howroyd, Aug 31 2020

a(n) = A317826(A033312(n+1)) = A317826((n+1)!-1) = A001055(A076954(n)).
a(n) = A001055(A006939(n)). - Gus Wiseman, Aug 21 2020
a(n) = A318284(A002110(n)). - Andrew Howroyd, Aug 31 2020

a(0)=1 prepended and a(7) added by Alois P. Heinz, Jul 26 2020

a(8)-a(13) from Andrew Howroyd, Aug 31 2020

A085082 Number of distinct prime signatures arising among the divisors of n.

Original entry on oeis.org

1, 2, 2, 3, 2, 3, 2, 4, 3, 3, 2, 5, 2, 3, 3, 5, 2, 5, 2, 5, 3, 3, 2, 7, 3, 3, 4, 5, 2, 4, 2, 6, 3, 3, 3, 6, 2, 3, 3, 7, 2, 4, 2, 5, 5, 3, 2, 9, 3, 5, 3, 5, 2, 7, 3, 7, 3, 3, 2, 7, 2, 3, 5, 7, 3, 4, 2, 5, 3, 4, 2, 9, 2, 3, 5, 5, 3, 4, 2, 9, 5, 3, 2, 7, 3, 3, 3, 7, 2, 7, 3, 5, 3, 3, 3, 11, 2, 5, 5, 6, 2, 4, 2, 7, 4

Offset: 1

Amarnath Murthy, Jul 01 2003

For a squarefree number n with k distinct prime divisors, a(n) = k+1.
If n = p^r then a(n) = tau(n) = r+1.
Question: Find a(n) in the following cases:
1. n = m^k where m is a squarefree number with r distinct prime divisors.
2. n = Product_{i=1..r} (p_i)^i, where p_i is the i-th distinct prime divisor of n.
Answers: 1. (r+k)!/(r!k!). 2. A000108(r+1). - David Wasserman, Jan 20 2005
I have submitted comments for A000108 and A016098 that each include a combinatorial statement equivalent to the second problem and its solution. - Matthew Vandermast, Nov 22 2010

			a(30) = 4 and the divisors with distinct prime signatures are 1, 2, 6 and 30. The divisors 3 and 5 with the same prime signature as of 2 and the divisors 10 and 15 with the same prime signature as that of 6 are not counted.
The divisors of 36 are 1, 2, 3, 4, 6, 9, 12 and 36. We can group them as (1), (2, 3), (6), (4, 9), (12, 18), (36) so that every group contains divisors with the same prime signature and we have a(36) = 6.

Alois P. Heinz, Table of n, a(n) for n = 1..10000

Cf. A000108.

The second problem describes A076954(i). See also A006939.

with(numtheory):
a:= n-> nops({seq(sort(map(x->x[2], ifactors(d)[2])), d=divisors(n))}):
seq(a(n), n=1..120);  # Alois P. Heinz, Jun 12 2012

ps[1] = {}; ps[n_] := FactorInteger[n][[All, 2]] // Sort; a[n_] := ps /@ Divisors[n] // Union // Length; Array[a, 120] (* Jean-François Alcover, Jun 10 2015 *)

a(n)=my(f=vecsort(factor(n)[,2]),v=[1],s); for(i=1,#f, s=0; v=vector(f[i]+1,i, if(i<=#v, s+=v[i]); s)); vecsum(v) \\ Charles R Greathouse IV, Feb 03 2017

More terms from David Wasserman, Jan 20 2005

A336420 Irregular triangle read by rows where T(n,k) is the number of divisors of the n-th superprimorial A006939(n) with distinct prime multiplicities and k prime factors counted with multiplicity.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 1, 3, 2, 5, 2, 1, 1, 1, 4, 3, 11, 7, 7, 10, 5, 2, 1, 1, 1, 5, 4, 19, 14, 18, 37, 25, 23, 15, 23, 10, 5, 2, 1, 1, 1, 6, 5, 29, 23, 33, 87, 70, 78, 74, 129, 84, 81, 49, 39, 47, 23, 10, 5, 2, 1, 1, 1, 7, 6, 41, 34, 52, 165, 144, 183, 196, 424, 317, 376, 325, 299, 431, 304, 261, 172, 129, 81, 103, 47, 23, 10, 5, 2, 1, 1

Offset: 0

Gus Wiseman, Jul 25 2020

A number's prime signature (row n of A124010) is the sequence of positive exponents in its prime factorization, so a number has distinct prime multiplicities iff all the exponents in its prime signature are distinct.
The n-th superprimorial or Chernoff number is A006939(n) = Product_{i = 1..n} prime(i)^(n - i + 1).
T(n,k) is also the number of length-n vectors 0 <= v_i <= i summing to k whose nonzero values are all distinct.

			Triangle begins:
  1
  1  1
  1  2  1  1
  1  3  2  5  2  1  1
  1  4  3 11  7  7 10  5  2  1  1
  1  5  4 19 14 18 37 25 23 15 23 10  5  2  1  1
The divisors counted in row n = 4 are:
  1  2  4     8   16   48   144   432  2160  10800  75600
     3  9    12   24   72   360   720  3024
     5  25   18   40   80   400  1008
     7       20   54  108   504  1200
             27   56  112   540  2800
             28  135  200   600
             45  189  675   756
             50            1350
             63            1400
             75            4725
            175

A000110 gives row sums.
A000124 gives row lengths.
A000142 counts divisors of superprimorials.
A006939 lists superprimorials or Chernoff numbers.
A008278 is the version counting only distinct prime factors.
A008302 counts divisors of superprimorials by bigomega.
A022915 counts permutations of prime indices of superprimorials.
A076954 can be used instead of A006939.
A130091 lists numbers with distinct prime multiplicities.
A146291 counts divisors by bigomega.
A181796 counts divisors with distinct prime multiplicities.
A181818 gives products of superprimorials.
A317829 counts factorizations of superprimorials.
A336417 counts perfect-power divisors of superprimorials.
A336498 counts divisors of factorials by bigomega.
A336499 uses factorials instead superprimorials.
Cf. A000005, A001222, A008278, A027423, A071625, A124010, A327498, A336419, A336421, A336426, A336500, A336568.

chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
Table[Length[Select[Divisors[chern[n]],PrimeOmega[#]==k&&UnsameQ@@Last/@FactorInteger[#]&]],{n,0,5},{k,0,n*(n+1)/2}]

A336417 Number of perfect-power divisors of superprimorials A006939.

Original entry on oeis.org

1, 1, 2, 5, 15, 44, 169, 652, 3106, 15286, 89933, 532476, 3698650, 25749335, 204947216, 1636097441, 14693641859, 132055603656, 1319433514898, 13186485900967, 144978145009105, 1594375302986404, 19128405558986057, 229508085926717076, 2983342885319348522

Offset: 0

Gus Wiseman, Jul 24 2020

nonn

A number is a perfect power iff it is 1 or its prime exponents (signature) are not relatively prime.

The n-th superprimorial number is A006939(n) = Product_{i = 1..n} prime(i)^(n - i + 1).

			The a(0) = 1 through a(4) = 15 divisors:
  1  2  12  360  75600
-------------------------
  1  1   1    1      1
         4    4      4
              8      8
              9      9
             36     16
                    25
                    27
                    36
                   100
                   144
                   216
                   225
                   400
                   900
                  3600

Andrew Howroyd, Table of n, a(n) for n = 0..200

A000325 is the uniform version.
A076954 can be used instead of A006939.
A336416 gives the same for factorials instead of superprimorials.
A000217 counts prime power divisors of superprimorials.
A000961 gives prime powers.
A001597 gives perfect powers, with complement A007916.
A006939 gives superprimorials or Chernoff numbers.
A022915 counts permutations of prime indices of superprimorials.
A091050 counts perfect power divisors.
A181818 gives products of superprimorials.
A294068 counts factorizations using perfect powers.
A317829 counts factorizations of superprimorials.
Cf. A000005, A027423, A090630, A118914, A124010, A203025, A251753, A327527, A336419, A336420, A336421, A336426.

chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
perpouQ[n_]:=Or[n==1,GCD@@FactorInteger[n][[All,2]]>1];
Table[Length[Select[Divisors[chern[n]],perpouQ]],{n,0,5}]

a(n) = {1 + sum(k=2, n, moebius(k)*(1 - prod(i=1, n, 1 + i\k)))} \\ Andrew Howroyd, Aug 30 2020

a(n) = A091050(A006939(n)).

a(n) = 1 + Sum_{k=2..n} mu(k)*(1 - Product_{i=1..n} 1 + floor(i/k)). - Andrew Howroyd, Aug 30 2020

Terms a(10) and beyond from Andrew Howroyd, Aug 30 2020

A336419 Number of divisors d of the n-th superprimorial A006939(n) with distinct prime exponents such that the quotient A006939(n)/d also has distinct prime exponents.

Original entry on oeis.org

1, 2, 4, 10, 24, 64, 184, 536, 1608, 5104, 16448, 55136, 187136, 658624, 2339648, 8618208, 31884640, 121733120, 468209408, 1849540416, 7342849216

Offset: 0

Gus Wiseman, Jul 25 2020

A number has distinct prime exponents iff its prime signature is strict.

The n-th superprimorial or Chernoff number is A006939(n) = Product_{i = 1..n} prime(i)^(n - i + 1).

			The a(0) = 1 through a(3) = 10 divisors:
  1  2  12  360
-----------------
  1  1   1    1
     2   3    5
         4    8
        12    9
             18
             20
             40
             45
             72
            360

A000110 shifted once to the left dominates this sequence.
A006939 lists superprimorials or Chernoff numbers.
A022915 counts permutations of prime indices of superprimorials.
A130091 lists numbers with distinct prime exponents.
A181796 counts divisors with distinct prime exponents.
A181818 gives products of superprimorials.
A317829 counts factorizations of superprimorials.
A336417 counts perfect-power divisors of superprimorials.
Cf. A000005, A000178, A008278, A071625, A076954, A118914, A124010, A327498, A327527, A336420, A336421, A336426, A336500, A336568.

chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
Table[Length[Select[Divisors[chern[n]],UnsameQ@@Last/@FactorInteger[#]&&UnsameQ@@Last/@FactorInteger[chern[n]/#]&]],{n,0,6}]

recurse(n,k,b,d)={if(k>n, 1, sum(i=0, k, if((i==0||!bittest(b,i)) && (i==k||!bittest(d,k-i)), self()(n, k+1, bitor(b, 1<Andrew Howroyd, Aug 30 2020

a(10)-a(20) from Andrew Howroyd, Aug 31 2020

A336426 Numbers that cannot be written as a product of superprimorials {2, 12, 360, 75600, ...}.

Original entry on oeis.org

3, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 25, 26, 27, 28, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76

Offset: 1

Gus Wiseman, Jul 26 2020

nonn

The n-th superprimorial is A006939(n) = Product_{i = 1..n} prime(i)^(n - i + 1).

			We have 288 = 2*12*12 so 288 is not in the sequence.

A181818 is the complement.
A336497 is the version for superfactorials.
A001055 counts factorizations.
A006939 lists superprimorials or Chernoff numbers.
A022915 counts permutations of prime indices of superprimorials.
A317829 counts factorizations of superprimorials.
A336417 counts perfect-power divisors of superprimorials.
Cf. A000325, A005117, A076954, A124010, A294068, A336419, A336420, A336421, A336496, A336500, A336568.

chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
facsusing[s_,n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[facsusing[Select[s,Divisible[n/d,#]&],n/d],Min@@#>=d&],{d,Select[s,Divisible[n,#]&]}]];
Select[Range[100],facsusing[Array[chern,30],#]=={}&]

A048742 a(n) = n! - (n-th Bell number).

Original entry on oeis.org

0, 0, 0, 1, 9, 68, 517, 4163, 36180, 341733, 3512825, 39238230, 474788003, 6199376363, 86987391878, 1306291409455, 20912309745853, 355604563226196, 6401691628921841, 121639267666626943, 2432850284018404628, 51090467301893283249, 1123996221061869232677

Offset: 0

N. J. A. Sloane

nonn

Number of permutations of [n] which have at least one cycle that has at least one inversion when written with its smallest element in the first position. Example: a(4)=9 because we have (1)(243), (1432), (142)(3), (132)(4), (1342), (1423), (1243), (143)(2) and (1324). - Emeric Deutsch, Apr 29 2008
Number of permutations of [n] having consecutive runs of increasing elements with initial elements in increasing order. a(4) = 9: `124`3, `13`24, `134`2, `14`23, `14`3`2, `2`14`3, `24`3`1, `3`14`2, `4`13`2. - Alois P. Heinz, Apr 27 2016
From Gus Wiseman, Aug 11 2020: (Start)
Also the number of divisors of the superfactorial A006939(n - 1) without distinct prime multiplicities. For example, the a(4) = 9 divisors together with their prime signatures are the following. Note that A076954 can be used here instead of A006939.
    6: (1,1)
   10: (1,1)
   15: (1,1)
   30: (1,1,1)
   36: (2,2)
   60: (2,1,1)
   90: (1,2,1)
  120: (3,1,1)
  180: (2,2,1)
(End)

Alois P. Heinz, Table of n, a(n) for n = 0..450
Index entries for sequences related to factorial numbers

A000110 lists Bell numbers.
A000142 lists factorial numbers.
A006939 lists superprimorials or Chernoff numbers.
A181796 counts divisors with distinct prime multiplicities.
A336414 counts divisors of n! with distinct prime multiplicities.
Cf. A000005, A022915, A027423, A076954, A095149, A130091, A336420, A336421.

with(combinat): seq(factorial(n)-bell(n),n=0..21); # Emeric Deutsch, Apr 29 2008

Table[n! - BellB[n], {n, 0, 40}] (* Vladimir Joseph Stephan Orlovsky, Jul 06 2011 *)

[factorial(m) - bell_number(m) for m in range(23)]  # Zerinvary Lajos, Jul 06 2008

a(n) = A000142(n) - A000110(n).

E.g.f.: 1/(1-x) - exp(exp(x)-1). - Alois P. Heinz, Apr 27 2016

A337070 Number of strict chains of divisors starting with the superprimorial A006939(n).

Original entry on oeis.org

1, 2, 16, 1208, 1383936, 32718467072, 20166949856488576, 391322675415566237681536

Offset: 0

Gus Wiseman, Aug 15 2020

The n-th superprimorial is A006939(n) = Product_{i = 1..n} prime(i)^(n - i + 1).

			The a(0) = 1 through a(2) = 16 chains:
  1  2    12
     2/1  12/1
          12/2
          12/3
          12/4
          12/6
          12/2/1
          12/3/1
          12/4/1
          12/4/2
          12/6/1
          12/6/2
          12/6/3
          12/4/2/1
          12/6/2/1
          12/6/3/1

A022915 is the maximal case.
A076954 can be used instead of A006939 (cf. A307895, A325337).
A336571 is the case with distinct prime multiplicities.
A336941 is the case ending with 1.
A337071 is the version for factorials.
A000005 counts divisors.
A000142 counts divisors of superprimorials.
A006939 lists superprimorials or Chernoff numbers.
A067824 counts chains of divisors starting with n.
A074206 counts chains of divisors from n to 1.
A253249 counts chains of divisors.
A317829 counts factorizations of superprimorials.
Cf. A001055, A002033, A167865, A181818, A336420, A336423, A336942, A337074.

chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
chnsc[n_]:=If[n==1,{{1}},Prepend[Join@@Table[Prepend[#,n]&/@chnsc[d],{d,Most[Divisors[n]]}],{n}]];
Table[Length[chnsc[chern[n]]],{n,0,3}]

a(n) = 2*A336941(n) for n > 0.

a(n) = A067824(A006939(n)).

A095149 Triangle read by rows: Aitken's array (A011971) but with a leading diagonal before it given by the Bell numbers (A000110), 1, 1, 2, 5, 15, 52, ...

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 5, 2, 3, 5, 15, 5, 7, 10, 15, 52, 15, 20, 27, 37, 52, 203, 52, 67, 87, 114, 151, 203, 877, 203, 255, 322, 409, 523, 674, 877, 4140, 877, 1080, 1335, 1657, 2066, 2589, 3263, 4140, 21147, 4140, 5017, 6097, 7432, 9089, 11155, 13744, 17007, 21147

Offset: 0

Gary W. Adamson, May 30 2004

Or, prefix Aitken's array (A011971) with a leading diagonal of 0's and take the differences of each row to get the new triangle.
With offset 1, triangle read by rows: T(n,k) is the number of partitions of the set {1,2,...,n} in which k is the largest entry in the block containing 1 (1 <= k <= n). - Emeric Deutsch, Oct 29 2006
Row term sums = the Bell numbers starting with A000110(1): 1, 2, 5, 15, ...
The k-th term in the n-th row is the number of permutations of length n starting with k and avoiding the dashed pattern 23-1. Equivalently, the number of permutations of length n ending with k and avoiding 1-32. - Andrew Baxter, Jun 13 2011
From Gus Wiseman, Aug 11 2020: (Start)
Conjecture: Also the number of divisors d with distinct prime multiplicities of the superprimorial A006939(n) that are of the form d = m * 2^k where m is odd. For example, row n = 4 counts the following divisors:
  1     2     4    8     16
  3     18    12   24    48
  5     50    20   40    80
  7     54    28   56    112
  9     1350  108  72    144
  25          540  200   400
  27          756  360   432
  45               504   720
  63               600   1008
  75               1400  1200
  135                    2160
  175                    2800
  189                    3024
  675                    10800
  4725                   75600
Equivalently, T(n,k) is the number of length-n vectors 0 <= v_i <= i whose nonzero values are distinct and such that v_n = k.
Crossrefs:
A008278 is the version counted by omega A001221.
A336420 is the version counted by Omega A001222.
A006939 lists superprimorials or Chernoff numbers.
A008302 counts divisors of superprimorials by Omega.
A022915 counts permutations of prime indices of superprimorials.
A098859 counts partitions with distinct multiplicities.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
Cf. A000005, A000142, A027423, A076954, A124010, A146291, A181818, A336417, A336419, A336421, A336499, A336942.
(End)

			Triangle starts:
   1;
   1,  1;
   2,  1,  2;
   5,  2,  3,  5;
  15,  5,  7, 10, 15;
  52, 15, 20, 27, 37, 52;
From _Gus Wiseman_, Aug 11 2020: (Start)
Row n = 3 counts the following set partitions (described in Emeric Deutsch's comment above):
  {1}{234}      {12}{34}    {123}{4}    {1234}
  {1}{2}{34}    {12}{3}{4}  {13}{24}    {124}{3}
  {1}{23}{4}                {13}{2}{4}  {134}{2}
  {1}{24}{3}                            {14}{23}
  {1}{2}{3}{4}                          {14}{2}{3}
(End)

Alois P. Heinz, Rows n = 0..150, flattened (first 51 rows from Chai Wah Wu)
Andrew M. Baxter and Lara K. Pudwell, Enumeration schemes for dashed patterns, arXiv:1108.2642 [math.CO], 2011.
Anders Claesson, Generalized pattern avoidance, Europ. J. Combin., 22 7 (2001), 961-971. See Proposition 3.
A. Bernini, M. Bouvel and L. Ferrari, Some statistics on permutations avoiding generalized patterns, PU.M.A. Vol. 18 (2007), No. 3-4, pp. 223-237 (see transposed array p. 227).

Cf. A000110, A005493, A008278, A011971, A188919, A271466.

T(2n,n) gives A020556.

with(combinat): T:=proc(n,k) if k=1 then bell(n-1) elif k=2 and n>=2 then bell(n-2) elif k<=n then add(binomial(k-2,i)*bell(n-2-i),i=0..k-2) else 0 fi end: matrix(8,8,T): for n from 1 to 11 do seq(T(n,k),k=1..n) od; # yields sequence in triangular form
Q[1]:=t*s: for n from 2 to 11 do Q[n]:=expand(t^n*subs(t=1,Q[n-1])+s*diff(Q[n-1],s)-Q[n-1]+s*Q[n-1]) od: for n from 1 to 11 do P[n]:=sort(subs(s=1,Q[n])) od: for n from 1 to 11 do seq(coeff(P[n],t,k),k=1..n) od; # yields sequence in triangular form - Emeric Deutsch, Oct 29 2006
A011971 := proc(n,k) option remember ; if k = 0 then if n=0 then 1; else A011971(n-1,n-1) ; fi ; else A011971(n,k-1)+A011971(n-1,k-1) ; fi ; end: A000110 := proc(n) option remember; if n<=1 then 1 ; else add( binomial(n-1,i)*A000110(n-1-i),i=0..n-1) ; fi ; end: A095149 := proc(n,k) option remember ; if k = 0 then A000110(n) ; else A011971(n-1,k-1) ; fi ; end: for n from 0 to 11 do for k from 0 to n do printf("%d, ",A095149(n,k)) ; od ; od ; # R. J. Mathar, Feb 05 2007
# alternative Maple program:
b:= proc(n, m, k) option remember; `if`(n=0, 1, add(
      b(n-1, max(j, m), max(k-1, -1)), j=`if`(k=0, m+1, 1..m+1)))
    end:
T:= (n, k)-> b(n, 0, n-k):
seq(seq(T(n, k), k=0..n), n=0..10);  # Alois P. Heinz, Dec 20 2018

nmax = 10; t[n_, 1] = t[n_, n_] = BellB[n-1]; t[n_, 2] = BellB[n-2]; t[n_, k_] /; n >= k >= 3 := t[n, k] = t[n, k-1] + t[n-1, k-1]; Flatten[ Table[ t[n, k], {n, 1, nmax}, {k, 1, n}]] (* Jean-François Alcover, Nov 15 2011, from formula with offset 1 *)

# requires Python 3.2 or higher.
from itertools import accumulate
A095149_list, blist = [1,1,1], [1]
for _ in range(2*10**2):
    b = blist[-1]
    blist = list(accumulate([b]+blist))
    A095149_list += [blist[-1]]+ blist
# Chai Wah Wu, Sep 02 2014, updated Chai Wah Wu, Sep 20 2014

With offset 1, T(n,1) = T(n,n) = T(n+1,2) = B(n-1) = A000110(n-1) (the Bell numbers). T(n,k) = T(n,k-1) + T(n-1,k-1) for n >= k >= 3. T(n,n-1) = B(n-1) - B(n-2) = A005493(n-3) for n >= 3 (B(q) are the Bell numbers A000110). T(n,k) = A011971(n-2,k-2) for n >= k >= 2. In other words, deleting the first row and first column we obtain Aitken's array A011971 (also called Bell or Pierce triangle; offset in A011971 is 0). - Emeric Deutsch, Oct 29 2006

T(n,1) = B(n-1); T(n,2) = B(n-2) for n >= 2; T(n,k) = Sum_{i=0..k-2} binomial(k-2,i)*B(n-2-i) for 3 <= k <= n, where B(q) are the Bell numbers (A000110). Generating polynomial of row n is P[n](t) = Q[n](t,1), where Q[n](t,s) = t^n*Q[n-1](1,s) + s*dQ[n-1](t,s)/ds + (s-1) Q[n-1](t,s); Q[1](t,s) = ts. - Emeric Deutsch, Oct 29 2006

Corrected and extended by R. J. Mathar, Feb 05 2007

Entry revised by N. J. A. Sloane, Jun 01 2005, Jun 16 2007

cp's OEIS Frontend

A008278 Reflected triangle of Stirling numbers of 2nd kind, S(n,n-k+1), n >= 1, 1 <= k <= n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Formula

Extensions

A317829 Number of set partitions of multiset {1, 2, 2, 3, 3, 3, ..., n X n}.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A085082 Number of distinct prime signatures arising among the divisors of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Extensions

A336420 Irregular triangle read by rows where T(n,k) is the number of divisors of the n-th superprimorial A006939(n) with distinct prime multiplicities and k prime factors counted with multiplicity.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A336417 Number of perfect-power divisors of superprimorials A006939.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A336419 Number of divisors d of the n-th superprimorial A006939(n) with distinct prime exponents such that the quotient A006939(n)/d also has distinct prime exponents.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

Extensions

A336426 Numbers that cannot be written as a product of superprimorials {2, 12, 360, 75600, ...}.

Views

Author

Keywords