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

A006939 Chernoff sequence: a(n) = Product_{k=1..n} prime(k)^(n-k+1).

Original entry on oeis.org

1, 2, 12, 360, 75600, 174636000, 5244319080000, 2677277333530800000, 25968760179275365452000000, 5793445238736255798985527240000000, 37481813439427687898244906452608585200000000, 7517370874372838151564668004911177464757864076000000000, 55784440720968513813368002533861454979548176771615744085560000000000
Offset: 0

Views

Author

Keywords

Comments

Product of first n primorials: a(n) = Product_{i=1..n} A002110(i).
Superprimorials, from primorials by analogy with superfactorials.
Smallest number k with n distinct exponents in its prime factorization, i.e., A071625(k) = n.
Subsequence of A130091. - Reinhard Zumkeller, May 06 2007
Hankel transform of A171448. - Paul Barry, Dec 09 2009
This might be a good place to explain the name "Chernoff sequence" since his name does not appear in the References or Links as of Mar 22 2014. - Jonathan Sondow, Mar 22 2014
Pickover (1992) named this sequence after Paul Chernoff of California, who contributed this sequence to his book. He was possibly referring to American mathematician Paul Robert Chernoff (1942 - 2017), a professor at the University of California. - Amiram Eldar, Jul 27 2020

Examples

			a(4) = 360 because 2^3 * 3^2 * 5 = 1 * 2 * 6 * 30 = 360.
a(5) = 75600 because 2^4 * 3^3 * 5^2 * 7 = 1 * 2 * 6 * 30 * 210 = 75600.
		

References

  • Clifford A. Pickover, Mazes for the Mind, St. Martin's Press, NY, 1992, p. 351.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • James K. Strayer, Elementary number theory, Waveland Press, Inc., Long Grove, IL, 1994. See p. 37.

Crossrefs

Cf. A000178 (product of first n factorials), A007489 (sum of first n factorials), A060389 (sum of first n primorials).
A000142 counts divisors of superprimorials.
A000325 counts uniform divisors of superprimorials.
A008302 counts divisors of superprimorials by bigomega.
A022915 counts permutations of prime indices of superprimorials.
A076954 is a sister-sequence.
A118914 has row a(n) equal to {1..n}.
A124010 has row a(n) equal to {n..1}.
A130091 lists numbers with distinct prime multiplicities.
A317829 counts factorizations of superprimorials.
A336417 counts perfect-power divisors of superprimorials.
A336426 gives non-products of superprimorials.

Programs

  • Haskell
    a006939 n = a006939_list !! n
    a006939_list = scanl1 (*) a002110_list -- Reinhard Zumkeller, Jul 21 2012
    
  • Magma
    [1] cat [(&*[NthPrime(k)^(n-k+1): k in [1..n]]): n in [1..15]]; // G. C. Greubel, Oct 14 2018
    
  • Maple
    a := []; printlevel := -1; for k from 0 to 20 do a := [op(a),product(ithprime(i)^(k-i+1),i=1..k)] od; print(a);
  • Mathematica
    Rest[FoldList[Times,1,FoldList[Times,1,Prime[Range[15]]]]] (* Harvey P. Dale, Jul 07 2011 *)
    Table[Times@@Table[Prime[i]^(n - i + 1), {i, n}], {n, 12}] (* Alonso del Arte, Sep 30 2011 *)
  • PARI
    a(n)=prod(k=1,n,prime(k)^(n-k+1)) \\ Charles R Greathouse IV, Jul 25 2011
    
  • Python
    from math import prod
    from sympy import prime
    def A006939(n): return prod(prime(k)**(n-k+1) for k in range(1,n+1)) # Chai Wah Wu, Aug 12 2025

Formula

a(n) = m(1)*m(2)*m(3)*...*m(n), where m(n) = n-th primorial number. - N. J. A. Sloane, Feb 20 2005
a(0) = 1, a(n) = a(n - 1)p(n)#, where p(n)# is the n-th primorial A002110(n) (the product of the first n primes). - Alonso del Arte, Sep 30 2011
log a(n) = n^2(log n + log log n - 3/2 + o(1))/2. - Charles R Greathouse IV, Mar 14 2011
A181796(a(n)) = A000110(n+1). It would be interesting to have a bijective proof of this theorem, which is stated at A181796 without proof. See also A336420. - Gus Wiseman, Aug 03 2020

