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-7 of 7 results.

A337105 Number of strict chains of divisors from n! to 1.

Original entry on oeis.org

1, 1, 1, 3, 20, 132, 1888, 20128, 584000, 17102016, 553895936, 11616690176, 743337949184, 19467186157568, 999551845713920, 66437400489711616, 10253161206302064640, 388089999627661557760, 53727789519052432998400, 2325767421950553303285760, 365546030278816140131041280
Offset: 0

Views

Author

Gus Wiseman, Aug 17 2020

Keywords

Examples

			The a(4) = 20 chains:
  24/1  24/2/1   24/4/2/1   24/8/4/2/1
        24/3/1   24/6/2/1   24/12/4/2/1
        24/4/1   24/6/3/1   24/12/6/2/1
        24/6/1   24/8/2/1   24/12/6/3/1
        24/8/1   24/8/4/1
        24/12/1  24/12/2/1
                 24/12/3/1
                 24/12/4/1
                 24/12/6/1
		

Crossrefs

A325617 is the maximal case.
A336941 is the version for superprimorials.
A337104 counts the case with distinct prime multiplicities.
A337071 is the case not necessarily ending with 1.
A000005 counts divisors.
A000142 lists factorial numbers.
A001055 counts factorizations.
A027423 counts divisors of factorial numbers.
A067824 counts chains of divisors starting with n.
A074206 counts chains of divisors from n to 1.
A076716 counts factorizations of factorial numbers.
A253249 counts chains of divisors.
A336423 counts chains using A130091, with maximal case A336569.
A336942 counts chains using A130091 from A006939(n) to 1.

