A070003 -id:A070003 - cp's OEIS frontend

A122111 Self-inverse permutation of the positive integers induced by partition enumeration in A112798 and partition conjugation.

1, 2, 4, 3, 8, 6, 16, 5, 9, 12, 32, 10, 64, 24, 18, 7, 128, 15, 256, 20, 36, 48, 512, 14, 27, 96, 25, 40, 1024, 30, 2048, 11, 72, 192, 54, 21, 4096, 384, 144, 28, 8192, 60, 16384, 80, 50, 768, 32768, 22, 81, 45, 288, 160, 65536, 35, 108, 56, 576, 1536, 131072, 42

Offset: 1

Franklin T. Adams-Watters, Oct 18 2006

Factor n; replace each prime(i) with i, take the conjugate partition, replace parts i with prime(i) and multiply out.
From Antti Karttunen, May 12-19 2014: (Start)
For all n >= 1, A001222(a(n)) = A061395(n), and vice versa, A061395(a(n)) = A001222(n).
Because the partition conjugation doesn't change the partition's total sum, this permutation preserves A056239, i.e., A056239(a(n)) = A056239(n) for all n.
(Similarly, for all n, A001221(a(n)) = A001221(n), because the number of steps in the Ferrers/Young-diagram stays invariant under the conjugation. - Note added Apr 29 2022).
Because this permutation commutes with A241909, in other words, as a(A241909(n)) = A241909(a(n)) for all n, from which follows, because both permutations are self-inverse, that a(n) = A241909(a(A241909(n))), it means that this is also induced when partitions are conjugated in the partition enumeration system A241918. (Not only in A112798.)
(End)
From Antti Karttunen, Jul 31 2014: (Start)
Rows in arrays A243060 and A243070 converge towards this sequence, and also, assuming no surprises at the rate of that convergence, this sequence occurs also as the central diagonal of both.
Each even number is mapped to a unique term of A102750 and vice versa.
Conversely, each odd number (larger than 1) is mapped to a unique term of A070003, and vice versa. The permutation pair A243287-A243288 has the same property. This is also used to induce the permutations A244981-A244984.
Taking the odd bisection and dividing out the largest prime factor results in the permutation A243505.
Shares with A245613 the property that each term of A028260 is mapped to a unique term of A244990 and each term of A026424 is mapped to a unique term of A244991.
Conversely, with A245614 (the inverse of above), shares the property that each term of A244990 is mapped to a unique term of A028260 and each term of A244991 is mapped to a unique term of A026424.
(End)
The Maple program follows the steps described in the first comment. The subprogram C yields the conjugate partition of a given partition. - Emeric Deutsch, May 09 2015
The Heinz number of the partition that is conjugate to the partition with Heinz number n. The Heinz number of a partition p = [p_1, p_2, ..., p_r] is defined as Product(p_j-th prime, j=1...r). Example: a(3) = 4. Indeed, the partition with Heinz number 3 is [2]; its conjugate is [1,1] having Heinz number 4. - Emeric Deutsch, May 19 2015

Alois P. Heinz, Table of n, a(n) for n = 1..10000 (first 1024 terms from Antti Karttunen)
Antti Karttunen, A few notes on A122111, A241909 & A241916.
Index entries for sequences that are permutations of the natural numbers

Cf. A088902 (fixed points).
Cf. A112798, A241918 (conjugates the partitions listed in these two tables).
Cf. A243060 and A243070. (Limit of rows in these arrays, and also their central diagonal).
Cf. also A080576, A036036, A001221, A001222, A061395, A243503, A056239, A052126, A003961, A064989, A071178, A241917, A241919, A242378, A242424, A006530, A105560, A070003, A102750, A066829, A028260, A026424, A244990, A244991, A244992, A108951, A181815, A181819, A181820, A238745, A238690, A242421.
Cf. A319988 (parity of this sequence for n > 1), A336124 (a(n) mod 4).
{A000027, A122111, A241909, A241916} form a 4-group.
{A000027, A122111, A153212, A242419} form also a 4-group.
Other related permutations: A048673-A064216, A244981-A244982, A244983-A244984, A243287-A243288, A243505-A243506, A245613-A245614, A075157, A075158, A129594, A069799, A242415, A245451, A245452, A245454, also A336321 & A336322 (composed with A225546).
Cf. also array A350066 [A(i, j) = a(a(i)*a(j))].

