A102750 -id:A102750 - cp's OEIS frontend

A359178 Numbers with a unique smallest prime exponent.

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, 61, 63, 64, 67, 68, 71, 72, 73, 75, 76, 79, 80, 81, 83, 88, 89, 92, 96, 97, 98, 99, 101, 103, 104, 107, 108, 109, 112, 113, 116, 117

Offset: 1

Jens Ahlström, Jan 08 2023

nonn

180 is the smallest number with a unique smallest prime exponent that is not a member of A130091.

			2 = 2^1 is a term since it has 1 as a unique smallest exponent.
6 = 2^1 * 3^1 is not a term since it has two primes with the same smallest exponent.
180 = 2^2 * 3^2 * 5^1 is a term since it has 1 as a unique smallest exponent.

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

For parts instead of multiplicities we have A247180, counted by A002865.
For greatest instead of smallest we have A356862, counted by A362608.
The complement is A362606, counted by A362609.
Partitions of this type are counted by A362610.
These are the positions of 1's in A362613, for modes A362611.
A001221 counts prime exponents and A001222 adds them up.
A027746 lists prime factors, A112798 indices, A124010 exponents.
Cf. A102750, A130091, A327473, A327476, A359908, A362605, A362615.

q[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, Count[e, Min[e]] == 1]; Select[Range[2, 200], q] (* Amiram Eldar, Jan 08 2023 *)

isok(n) = if (n>1, my(f=factor(n), e = vecmin(f[,2])); #select(x->(x==e), f[,2], 1) == 1); \\ Michel Marcus, Jan 27 2023

from sympy import factorint
def ok(k):
  c = sorted(factorint(k).values())
  return len(c) == 1 or c[0] != c[1]
print([k for k in range(2, 118) if ok(k)])

from itertools import count, islice
from sympy import factorint
def A359178_gen(startvalue=2): # generator of terms >= startvalue
    return filter(lambda n:(f:=list(factorint(n).values())).count(min(f))==1,count(max(startvalue,2)))
A359178_list = list(islice(A359178_gen(),20)) # Chai Wah Wu, Feb 08 2023

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

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

A253550 Shift one instance of the largest prime one step towards larger primes: a(1) = 1, for n>1: a(n) = (n / prime(g)) * prime(g+1), where g = A061395(n), index of the greatest prime dividing n.

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Jan 03 2015

nonn

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

Inverse: A252462.
Cf. A102750 (same terms, but with 2 instead of 1, sorted into ascending order).
Cf. A003961, A052126, A065091, A061395, A122111, A253553, A253554, A253560, A253563, A253565.

(define (A253550 n) (if (= 1 n) n (* (A065091 (A061395 n)) (A052126 n))))

a(1) = 1; for n>1: a(n) = A065091(A061395(n)) * A052126(n).
Other identities. For all n >= 1:
A252462(a(n)) = n. [A252462 works as an inverse function for this injection.]
a(n) <= A253560(n).

A348041 Square array read by antidiagonals. A(n,k) is the nearest common ancestor of n and k in the Doudna tree (A005940).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 2, 3, 2, 1, 1, 2, 2, 2, 2, 1, 1, 2, 3, 4, 3, 2, 1, 1, 2, 3, 2, 2, 3, 2, 1, 1, 2, 3, 2, 5, 2, 3, 2, 1, 1, 2, 2, 2, 3, 3, 2, 2, 2, 1, 1, 2, 2, 4, 5, 6, 5, 4, 2, 2, 1, 1, 2, 3, 4, 2, 3, 3, 2, 4, 3, 2, 1, 1, 2, 3, 2, 2, 2, 7, 2, 2, 2, 3, 2, 1, 1, 2, 3, 2, 5, 2, 2, 2, 2, 5, 2, 3, 2, 1

Offset: 1

Antti Karttunen and Peter Munn, Sep 27 2021

Array is symmetric and is read by antidiagonals as A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), ... .
Also the nearest common ancestor of n and k in the tree depicted in A163511 (the mirror image of the Doudna tree).
The first fork in the Doudna tree separates numbers divisible by the square of their largest prime factor (on one main branch) from other numbers greater than 2 (on the other main branch). If n and m are on different main branches then A(n, m) = 2.
In more general terms A(.,.) can be considered as a binary operation that evaluates certain differences between the prime factors of its operands. To start, compare the largest prime factor of each operand with the 2nd largest prime factor. As described above, 2 is the result if these 2 factors are the same in one operand, but are different in the other operand; otherwise 3 is the result if these 2 factors are consecutive primes in one operand, but are nonconsecutive primes in the other operand. Further cases are covered in the examples, but note it is the difference between the indices of the prime numbers that is significant.

			The top left 17x17 corner of the array:
  n/k |  1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17
