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

A368682 Products of primorials that are perfect powers but not prime powers.

Original entry on oeis.org

36, 144, 216, 576, 900, 1296, 1728, 2304, 3600, 5184, 7776, 9216, 13824, 14400, 20736, 27000, 32400, 36864, 44100, 46656, 57600, 82944, 110592, 129600, 147456, 176400, 186624, 216000, 230400, 248832, 279936, 331776, 373248, 518400, 589824, 705600, 746496, 810000
Offset: 1

Views

Author

Michael De Vlieger, Jan 02 2024

Keywords

Comments

Intersection of A025487 and A131605.
Proper subset of A286708.
Contains A365308 (perfect powers of composite primorials) and A368508 (perfect powers of composite superprimorials).
These numbers are perfect powers of some smaller product of primorials.

Examples

			b(n) = A025487(n).
a(1) = b(11) = 36 = 6^2 = b(4)^2,
a(2) = b(19) = 144 = 12^2 = b(6)^2,
a(3) = b(23) = 216 = 6^3 = b(4)^3,
a(4) = b(33) = 576 = 24^2 = b(8)^2,
a(5) = b(38) = 900 = 30^2 = b(9)^2, etc.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[36, 2^18, 2], And[Union@ Differences@ PrimePi@ #1 == {1}, AllTrue[Union@ Differences@ #2, # <= 0 &], GCD @@ #2 > 1] & @@ Transpose@ FactorInteger[#] &]

Formula

This sequence is { A368681 \ A000079 }.

A380446 Perfect powers k^m, m > 1, omega(k) > 1, such that A053669(k) > A006530(k), where omega = A001221.

Original entry on oeis.org

36, 144, 216, 324, 576, 900, 1296, 1728, 2304, 2916, 3600, 5184, 5832, 7776, 8100, 9216, 11664, 13824, 14400, 20736, 22500, 26244, 27000, 32400, 36864, 44100, 46656, 57600, 72900, 82944, 90000, 104976, 110592, 129600, 147456, 157464, 176400, 186624, 202500, 216000
Offset: 1

Views

Author

Michael De Vlieger, Jul 25 2025

Keywords

Comments

Perfect powers k^m, m > 1, for k in A055932.
Union of {k^m : rad(k) | P(i), m >= 2}, rad = A007947, P = A002110. Therefore perfect powers in A033845, A143207, A147571, A147572, etc. are proper subsets.
Terms are even. For a(n) such that omega(a(n)) > 2, a(n) mod 10 = 0, where omega = A001221.

Examples

			Table of n, a(n) for select n, showing exponents m of prime power factors p^m | a(n) for primes p listed in the heading. Terms that also appear in A368682 are marked by "#":
                         Exponents
 n      a(n)             2.3.5.7.11
-----------------------------------
 1       36 =    6^2  #  2.2
 2      144 =   12^2  #  4.2
 3      216 =    6^3  #  3.3
 4      324 =   18^2     2.4
 5      576 =   24^2  #  6.2
 6      900 =   30^2  #  2.2.2
 7     1296 =    6^4  #  4.4
 8     1728 =   12^3  #  6.3
 9     2304 =   48^2  #  8.2
10     2916 =   54^2     2.6
11     3600 =   60^2  #  4.2.2
12     5184 =   72^2  #  6.4
26    44100 =  210^2  #  2.2.2.2
90  5336100 = 2310^2  #  2.2.2.2.2
		

Crossrefs

Programs

  • Mathematica
    (* Load linked Mathematica algorithm, then: *)
    Select[Union@ Flatten[a055932[7][[3 ;; -1, 2 ;; -1]] ], And[Divisible[#1, Apply[Times, #2[[All, 1]] ]^2], GCD @@ #2[[All, -1]] > 1] & @@ {#, FactorInteger[#]} &]

Formula

Intersection of A131605 and A055932 = A304250 \ A246547.

A368681 Products of primorials that are perfect powers.

Original entry on oeis.org

1, 4, 8, 16, 32, 36, 64, 128, 144, 216, 256, 512, 576, 900, 1024, 1296, 1728, 2048, 2304, 3600, 4096, 5184, 7776, 8192, 9216, 13824, 14400, 16384, 20736, 27000, 32400, 32768, 36864, 44100, 46656, 57600, 65536, 82944, 110592, 129600, 131072, 147456, 176400, 186624
Offset: 1

Views

Author

Michael De Vlieger, Jan 02 2024

Keywords

Comments

Intersection of A025487 and A001597.
Contains A365308 (perfect powers of composite primorials), A368508 (perfect powers of composite superprimorials), and A368682.
These numbers are perfect powers of some smaller product of primorials.

Examples

			Let b(n) = A025487(n).
a(1) = b(1) = 1 = 1^k = b(1)^k, k >= 2,
a(2) = b(3) = 4 = 2^2 = b(2)^2,
a(3) = b(5) = 8 = 2^3 = b(2)^3,
a(6) = b(11) = 36 = 6^2 = b(4)^2,
a(9) = b(19) = 144 = 12^2 = b(6)^2, etc.
2 is not in the sequence since 2 is squarefree and not in A001597.
		

Crossrefs

Programs

  • Mathematica
    {1}~Join~Select[Range[4, 200000, 2], Or[PrimePowerQ[#], And[Union@ Differences@ PrimePi@ #1 == {1}, AllTrue[Union@ Differences@ #2, # <= 0 &], GCD @@ #2 > 1] & @@ Transpose@ FactorInteger[#]] &]

A365745 Powers s^m, m > 1, where s is a composite squarefree number that is not a primorial.

Original entry on oeis.org

100, 196, 225, 441, 484, 676, 1000, 1089, 1156, 1225, 1444, 1521, 1764, 2116, 2601, 2744, 3025, 3249, 3364, 3375, 3844, 4225, 4356, 4761, 4900, 5476, 5929, 6084, 6724, 7225, 7396, 7569, 8281, 8649, 8836, 9025, 9261, 10000, 10404, 10648, 11025, 11236, 12100, 12321
Offset: 1

Views

Author

Michael De Vlieger, Dec 10 2023

Keywords

Examples

			100 is the first term in this sequence since it is 10^2; 10 is a squarefree composite number. All powers 10^m, m > 1 are in the sequence.
36 is not in the sequence since it is 6^2, where 6 is a product of the smallest 2 primes; none of the powers 6^m, m > 1 are in the sequence.
		

Crossrefs

Programs

  • Mathematica
    nn = 2^30; ss = Floor@ Sqrt[nn]; p = 3;
    s = Complement[
      Select[Range[ss], And[SquareFreeQ[#], CompositeQ[#]] &],
      NestWhileList[(Set[p, NextPrime[p]]; # p) &, 6, # <= ss &] ];
    Union@ Reap[Do[k = 2; While[s[[i]]^k <= nn, Sow[s[[i]]^k]; k++],
      {i, Length[s]}] ][[-1, 1]]

Formula

This sequence is {A303606 \ A365308}.
This sequence contains powers s^m, m > 1, for s in {A120944 \ A002110}.

A367511 Highly composite numbers h(k) = A002182(k) such that h >= rad(h)^2, where rad() = A007947().

Original entry on oeis.org

1, 4, 36, 48, 45360, 50400
Offset: 1

Views

Author

Michael De Vlieger, Feb 08 2024

Keywords

Comments

Alternatively, this sequence lists h(k) such that A301413(k) >= A002110(A108602(k)), where A301413 is the "variable part" v described on page 5 of 12 of the Siano paper.
This sequence is likely finite and full. See Chapter III regarding the structure of "Highly Composite Numbers".
Terms larger than 36 are in A366250; A366250 is in A364702, which is in turn a proper subset of A332785, itself contained in A126706.
36 is in A365308, a proper subset of A303606, contained in A131605, in turn contained in A286708.

Examples

			Let P(n) = A002110(n).
a(1) = h(1) = 1 since 1 >= 1^2.
a(2) = h(3) = 4 since 4 >= P(1)^2, 4 >= 2^2.
a(3) = h(7) = 36 since 36 >= P(2)^2, 36 >= 6^2.
a(4) = h(8) = 48 since 48 >= P(2)^2, 48 >= 6^2.
a(5) = h(26) = 43560 since 43560 >= P(4)^2, where P(4) = 210, and 210^2 = 44100.
a(6) = h(27) = 50400 since 50400 >= P(4)^2.
Let V(i) = A301414(i) and let P(j) = A002110(j).
Plot of highly composite h = V(i)*P(j) at (x,y) = (j,i), i = 1..16, j = 1..7, showing h in this sequence in parentheses, and h in A168263 marked with an asterisk (*):
V(i)\P(j) 1   2    6   30   210    2310    30030 ...
        +---------------------------------------
      1 |(1*) 2*   6*
      2 |    (4*) 12*  60*
      4 |         24* 120*  840*
      6 |        (36) 180* 1260*
      8 |        (48) 240  1680*
     12 |             360  2520   27720*
     24 |             720  5040   55440   720720
     36 |                  7560   83160  1081080
     48 |                 10080  110880  1441440
     72 |                 15120  166320  2162160
     96 |                 20160  221760  2882880
    120 |                 25200  277200  3603600
    144 |                        332640  4324320
    216 |                (45360) 498960  6486480
    240 |                (50400) 554400  7207200
    ...
		

Crossrefs

Programs

  • Mathematica
    (* First load function f at A025487, then run the following: *)
    s = Union@ Flatten@ f[12];
    t = Map[DivisorSigma[0, #] &, s];
    h = Map[s[[FirstPosition[t, #][[1]]]] &, Union@ FoldList[Max, t]];
    Reap[Do[If[# >= Product[Prime[j], {j, PrimeNu[#]}]^2, Sow[#]] &[ h[[i]] ],
      {i, Length[h]}] ][[-1, 1]]
Showing 1-5 of 5 results.