with(numtheory): c := proc (n) local B, C: B := proc (n) local pf: pf := op(2, ifactors(n)): [seq(seq(pi(op(1, op(i, pf))), j = 1 .. op(2, op(i, pf))), i = 1 .. nops(pf))] end proc: C := proc (P) local a: a := proc (j) local c, i: c := 0; for i to nops(P) do if j <= P[i] then c := c+1 else  end if end do: c end proc: [seq(a(k), k = 1 .. max(P))] end proc: mul(ithprime(C(B(n))[q]), q = 1 .. nops(C(B(n)))) end proc: seq(c(n), n = 1 .. 59); # Emeric Deutsch, May 09 2015
# second Maple program:
a:= n-> (l-> mul(ithprime(add(`if`(jAlois P. Heinz, Sep 30 2017

A122111[1] = 1; A122111[n_] := Module[{l = #, m = 0}, Times @@ Power @@@ Table[l -= m; l = DeleteCases[l, 0]; {Prime@Length@l, m = Min@l}, Length@Union@l]] &@Catenate[ConstantArray[PrimePi[#1], #2] & @@@ FactorInteger@n]; Array[A122111, 60] (* JungHwan Min, Aug 22 2016 *)
a[n_] := Function[l, Product[Prime[Sum[If[jJean-François Alcover, Sep 23 2020, after Alois P. Heinz *)

A122111(n) = if(1==n,n,my(f=factor(n), es=Vecrev(f[,2]),is=concat(apply(primepi,Vecrev(f[,1])),[0]),pri=0,m=1); for(i=1, #es, pri += es[i]; m *= prime(pri)^(is[i]-is[1+i])); (m)); \\ Antti Karttunen, Jul 20 2020

from sympy import factorint, prevprime, prime, primefactors
from operator import mul
def a001222(n): return 0 if n==1 else a001222(n/primefactors(n)[0]) + 1
def a064989(n):
    f=factorint(n)
    return 1 if n==1 else reduce(mul, [1 if i==2 else prevprime(i)**f[i] for i in f])
def a105560(n): return 1 if n==1 else prime(a001222(n))
def a(n): return 1 if n==1 else a105560(n)*a(a064989(n))
[a(n) for n in range(1, 101)] # Indranil Ghosh, Jun 15 2017

;; Uses Antti Karttunen's IntSeq-library.
(definec (A122111 n) (if (<= n 1) n (* (A000040 (A001222 n)) (A122111 (A064989 n)))))
;; Antti Karttunen, May 12 2014

;; Uses Antti Karttunen's IntSeq-library.
(definec (A122111 n) (if (<= n 1) n (* (A000079 (A241917 n)) (A003961 (A122111 (A052126 n))))))
;; Antti Karttunen, May 12 2014

;; Uses Antti Karttunen's IntSeq-library.
(definec (A122111 n) (if (<= n 1) n (* (expt (A000040 (A071178 n)) (A241919 n)) (A242378bi (A071178 n) (A122111 (A051119 n))))))
;; Antti Karttunen, May 12 2014

From Antti Karttunen, May 12-19 2014: (Start)
a(1) = 1, a(p_i) = 2^i, and for other cases, if n = p_i1 * p_i2 * p_i3 * ... * p_{k-1} * p_k, where p's are primes, not necessarily distinct, sorted into nondescending order so that i1 <= i2 <= i3 <= ... <= i_{k-1} <= ik, then a(n) = 2^(ik-i_{k-1}) * 3^(i_{k-1}-i_{k-2}) * ... * p_{i_{k-1}}^(i2-i1) * p_ik^(i1).
This can be implemented as a recurrence, with base case a(1) = 1,
and then using any of the following three alternative formulas:
  a(n) = A105560(n) * a(A064989(n)) = A000040(A001222(n)) * a(A064989(n)). [Cf. the formula for A242424.]
  a(n) = A000079(A241917(n)) * A003961(a(A052126(n))).
  a(n) = (A000040(A071178(n))^A241919(n)) * A242378(A071178(n), a(A051119(n))). [Here ^ stands for the ordinary exponentiation, and the bivariate function A242378(k,n) changes each prime p(i) in the prime factorization of n to p(i+k), i.e., it's the result of A003961 iterated k times starting from n.]
a(n) = 1 + A075157(A129594(A075158(n-1))). [Follows from the commutativity with A241909, please see the comments section.]
(End)
From Antti Karttunen, Jul 31 2014: (Start)
As a composition of related permutations:
a(n) = A153212(A242419(n)) = A242419(A153212(n)).
a(n) = A241909(A241916(n)) = A241916(A241909(n)).
a(n) = A243505(A048673(n)).
a(n) = A064216(A243506(n)).
Other identities. For all n >= 1, the following holds:
A006530(a(n)) = A105560(n). [The latter sequence gives greatest prime factor of the n-th term].
a(2n)/a(n) = A105560(2n)/A105560(n), which is equal to A003961(A105560(n))/A105560(n) when n > 1.
A243505(n) = A052126(a(2n-1)) = A052126(a(4n-2)).
A066829(n) = A244992(a(n)) and vice versa, A244992(n) = A066829(a(n)).
A243503(a(n)) = A243503(n). [Because partition conjugation does not change the partition size.]
A238690(a(n)) = A238690(n). - per Matthew Vandermast's note in that sequence.
A238745(n) = a(A181819(n)) and a(A238745(n)) = A181819(n). - per Matthew Vandermast's note in A238745.
A181815(n) = a(A181820(n)) and a(A181815(n)) = A181820(n). - per Matthew Vandermast's note in A181815.
(End)
a(n) = A181819(A108951(n)). [Prime shadow of the primorial inflation of n] - Antti Karttunen, Apr 29 2022

A071178 Exponent of the largest prime factor of n.

Original entry on oeis.org

0, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 1, 1, 1, 1, 4, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 3, 1, 1, 1, 1, 5, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 6, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1

Offset: 1

Benoit Cloitre, Jun 10 2002

a(n) = the multiplicity of the largest part in the partition having Heinz number n. We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product(p_j-th prime, j=1...r) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] we get 2*2*3*7*29 = 2436. Example: a(18) = 2; indeed, the partition having Heinz number 18 = 2*3*3 is [1,2,2]. - Emeric Deutsch, Jun 04 2015

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A001221, A001222, A006530, A053585, A067029, A067255, A070003, A124010, A215366.

a071178 = last . a124010_row -- Reinhard Zumkeller, Aug 27 2011

with(numtheory): with(padic):
a:= n-> `if`(n=1, 0, ordp(n, max(factorset(n)[]))):
seq(a(n), n=1..120);  # Alois P. Heinz, Jun 04 2015

a[n_] := FactorInteger[n] // Last // Last; Table[a[n], {n, 1, 120}] (* Jean-François Alcover, Jun 12 2015 *)
Join[{0},Table[FactorInteger[n][[-1,2]],{n,2,120}]] (* Harvey P. Dale, Aug 02 2025 *)

a(n) = if(n == 1, 0, my(e = factor(n)[, 2]); e[#e]); \\ Amiram Eldar, Oct 02 2024

from sympy import factorint
def A071178(n): return max(factorint(n).items())[1] if n>1 else 0 # Chai Wah Wu, Oct 10 2023

a(n) = A124010(n, A001221(n)); A053585(n) = A006530(n)^a(n). - Reinhard Zumkeller, Aug 27 2011
a(n) = A067255(n, A001222(n)). - Reinhard Zumkeller, Jun 11 2013
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1 (since the asymptotic density of A070003 is 0). - Amiram Eldar, Oct 02 2024

A356862 Numbers with a unique largest prime exponent.

Original entry on oeis.org

2, 3, 4, 5, 7, 8, 9, 11, 12, 13, 16, 17, 18, 19, 20, 23, 24, 25, 27, 28, 29, 31, 32, 37, 40, 41, 43, 44, 45, 47, 48, 49, 50, 52, 53, 54, 56, 59, 60, 61, 63, 64, 67, 68, 71, 72, 73, 75, 76, 79, 80, 81, 83, 84, 88, 89, 90, 92, 96, 97, 98, 99, 101, 103, 104

Offset: 1

Jens Ahlström, Sep 01 2022

If the prime factorization of k has a unique largest exponent, then k is a term.
Numbers whose multiset of prime factors (with multiplicity) has a unique mode. - Gus Wiseman, May 12 2023
Disjoint union of A246655 and A376250. The asymptotic density of this sequence, 0.3660366524547281232052..., is equal to the density of A376250 since the prime powers have a zero density. - Amiram Eldar, Sep 17 2024

			Prime powers (A246655) are in the sequence, since they have only one prime exponent in their prime factorization, hence a unique largest exponent.
144 is in the sequence, since 144 = 2^4 * 3^2 and there is the unique largest exponent 4.
225 is not in the sequence, since 225 = 3^2 * 5^2 and the largest exponent 2 is not unique, but rather it is the exponent of both the prime factor 3 and of the prime factor 5.

Jens Ahlström, Table of n, a(n) for n = 1..10000
Wikipedia, Mode (statistics).

Subsequence of A319161 (which has additional terms 1, 180, 252, 300, 396, 450, 468, ...).
Cf. A246655, A356838, A356840, A376250.
For factors instead of exponents we have A102750.
For smallest instead of largest we have A359178, counted by A362610.
The complement is A362605, counted by A362607.
The complement for co-mode is A362606, counted by A362609.
Partitions of this type are counted by A362608.
These are the positions of 1's in A362611, for co-modes A362613.
A001221 is the number of prime exponents, sum A001222.
A027746 lists prime factors, A112798 indices, A124010 exponents.
A362614 counts partitions by number of modes, A362615 co-modes.
Cf. A000041, A002865, A051903, A056239, A070003, A247180, A283050, A327473.

Select[Range[2, 100], Count[(e = FactorInteger[#][[;; , 2]]), Max[e]] == 1 &] (* Amiram Eldar, Sep 01 2022 *)

isok(k) = if (k>1, my(f=factor(k), m=vecmax(f[,2]), w=select(x->(f[x,2] == m), [1..#f~])); #w == 1); \\ Michel Marcus, Sep 01 2022

from sympy import factorint
from collections import Counter
def ok(k):
    c = Counter(factorint(k)).most_common(2)
    return not (len(c) > 1 and c[0][1] == c[1][1])
print([k for k in range(2, 105) if ok(k)])

from sympy import factorint
from itertools import count, islice
def A356862_gen(startvalue=2): # generator of terms >= startvalue
    return filter(lambda n:len(f:=sorted(factorint(n).values(),reverse=True))==1 or f[0]!=f[1],count(max(startvalue,2)))
A356862_list = list(islice(A356862_gen(),30)) # Chai Wah Wu, Sep 10 2022

A362612 Number of integer partitions of n such that the greatest part is the unique mode.

Original entry on oeis.org

0, 1, 2, 2, 3, 3, 4, 4, 6, 6, 7, 9, 10, 12, 15, 16, 19, 23, 26, 32, 37, 41, 48, 58, 65, 75, 88, 101, 115, 135, 151, 176, 200, 228, 261, 300, 336, 385, 439, 498, 561, 641, 717, 818, 921, 1036, 1166, 1321, 1477, 1667, 1867, 2099, 2346, 2640, 2944, 3303, 3684

Offset: 0

Gus Wiseman, May 03 2023

nonn

A mode in a multiset is an element that appears at least as many times as each of the others. For example, the modes of {a,a,b,b,b,c,d,d,d} are {b,d}.

			The a(1) = 1 through a(10) = 7 partitions (A = 10):
  1  2   3    4     5      6       7        8         9          A
     11  111  22    221    33      331      44        333        55
              1111  11111  222     2221     332       441        442
                           111111  1111111  2222      3321       3331
                                            22211     22221      22222
                                            11111111  111111111  222211
                                                                 1111111111

John Tyler Rascoe, Table of n, a(n) for n = 0..500

For median instead of mode we have A053263, complement A237821.
These partitions have ranks A362616.
A000041 counts integer partitions.
A275870 counts collapsible partitions.
A359893 counts partitions by median.
A362607 counts partitions with more than one mode, ranks A362605.
A362608 counts partitions with a unique mode, ranks A356862.
A362611 counts modes in prime factorization.
A362614 counts partitions by number of modes, co-modes A362615.
Cf. A002865, A008284, A070003, A098859, A102750, A237984, A238478, A238479, A327472, A362609, A362610.

Table[Length[Select[IntegerPartitions[n],Commonest[#]=={Max[#]}&]],{n,0,30}]

A_x(N)={my(x='x+O('x^N), g=sum(i=1, N, sum(j=1, N/i, x^(i*j)*prod(k=1,i-1,(1-x^(j*k))/(1-x^k))))); concat([0],Vec(g))}
A_x(60) \\ John Tyler Rascoe, Apr 03 2024

G.f.: Sum_{i, j>0} x^(i*j) * Product_{k=1,i-1} ((1-x^(j*k))/(1-x^k)). - John Tyler Rascoe, Apr 03 2024

A102750 Numbers n such that square of largest prime dividing n does not divide n.

Original entry on oeis.org

2, 3, 5, 6, 7, 10, 11, 12, 13, 14, 15, 17, 19, 20, 21, 22, 23, 24, 26, 28, 29, 30, 31, 33, 34, 35, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 51, 52, 53, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 76, 77, 78, 79, 80, 82, 83, 84, 85

Offset: 1

Leroy Quet, Feb 09 2005

nonn

Numbers n such that the exponent of the largest prime dividing n is one. - Harvey P. Dale, May 02 2019
From Peter Munn, Sep 30 2020: (Start)
2 together with numbers on the left half of the Doudna sequence tree depicted in Antti Karttunen's 2014 comment in A005940.
This sequence and A335738, considered as sets, are related by the self-inverse function A225546(.), which maps the members of either set 1:1 onto the other set.
(End)

			63 is included because 63 = 3^2 *7 and 7 (the largest prime dividing 63) only divides 63 once.

Rémy Sigrist, Table of n, a(n) for n = 1..25000

Cf. A070003 (complement, apart from the term 1 that is in neither sequence).
Related to A335738 via A225546.
Cf. A005940.

Select[Range[2,100],FactorInteger[#][[-1,2]]==1&] (* Harvey P. Dale, May 02 2019 *)

isok(n) = my(f = factor(n)); n % f[#f~, 1]^2; \\ Michel Marcus, May 20 2014

More terms from Erich Friedman, Aug 08 2005

A241916 a(2^k) = 2^k, and for other numbers, if n = 2^e1 * 3^e2 * 5^e3 * ... p_k^e_k, then a(n) = 2^(e_k - 1) * 3^(e_{k-1}) * ... * p_{k-1}^e2 * p_k^(e1+1). Here p_k is the greatest prime factor of n (A006530), and e_k is its exponent (A071178), and the exponents e1, ..., e_{k-1} >= 0.

Original entry on oeis.org

1, 2, 3, 4, 5, 9, 7, 8, 6, 25, 11, 27, 13, 49, 15, 16, 17, 18, 19, 125, 35, 121, 23, 81, 10, 169, 12, 343, 29, 75, 31, 32, 77, 289, 21, 54, 37, 361, 143, 625, 41, 245, 43, 1331, 45, 529, 47, 243, 14, 50, 221, 2197, 53, 36, 55, 2401, 323, 841, 59, 375, 61, 961, 175, 64

Offset: 1

Antti Karttunen, May 03 2014

nonn

For other numbers than the powers of 2 (that are fixed), this permutation reverses the sequence of exponents in the prime factorization of n from the exponent of 2 to that of the largest prime factor, except that the exponents of 2 and the greatest prime factor present are adjusted by one. Note that some of the exponents might be zeros.
Self-inverse permutation of natural numbers, composition of A122111 & A241909 in either order: a(n) = A122111(A241909(n)) = A241909(A122111(n)).
This permutation preserves both bigomega and the (index of) largest prime factor: for all n it holds that A001222(a(n)) = A001222(n) and A006530(a(n)) = A006530(n) [equally: A061395(a(n)) = A061395(n)].
From the above it follows, that this fixes both primes (A000040) and powers of two (A000079), among other numbers.
Even positions from n=4 onward contain only terms of A070003, and the odd positions only the terms of A102750, apart from 1 which is at a(1), and 2 which is at a(2).

Antti Karttunen, Table of n, a(n) for n = 1..8192
A. Karttunen, A few notes on A122111, A241909 & A241916.
Index entries for sequences that are permutations of the natural numbers

A241912 gives the fixed points; A241913 their complement.
Cf. A006530, A137502, A070003, A102750, A278525.
{A000027, A122111, A241909, A241916} form a 4-group.
The sum of prime indices of a(n) is A243503(n).
Even bisection of A358195 = Heinz numbers of rows of A358172.
A112798 lists prime indices, length A001222, sum A056239.
Cf. A005940, A019565, A161511, A246277, A253565, A326844, A356958, A358137.

nn = 65; f[n_] := If[n == 1, {0}, Function[f, ReplacePart[Table[0, {PrimePi[f[[-1, 1]]]}], #] &@ Map[PrimePi@ First@ # -> Last@ # &, f]]@ FactorInteger@ n]; g[w_List] := Times @@ Flatten@ MapIndexed[Prime[#2]^#1 &, w]; Table[If[IntegerQ@ #, n/4, g@ Reverse@(# - Join[{1}, ConstantArray[0, Length@ # - 2], {1}] &@ f@ n)] &@ Log2@ n, {n, 4, 4 nn, 4}] (* Michael De Vlieger, Aug 27 2016 *)

A209229(n) = (n && !bitand(n,n-1));
A241916(n) = if(1==A209229(n), n, my(f = factor(2*n), nbf = #f~, igp = primepi(f[nbf,1]), g = f); for(i=1,nbf,g[i,1] = prime(1+igp-primepi(f[i,1]))); factorback(g)/2); \\ Antti Karttunen, Jul 02 2018

(define (A241916 n) (A122111 (A241909 n)))

a(1)=1, and for n>1, a(n) = A006530(n) * A137502(n)/2.
a(n) = A122111(A241909(n)) = A241909(A122111(n)).
If 2n has prime factorization Product_{i=1..k} prime(x_i), then a(n) = Product_{i=1..k-1} prime(x_k-x_i+1). The opposite version is A000027, even bisection of A246277. - Gus Wiseman, Dec 28 2022

Description clarified by Antti Karttunen, Jul 02 2018

A253560 Multiply n by its largest prime factor: a(n) = A006530(n) * n.

Original entry on oeis.org

1, 4, 9, 8, 25, 18, 49, 16, 27, 50, 121, 36, 169, 98, 75, 32, 289, 54, 361, 100, 147, 242, 529, 72, 125, 338, 81, 196, 841, 150, 961, 64, 363, 578, 245, 108, 1369, 722, 507, 200, 1681, 294, 1849, 484, 225, 1058, 2209, 144, 343, 250, 867, 676, 2809, 162, 605, 392, 1083, 1682, 3481, 300, 3721, 1922, 441

Offset: 1

Antti Karttunen, Jan 03 2015

nonn

Antti Karttunen, Table of n, a(n) for n = 1..10000

Essentially the same as A129598, except that here we have a(1) = 1.
Cf. A070003 (same sequence without 1, sorted into ascending order).
Cf. A000040, A006530, A052126, A061395, A122111, A253550, A253561, A253563, A253565.
Differs from A072995 for the first time at n=15, where a(15) = 75, while A072995(15) = 225.

A253560 := proc(n)
    n*A006530(n) # code re-use
end proc:
seq( A253560(n),n=1..80) ; # R. J. Mathar, Jun 27 2024

a[n_] := n FactorInteger[n][[-1, 1]];
Array[a, 100] (* Jean-François Alcover, Feb 15 2021 *)

a(n) = if (n==1, 1, n*vecmax(factor(n)[,1])); \\ Michel Marcus, Feb 15 2021

Scheme
```
(define (A253560 n) (* (A006530 n) n))
```

a(1) = 1; for n > 1, a(n) = A006530(n) * n = A000040(A061395(n)) * n.
Other identities:
a(n) >= A253550(n) for all n >= 1.
a(n) = A129598(n) for all n >= 2.
A052126(a(n)) = n. [A052126 works as an inverse function for this injection.]

A243505 Permutation of natural numbers, take the odd bisection of A122111 and divide the largest prime factor out: a(n) = A052126(A122111(2n-1)).

Original entry on oeis.org

1, 2, 4, 8, 3, 16, 32, 6, 64, 128, 12, 256, 9, 5, 512, 1024, 24, 18, 2048, 48, 4096, 8192, 10, 16384, 27, 96, 32768, 36, 192, 65536, 131072, 20, 72, 262144, 384, 524288, 1048576, 15, 54, 2097152, 7, 4194304, 144, 768, 8388608, 108, 1536, 288, 16777216, 40, 33554432, 67108864, 30

Offset: 1

Antti Karttunen, Jun 25 2014

nonn

Antti Karttunen, Table of n, a(n) for n = 1..1024
Indranil Ghosh, Python program for computing the sequence
Index entries for sequences that are permutations of the natural numbers

Inverse: A243506.
Cf. A000040, A000079, A006254, A007051, A052126, A105560.
Related or similar permutations: A122111, A064216, A241916, A243065-A243066, A244981-A244982, A244983-A244984, A244153-A244154.

A052126[n_] := n/FactorInteger[n][[-1, 1]];
A122111[n_] := Product[Prime[Sum[If[j < i, 0, 1], {j, #}]], {i, Max[#]}]&[ Flatten[Table[Table[PrimePi[f[[1]]], {f[[2]]}], {f, FactorInteger[n]}]]];
a[n_] := A052126[A122111[2n-1]];
Array[a, 60] (* Jean-François Alcover, Sep 23 2020 *)

(define (A243505 n) (A122111 (A064216 n)))

a(n) = A052126(A122111((2*n)-1)).
a(n) = A122111((2*n)-1) / A105560((2*n)-1).
As a composition of related permutations:
a(n) = A122111(A064216(n)).
a(n) = A241916(A243065(n)).
Other identities:
For all n >= 2, a(n) = A070003(A244984(n)-1) / A105560((2*n)-1).
For all n >= 1, a(A006254(n)) = A000079(n) and a(A007051(n)) = A000040(n).
For all n >= 1, A105560(2n-1) divides a(n).

A241917 If n is a prime with index i, p_i, a(n) = i, (with a(1)=0), otherwise difference (i-j) of the indices of the two largest primes p_i, p_j, i >= j in the prime factorization of n: a(n) = A061395(n) - A061395(A052126(n)).

Original entry on oeis.org

0, 1, 2, 0, 3, 1, 4, 0, 0, 2, 5, 1, 6, 3, 1, 0, 7, 0, 8, 2, 2, 4, 9, 1, 0, 5, 0, 3, 10, 1, 11, 0, 3, 6, 1, 0, 12, 7, 4, 2, 13, 2, 14, 4, 1, 8, 15, 1, 0, 0, 5, 5, 16, 0, 2, 3, 6, 9, 17, 1, 18, 10, 2, 0, 3, 3, 19, 6, 7, 1, 20, 0, 21, 11, 0, 7, 1, 4, 22, 2, 0, 12, 23

Offset: 1

Antti Karttunen, May 13 2014

nonn

Note: the two largest primes in the multiset of prime divisors of n are equal for all numbers that are in A070003, thus, after a(1)=0, A070003 gives the positions of the other zeros in this sequence.

Cf. A052126, A061395, A070003, A122111, A241909.
Cf. A049084, A027746.
Cf. A241919, A242411, A243055 for other variants.

a241917 n = i - j where
            (i:j:_) = map a049084 $ reverse (1 : a027746_row n)
-- Reinhard Zumkeller, May 15 2014

A241917(n) = if(isprime(n), primepi(n), if(1>=omega(n), 0, my(f=factor(n)); if(f[#f~,2]>1, 0, primepi(f[#f~,1])-primepi(f[(#f~)-1,1])))); \\ Antti Karttunen, Jul 10 2024

from sympy import primefactors, primepi
def a061395(n): return 0 if n==1 else primepi(primefactors(n)[-1])
def a052126(n): return 1 if n==1 else n/primefactors(n)[-1]
def a(n): return 0 if n==1 else a061395(n) - a061395(a052126(n)) # Indranil Ghosh, May 19 2017

(define (A241917 n) (- (A061395 n) (A061395 (A052126 n))))

a(n) = A061395(n) - A061395(A052126(n)).

A070812 a(n) = phi(gpf(n)) - gpf(phi(n)) = A000010(A006530(n)) - A006530(A000010(n)).

Original entry on oeis.org

0, -1, 2, 0, 3, -1, -1, 2, 5, 0, 9, 3, 2, -1, 14, -1, 15, 2, 3, 5, 11, 0, -1, 9, -1, 3, 21, 2, 25, -1, 5, 14, 3, -1, 33, 15, 9, 2, 35, 3, 35, 5, 1, 11, 23, 0, -1, -1, 14, 9, 39, -1, 5, 3, 15, 21, 29, 2, 55, 25, 3, -1, 9, 5, 55, 14, 11, 3, 63, -1, 69, 33, -1, 15, 5, 9, 65, 2, -1, 35, 41, 3, 14, 35, 21, 5, 77, 1, 9, 11, 25, 23, 15, 0, 93, -1, 5

Offset: 3

Labos Elemer, May 09 2002

sign

Value of commutator[A000010, A006530] at n.

			Cases of n when a(n) = 1, -1, 2 or 0 are listed in A070002, A070003, A070004, A007283 respectively. Further regular solutions: if a(n)=3, then n=7k, where k has prime divisors < 7; if a(n)=5, then n=11k, where k has no prime divisors >=11; if a(n)=25, then mostly (not always!) n=31k ...

Michael De Vlieger, Table of n, a(n) for n = 3..10000
Michael De Vlieger, Log log scatterplot of a(n)+2, n = 3..10^5.

Cf. A000010, A006530, A070777, A068211, A070002, A070003, A070004, A007283.

pf[x_] := Part[Reverse[Flatten[FactorInteger[x]]], 2] Table[EulerPhi[pf[u]]-pf[EulerPhi[u]], {u, 3, 128}]

gpf(n)=my(f=factor(n)[,1]);f[#f]
a(n)=gpf(n)-gpf(eulerphi(n))-1 \\ Charles R Greathouse IV, Feb 19 2013

a(n) = A070777(n) - A068211(n).

cp's OEIS Frontend

A122111 Self-inverse permutation of the positive integers induced by partition enumeration in A112798 and partition conjugation.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Scheme

Scheme

Scheme

Formula

A071178 Exponent of the largest prime factor of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

A356862 Numbers with a unique largest prime exponent.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Python

A362612 Number of integer partitions of n such that the greatest part is the unique mode.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A102750 Numbers n such that square of largest prime dividing n does not divide n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A241916 a(2^k) = 2^k, and for other numbers, if n = 2^e1 * 3^e2 * 5^e3 * ... p_k^e_k, then a(n) = 2^(e_k - 1) * 3^(e_{k-1}) * ... * p_{k-1}^e2 * p_k^(e1+1). Here p_k is the greatest prime factor of n (A006530), and e_k is its exponent (A071178), and the exponents e1, ..., e_{k-1} >= 0.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Scheme

Formula