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

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

A336415 Number of divisors of n! with equal prime multiplicities.

Original entry on oeis.org

1, 1, 2, 4, 6, 10, 13, 21, 24, 28, 33, 49, 53, 85, 94, 100, 104, 168, 173, 301, 307, 317, 334, 590, 595, 603, 636, 642, 652, 1164, 1171, 2195, 2200, 2218, 2283, 2295, 2301, 4349, 4478, 4512, 4519, 8615, 8626, 16818, 16836, 16844, 17101, 33485, 33491, 33507, 33516, 33582
Offset: 0

Views

Author

Gus Wiseman, Jul 22 2020

Keywords

Comments

A number k has "equal prime multiplicities" (or is "uniform") iff its prime signature is constant, meaning that k is a power of a squarefree number.

Examples

			The a(n) uniform divisors of n for n = 1, 2, 6, 8, 30, 36 are the columns:
  1  2  6  8  30  36
     1  3  6  15  30
        2  4  10  16
        1  3   8  15
           2   6  10
           1   5   9
               4   8
               3   6
               2   5
               1   4
                   3
                   2
                   1
In 20!, the multiplicity of the third prime (5) is 4 but the multiplicity of the fourth prime (7) is 2. Hence there are 2^3 - 1 = 3 divisors with all exponents 3 (we subtract |{1}| = 1 from that count as 1 has no exponent 3). - _David A. Corneth_, Jul 27 2020
		

Crossrefs

The version for distinct prime multiplicities is A336414.
The version for nonprime perfect powers is A336416.
Uniform partitions are counted by A047966.
Uniform numbers are A072774, with nonprime terms A182853.
Numbers with distinct prime multiplicities are A130091.
Divisors with distinct prime multiplicities are counted by A181796.
Maximum divisor with distinct prime multiplicities is A327498.
Uniform divisors are counted by A327527.
Maximum uniform divisor is A336618.
1st differences are given by A048675.

