A363098 -id:A363098 - cp's OEIS frontend

A363063 Positive integers k such that the largest power of p dividing k is larger than or equal to the largest power of q dividing k (i.e., A305720(k,p) >= A305720(k,q)) for all primes p and q with p < q.

1, 2, 4, 8, 12, 16, 24, 32, 48, 64, 96, 128, 144, 192, 256, 288, 384, 512, 576, 720, 768, 864, 1024, 1152, 1440, 1536, 1728, 2048, 2304, 2880, 3072, 3456, 4096, 4320, 4608, 5760, 6144, 6912, 8192, 8640, 9216, 10368, 11520, 12288, 13824, 16384, 17280, 18432

Offset: 1

Pontus von Brömssen, May 16 2023

nonn

Includes all products of terms in A347284, but there are also other terms such as 4320.

Closed under multiplication. - Peter Munn, May 21 2023

			151200 = 2^5 * 3^3 * 5^2 * 7 is a term, because 2^5 >= 3^3 >= 5^2 >= 7.
72 = 2^3 * 3^2 is not a term, because 2^3 < 3^2.
40 = 2^3 * 3^0 * 5 is not a term, because 3^0 < 5.
From _Michael De Vlieger_, May 19 2023: (Start)
Sequence read as an irregular triangle delimited by appearance of 2^m:
     1
     2
     4
     8     12
    16     24
    32     48
    64     96
   128    144    192
   256    288    384
   512    576    720    768    864
  1024   1152   1440   1536   1728
  2048   2304   2880   3072   3456
  4096   4320   4608   5760   6144   6912
  8192   8640   9216  10368  11520  12288  13824
  ... (End)

Pontus von Brömssen, Table of n, a(n) for n = 1..10000
Michael De Vlieger, Plot p^e | a(n) at (x,y) = (n, pi(p)), n = 1..1024, showing multiplicity e with a color function such that e = 1 is black, e = 2 is red, e = 3 is orange, etc., 12X vertical exaggeration. On the bottom, a color code represents a(n) is empty product (black), prime (red), composite prime power (gold), neither squarefree nor prime power (blue).
Michael De Vlieger, Plot multiplicities e in a(n) = Product p^e at (x,y) = (e, -n) for n = 1..1024, 8X horizontal exaggeration.

Cf. A181818, A305720, A347284, A363098 (primitive terms).

Subsequence of: A087980, {1} U A363122.

Select[Range[20000], # == 1 || PrimePi[(f = FactorInteger[#])[[-1, 1]]] == Length[f] && Greater @@ (Power @@@ f) &] (* Amiram Eldar, May 16 2023 *)

from sympy import nextprime
primes = [2] # global list of first primes
def f(kmax,pi,ppmax):
    # Generate numbers up to kmax with nonincreasing prime-powers <= ppmax, starting at the (pi+1)-st prime.
    if len(primes) <= pi: primes.append(nextprime(primes[-1]))
    p0 = primes[pi]
    ppmax = min(ppmax,kmax)
    if ppmax < p0:
        yield 1
        return
    pp = 1
    while pp <= ppmax:
        for x in f(kmax//pp,pi+1,pp):
            yield pp*x
        pp *= p0
def A363063_list(kmax):
    return sorted(f(kmax,0,kmax))

A363234 Least number divisible by the first n primes whose factorization into maximal prime powers, if ordered by increasing prime divisor, then has these prime power factors in decreasing order.

Original entry on oeis.org

1, 2, 12, 720, 151200, 4191264000, 251727315840000, 1542111744113740800000, 10769764221549079560253440000000, 12109394351419848024974600399142912000000000, 78344066654781231654807043124290195568885760000000000, 188552692884723759943358058475004257579791386442930585600000000000

Offset: 0

Michael De Vlieger and David A. Corneth, May 22 2023

a(n) is the least number in A347284 divisible by prime(n).
Also a(n) is the smallest positive integer divisible by prime(n) and prime(i)^e(i) > prime(i + 1)^e(i + 1) where e(k) is the valuation of prime(k) in a(n) and 1 <= i < n. - David A. Corneth, May 24 2023
Equivalently, we can say a(n) is the least number divisible by prime(n) in A363063. This is true also of A363098, the primitive terms of A363063. {a(n)} is the intersection of A347284 and A363098. - Peter Munn, May 29 2023
If we change the end of the sequence name from "decreasing order" to "increasing order", we get the primorial numbers (A002110). - Peter Munn, Jun 04 2023

			Table shows a(n) = A347284(j) = Product p(i)^m(i), m(i) is the i-th term read from left to right, delimited by ".", in row a(n) of A067255. Example: "4.2.1" signifies 2^4 * 3^2 * 5^1 = 720.
n    j   A067255(a(n))                                   a(n)
-------------------------------------------------------------
0    0                                                      1
1    1   1                                                  2
2    2   2.1                                               12
3    4   4.2.1                                            720
4    5   5.3.2.1                                       151200
5    8   8.5.3.2.1                                 4191264000
6   10   10.6.4.3.2.1                         251727315840000
7   13   13.8.5.4.3.2.1                1542111744113740800000
8   18   18.11.7.5.4.3.2.1   10769764221549079560253440000000
...

Michael De Vlieger, Plot prime(i)^k | a(n) at (x,y) = (k,-n) for n = 1..503.

Cf. A002110, A067255, A347355.

Subsequence of A347284, A363063, A363098.

nn = 120; a[0] = {0}; Do[b = {2^k}; Do[If[Last[b] == 1, Break[], i = 1; p = Prime[j]; While[p^i < b[[j - 1]], i++]; AppendTo[b, p^(i - 1)]], {j, 2, k}]; Set[a[k], b], {k, nn}]; s = DeleteCases[Array[a, nn], 1, {2}]; {1}~Join~Table[Times @@ s[[FirstPosition[s, _?(Length[#] == k &)][[1]]]], {k, Max[Length /@ s]}]
(* Generate terms from the linked image. Caution, terms become very large. *)
img = Import["https://oeis.org/A363234/a363234.png", "Image"]; Map[Times @@ MapIndexed[Prime[First[#2]]^#1 &, Reverse@ #] &, SplitBy[Position[ImageData[img][[1 ;; 12]], 0.], First][[All, All, -1]] ]

a(n) = {resf = matrix(n, 2); resf[,1] = primes(n)~; resf[n, 2] = 1; forstep(j = n-1, 1, -1, resf[j, 2] = logint(resf[j+1, 1]^resf[j+1, 2], resf[j, 1]) + 1); factorback(resf)} \\ David A. Corneth, May 24 2023

a(n) = A347284(A347355(n)).

cp's OEIS Frontend

A363063 Positive integers k such that the largest power of p dividing k is larger than or equal to the largest power of q dividing k (i.e., A305720(k,p) >= A305720(k,q)) for all primes p and q with p < q.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Python