------+-------------------------------------------------------------
    1 |  1, 1, 1, 1, 1, 1, 1, 1, 1,  1,  1,  1,  1,  1,  1,  1,  1,
    2 |  1, 2, 2, 2, 2, 2, 2, 2, 2,  2,  2,  2,  2,  2,  2,  2,  2,
    3 |  1, 2, 3, 2, 3, 3, 3, 2, 2,  3,  3,  3,  3,  3,  3,  2,  3,
    4 |  1, 2, 2, 4, 2, 2, 2, 4, 4,  2,  2,  2,  2,  2,  2,  4,  2,
    5 |  1, 2, 3, 2, 5, 3, 5, 2, 2,  5,  5,  3,  5,  5,  3,  2,  5,
    6 |  1, 2, 3, 2, 3, 6, 3, 2, 2,  3,  3,  6,  3,  3,  6,  2,  3,
    7 |  1, 2, 3, 2, 5, 3, 7, 2, 2,  5,  7,  3,  7,  7,  3,  2,  7,
    8 |  1, 2, 2, 4, 2, 2, 2, 8, 4,  2,  2,  2,  2,  2,  2,  8,  2,
    9 |  1, 2, 2, 4, 2, 2, 2, 4, 9,  2,  2,  2,  2,  2,  2,  4,  2,
   10 |  1, 2, 3, 2, 5, 3, 5, 2, 2, 10,  5,  3,  5,  5,  3,  2,  5,
   11 |  1, 2, 3, 2, 5, 3, 7, 2, 2,  5, 11,  3, 11,  7,  3,  2, 11,
   12 |  1, 2, 3, 2, 3, 6, 3, 2, 2,  3,  3, 12,  3,  3,  6,  2,  3,
   13 |  1, 2, 3, 2, 5, 3, 7, 2, 2,  5, 11,  3, 13,  7,  3,  2, 13,
   14 |  1, 2, 3, 2, 5, 3, 7, 2, 2,  5,  7,  3,  7, 14,  3,  2,  7,
   15 |  1, 2, 3, 2, 3, 6, 3, 2, 2,  3,  3,  6,  3,  3, 15,  2,  3,
   16 |  1, 2, 2, 4, 2, 2, 2, 8, 4,  2,  2,  2,  2,  2,  2, 16,  2,
   17 |  1, 2, 3, 2, 5, 3, 7, 2, 2,  5, 11,  3, 13,  7,  3,  2, 17,
.
The nearest common ancestor of 7 and 15 in the Doudna tree (see diagram in the links and A005940) is 3, thus A(7,15) = A(15,7) = 3.
The nearest common ancestor of 12 and 15 in the Doudna tree is 6, thus A(12,15) = A(15,12) = 6.
The nearest common ancestor of 4 and 27 is 4 because 27 is a descendant of 4 in the Doudna tree, thus A(4,27) = A(27,4) = 4.
Example without reference to the Doudna tree: (Start)
The method below works in general for A(.,.) considered as a binary operation, but we use A(20, 42) as our example.
(1) Write each operand as a product of primes in nondecreasing order, convert to a tuple of prime indices, decrement each index, take first differences, then reverse the order:
  20 = 2*2*5 = prime(1) * prime(1) * prime(3) -> (1,1,3) -> (0,0,2) -> (0,0,2) -> (2,0,0);
  42 = 2*3*7 = prime(1) * prime(2) * prime(4) -> (1,2,4) -> (0,1,3) -> (0,1,2) -> (2,1,0).