Programs

  • Mathematica
    Table[Length[Select[Divisors[n!],SameQ@@Last/@FactorInteger[#]&]],{n,0,15}]
  • PARI
    a(n) = sumdiv(n!, d, my(ex=factor(d)[,2]); (#ex==0) || (vecmin(ex) == vecmax(ex))); \\ Michel Marcus, Jul 24 2020
    
  • PARI
    a(n) = {if(n<2, return(1)); my(f = primes(primepi(n)), res = 1, t = #f); f = vector(#f, i, val(n, f[i])); for(i = 1, f[1], while(f[t] < i, t--; ); res+=(1<David A. Corneth, Jul 27 2020

Formula

a(n) = A327527(n!).

Extensions

Terms a(31) and onwards from David A. Corneth, Jul 27 2020

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

Views

Author

Keywords

Comments

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)

Crossrefs

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.

Programs

Formula

a(n) = A000142(n) - A000110(n).
E.g.f.: 1/(1-x) - exp(exp(x)-1). - Alois P. Heinz, Apr 27 2016

A336569 Number of maximal strict chains of divisors from n to 1 using elements of A130091 (numbers with distinct prime multiplicities).

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jul 29 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.

Examples

			The a(n) chains for n = 12, 72, 144, 192 (ones not shown):
  12/3    72/18/2       144/72/18/2       192/96/48/24/12/3
  12/4/2  72/18/9/3     144/72/18/9/3     192/64/32/16/8/4/2
          72/24/12/3    144/48/24/12/3    192/96/32/16/8/4/2
          72/24/8/4/2   144/72/24/12/3    192/96/48/16/8/4/2
          72/24/12/4/2  144/48/16/8/4/2   192/96/48/24/8/4/2
                        144/48/24/8/4/2   192/96/48/24/12/4/2
                        144/72/24/8/4/2
                        144/48/24/12/4/2
                        144/72/24/12/4/2
		

Crossrefs

A336423 is the non-maximal version.
A336570 is the version for chains not necessarily containing n.
A000005 counts divisors.
A001055 counts factorizations.
A001222 counts prime factors with multiplicity.
A007425 counts divisors of divisors.
A032741 counts proper divisors.
A045778 counts strict factorizations.
A071625 counts distinct prime multiplicities.
A074206 counts strict 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.
A336422 counts divisible pairs of divisors, both in A130091.
A336424 counts factorizations using A130091.
A336571 counts divisor sets of elements of A130091.

Programs

  • Mathematica
    strsigQ[n_]:=UnsameQ@@Last/@FactorInteger[n];
    fasmax[y_]:=Complement[y,Union@@(Most[Subsets[#]]&/@y)];
    strchs[n_]:=If[n==1,{{}},If[!strsigQ[n],{},Join@@Table[Prepend[#,d]&/@strchs[d],{d,Select[Most[Divisors[n]],strsigQ]}]]];
    Table[Length[fasmax[strchs[n]]],{n,100}]

A336570 Number of maximal sets of proper divisors d|n, d < n, all belonging to A130091 (numbers with distinct prime multiplicities) and forming a divisibility chain.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jul 29 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.

Examples

			The a(n) sets for n = 36, 120, 144, 180 (ones not shown):
  {2,18}    {3,12,24}    {2,18,72}       {2,18}
  {3,12}    {5,20,40}    {3,9,18,72}     {3,12}
  {2,4,12}  {2,4,8,24}   {3,12,24,48}    {5,20}
  {3,9,18}  {2,4,8,40}   {3,12,24,72}    {5,45}
            {2,4,12,24}  {2,4,8,16,48}   {2,4,12}
            {2,4,20,40}  {2,4,8,24,48}   {2,4,20}
                         {2,4,8,24,72}   {3,9,18}
                         {2,4,12,24,48}  {3,9,45}
                         {2,4,12,24,72}
		

Crossrefs

A336569 is the version for chains containing n.
A336571 is the non-maximal version.
A000005 counts divisors.
A001055 counts factorizations.
A007425 counts divisors of divisors.
A032741 counts proper divisors.
A045778 counts strict factorizations.
A071625 counts distinct prime multiplicities.
A074206 counts strict 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.
A336422 counts divisible pairs of divisors, both in A130091.
A336424 counts factorizations using A130091.
A336500 counts divisors of n in A130091 with quotient also in A130091.

Programs

  • Mathematica
    strsigQ[n_]:=UnsameQ@@Last/@FactorInteger[n];
    fasmax[y_]:=Complement[y,Union@@(Most[Subsets[#]]&/@y)];
    strses[n_]:=If[n==1,{{}},Join@@Table[Append[#,d]&/@strses[d],{d,Select[Most[Divisors[n]],strsigQ]}]];
    Table[Length[fasmax[strses[n]]],{n,100}]

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.

A336425 Number of ways to choose a divisor with distinct prime exponents of a divisor with distinct prime exponents of n!.

Original entry on oeis.org

1, 1, 3, 5, 24, 38, 132, 195, 570, 1588, 4193, 6086, 14561, 19232, 37142, 106479, 207291, 266871, 549726, 674330, 1465399, 3086598, 5939574, 7182133, 12324512, 28968994, 46819193, 82873443, 165205159, 196666406, 350397910, 406894074, 593725529, 1229814478, 1853300600, 4024414209, 6049714096, 6968090487, 9700557121, 16810076542, 26339337285
Offset: 0

Views

Author

Gus Wiseman, Aug 06 2020

Keywords

Examples

			The a(4) = 24 divisors of divisors:
  1/1  2/1  3/1  4/1  8/1  12/1   24/1
       2/2  3/3  4/2  8/2  12/2   24/2
                 4/4  8/4  12/3   24/3
                      8/8  12/4   24/4
                           12/12  24/8
                                  24/12
                                  24/24
		

Crossrefs

A336422 is the non-factorial generalization.
A130091 lists numbers with distinct prime exponents.
A181796 counts divisors with distinct prime exponents.
A327526 gives the maximum divisor of n with equal prime exponents.
A327498 gives the maximum divisor of n with distinct prime exponents.
A336414 counts divisors of n! with distinct prime exponents.
A336415 counts divisors of n! with equal prime exponents.
A336423 counts chains in A130091, with maximal version A336569.

Programs

  • Mathematica
    strsigQ[n_]:=UnsameQ@@Last/@FactorInteger[n];
    Table[Total[Cases[Divisors[n!],d_?strsigQ:>Count[Divisors[d],e_?strsigQ]]],{n,0,20}]

Extensions

Terms a(21) onward from Max Alekseyev, Nov 07 2024

A336498 Irregular triangle read by rows where T(n,k) is the number of divisors of n! with k prime factors, counted with multiplicity.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 2, 2, 2, 1, 1, 3, 4, 4, 3, 1, 1, 3, 5, 6, 6, 5, 3, 1, 1, 4, 8, 11, 12, 11, 8, 4, 1, 1, 4, 8, 11, 12, 12, 12, 12, 11, 8, 4, 1, 1, 4, 8, 12, 16, 19, 20, 20, 19, 16, 12, 8, 4, 1, 1, 4, 9, 15, 21, 26, 29, 30, 30, 29, 26, 21, 15, 9, 4, 1
Offset: 0

Views

Author

Gus Wiseman, Aug 03 2020

Keywords

Comments

Row n is row n! of A146291. Row lengths are A022559(n) + 1.

Examples

			Triangle begins:
  1
  1
  1  1
  1  2  1
  1  2  2  2  1
  1  3  4  4  3  1
  1  3  5  6  6  5  3  1
  1  4  8 11 12 11  8  4  1
  1  4  8 11 12 12 12 12 11  8  4  1
  1  4  8 12 16 19 20 20 19 16 12  8  4  1
Row n = 6 counts the following divisors:
  1  2   4   8  16   48  144  720
     3   6  12  24   72  240
     5   9  18  36   80  360
        10  20  40  120
        15  30  60  180
            45  90
Row n = 7 counts the following divisors:
  1  2   4    8   16   48   144   720  5040
     3   6   12   24   72   240  1008
     5   9   18   36   80   336  1680
     7  10   20   40  112   360  2520
        14   28   56  120   504
        15   30   60  168   560
        21   42   84  180   840
        35   45   90  252  1260
             63  126  280
             70  140  420
            105  210  630
                 315
		

Crossrefs

A000720 is column k = 1.
A008302 is the version for superprimorials.
A022559 gives row lengths minus one.
A027423 gives row sums.
A146291 is the generalization to non-factorials.
A336499 is the restriction to divisors in A130091.
A000142 lists factorial numbers.
A336415 counts uniform divisors of n!.

Programs

  • Mathematica
    Table[Length[Select[Divisors[n!],PrimeOmega[#]==k&]],{n,0,10},{k,0,PrimeOmega[n!]}]

A336616 Maximum divisor of n! with distinct prime multiplicities.

Original entry on oeis.org

1, 1, 2, 3, 24, 40, 720, 1008, 8064, 72576, 3628800, 5702400, 68428800, 80870400, 317011968, 118879488000, 1902071808000, 2487324672000, 44771844096000, 50039119872000, 1000782397440000, 21016430346240000, 5085976143790080000, 6156707963535360000
Offset: 0

Views

Author

Gus Wiseman, Jul 29 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.

Examples

			The sequence of terms together with their prime signatures begins:
             1: ()
             1: ()
             2: (1)
             3: (1)
            24: (3,1)
            40: (3,1)
           720: (4,2,1)
          1008: (4,2,1)
          8064: (7,2,1)
         72576: (7,4,1)
       3628800: (8,4,2,1)
       5702400: (8,4,2,1)
      68428800: (10,5,2,1)
      80870400: (10,5,2,1)
     317011968: (11,5,2,1)
  118879488000: (11,6,3,2,1)
		

Crossrefs

A327498 is the version not restricted to factorials, with quotient A327499.
A336414 counts these divisors.
A336617 is the quotient n!/a(n).
A336618 is the version for equal prime multiplicities.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A327526 gives the maximum divisor of n with equal prime multiplicities.
A336415 counts divisors of n! with equal prime multiplicities.

Programs

  • Mathematica
    Table[Max@@Select[Divisors[n!],UnsameQ@@Last/@If[#==1,{},FactorInteger[#]]&],{n,0,15}]
  • PARI
    a(n) = { if(n < 2, return(1)); my(pr = primes(primepi(n)), res = pr[#pr]); for(i = 1, #pr, pr[i] = [pr[i], val(n, pr[i])] ); forstep(i = #pr, 2, -1, if(pr[i][2] < pr[i-1][2], res*=pr[i-1][1]^pr[i-1][2] ) ); res }
    val(n, p) = my(r=0); while(n, r+=n\=p); r \\ David A. Corneth, Aug 25 2020

Formula

a(n) = A327498(n!).

A336617 a(n) = n!/d where d = A336616(n) is the maximum divisor of n! with distinct prime multiplicities.

Original entry on oeis.org

1, 1, 1, 2, 1, 3, 1, 5, 5, 5, 1, 7, 7, 77, 275, 11, 11, 143, 143, 2431, 2431, 2431, 221, 4199, 4199, 4199, 39083, 39083, 39083, 898909, 898909, 26068361, 26068361, 215441, 2141737, 2141737, 2141737, 66393847, 1009885357, 7953594143, 7953594143, 294282983291
Offset: 0

Views

Author

Gus Wiseman, Jul 29 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.

Examples

			The maximum divisor of 13! with distinct prime multiplicities is 80870400, so a(13) = 13!/80870400 = 77.
		

Crossrefs

A327499 is the non-factorial generalization, with quotient A327498.
A336414 counts these divisors.
A336616 is the maximum divisor d.
A336619 is the version for equal prime multiplicities.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A336415 counts divisors of n! with equal prime multiplicities.

Programs

  • Mathematica
    Table[n!/Max@@Select[Divisors[n!],UnsameQ@@Last/@If[#==1,{},FactorInteger[#]]&],{n,0,15}]

Formula

a(n) = A327499(n!).

Extensions

More terms from Jinyuan Wang, Jul 31 2020
Showing 1-10 of 23 results. Next