A335740 -id:A335740 - cp's OEIS frontend

A070003 Numbers divisible by the square of their largest prime factor.

4, 8, 9, 16, 18, 25, 27, 32, 36, 49, 50, 54, 64, 72, 75, 81, 98, 100, 108, 121, 125, 128, 144, 147, 150, 162, 169, 196, 200, 216, 225, 242, 243, 245, 250, 256, 288, 289, 294, 300, 324, 338, 343, 361, 363, 375, 392, 400, 432, 441, 450, 484, 486, 490, 500, 507

Offset: 1

Labos Elemer, May 07 2002

nonn

Numbers n such that P(phi(n)) - phi(P(n)) = 1, where P(x) is the largest prime factor of x. P(phi(n)) - phi(P(n)) = A006530(A000010(n)) - A000010(A006530(n)).
Numbers n such that the value of the commutator of phi and P functions at n is -1.
Equivalently, n such that n and phi(n) have the same largest prime factor since Phi(p) = p-1 if p is prime. - Benoit Cloitre, Jun 08 2002
Since n is divisible by P(n)^2, n cannot divide P(n)! and so A057109 is a supersequence. Hence all A002034(a(n)) are composite. - Jonathan Sondow, Dec 28 2004
A225546 defines a self-inverse bijection between this sequence and A335740, considered as sets. - Peter Munn, Jul 19 2020

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
Paul Erdős and Ron L. Graham, On products of factorials, Bull. Inst. Math. Acad. Sinica 4:2 (1976), pp. 337-355. [alternate link]
Paul Erdős and Ilias Kastanas, Solution 6674:The smallest factorial that is a multiple of n, Amer. Math. Monthly 101 (1994) 179.
A. J. Kempner, Miscellanea, Amer. Math. Monthly, 25 (1918), 201-210. See Section II, "Concerning the smallest integer m! divisible by a given integer n."
Eric Weisstein's World of Mathematics, Greatest Prime Factor
Eric Weisstein's World of Mathematics, Totient Function

Subsequence of A057109, A122145.
Complement within A020725 of A102750.
Cf. A000010, A006530, A068211, A070777, A070812, A070002, A070004, A007283, A070813, A070814, A070815, A070816, A002034, A102067, A102068.
Related to A335740 via A225546.
A195212 is a subsequence.
Cf. A319988 (characteristic function). Positions of odd terms > 1 in A122111.

isA070003 := proc(n)
    if modp(n,A006530(n)^2) = 0 then # code re-use
        true;
    else
        false;
    end if;
end proc:
A070003 := proc(n)
    option remember ;
    if n =1 then
        4;
    else
        for a from procname(n-1)+1 do
            if isA070003(a) then
                return a
            end if;
        end do:
    end if;
end proc:
seq( A070003(n),n=1..80) ; # R. J. Mathar, Jun 27 2024