(2) Truncate each tuple after the first elements that differ between them (or at the length of the shorter tuple):
  (2,0,0) -> (2,0); (2,1,0) -> (2,1).
(3) Choose the lesser tuple: (2,0).
(4) Determine which number would generate this tuple by the process from step (1):
  10 = 2*5 = prime(1) * prime(3) -> (1,3) -> (0,2) -> (0,2) -> (2,0).
This gives A(20, 42) = 10.
(End)

Michael De Vlieger, Doudna Tree diagram showing 8 rows.
Index entries for sequences computed from indices in prime factorization

Cf. A000027 (main diagonal).
Cf. A000040, A005940, A052126, A061395, A064989, A070003, A102750, A156552, A163511, A241917, A347879 [= a(n,sigma(n))], A348040, A348041, A348042, A348043, A348044 [= a(n,n^2)].
Cf. also A341510, A347380, A347381.

up_to = 105;
Abincompreflen(n, m) = { my(x=binary(n),y=binary(m),u=min(#x,#y)); for(i=1,u,if(x[i]!=y[i],return(i-1))); (u);};
Abinprefix(n,k) = { my(digs=binary(n)); fromdigits(vector(k,i,digs[i]),2); };
A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); (t); };
A156552(n) = {my(f = factor(n), p, p2 = 1, res = 0); for(i = 1, #f~, p = 1 << (primepi(f[i, 1]) - 1); res += (p * p2 * (2^(f[i, 2]) - 1)); p2 <<= f[i, 2]); res}; \\ From A156552
A348040sq(x,y) = Abincompreflen(A156552(x), A156552(y));
A348041sq(x,y) = A005940(1+Abinprefix(A156552(x),A348040sq(x,y)));
A348041list(up_to) = { my(v = vector(up_to), i=0); for(a=1,oo, for(col=1,a, i++; if(i > up_to, return(v)); v[i] = A348041sq(col,(a-(col-1))))); (v); };
v348041 = A348041list(up_to);
A348041(n) = v348041[n];

\\ A348041sq can be defined also as:
A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
A252463(n) = if(!(n%2),n/2,A064989(n));
A348041sq(x,y) = if(1==x||1==y,1, my(lista=List([]), i, k=x, stemvec, h=y); while(k>1, listput(lista,k); k = A252463(k)); stemvec = Vecrev(Vec(lista)); while(1, if((i=vecsearch(stemvec,h))>0, return(stemvec[i])); h = A252463(h)));

A(n, 1) = A(1, n) = 1; otherwise if A241917(n) <> A241917(m) then A(n, m) = A000040(1 + min(A241917(2*n), A241917(2*m))); otherwise A(n, m) = x * A000040(A061395(x)+A241917(n)), where x = A(A052126(n), A052126(m)).
A(i, j) = A(j, i).
A(n, n) = n.
A(2, n) = 2 for all n > 1.
A(p, q) = min(p, q) for any primes p and q.
A(A070003(n), A102750(m)) = 2.
A(u^2, v^2) = A(u, v)^2.
A(4k+2, 6k+3) = A064989(2k+1) for all k >= 1.

A329888 a(n) = A329900(A329602(n)); Heinz number of the even bisection (even-indexed parts) of the integer partition with Heinz number n.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 2, 3, 2, 1, 2, 1, 2, 3, 4, 1, 3, 1, 2, 3, 2, 1, 4, 5, 2, 3, 2, 1, 3, 1, 4, 3, 2, 5, 6, 1, 2, 3, 4, 1, 3, 1, 2, 3, 2, 1, 4, 7, 5, 3, 2, 1, 6, 5, 4, 3, 2, 1, 6, 1, 2, 3, 8, 5, 3, 1, 2, 3, 5, 1, 6, 1, 2, 5, 2, 7, 3, 1, 4, 9, 2, 1, 6, 5, 2, 3, 4, 1, 6, 7, 2, 3, 2, 5, 8, 1, 7, 3, 10, 1, 3, 1, 4, 5

Offset: 1

Antti Karttunen, Dec 22 2019

nonn

From Gus Wiseman, Aug 05 2021 and Antti Karttunen, Oct 13 2021: (Start)
Also the product of primes at even positions in the weakly decreasing list (with multiplicity) of prime factors of n. For example, the prime factors of 108 are (3,3,3,2,2), with even bisection (3,2), with product 6, so a(108) = 6.
Proof: A108951(n) gives a number with the same largest prime factor (A006530) and its exponent (A071178) as in n, and with each smaller prime p = 2, 3, 5, 7, ... < A006530(n) having as its exponent the partial sum of the exponents of all prime factors >= p present in n (with primes not present in n having the exponent 0). Then applying A000188 replaces each such "partial sum exponent" k with floor(k/2). Finally, A319626 replaces those halved exponents with their first differences (here the exponent of the largest prime present stays intact, because the next larger prime's exponent is 0 in n). It should be easy to see that if prime q is not present in n (i.e., does not divide it), then neither it is present in a(n). Moreover, if the partial sum exponent of q is odd and only one larger than the partial sum exponent of the next larger prime factor of n, then q will not be present in a(n), while in all other cases q is present in a(n). See also the last example.
(End)

			From _Gus Wiseman_, Aug 15 2021: (Start)
The list of all numbers with image 12 and their corresponding prime factors begins:
  144: (3,3,2,2,2,2)
  216: (3,3,3,2,2,2)
  240: (5,3,2,2,2,2)
  288: (3,3,2,2,2,2,2)
  336: (7,3,2,2,2,2)
  360: (5,3,3,2,2,2)
(End)
The positions from the left are indexed as 1, 2, 3, ..., etc, so e.g., for 240 we pick the second, the fourth and the sixth prime factor, 3, 2 and 2, to obtain a(240) = 3*2*2 = 12. For 288, we similarly pick the second (3), the fourth (2) and the sixth (2) to obtain a(288) = 3*2*2 = 12. - _Antti Karttunen_, Oct 13 2021
Consider n = 11945934 = 2*3*3*3*7*11*13*13*17. Its primorial inflation is A108951(11945934) = 96478365991115908800000 = 2^9 * 3^8 * 5^5 * 7^5 * 11^4 * 13^3 * 17^1. Applying A000188 to this halves each exponent (floored down if the exponent is odd), leaving the factors 2^4 * 3^4 * 5^2 * 7^2 * 11^2 * 13^1 = 2497294800. Then applying A319626 to this number retains the largest prime factor (and its exponent), and subtracts from the exponent of each of the rest of primes the exponent of the next larger prime, so from 2^4 * 3^4 * 5^2 * 7^2 * 11^2 * 13^1 we get 2^(4-4) * 3^(4-2) * 5^(2-2) * 7^(2-2) * 11^(2-1) * 13^1 = 3^2 * 11^1 * 13^1 = 1287 = a(11945934), which is obtained also by selecting every second prime from the list [17, 13, 13, 11, 7, 3, 3, 3, 2] and taking their product. - _Antti Karttunen_, Oct 15 2021

Antti Karttunen, Table of n, a(n) for n = 1..16384
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537
Index entries for sequences related to Heinz numbers
Index entries for sequences related to primorial numbers

Cf. A000188, A108951, A319626, A329602, A329900.
A left inverse of A000290.
Positions of 1's are A008578.
Positions of primes are A168645.
The sum of prime indices of a(n) is A346700(n).
The odd version is A346701.
The odd non-reverse version is A346703.
The non-reverse version is A346704.
The version for standard compositions is A346705, odd A346702.
A001221 counts distinct prime factors.
A001222 counts all prime factors.
A001414 adds up prime factors, row sums of A027746.
A027187 counts partitions of even length, ranked by A028260.
A056239 adds up prime indices, row sums of A112798.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A346633 adds up the even bisection of standard compositions.
A346698 adds up the even bisection of prime indices.
Cf. A000097, A035363, A102750, A236913, A253553, A344606, A344617, A344653, A345957, A345958, A345959.

Table[Times@@Last/@Partition[Reverse[Flatten[Apply[ConstantArray,FactorInteger[n],{1}]]],2],{n,100}] (* Gus Wiseman, Oct 13 2021 *)

PARI
```
A329888(n) = A329900(A329602(n));
```

A329888(n) = if(1==n,n,my(f=factor(n),m=1,p=0); forstep(k=#f~,1,-1,while(f[k,2], m *= f[k,1]^(p%2); f[k,2]--; p++)); (m)); \\ (After Wiseman's new interpretation) - Antti Karttunen, Sep 21 2021

a(n) = A329900(A329602(n)) = A329900(A000188(A108951(n))).
A108951(a(n)) = A329602(n).
a(n^2) = n for all n >= 1.
a(n) * A346701(n) = n. - Gus Wiseman, Aug 07 2021
A056239(a(n)) = A346700(n). - Gus Wiseman, Aug 07 2021
Antti Karttunen, Sep 21 2021
From Antti Karttunen, Oct 13 2021: (Start)
a(n) = A319626(A329602(n)) = A319626(A000188(A108951(n))).
For all x in A102750, a(x) = a(A253553(x)). (End)

Name amended with Gus Wiseman's new interpretation - Antti Karttunen, Oct 13 2021

A034876 Number of ways to write n! as a product of smaller factorials each greater than 1.

Original entry on oeis.org

0, 0, 0, 1, 0, 1, 0, 1, 1, 2, 0, 1, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0

Offset: 1

Erich Friedman

By definition, a(n) > 0 if and only if n is a member of A034878. If n > 2, then a(n!) > max(a(n), a(n!-1)), as (n!)! = n!*(n!-1)!. Similarly, a(A001013(n)) > 0 for n > 2. Clearly a(n)=0 if n is a prime A000040. So a(n+1)=1 if n=2^p-1 is a Mersenne prime A000668, as (n+1)!=(2!)^p*n! and n is prime. - Jonathan Sondow, Dec 15 2004
From Antti Karttunen, Dec 25 2018: (Start)
If n! = a! * x! * y! * ... * z!, with a > x >= y >= z, then A006530(n!) = A006530(a!) > A006530(x!). This follows because all rows in A115627 end with 1, that is, because all factorials >= 2 are in A102750.
If all the two-term solutions are of the form n! = a! * x! = b! * y! = ... = c! * z! (that is, all are products of two factorials larger than one), with a > x, b > y, ..., c > z, then a(n) = (a(x)+1 + a(y)+1 + ... + a(z)+1).
Values 0..5 occur for the first time at n = 1, 4, 10, 576, 13824, 69120.
In range 1..69120 differs from A322583 only at positions n = 1, 2, 9, 10 and 16.
(End)

			a(10) = 2 because 10! = 3! * 5! * 7! = 6! * 7! are the only two ways to write 10! as a product of smaller factorials > 1.
From _Antti Karttunen_, Dec 25 2018: (Start)
a(8) = 1 because 8! = 7! * (2!)^3.
a(9) = 1 because 9! = 7! * 3! * 3! * 2!.
a(16) = 2 because 16! = 15! * (2!)^4 = 14! * 5! * 2!.
a(144) = 2 because 144! = 143! * 4! * 3! = 143! * 3! * 3! * 2! * 2!.
a(576) = 3 because 576! = 575! * 4! * 4! = 575! * 4! * 3! * 2! * 2! = 575! * 3! * 3! * 2! * 2! * 2! * 2!.
a(720) = 2 because 720! = 719! * 6! = 719! * 5! * 3!.
a(3456) = 3 because 3456! = 3455! * 4! * 4! * 3! = 3455! * 4! * 3! * 3! * 2! * 2! = 3455! * 3! * 3! * 3! * 2! * 2! * 2! * 2!.
(End)

R. K. Guy, Unsolved Problems in Number Theory, B23.

Antti Karttunen, Table of n, a(n) for n = 1..69120
Eric Weisstein's World of Mathematics, Factorial Products
Index entries for sequences related to factorial numbers

Cf. A034878, A001013, A075082, A085604, A115627, A322583.

A034876aux(n, m, p) = if(1==n, 1, my(s=0); forstep(i=m, p, -1, my(f=i!); if(!(n%f), s += A034876aux(n/f, i, 2))); (s));
A034876(n) = if(1==n,0,A034876aux(n!, n-1, precprime(n))); \\ (Slow) - Antti Karttunen, Dec 24 2018

A322583aux(n, m) = if(1==n, 1, my(s=0); for(i=2, oo, my(f=i!); if(f>m, return(s)); if(!(n%f), s += A322583aux(n/f, f))));
memoA322583 = Map();
A322583(n) = { my(c); if(mapisdefined(memoA322583,n,&c), c, c = A322583aux(n,n); mapput(memoA322583,n,c); (c)); };
A034876aux(n, m, p) = if(1==n, 1, my(s=0); forstep(i=m, p, -1, my(f=i!); s += A322583(n/f)); (s));
A034876(n) = if(1==n, 0, A034876aux(n!, n-1, precprime(n))); \\ Antti Karttunen, Dec 25 2018

a(1) = 0; for n > 1, a(n) = Sum_{x=A007917(n)..n-1} A322583(n!/x!) when n is a composite, and a(n) = 0 when n is a prime. - Antti Karttunen, Dec 25 2018

Corrected by Jonathan Sondow, Dec 18 2004

A243057 If n = p_a * p_b * ... * p_h * p_i * p_j * p_k, with p_a <= p_b <= ... <= p_k, then a(n) = p_{k-j} * p_{k-i} * p_{k-h} * ... * p_{k-a} * p_k, where p_{0} = 1 and for k>=1, p_{k} = A000040(k).

Original entry on oeis.org

1, 2, 3, 2, 5, 6, 7, 2, 3, 15, 11, 12, 13, 35, 10, 2, 17, 6, 19, 45, 21, 77, 23, 24, 5, 143, 3, 175, 29, 30, 31, 2, 55, 221, 14, 12, 37, 323, 91, 135, 41, 105, 43, 539, 20, 437, 47, 48, 7, 15, 187, 1573, 53, 6, 33, 875, 247, 667, 59, 90, 61, 899, 63, 2, 65, 385

Offset: 1

Antti Karttunen, May 31 2014

nonn

A243058 gives all n such that a(n) = n (the fixed points of this sequence, which include primes).
A102750 gives such n that a(a(n)) = n. A243286 is the self-inverse permutation induced when the domain is restricted to A102750. Cf. also A242420.
A070003 gives all n such that a(a(n)) <> n. Another variant, A243059, is defined to be zero when n is one of the terms of A070003.

			For n = 9 = 3*3 = p_2 * p_2, we have a(n) = p_{3-3} * p_3 = 1*3 = 3. [Like all terms in A070003 this is an example of "degenerate case", where some p's in the product get index 0, and thus are set to 1 by convention used here.]
For n = 10 = 2*5 = p_1 * p_3, we have a(n) = p_{3-1} * p_3 = 3*5 = 15.
For n = 12 = 2*2*3 = p_1 * p_1 * p_2, we have a(n) = p_{2-1} * p{2-1} * p_2 = p_1^2 * p_2 = 12.
For n = 15 = 3*5 = p_2 * p_3, we have a(n) = p_{3-2} * p_3 = 2*5 = 10.
For n = 2200 = 2*2*2*5*5*11 = p_1 * p_1 * p_1 * p_3 * p_3 * p_5, we have a(n) = p_{5-3} * p_{5-3} * p_{5-1} * p_{5-1} * p_{5-1} * p_5 = 3*3*7*7*7*11 = 33957.
For n = 33957 = 3*3*7*7*7*11 = p_2 * p_2 * p_4 * p_4 * p_4 * p_5, we have a(n) = p_{5-4} * p_{5-4} * p_{5-4} * p_{5-2} * p_{5-2} * p_5 = 2*2*2*5*5*11 = 2200.

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

Fixed points: A243058 (includes primes).

Cf. A243059, A000040, A032742, A243056, A242419, A242420, A243286, A070003, A102750, A243074.

If n = p_a * p_b * ... * p_h * p_i * p_j * p_k, where p_a <= p_b <= ... <= p_k are (not necessarily distinct) primes (sorted into nondescending order) in the prime factorization of n, then a(n) = p_{k-j} * p_{k-i} * p_{k-h} * ... * p_{k-a} * p_k, where p_{0} = 1 and for k>=1, p_{k} = A000040(k).
a(1)=1, and for n>1, a(n) = p_{A243056(n)} * a(A032742(n)). Here p_{k} stands for 1 when k=0, and otherwise for the k-th prime, A000040(k).
For all n, a(n) = a(A243074(n)).
For all k in A102750, a(k) = A242420(k).

A243285 Number of integers 1 <= k <= n which are not divisible by the square of their largest noncomposite divisor.

Original entry on oeis.org

0, 1, 2, 2, 3, 4, 5, 5, 5, 6, 7, 8, 9, 10, 11, 11, 12, 12, 13, 14, 15, 16, 17, 18, 18, 19, 19, 20, 21, 22, 23, 23, 24, 25, 26, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 38, 38, 39, 40, 41, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 50, 51, 52, 53, 54, 55, 56, 57

Offset: 1

Antti Karttunen, Jun 02 2014

nonn

a(n) tells how many natural numbers <= n there are which are not divisible by the square of their largest noncomposite divisor.
The largest noncomposite divisor of 1 is 1 itself, and 1 is divisible by 1^2, thus 1 is not included in the count, and a(1)=0.
The "largest noncomposite divisor" for any integer > 1 means the same thing as the largest prime divisor, and thus we are counting the terms of A102750 (Numbers n such that square of largest prime dividing n does not divide n).
Thus this is the partial sums of the characteric function for A102750.

			For n = 9, there are numbers 2, 3, 5, 6 and 7 which are not divisible by the square of their largest prime factor, while 1 is excluded (no prime factors) and 4 and 8 are divisible both by 2^2 and 9 is divisible by 3^2. Thus a(9) = 5.

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

Cf. A243282, A243283, A102750, A070003, A013928, A057627.

ndsQ[n_]:=Mod[n,Max[Select[Divisors[n],!CompositeQ[#]&]]^2]!=0; Accumulate[Table[If[ ndsQ[n],1,0],{n,80}]] (* Harvey P. Dale, Oct 14 2023 *)

from sympy import primefactors
def a243285(n): return 0 if n==1 else sum([1 for k in range(2, n + 1) if k%(primefactors(k)[-1]**2)!=0]) # Indranil Ghosh, Jun 15 2017

Scheme
```
(define (A243285 n) (- n (A243283 n)))
```

a(n) = n - A243283(n).

For all n, a(A102750(n)) = n, thus this sequence works also as an inverse function for the injection A102750.

A242420 Self-inverse permutation of positive integers: a(n) = (A006530(n)^(A071178(n)-1)) * A243057(n).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 15, 11, 12, 13, 35, 10, 16, 17, 18, 19, 45, 21, 77, 23, 24, 25, 143, 27, 175, 29, 30, 31, 32, 55, 221, 14, 36, 37, 323, 91, 135, 41, 105, 43, 539, 20, 437, 47, 48, 49, 75, 187, 1573, 53, 54, 33, 875, 247, 667, 59, 90, 61, 899, 63, 64, 65

Offset: 1

Antti Karttunen, May 31 2014

nonn

This self-inverse permutation (involution) of positive integers preserves both the total number of prime divisors and the (index of) largest prime factor of n, i.e., for all n it holds that A001222(a(n)) = A001222(n) and A006530(a(n)) = A006530(n) [equally: A061395(a(n)) = A061395(n)].
It also preserves the exponent of the largest prime factor (A071178), from which follows that the sequence A102750 is closed with respect to this permutation, i.e., for all n in A102750, a(n) is either same n or some other term of A102750.
Considered as an operation on partitions encoded by the indices of primes in the prime factorization of n (as in table A112798), this implements a self-inverse bijection which is a composition of the effects of A242419 and A225891. (Or equally: A105119 and A242419). For details, please see the respective Comments sections and/or Example section of this entry.

			For n = 2200, we see that it encodes the partition (1,1,1,3,3,5) in A112798 as 2200 = p_1 * p_1 * p_1 * p_3 * p_3 * p_5 = 2^3 * 5^2 * 11. This in turn corresponds to the following Young diagram in French notation:
   _
  | |
  | |
  | |_ _
  |     |
  |     |_ _
  |_ _ _ _ _|
First we apply A242419, which reverses the order of "steps", so that each horizontal and vertical line segment centered around a "convex corner" moves as a whole, so that the first stair from the top (one unit wide and three units high) is moved to the last position, the second one (two units wide and two units high) stays in the middle, and the original bottom step (two units wide and one unit high) will be the new topmost step, thus we get the following Young diagram:
   _ _
  |   |_ _
  |       |
  |       |_
  |         |
  |         |
  |_ _ _ _ _|
which represents the partition (2,4,4,5,5,5), encoded in A112798 by p_2 * p_4^2 * p_5^3 = 3 * 7^2 * 11^3 = 195657.
Then we apply A225891, which rotates the exponents of distinct primes in the factorization of n one left, in this context the vertical line segments one step up, with the top-one going to the bottomost, and so we get:
   _ _
  |   |
  |   |_ _
  |       |
  |       |
  |       |_
  |_ _ _ _ _|
which represents the partition (2,2,4,4,4,5), encoded in A112798 by p_2^2 * p_4^3 * p_5 = 3^2 * 7^3 * 11 = 33957, thus a(2200) = 33957.

Antti Karttunen, Table of n, a(n) for n = 1..8192
Wikipedia, Young diagram
Index entries for sequences that are permutations of the natural numbers

Fixed points: A243068.

Cf. A112798, A242419, A105119, A225891, A006530, A071178, A102750, A243057, A243059.

(define (A242420 n) (* (expt (A006530 n) (- (A071178 n) 1)) (A243057 n)))

a(n) = (A006530(n)^(A071178(n)-1)) * A243057(n).
For all k in A102750, a(k) = A243057(k) = A243059(k).
By composing related permutations:
a(n) = A225891(A242419(n)) = A242419(A105119(n)).

cp's OEIS Frontend

A359178 Numbers with a unique smallest 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

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

Extensions

A253550 Shift one instance of the largest prime one step towards larger primes: a(1) = 1, for n>1: a(n) = (n / prime(g)) * prime(g+1), where g = A061395(n), index of the greatest prime dividing n.

Views

Author

Keywords

Links

Crossrefs

Programs

Scheme

Formula

A348041 Square array read by antidiagonals. A(n,k) is the nearest common ancestor of n and k in the Doudna tree (A005940).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

PARI

Formula

A329888 a(n) = A329900(A329602(n)); Heinz number of the even bisection (even-indexed parts) of the integer partition with Heinz number n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

Extensions

A034876 Number of ways to write n! as a product of smaller factorials each greater than 1.

Views

Author

Keywords

Comments

Examples