Programs

  • Maple
    b:= proc(n) option remember; 1 +
          add(b(d), d=numtheory[divisors](n) minus {n})
        end:
    a:= n-> ceil(b(n!)/2):
    seq(a(n), n=0..14);  # Alois P. Heinz, Aug 23 2020
  • Mathematica
    chnsc[n_]:=Prepend[Join@@Table[Prepend[#,n]&/@chnsc[d],{d,DeleteCases[Divisors[n],1|n]}],{n}];
    Table[Length[chnsc[n!]],{n,0,5}]

Formula

a(n) = A337071(n)/2 for n > 1.
a(n) = A074206(n!).

Extensions

a(19)-a(20) from Alois P. Heinz, Aug 22 2020

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

Views

Author

Gus Wiseman, Aug 15 2020

Keywords

Comments

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

Examples

			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
		

Crossrefs

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.

Programs

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

Formula

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

A337071 Number of strict chains of divisors starting with n!.

Original entry on oeis.org

1, 1, 2, 6, 40, 264, 3776, 40256, 1168000, 34204032, 1107791872, 23233380352, 1486675898368, 38934372315136, 1999103691427840, 132874800979423232, 20506322412604129280, 776179999255323115520, 107455579038104865996800, 4651534843901106606571520, 731092060557632280262082560
Offset: 0

Views

Author

Gus Wiseman, Aug 16 2020

Keywords

Examples

			The a(1) = 1 through a(3) = 6 chains:
  1  2    6
     2/1  6/1
          6/2
          6/3
          6/2/1
          6/3/1
The a(4) = 40 chains:
  24  24/1   24/2/1   24/4/2/1   24/8/4/2/1
      24/2   24/3/1   24/6/2/1   24/12/4/2/1
      24/3   24/4/1   24/6/3/1   24/12/6/2/1
      24/4   24/4/2   24/8/2/1   24/12/6/3/1
      24/6   24/6/1   24/8/4/1
      24/8   24/6/2   24/8/4/2
      24/12  24/6/3   24/12/2/1
             24/8/1   24/12/3/1
             24/8/2   24/12/4/1
             24/8/4   24/12/4/2
             24/12/1  24/12/6/1
             24/12/2  24/12/6/2
             24/12/3  24/12/6/3
             24/12/4
             24/12/6
		

Crossrefs

A325617 is the maximal case.
A337070 is the version for superprimorials.
A337074 counts the case with distinct prime multiplicities.
A337105 is the case ending with one.
A000005 counts divisors.
A000142 lists factorial numbers.
A027423 counts divisors of factorial numbers.
A067824 counts chains of divisors starting with n.
A074206 counts chains of divisors from n to 1.
A076716 counts factorizations of factorial numbers.
A253249 counts chains of divisors.

Programs

  • Mathematica
    chnsc[n_]:=Prepend[Join@@Table[Prepend[#,n]&/@chnsc[d],{d,Most[Divisors[n]]}],{n}];
    Table[Length[chnsc[n!]],{n,0,5}]

Formula

a(n) = 2*A337105(n) for n > 1.
a(n) = A067824(n!).

Extensions

a(19)-a(20) from Alois P. Heinz, Aug 23 2020

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

A336941 Number of strict chains of divisors starting with the superprimorial A006939(n) and ending with 1.

Original entry on oeis.org

1, 1, 8, 604, 691968, 16359233536, 10083474928244288, 195661337707783118840768, 139988400203593571474134024847360, 4231553868972506381329450624389969130848256, 6090860257621637852755610879241895108657182173073604608, 464479854191019594417264488167571483344961210693790188774166838214656
Offset: 0

Views

Author

Gus Wiseman, Aug 13 2020

Keywords

Examples

			The a(2) = 8 chains:
  12/1
  12/2/1
  12/3/1
  12/4/1
  12/6/1
  12/4/2/1
  12/6/2/1
  12/6/3/1
		

Crossrefs

A022915 is the maximal case.
A076954 can be used instead of A006939.
A336571 is the case with distinct prime multiplicities.
A336942 is the case using members of A130091.
A337070 is the version ending with any divisor of A006939(n).
A000005 counts divisors.
A074206 counts chains of divisors from n to 1.
A006939 lists superprimorials or Chernoff numbers.
A067824 counts divisor chains starting with n.
A181818 gives products of superprimorials, with complement A336426.
A253249 counts chains of divisors.
A317829 counts factorizations of superprimorials.
A336423 counts chains using A130091, with maximal case A336569.

Programs

  • Mathematica
    chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
    chns[n_]:=If[n==1,1,Sum[chns[d],{d,Most[Divisors[n]]}]];
    Table[chns[chern[n]],{n,0,3}]
  • PARI
    a(n)={my(sig=vector(n,i,i), m=vecsum(sig)); sum(k=0, m, prod(i=1, #sig, binomial(sig[i]+k-1, k-1))*sum(r=k, m, binomial(r,k)*(-1)^(r-k)))} \\ Andrew Howroyd, Aug 30 2020

Formula

a(n) = A337070(n)/2 for n > 0.
a(n) = A074206(A006939(n)).

Extensions

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

A337074 Number of strict chains of divisors in A130091 (numbers with distinct prime multiplicities), starting with n!.

Original entry on oeis.org

1, 1, 2, 0, 28, 0, 768, 0, 0, 0, 42155360, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Gus Wiseman, Aug 16 2020

Keywords

Comments

Support appears to be {0, 1, 2, 4, 6, 10}.

Examples

			The a(4) = 28 chains:
  24  24/1   24/2/1   24/4/2/1   24/8/4/2/1
      24/2   24/3/1   24/8/2/1   24/12/4/2/1
      24/3   24/4/1   24/8/4/1
      24/4   24/4/2   24/8/4/2
      24/8   24/8/1   24/12/2/1
      24/12  24/8/2   24/12/3/1
             24/8/4   24/12/4/1
             24/12/1  24/12/4/2
             24/12/2
             24/12/3
             24/12/4
		

Crossrefs

A336867 is the complement of the support.
A336868 is the characteristic function (image under A057427).
A336942 is half the version for superprimorials (n > 1).
A337071 does not require distinct prime multiplicities.
A337104 is the case of chains ending with 1.
A000005 counts divisors.
A000142 lists factorial numbers.
A027423 counts divisors of factorial numbers.
A067824 counts chains of divisors starting with n.
A074206 counts chains of divisors from n to 1.
A076716 counts factorizations of factorial numbers.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A253249 counts chains of divisors.
A327498 gives the maximum divisor with distinct prime multiplicities.
A336414 counts divisors of n! with distinct prime multiplicities.
A336415 counts divisors of n! with equal prime multiplicities.
A336423 counts chains using A130091, with maximal case A336569.
A336571 counts chains of divisors 1 < d < n using A130091.

Programs

  • Mathematica
    chnsc[n_]:=If[!UnsameQ@@Last/@FactorInteger[n],{},If[n==1,{{1}},Prepend[Join@@Table[Prepend[#,n]&/@chnsc[d],{d,Most[Divisors[n]]}],{n}]]];
    Table[Length[chnsc[n!]],{n,0,6}]

Formula

a(n) = 2*A337104(n) = 2*A336423(n!) for n > 1.

A337104 Number of strict chains of divisors from n! to 1 using terms of A130091 (numbers with distinct prime multiplicities).

Original entry on oeis.org

1, 1, 1, 0, 14, 0, 384, 0, 0, 0, 21077680, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Gus Wiseman, Aug 17 2020

Keywords

Comments

The support appears to be {0, 1, 2, 4, 6, 10}.

Examples

			The a(4) = 14 chains:
  24/1
  24/2/1
  24/3/1
  24/4/1
  24/8/1
  24/12/1
  24/4/2/1
  24/8/2/1
  24/8/4/1
  24/12/2/1
  24/12/3/1
  24/12/4/1
  24/8/4/2/1
  24/12/4/2/1
		

Crossrefs

A336867 appears to be the positions of zeros.
A336868 is the characteristic function (image under A057427).
A336942 is the version for superprimorials (n > 1).
A337105 does not require distinct prime multiplicities.
A337074 does not require chains to end with 1.
A337075 is the version for chains not containing n!.
A000005 counts divisors.
A000142 lists factorial numbers.
A001055 counts factorizations.
A027423 counts divisors of factorial numbers.
A067824 counts chains of divisors starting with n.
A074206 counts chains of divisors from n to 1.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A253249 counts chains of divisors.
A327498 gives the maximum divisor with distinct prime multiplicities.
A336414 counts divisors of n! with distinct prime multiplicities.
A336423 counts chains using A130091, with maximal case A336569.
A336425 counts divisible pairs of divisors of n!, both in A130091.
A336571 counts chains of divisors 1 < d < n using A130091.
A337071 counts chains of divisors starting with n!.

Programs

  • Mathematica
    strchns[n_]:=If[n==1,1,If[!UnsameQ@@Last/@FactorInteger[n],0,Sum[strchns[d],{d,Select[DeleteCases[Divisors[n],n],UnsameQ@@Last/@FactorInteger[#]&]}]]];
    Table[strchns[n!],{n,0,8}]

Formula

a(n) = A337075(n) whenever A337075(n) != 0.
a(n) = A337074(n)/2 for n > 1.
a(n) = A336423(n!).
Showing 1-7 of 7 results.