Extensions

Corrected and extended by Labos Elemer, May 30 2001

A022915 Multinomial coefficients (0, 1, ..., n)! = C(n+1,2)!/(0!*1!*2!*...*n!).

Original entry on oeis.org

1, 1, 3, 60, 12600, 37837800, 2053230379200, 2431106898187968000, 73566121315513295589120000, 65191584694745586153436251091200000, 1906765806522767212441719098019963758016000000, 2048024348726152339387799085049745725891853852479488000000
Offset: 0

Views

Author

Keywords

Comments

Number of ways to put numbers 1, 2, ..., n*(n+1)/2 in a triangular array of n rows in such a way that each row is increasing. Also number of ways to choose groups of 1, 2, 3, ..., n-1 and n objects out of n*(n+1)/2 objects. - Floor van Lamoen, Jul 16 2001
a(n) is the number of ways to linearly order the multiset {1,2,2,3,3,3,...n,n,...n}. - Geoffrey Critzer, Mar 08 2009
Also the number of distinct adjacency matrices in the n-triangular honeycomb rook graph. - Eric W. Weisstein, Jul 14 2017

Examples

			From _Gus Wiseman_, Aug 12 2020: (Start)
The a(3) = 60 permutations of the prime indices of A006939(3) = 360:
  (111223)  (121123)  (131122)  (212113)  (231211)
  (111232)  (121132)  (131212)  (212131)  (232111)
  (111322)  (121213)  (131221)  (212311)  (311122)
  (112123)  (121231)  (132112)  (213112)  (311212)
  (112132)  (121312)  (132121)  (213121)  (311221)
  (112213)  (121321)  (132211)  (213211)  (312112)
  (112231)  (122113)  (211123)  (221113)  (312121)
  (112312)  (122131)  (211132)  (221131)  (312211)
  (112321)  (122311)  (211213)  (221311)  (321112)
  (113122)  (123112)  (211231)  (223111)  (321121)
  (113212)  (123121)  (211312)  (231112)  (321211)
  (113221)  (123211)  (211321)  (231121)  (322111)
(End)
		

Crossrefs

A190945 counts the case of anti-run permutations.
A317829 counts partitions of this multiset.
A325617 is the version for factorials instead of superprimorials.
A006939 lists superprimorials or Chernoff numbers.
A008480 counts permutations of prime indices.
A181818 gives products of superprimorials, with complement A336426.

Programs

  • Maple
    with(combinat):
    a:= n-> multinomial(binomial(n+1, 2), $0..n):
    seq(a(n), n=0..12);  # Alois P. Heinz, May 18 2013
  • Mathematica
    Table[Apply[Multinomial ,Range[n]], {n, 0, 20}]  (* Geoffrey Critzer, Dec 09 2012 *)
    Table[Multinomial @@ Range[n], {n, 0, 20}] (* Eric W. Weisstein, Jul 14 2017 *)
    Table[Binomial[n + 1, 2]!/BarnesG[n + 2], {n, 0, 20}] (* Eric W. Weisstein, Jul 14 2017 *)
    Table[Length[Permutations[Join@@Table[i,{i,n},{i}]]],{n,0,4}] (* Gus Wiseman, Aug 12 2020 *)
  • PARI
    a(n) = binomial(n+1,2)!/prod(k=1, n, k^(n+1-k)); \\ Michel Marcus, May 02 2019

Formula

a(n) = (n*(n+1)/2)!/(0!*1!*2!*...*n!).
a(n) = a(n-1) * A014068(n). - Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Apr 08 2001.
a(n) = A052295(n)/A000178(n). - Lekraj Beedassy, Feb 19 2004
a(n) = A208437(n*(n+1)/2,n). - Alois P. Heinz, Apr 08 2016
a(n) ~ A * exp(n^2/4 + n + 1/6) * n^(n^2/2 + 7/12) / (2^((n+1)^2/2) * Pi^(n/2)), where A is the Glaisher-Kinkelin constant A074962. - Vaclav Kotesovec, May 02 2019
a(n) = A327803(n*(n+1)/2,n). - Alois P. Heinz, Sep 25 2019
a(n) = A008480(A006939(n)). - Gus Wiseman, Aug 12 2020

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Apr 11 2001
More terms from Michel ten Voorde, Apr 12 2001
Better definition from L. Edson Jeffery, May 18 2013