p[n_] := FactorInteger[n][[-1, 1]]; ep[n_] := EulerPhi[n]; fQ[n_] := p[ep[n]] == 1 + ep[p[n]]; Select[ Range[ 510], fQ] (* Robert G. Wilson v, Mar 26 2012 *)
Select[Range[500], FactorInteger[#][[-1,2]] > 1 &] (* T. D. Noe, Dec 06 2012 *)

for(n=3,1000,if(component(component(factor(n),1),omega(n))==component(component(factor(eulerphi(n)),1),omega(eulerphi(n))),print1(n,",")))

is(n)=my(f=factor(n)[,2]);f[#f]>1 \\ Charles R Greathouse IV, Mar 21 2012

sm(lim,mx)=if(mx==2,return(vector(log(lim+.5)\log(2)+1,i,1<<(i-1))));my(v=[1]);forprime(p=2,min(mx,lim),v=concat(v,p*sm(lim\p,p)));vecsort(v)
list(lim)=my(v=[]);forprime(p=2,sqrt(lim),v=concat(v,p^2*sm(lim\p^2,p)));vecsort(v) \\ Charles R Greathouse IV, Mar 27 2012

from sympy import factorint
def ok(n): f = factorint(n); return f[max(f)] >= 2
print(list(filter(ok, range(4, 508)))) # Michael S. Branicky, Apr 08 2021

Erdős proved that there are x * exp(-(1 + o(1))sqrt(log x log log x)) members of this sequence up to x. - Charles R Greathouse IV, Mar 26 2012

New name from Jonathan Sondow and Charles R Greathouse IV, Mar 27 2012

A335738 Factorize each integer m >= 2 as the product of powers of nonunit squarefree numbers with distinct exponents that are powers of 2. The sequence lists m such that the factor with the largest exponent is a power of 2.

Original entry on oeis.org

2, 4, 8, 12, 16, 20, 24, 28, 32, 40, 44, 48, 52, 56, 60, 64, 68, 76, 80, 84, 88, 92, 96, 104, 112, 116, 120, 124, 128, 132, 136, 140, 144, 148, 152, 156, 160, 164, 168, 172, 176, 184, 188, 192, 204, 208, 212, 220, 224, 228, 232, 236, 240, 244, 248, 256, 260, 264, 268, 272

Offset: 1

Peter Munn, Jun 20 2020

nonn

2 is the only term not divisible by 4. All powers of 2 are present. Every term divisible by an odd square is divisible by 16, the first such being 144.
The defined factorization is unique. Every positive number is a product of at most one squarefree number (A005117), at most one square of a squarefree number (A062503), at most one 4th power of a squarefree number (A113849), at most one 8th power of a squarefree number, and so on.
Iteratively map m using A000188, until 1 is reached, as A000188^k(m), for some k >= 1. m is in the sequence if and only if the preceding number, A000188^(k-1)(m), is 2. k can be shown to be A299090(m).
Closed under squaring, but not closed under multiplication: 12 = 3^1 * 2^2 and 432 = 3^1 * 3^2 * 2^4 are in the sequence, but 12 * 432 = 5184 = 3^4 * 2^6 = 2^2 * 6^4 is not.
The asymptotic density of this sequence is Sum_{k>=0} (d(2^(k+1)) - d(2^k))/2^(2^(k+1)-1) = 0.21363357193921052068..., where d(k) = 2^(k-1)/((2^k-1)*zeta(k)) is the asymptotic density of odd k-free numbers for k >= 2, and d(1) = 0. - Amiram Eldar, Feb 10 2024

			6 is a squarefree number, so its factorization for the definition (into powers of nonunit squarefree numbers with distinct exponents that are powers of 2) is the trivial "6^1". 6^1 is therefore the factor with the largest exponent, and is not a power of 2, so 6 is not in the sequence.
48 factorizes for the definition as 3^1 * 2^4. The factor with the largest exponent is 2^4, which is a power of 2, so 48 is in the sequence.
10^100 (a googol) factorizes in this way as 10^4 * 10^32 * 10^64. The factor with the largest exponent, 10^64, is a power of 10, not a power of 2, so 10^100 is not in the sequence.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Complement within A020725 of A335740.
A000188, A299090 are used in a formula defining this sequence.
Powers of squarefree numbers: A005117(1), A144338(1), A062503(2), A113849(4).
Subsequences: A000079\{1}, A001749, A181818\{1}, A273798.
Numbers in the even bisection of A336322.
Row m of A352780 essentially gives the defined factorization of m.

f[p_, e_] := p^Floor[e/2]; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; Select[Range[2, 300], FixedPointList[s, #] [[-3]] == 2 &] (* Amiram Eldar, Nov 27 2020 *)

is(n) = {my(e = valuation(n, 2), o = n >> e); if(e == 0, 0, if(o == 1, n > 1, floor(logint(e, 2)) > floor(logint(vecmax(factor(o)[,2]), 2))));} \\ Amiram Eldar, Feb 10 2024

{a(n)} = {m : m >= 2 and A000188^(k-1)(m) = 2, where k = A299090(m)}.

{a(n)} = {m : m >= 2 and A352780(m,e) = 2^(2^e), where e = A299090(m)-1}. - Peter Munn, Jun 24 2022

A336321 a(n) = A122111(A225546(n)).

Original entry on oeis.org

1, 2, 3, 4, 7, 5, 19, 6, 9, 11, 53, 10, 131, 23, 13, 8, 311, 15, 719, 22, 29, 59, 1619, 14, 49, 137, 21, 46, 3671, 17, 8161, 12, 61, 313, 37, 25, 17863, 727, 139, 26, 38873, 31, 84017, 118, 39, 1621, 180503, 20, 361, 77, 317, 274, 386093, 33, 71, 58, 733, 3673, 821641, 34, 1742537, 8167, 87, 18, 151, 67, 3681131, 626, 1627, 41, 7754077, 35, 16290047

Offset: 1

Antti Karttunen and Peter Munn, Jul 17 2020

nonn

A122111 and A225546 are both self-inverse permutations of the positive integers based on prime factorizations, and they share further common properties. For instance, they map the prime numbers to powers of 2: A122111 maps the k-th prime to 2^k, whereas A225546 maps it to 2^2^(k-1).

In composing these permutations, this sequence maps the squarefree numbers, as listed in A019565, to the prime numbers in increasing order; and the list of powers of 2 to the "normal" numbers (A055932), as listed in A057335.

			From _Peter Munn_, Jan 04 2021: (Start)
In this set of examples we consider [a(n)] as a function a(.) with an inverse, a^-1(.).
First, a table showing mapping of the powers of 2:
  n     a^-1(2^n) =    2^n =        a(2^n) =
        A001146(n-1)   A000079(n)   A057335(n)
  0             (1)         1            1
  1               2         2            2
  2               4         4            4
  3              16         8            6
  4             256        16            8
  5           65536        32           12
  6      4294967296        64           18
  ...
Next, a table showing mapping of the squarefree numbers, as listed in A019565 (a lexicographic ordering by prime factors):
  n   a^-1(A019565(n))   A019565(n)      a(A019565(n))   a^2(A019565(n))
      Cf. {A337533}      Cf. {A005117}   = prime(n)      = A033844(n-1)
  0              1               1             (1)               (1)
  1              2               2               2                 2
  2              3               3               3                 3
  3              8               6               5                 7
  4              6               5               7                19
  5             12              10              11                53
  6             18              15              13               131
  7            128              30              17               311
  8              5               7              19               719
  9             24              14              23              1619
  ...
As sets, the above columns are A337533, A005117, A008578, {1} U A033844.
Similarly, we get bijections between sets A000290\{0} -> {1} U A070003; and {1} U A335740 -> A005408 -> A066207.
(End)

A122111 composed with A225546.
Cf. A336322 (inverse permutation).
Other sequences used in a definition of this sequence: A000040, A000188, A019565, A248663, A253550, A253560.
Sequences used to express relationship between terms of this sequence: A003159, A003961, A297002, A334747.
Cf. A057335.
A mapping between the binary tree sequences A334866 and A253563.
Lists of sets (S_1, S_2, ... S_j) related by the bijection defined by the sequence: (A000290\{0}, {1} U A070003), ({1} U A001146, A000079, A055932), ({1} U A335740, A005408, A066207), (A337533, A005117, A008578, {1} U A033844).

a(n) = A122111(A225546(n)).
Alternative definition: (Start)
Write n = m^2 * A019565(j), where m = A000188(n), j = A248663(n).
a(1) = 1; otherwise for m = 1, a(n) = A000040(j), for m > 1, a(n) = A253550^j(A253560(a(m))).
(End)
a(A000040(m)) = A033844(m-1).
a(A001146(m)) = 2^(m+1).
a(2^n) = A057335(n).
a(n^2) = A253560(a(n)).
For n in A003159, a(2n) = b(a(n)), where b(1) = 2, b(n) = A253550(n), n >= 2.
More generally, a(A334747(n)) = b(a(n)).
a(A003961(n)) = A297002(a(n)).
a(A334866(m)) = A253563(m).

A336322 a(n) = A225546(A122111(n)).

Original entry on oeis.org

1, 2, 3, 4, 6, 8, 5, 16, 9, 12, 10, 32, 15, 24, 18, 256, 30, 64, 7, 48, 27, 20, 14, 512, 36, 40, 81, 96, 21, 128, 42, 65536, 54, 60, 72, 1024, 35, 120, 45, 768, 70, 192, 105, 80, 162, 28, 210, 131072, 25, 144, 90, 160, 11, 4096, 108, 1536, 135, 56, 22, 2048, 33, 84, 243, 4294967296, 216, 384, 66, 240, 270, 288, 55, 262144, 110, 168, 324, 480, 50

Offset: 1

Antti Karttunen and Peter Munn, Jul 17 2020

nonn

A225546 and A122111 are both self-inverse permutations of the positive integers based on prime factorizations, and they share further common properties. For instance, they map the prime numbers to powers of 2: A225546 maps the k-th prime to 2^2^(k-1), whereas A122111 maps it to 2^k.

In composing these permutations, this sequence maps the list of prime numbers to the squarefree numbers, as listed in A019565; and the "normal" numbers (A055932), as listed in A057335, to ascending powers of 2.

Index entries for sequences that are permutations of the natural numbers

A225546 composed with A122111.
Sorted even bisection: A335738.
Sorted odd bisection (excluding 1): A335740.
Sequences used to express relationship between terms of this sequence: A001222, A003961, A253560, A331590, A350066.
Sequences of sequences (S_1, S_2, ... S_j) with the property a(S_i) = S_{i+1}, or essentially so: (A033844, A000040, A019565), (A057335, A000079, A001146), (A000244, A011764), (A001248, A334110), (A253563, A334866).
The inverse permutation, A336321, lists sequences where the property is weaker (between the sets of terms).

a(A033844(m)) = A000040(m+1). [Offset corrected Peter Munn, Feb 14 2022]
a(A000040(m)) = A019565(m).
a(A057335(m)) = 2^m.
For m >= 1, a(2^m) = A001146(m-1).
a(A253563(m)) = A334866(m).
From Peter Munn, Feb 14 2022: (Start)
a(A253560(n)) = a(n)^2.
For n >= 2, a(A003961(n)) = A331590(a(n), 2^2^(A001222(n)-1)).
a(A350066(n, k)) = A331590(a(n), a(k)).
(End)

A381439 Numbers whose exponent of 2 in their canonical prime factorization is not larger than all the other exponents.

Original entry on oeis.org

3, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 18, 19, 21, 22, 23, 25, 26, 27, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 41, 42, 43, 45, 46, 47, 49, 50, 51, 53, 54, 55, 57, 58, 59, 61, 62, 63, 65, 66, 67, 69, 70, 71, 73, 74, 75, 77, 78, 79, 81, 82, 83, 85, 86, 87, 89

Offset: 1

Gus Wiseman, Mar 02 2025

nonn

First differs from A335740 in lacking 72, which has prime indices {1,1,1,2,2} and section-sum partition (3,3,1).
Also numbers whose section-sum partition of prime indices (A381436) ends with a number > 1.
Includes all squarefree numbers (A005117) except 2.

			The terms together with their prime indices begin:
     3: {2}        25: {3,3}        45: {2,2,3}
     5: {3}        26: {1,6}        46: {1,9}
     6: {1,2}      27: {2,2,2}      47: {15}
     7: {4}        29: {10}         49: {4,4}
     9: {2,2}      30: {1,2,3}      50: {1,3,3}
    10: {1,3}      31: {11}         51: {2,7}
    11: {5}        33: {2,5}        53: {16}
    13: {6}        34: {1,7}        54: {1,2,2,2}
    14: {1,4}      35: {3,4}        55: {3,5}
    15: {2,3}      36: {1,1,2,2}    57: {2,8}
    17: {7}        37: {12}         58: {1,10}
    18: {1,2,2}    38: {1,8}        59: {17}
    19: {8}        39: {2,6}        61: {18}
    21: {2,4}      41: {13}         62: {1,11}
    22: {1,5}      42: {1,2,4}      63: {2,2,4}
    23: {9}        43: {14}         65: {3,6}

The LHS (exponent of 2) is A007814.
The complement is A360013 = 2*A360015 (if we prepend 1), counted by A241131 (shifted right and starting with 1 instead of 0).
The case of equality is A360014, inclusive A360015.
The RHS (greatest exponent of an odd prime factor) is A375669.
These are positions of terms > 1 in A381437.
Partitions of this type are counted by A381544.
A000040 lists the primes, differences A001223.
A051903 gives greatest prime exponent, least A051904.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A122111 represents conjugation in terms of Heinz numbers.
A239455 counts Look-and-Say partitions, complement A351293.
A381436 gives section-sum partition of prime indices, Heinz number A381431.
A381438 counts partitions by last part part of section-sum partition.
Cf. A000720, A001221, A001222, A001694, A003557, A005117, A066328, A130091, A181819, A212166, A380955.

Select[Range[100],FactorInteger[2*#][[1,2]]-1<=Max@@Last/@Rest[FactorInteger[2*#]]&]

Positive integers n such that A007814(n) <= A375669(n).

cp's OEIS Frontend

A070003 Numbers divisible by the square of their largest prime factor.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

PARI

Python

Formula

Extensions

A335738 Factorize each integer m >= 2 as the product of powers of nonunit squarefree numbers with distinct exponents that are powers of 2. The sequence lists m such that the factor with the largest exponent is a power of 2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A336321 a(n) = A122111(A225546(n)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

A336322 a(n) = A225546(A122111(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Formula

A381439 Numbers whose exponent of 2 in their canonical prime factorization is not larger than all the other exponents.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula