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 31 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

A056171 a(n) = pi(n) - pi(floor(n/2)), where pi is A000720.

Original entry on oeis.org

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

Views

Author

Labos Elemer, Jul 27 2000

Keywords

Comments

Also, the number of unitary prime divisors of n!. A prime divisor of n is unitary iff its exponent is 1 in the prime power factorization of n. In general, gcd(p, n/p) = 1 or p. Here we count the cases when gcd(p, n/p) = 1.
A unitary prime divisor of n! is >= n/2, hence their number is pi(n) - pi(n/2). - Peter Luschny, Mar 13 2011
See also the references and links mentioned in A143227. - Jonathan Sondow, Aug 03 2008
From Robert G. Wilson v, Mar 20 2017: (Start)
First occurrence of k is at n = A080359(k).
The last occurrence of k is at n = A080360(k).
The number of times k appears is A080362(k). (End)
Lev Schnirelmann proved that for every n, a(n) > (1/log_2(n))*(n/3 - 4*sqrt(n)) - 1 - (3/2)*log_2(n). - Arkadiusz Wesolowski, Nov 03 2017

Examples

			10! = 2^8 * 3^2 * 5^2 * 7. The only unitary prime divisor is 7, so a(10) = 1.
		

References

  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 214.

Crossrefs

Programs

  • Maple
    A056171 := proc(x)
         numtheory[pi](x)-numtheory[pi](floor(x/2)) ;
    end proc:
    seq(A056171(n),n=1..130) ; # N. J. A. Sloane, Sep 01 2015
    A056171 := n -> nops(select(isprime,[$iquo(n,2)+1..n])):
    seq(A056171(i),i=1..98); # Peter Luschny, Mar 13 2011
  • Mathematica
    s=0; Table[If[PrimeQ[k], s++]; If[PrimeQ[k/2], s--]; s, {k,100}]
    Table[PrimePi[n]-PrimePi[Floor[n/2]],{n,100}] (* Harvey P. Dale, Sep 01 2015 *)
  • PARI
    A056171=n->primepi(n)-primepi(n\2) \\ M. F. Hasler, Dec 31 2016
    
  • Python
    from sympy import primepi
    [primepi(n) - primepi(n//2) for n in range(1,151)] # Indranil Ghosh, Mar 22 2017
    
  • Sage
    [prime_pi(n)-prime_pi(floor(n/2)) for n in range(1,99)] # Stefano Spezia, Apr 22 2025

Formula

a(n) = A000720(n) - A056172(n). - Robert G. Wilson v, Apr 09 2017
a(n) = A056169(n!). - Amiram Eldar, Jul 24 2024

Extensions

Definition simplified by N. J. A. Sloane, Sep 01 2015

A336414 Number of divisors of n! with distinct prime multiplicities.

Original entry on oeis.org

1, 1, 2, 3, 7, 10, 20, 27, 48, 86, 147, 195, 311, 390, 595, 1031, 1459, 1791, 2637, 3134, 4747, 7312, 10766, 12633, 16785, 26377, 36142, 48931, 71144, 82591, 112308, 128023, 155523, 231049, 304326, 459203, 568095, 642446, 812245, 1137063, 1441067, 1612998, 2193307, 2429362
Offset: 0

Views

Author

Gus Wiseman, Jul 22 2020

Keywords

Comments

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

Examples

			The first and second columns below are the a(6) = 20 counted divisors of 6! together with their prime signatures. The third column shows the A000005(6!) - a(6) = 10 remaining divisors.
      1: ()      20: (2,1)    |    6: (1,1)
      2: (1)     24: (3,1)    |   10: (1,1)
      3: (1)     40: (3,1)    |   15: (1,1)
      4: (2)     45: (2,1)    |   30: (1,1,1)
      5: (1)     48: (4,1)    |   36: (2,2)
      8: (3)     72: (3,2)    |   60: (2,1,1)
      9: (2)     80: (4,1)    |   90: (1,2,1)
     12: (2,1)  144: (4,2)    |  120: (3,1,1)
     16: (4)    360: (3,2,1)  |  180: (2,2,1)
     18: (1,2)  720: (4,2,1)  |  240: (4,1,1)
		

Crossrefs

Perfect-powers are A001597, with complement A007916.
Numbers with distinct prime multiplicities are A130091.
Divisors with distinct prime multiplicities are counted by A181796.
The maximum divisor with distinct prime multiplicities is A327498.
Divisors of n! with equal prime multiplicities are counted by A336415.

Programs

  • Mathematica
    Table[Length[Select[Divisors[n!],UnsameQ@@Last/@FactorInteger[#]&]],{n,0,15}]
  • PARI
    a(n) = sumdiv(n!, d, my(ex=factor(d)[,2]); #vecsort(ex,,8) == #ex); \\ Michel Marcus, Jul 24 2020

Formula

a(n) = A181796(n!).

Extensions

a(21)-a(41) from Alois P. Heinz, Jul 24 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

A056172 Number of non-unitary prime divisors of n!.

Original entry on oeis.org

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

Views

Author

Labos Elemer, Jul 27 2000

Keywords

Comments

A non-unitary prime divisor for n! cannot exceed n/2.

Examples

			10! = 2^8 * 3^4 * 5^2 * 7. The non-unitary prime divisors are 2, 3, and 5 because their exponents exceed 1, so a(10) = 3.  The only unitary prime divisor of 10! is 7.
		

Crossrefs

Programs

Formula

a(n) = pi(n/2).
A prime divisor of x is non-unitary iff its exponent is at least 2 in the prime power factorization of x. In general, GCD(p, x/p) = 1 or p. Cases are counted when GCD(p, n/p) > 1.
a(n) = A000720(n) - A056171(n). - Robert G. Wilson v, Apr 09 2017
a(n) = A056170(n!). - Amiram Eldar, Jul 24 2024

Extensions

Example corrected by Jon E. Schoenfield, Sep 30 2013

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