A376562 Second differences of consecutive non-perfect-powers (A007916). First differences of A375706.

Original entry on oeis.org

1, -1, 0, 2, -2, 0, 0, 0, 0, 1, -1, 0, 0, 0, 0, 0, 0, 1, 0, -1, 0, 0, 1, -1, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Oct 01 2024

Keywords

Comments

Non-perfect-powers (A007916) are numbers without a proper integer root.

Examples

			The non-perfect powers (A007916) are:
  2, 3, 5, 6, 7, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 24, 26, 28, ...
with first differences (A375706):
  1, 2, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, ...
with first differences (A376562):
  1, -1, 0, 2, -2, 0, 0, 0, 0, 1, -1, 0, 0, 0, 0, 0, 0, 1, 0, -1, 0, 0, 1, -1, 0, ...
		

Crossrefs

The version for A000002 is A376604, first differences of A054354.
For first differences we had A375706, ones A375740, complement A375714.
Positions of zeros are A376588, complement A376589.
Runs of non-perfect-powers:
- length: A375702 = A053289(n+1) - 1
- first: A375703 (same as A216765 with 2 exceptions)
- last: A375704 (same as A045542 with 8 removed)
- sum: A375705
A000961 lists prime-powers inclusive, exclusive A246655.
A007916 lists non-perfect-powers, complement A001597.
A112344 counts integer partitions into perfect-powers, factorizations A294068.
A333254 gives run-lengths of differences between consecutive primes.
For non-perfect-powers: A375706 (first differences), A376588 (inflections and undulations), A376589 (nonzero curvature).
For second differences: A036263 (prime), A073445 (composite), A376559 (perfect-power), A376590 (squarefree), A376593 (nonsquarefree), A376596 (prime-power inclusive), A376599 (non-prime-power inclusive).

Programs

  • Mathematica
    radQ[n_]:=n>1&&GCD@@Last/@FactorInteger[n]==1;
    Differences[Select[Range[100],radQ],2]
  • Python
    from itertools import count
    from sympy import mobius, integer_nthroot, perfect_power
    def A376562(n):
        def f(x): return int(n+1-sum(mobius(k)*(integer_nthroot(x, k)[0]-1) for k in range(2, x.bit_length())))
        m, k = n, f(n)
        while m != k: m, k = k, f(k)
        r = m+((k:=next(i for i in count(1) if not perfect_power(m+i)))<<1)
        return next(i for i in count(1-k) if not perfect_power(r+i)) # Chai Wah Wu, Oct 02 2024

A336416 Number of perfect-power divisors of n!.

Original entry on oeis.org

1, 1, 1, 1, 3, 3, 7, 7, 11, 18, 36, 36, 47, 47, 84, 122, 166, 166, 221, 221, 346, 416, 717, 717, 1001, 1360, 2513, 2942, 4652, 4652, 5675, 5675, 6507, 6980, 13892, 17212, 20408, 20408, 39869, 45329, 51018, 51018, 68758, 68758, 105573, 138617, 284718, 284718, 338126, 421126
Offset: 0

Views

Author

Gus Wiseman, Jul 22 2020

Keywords

Comments

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

Examples

			The a(1) = 0 through a(9) = 18 divisors:
       1: 1
       2: 1
       6: 1
      24: 1,4,8
     120: 1,4,8
     720: 1,4,8,9,16,36,144
    5040: 1,4,8,9,16,36,144
   40320: 1,4,8,9,16,32,36,64,128,144,576
  362880: 1,4,8,9,16,27,32,36,64,81,128,144,216,324,576,1296,1728,5184
		

Crossrefs

The maximum among these divisors is A090630, with quotient A251753.
The version for distinct prime exponents is A336414.
The uniform version is A336415.
Replacing factorials with Chernoff numbers (A006939) gives A336417.
Prime powers are A000961.
Perfect powers are A001597, with complement A007916.
Prime power divisors are counted by A022559.

Programs

  • Mathematica
    perpouQ[n_]:=Or[n==1,GCD@@FactorInteger[n][[All,2]]>1];
    Table[Length[Select[Divisors[n!],perpouQ]],{n,0,15}]
  • PARI
    a(n) = sumdiv(n!, d, (d==1) || ispower(d)); \\ Michel Marcus, Aug 19 2020
    
  • PARI
    addhelp(val, "exponent of prime p in n!")
    val(n, p) = my(r=0); while(n, r+=n\=p);r
    a(n) = {if(n<=3, return(1)); my(pr = primes(primepi(n\2)), v = vector(#pr, i, val(n, pr[i])), res = 1, cv); for(i = 2, v[1], if(issquarefree(i), cv = v\i; res-=(prod(i = 1, #cv, cv[i]+1)-1)*(-1)^omega(i) ) ); res } \\ David A. Corneth, Aug 19 2020

Formula

a(p) = a(p-1) for prime p. - David A. Corneth, Aug 19 2020

Extensions

a(26)-a(34) from Jinyuan Wang, Aug 19 2020
a(35)-a(49) from David A. Corneth, Aug 19 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

Views

Author

Antti Karttunen, Aug 10 2018

Keywords

Comments

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

Examples

			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.
		

Crossrefs

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.

Programs

  • Maple
    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
  • Mathematica
    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 *)
  • PARI
    \\ See A318284 for count.
    a(n) = {if(n==0, 1, count(vector(n,i,i)))} \\ Andrew Howroyd, Aug 31 2020

Formula

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

Extensions

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

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

Views

Author

Gus Wiseman, Jul 25 2020

Keywords

Comments

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.

Examples

			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
		

Crossrefs

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.

Programs

  • Mathematica
    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}]

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

Views

Author

Gus Wiseman, Jul 25 2020

Keywords

Comments

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).

Examples

			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
		

Crossrefs

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.

Programs

  • Mathematica
    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}]
  • PARI
    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

Extensions

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

Views

Author

Gus Wiseman, Jul 26 2020

Keywords

Comments

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

Examples

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

Crossrefs

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.

Programs

  • Mathematica
    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],#]=={}&]

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

Views

Author

Gary W. Adamson, May 30 2004

Keywords

Comments

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.
(End)

Examples

			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)
		

Crossrefs

Programs

  • Maple
    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
  • Mathematica
    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 *)
  • Python
    # 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

Formula

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

Extensions

Corrected and extended by R. J. Mathar, Feb 05 2007
Entry revised by N. J. A. Sloane, Jun 01 2005, Jun 16 2007

A336421 Number of ways to choose a divisor of a divisor, both having distinct prime exponents, of the n-th superprimorial number A006939(n).

Original entry on oeis.org

1, 3, 13, 76, 571, 5309, 59341, 780149
Offset: 0

Views

Author

Gus Wiseman, Jul 25 2020

Keywords

Comments

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).

Examples

			The a(2) = 13 ways:
  12/1/1  12/2/1  12/3/1  12/4/1  12/12/1
          12/2/2  12/3/3  12/4/2  12/12/2
                          12/4/4  12/12/3
                                  12/12/4
                                  12/12/12
		

Crossrefs

A000258 shifted once to the left is dominated by this sequence.
A336422 is the generalization to non-superprimorials.
A000110 counts divisors of superprimorials with distinct prime exponents.
A006939 lists superprimorials or Chernoff numbers.
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 exponents.
A181796 counts divisors with distinct prime exponents.
A181818 gives products of superprimorials.
A317829 counts factorizations of superprimorials.

Programs

  • Mathematica
    chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
    strsig[n_]:=UnsameQ@@Last/@FactorInteger[n];
    Table[Total[Cases[Divisors[chern[n]],d_?strsig:>Count[Divisors[d],e_?strsig]]],{n,0,5}]
Showing 1-10 of 14 results. Next