A307035 -id:A307035 - cp's OEIS frontend

A064989 Multiplicative with a(2^e) = 1 and a(p^e) = prevprime(p)^e for odd primes p.

1, 1, 2, 1, 3, 2, 5, 1, 4, 3, 7, 2, 11, 5, 6, 1, 13, 4, 17, 3, 10, 7, 19, 2, 9, 11, 8, 5, 23, 6, 29, 1, 14, 13, 15, 4, 31, 17, 22, 3, 37, 10, 41, 7, 12, 19, 43, 2, 25, 9, 26, 11, 47, 8, 21, 5, 34, 23, 53, 6, 59, 29, 20, 1, 33, 14, 61, 13, 38, 15, 67, 4, 71, 31, 18, 17, 35, 22, 73, 3, 16

Offset: 1

Vladeta Jovovic, Oct 30 2001

From Antti Karttunen, May 12 2014: (Start)
a(A003961(n)) = n for all n. [This is a left inverse function for the injection A003961.]
Bisections are A064216 (the terms at odd indices) and A064989 itself (the terms at even indices), i.e., a(2n) = a(n) for all n.
(End)
From Antti Karttunen, Dec 18-21 2014: (Start)
When n represents an unordered integer partition via the indices of primes present in its prime factorization (for n >= 2, n corresponds to the partition given as the n-th row of A112798) this operation subtracts one from each part. If n is of the form 2^k (a partition having just k 1's as its parts) the result is an empty partition (which is encoded by 1, having an "empty" factorization).
For all odd numbers n >= 3, a(n) tells which number is located immediately above n in square array A246278. Cf. also A246277.
(End)
Alternatively, if numbers are represented as the multiset of indices of prime factors with multiplicity, this operation subtracts 1 from each element and discards the 0's. - M. F. Hasler, Dec 29 2014

			a(20) = a(2^2*5) = a(2^2)*a(5) = prevprime(5) = 3.

Harry J. Smith (first 1000 terms) & Antti Karttunen, Table of n, a(n) for n = 1..10000

Cf. A064216 (odd bisection), A003961 (inverse), A151799.
Cf. A000040, A008578, A027746, A032742, A049084, A052126, A055396, A061395, A112798, A249817, A249818.
Other sequences whose definition involve or are some other way related with this sequence: A105560, A108951, A118306, A122111, A156552, A163511, A200746, A241909, A243070, A243071, A243072, A243073, A244319, A245605, A245607, A246165, A246266, A246268, A246277, A246278, A246361, A246362, A246371, A246372, A246373, A246374, A246376, A246380, A246675, A246682, A249745, A250470.
Similar prime-shifts towards smaller numbers: A252461, A252462, A252463.

a064989 1 = 1
a064989 n = product $ map (a008578 . a049084) $ a027746_row n
-- Reinhard Zumkeller, Apr 09 2012
(MIT/GNU Scheme, with Aubrey Jaffer's SLIB Scheme library)
(require 'factor)
(define (A064989 n) (if (= 1 n) n (apply * (map (lambda (k) (if (zero? k) 1 (A000040 k))) (map -1+ (map A049084 (factor n)))))))
;; Antti Karttunen, May 12 2014
(definec (A064989 n) (if (= 1 n) n (* (A008578 (A055396 n)) (A064989 (A032742 n))))) ;; One based on given recurrence and utilizing memoizing definec-macro.
(definec (A064989 n) (cond ((= 1 n) n) ((even? n) (A064989 (/ n 2))) (else (A163511 (/ (- (A243071 n) 1) 2))))) ;; Corresponds to one of the alternative formulas, but is very unpractical way to compute this sequence. - Antti Karttunen, Dec 18 2014

q:= proc(p) prevprime(p) end proc: q(2):= 1:
[seq(mul(q(f[1])^f[2], f = ifactors(n)[2]), n = 1 .. 1000)]; # Robert Israel, Dec 21 2014

Table[Times @@ Power[Which[# == 1, 1, # == 2, 1, True, NextPrime[#, -1]] & /@ First@ #, Last@ #] &@ Transpose@ FactorInteger@ n, {n, 81}] (* Michael De Vlieger, Jan 04 2016 *)

{ for (n=1, 1000, f=factor(n)~; a=1; j=1; if (n>1 && f[1, 1]==2, j=2); for (i=j, length(f), a*=precprime(f[1, i] - 1)^f[2, i]); write("b064989.txt", n, " ", a) ) } \\ Harry J. Smith, Oct 02 2009

a(n) = {my(f = factor(n)); for (i=1, #f~, if ((p=f[i,1]) % 2, f[i,1] = precprime(p-1), f[i,1] = 1);); factorback(f);} \\ Michel Marcus, Dec 18 2014

A064989(n)=factorback(Mat(apply(t->[max(precprime(t[1]-1),1),t[2]],Vec(factor(n)~))~)) \\ M. F. Hasler, Dec 29 2014

from sympy import factorint, prevprime
from operator import mul
from functools import reduce
def a(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])
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jun 15 2017

from math import prod
from sympy import prevprime, factorint
def A064989(n): return prod(prevprime(p)**e for p, e in  factorint(n>>(~n&n-1).bit_length()).items()) # Chai Wah Wu, Jan 05 2023

From Antti Karttunen, Dec 18 2014: (Start)
If n = product A000040(k)^e(k) then a(n) = product A008578(k)^e(k) [where A000040(n) gives the n-th prime, and A008578(n) gives 1 for 1 and otherwise the (n-1)-th prime].
a(1) = 1; for n > 1, a(n) = A008578(A055396(n)) * a(A032742(n)). [Above formula represented as a recurrence. Cf. A252461.]
a(1) = 1; for n > 1, a(n) = A008578(A061395(n)) * a(A052126(n)). [Compare to the formula of A252462.]
This prime-shift operation is used in the definitions of many other sequences, thus it can be expressed in many alternative ways:
a(n) = A200746(n) / n.
a(n) = A242424(n) / A105560(n).
a(n) = A122111(A122111(n)/A105560(n)) = A122111(A052126(A122111(n))). [In A112798-partition context: conjugate, remove the largest part (the largest prime factor), and conjugate again.]
a(1) = 1; for n > 1, a(2n) = a(n), a(2n+1) = A163511((A243071(2n+1)-1) / 2).
a(n) = A249818(A250470(A249817(n))). [A250470 is an analogous operation for "going one step up" in the square array A083221 (A083140).]
(End)
Product_{k=1..n} a(k) = n! / A307035(n). - Vaclav Kotesovec, Mar 21 2019
Sum_{k=1..n} a(k) ~ c * n^2, where c = (1/2) * Product_{p prime} ((p^2-p)/(p^2-q(p))) = 0.220703928... , where q(p) = prevprime(p) (A151799) if p > 2 and q(2) = 1. - Amiram Eldar, Nov 18 2022

A108951 Primorial inflation of n: Fully multiplicative with a(p) = p# for prime p, where x# is the primorial A034386(x).

Original entry on oeis.org

1, 2, 6, 4, 30, 12, 210, 8, 36, 60, 2310, 24, 30030, 420, 180, 16, 510510, 72, 9699690, 120, 1260, 4620, 223092870, 48, 900, 60060, 216, 840, 6469693230, 360, 200560490130, 32, 13860, 1021020, 6300, 144, 7420738134810, 19399380, 180180, 240, 304250263527210, 2520

Offset: 1

Paul Boddington, Jul 21 2005

This sequence is a permutation of A025487.
And thus also a permutation of A181812, see the formula section. - Antti Karttunen, Jul 21 2014
A previous description of this sequence was: "Multiplicative with a(p^e) equal to the product of the e-th powers of all primes at most p" (see extensions), Giuseppe Coppoletta, Feb 28 2015

			a(12) = a(2^2) * a(3) = (2#)^2 * (3#) = 2^2 * 6 = 24
a(45) = (3#)^2 * (5#) = (2*3)^2 * (2*3*5) = 1080 (as 45 = 3^2 * 5).

Amiram Eldar, Table of n, a(n) for n = 1..2370 (terms 1..256 from Antti Karttunen)
Index to divisibility sequences.
Index entries for sequences computed from indices in prime factorization.
Index entries for sequences related to primorial numbers.

Cf. A319626, A329900 (left inverses).

Cf. A034386, A002110, A025487, A048673, A064216, A064989, A085082, A122111, A124859, A161360, A181811, A181812, A181814, A181815, A181817, A181819, A181822, A238690, A283477, A283478, A307035, A324886, A324887, A324888, A324896, A325226, A329040, A329046, A329047, A329344, A329348, A329349, A329378, A329382, A329600, A329602, A329605, A329607, A329615, A329616, A329617, A329619, A329622, A319627, A329647, A331292, A337474, A346108, A346109, A344698, A344699.

a[n_] := a[n] = Module[{f = FactorInteger[n], p, e}, If[Length[f]>1, Times @@ a /@ Power @@@ f, {{p, e}} = f; Times @@ (Prime[Range[PrimePi[p]]]^e)]]; a[1] = 1; Table[a[n], {n, 1, 42}] (* Jean-François Alcover, Feb 24 2015 *)
Table[Times @@ Map[#1^#2 & @@ # &, FactorInteger[n] /. {p_, e_} /; e > 0 :> {Times @@ Prime@ Range@ PrimePi@ p, e}], {n, 42}] (* Michael De Vlieger, Mar 18 2017 *)

primorial(n)=prod(i=1,primepi(n),prime(i))
a(n)=my(f=factor(n)); prod(i=1,#f~, primorial(f[i,1])^f[i,2]) \\ Charles R Greathouse IV, Jun 28 2015

from sympy import primerange, factorint
from operator import mul
def P(n): return reduce(mul, [i for i in primerange(2, n + 1)])
def a(n):
    f = factorint(n)
    return 1 if n==1 else reduce(mul, [P(i)**f[i] for i in f])
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, May 14 2017

def sharp_primorial(n): return sloane.A002110(prime_pi(n))
def p(f):
    return sharp_primorial(f[0])^f[1]
[prod(p(f) for f in factor(n)) for n in range (1,51)]
# Giuseppe Coppoletta, Feb 07 2015

Dirichlet g.f.: 1/(1-2*2^(-s))/(1-6*3^(-s))/(1-30*5^(-s))...
Completely multiplicative with a(p_i) = A002110(i) = prime(i)#. [Franklin T. Adams-Watters, Jun 24 2009; typos corrected by Antti Karttunen, Jul 21 2014]
From Antti Karttunen, Jul 21 2014: (Start)
a(1) = 1, and for n > 1, a(n) = n * a(A064989(n)).
a(n) = n * A181811(n).
a(n) = A002110(A061395(n)) * A331188(n). - [added Jan 14 2020]
a(n) = A181812(A048673(n)).
Other identities:
A006530(a(n)) = A006530(n). [Preserves the largest prime factor of n.]
A071178(a(n)) = A071178(n). [And also its exponent.]
a(2^n) = 2^n. [Fixes the powers of two.]
A067029(a(n)) = A007814(a(n)) = A001222(n). [The exponent of the least prime of a(n), that prime always being 2 for n>1, is equal to the total number of prime factors in n.]
(End)
From Antti Karttunen, Nov 19 2019: (Start)
Further identities:
a(A307035(n)) = A000142(n).
a(A003418(n)) = A181814(n).
a(A025487(n)) = A181817(n).
a(A181820(n)) = A181822(n).
a(A019565(n)) = A283477(n).
A001221(a(n)) = A061395(n).
A001222(a(n)) = A056239(n).
A181819(a(n)) = A122111(n).
A124859(a(n)) = A181821(n).
A085082(a(n)) = A238690(n).
A328400(a(n)) = A329600(n). (smallest number with the same set of distinct prime exponents)
A000188(a(n)) = A329602(n). (square root of the greatest square divisor)
A072411(a(n)) = A329378(n). (LCM of exponents of prime factors)
A005361(a(n)) = A329382(n). (product of exponents of prime factors)
A290107(a(n)) = A329617(n). (product of distinct exponents of prime factors)
A000005(a(n)) = A329605(n). (number of divisors)
A071187(a(n)) = A329614(n). (smallest prime factor of number of divisors)
A267115(a(n)) = A329615(n). (bitwise-AND of exponents of prime factors)
A267116(a(n)) = A329616(n). (bitwise-OR of exponents of prime factors)
A268387(a(n)) = A329647(n). (bitwise-XOR of exponents of prime factors)
A276086(a(n)) = A324886(n). (prime product form of primorial base expansion)
A324580(a(n)) = A324887(n).
A276150(a(n)) = A324888(n). (digit sum in primorial base)
A267263(a(n)) = A329040(n). (number of distinct nonzero digits in primorial base)
A243055(a(n)) = A329343(n).
A276088(a(n)) = A329348(n). (least significant nonzero digit in primorial base)
A276153(a(n)) = A329349(n). (most significant nonzero digit in primorial base)
A328114(a(n)) = A329344(n). (maximal digit in primorial base)
A062977(a(n)) = A325226(n).
A097248(a(n)) = A283478(n).
A324895(a(n)) = A324896(n).
A324655(a(n)) = A329046(n).
A327860(a(n)) = A329047(n).
A329601(a(n)) = A329607(n).
(End)
a(A181815(n)) = A025487(n), and A319626(a(n)) = A329900(a(n)) = n. - Antti Karttunen, Dec 29 2019
From Antti Karttunen, Jul 09 2021: (Start)
a(n) = A346092(n) + A346093(n).
a(n) = A346108(n) - A346109(n).
a(A342012(n)) = A004490(n).
a(A337478(n)) = A336389(n).
A336835(a(n)) = A337474(n).
A342002(a(n)) = A342920(n).
A328571(a(n)) = A346091(n).
A328572(a(n)) = A344592(n).
(End)
Sum_{n>=1} 1/a(n) = A161360. - Amiram Eldar, Aug 04 2022

More terms computed by Antti Karttunen, Jul 21 2014
The name of the sequence was changed for more clarity, in accordance with the above remark of Franklin T. Adams-Watters (dated Jun 24 2009). It is implicitly understood that a(n) is then uniquely defined by completely multiplicative extension. - Giuseppe Coppoletta, Feb 28 2015
Name "Primorial inflation" (coined by Matthew Vandermast in A181815) prefixed to the name by Antti Karttunen, Jan 14 2020

A319626 Primorial deflation of n (numerator): Let f be the completely multiplicative function over the positive rational numbers defined by f(p) = A034386(p) for any prime number p; f constitutes a permutation of the positive rational numbers; let g be the inverse of f; for any n > 0, a(n) is the numerator of g(n).

Original entry on oeis.org

1, 2, 3, 4, 5, 3, 7, 8, 9, 10, 11, 6, 13, 14, 5, 16, 17, 9, 19, 20, 21, 22, 23, 12, 25, 26, 27, 28, 29, 5, 31, 32, 33, 34, 7, 9, 37, 38, 39, 40, 41, 21, 43, 44, 15, 46, 47, 24, 49, 50, 51, 52, 53, 27, 55, 56, 57, 58, 59, 10, 61, 62, 63, 64, 65, 33, 67, 68, 69

Offset: 1

Rémy Sigrist, Sep 25 2018

See A319627 for the corresponding denominators.
The restriction of f to the natural numbers corresponds to A108951.
The function g is completely multiplicative over the positive rational numbers with g(2) = 2 and g(q) = q/p for any pair (p, q) of consecutive prime numbers.
The ratio A319626(n)/A319627(n) can be viewed as a "primorial deflation" of n (see also A329900), with the inverse operation being n = A108951(A319626(n)) / A108951(A319627(n)), where A319627(k) = 1 for all k in A025487. - Daniel Suteu, Dec 29 2019

			f(21/5) = (2*3) * (2*3*5*7) / (2*3*5) = 42, hence g(42) = 21/5 and a(42) = 21.

Daniel Suteu, Table of n, a(n) for n = 1..10000

A left inverse of A108951. Coincides with A329900 on A025487.

Cf. A006530, A053585, A064989, A181815, A307035, A319627, A319630, A329902, A330749, A330750 (rgs-transform), A330751 (ordinal transform).

Array[#1/GCD[#1, #2] & @@ {#, Apply[Times, Map[If[#1 <= 2, 1, NextPrime[#1, -1]]^#2 & @@ # &, FactorInteger[#]]]} &, 120] (* Michael De Vlieger, Aug 27 2020 *)

a(n) = my (f=factor(n)); numerator(prod(i=1, #f~, my (p=f[i,1]); (p/if (p>2, precprime(p-1), 1))^f[i,2]))

a(n) = n / gcd(n, A064989(n)) = n / A330749(n).
a(n) <= n with equality iff n belongs to A319630.
A006530(a(n)) = A006530(n).
A053585(a(n)) = A053585(n).
From Antti Karttunen, Dec 29 2019: (Start)
a(A108951(n)) = n.
a(A025487(n)) = A329900(A025487(n)) = A181815(n).
Many of the formulas given in A329900 apply here as well:
a(n!) = A307035(n), a(A002182(n)) = A329902(n), and so on.
(End)

"Primorial deflation" prefixed to the name by Antti Karttunen, Dec 29 2019

A329900 Primorial deflation of n: starting from x = n, repeatedly divide x by the largest primorial A002110(k) that divides it, until x is an odd number. Then a(n) = Product prime(k_i), for primorial indices k_1 >= k_2 >= ..., encountered in the process.

Original entry on oeis.org

1, 2, 1, 4, 1, 3, 1, 8, 1, 2, 1, 6, 1, 2, 1, 16, 1, 3, 1, 4, 1, 2, 1, 12, 1, 2, 1, 4, 1, 5, 1, 32, 1, 2, 1, 9, 1, 2, 1, 8, 1, 3, 1, 4, 1, 2, 1, 24, 1, 2, 1, 4, 1, 3, 1, 8, 1, 2, 1, 10, 1, 2, 1, 64, 1, 3, 1, 4, 1, 2, 1, 18, 1, 2, 1, 4, 1, 3, 1, 16, 1, 2, 1, 6, 1, 2, 1, 8, 1, 5, 1, 4, 1, 2, 1, 48, 1, 2, 1, 4, 1, 3, 1, 8, 1

Offset: 1

Antti Karttunen, Dec 22 2019

nonn

When applied to arbitrary n, the "primorial deflation" (term coined by Matthew Vandermast in A181815) induces the splitting of n to two factors A328478(n)*A328479(n) = n, where we call A328478(n) the non-deflatable component of n (which is essentially discarded), while A328479(n) is the deflatable component. Only if n is in A025487, then the entire n is deflatable, i.e., A328478(n) = 1 and A328479(n) = n.

According to Daniel Suteu, also the ratio (A319626(n) / A319627(n)) can be viewed as a "primorial deflation". That definition coincides with this one when restricted to terms of A025487, as for all k in A025487, A319626(k) = a(k), and A319627(k) = 1. - Antti Karttunen, Dec 29 2019

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

A left inverse of A108951. Coincides with A319626 on A025487.

Cf. A002110, A002182, A111701, A181815, A181817, A181819, A181821, A276084, A304886, A319626, A319627, A328478, A328479, A329889, A329902, A330685, A330686, A330689.

Array[If[OddQ@ #, 1, Times @@ Prime@ # &@ Rest@ NestWhile[Append[#1, {#3, Drop[#, -LengthWhile[Reverse@ #, # == 0 &]] &[#2 - PadRight[ConstantArray[1, #3], Length@ #2]]}] & @@ {#1, #2, LengthWhile[#2, # > 0 &]} & @@ {#, #[[-1, -1]]} &, {{0, TakeWhile[If[# == 1, {0}, Function[f, ReplacePart[Table[0, {PrimePi[f[[-1, 1]]]}], #] &@ Map[PrimePi@ First@ # -> Last@ # &, f]]@ FactorInteger@ #], # > 0 &]}}, And[FreeQ[#[[-1, -1]], 0], Length[#[[-1, -1]] ] != 0] &][[All, 1]] ] &, 105] (* Michael De Vlieger, Dec 28 2019 *)
Array[Times @@ Prime@(TakeWhile[Reap[FixedPointList[Block[{k = 1}, While[Mod[#, Prime@ k] == 0, k++]; Sow[k - 1]; #/Product[Prime@ i, {i, k - 1}]] &, #]][[-1, 1]], # > 0 &]) &, 105] (* Michael De Vlieger, Jan 11 2020 *)

A329900(n) = { my(m=1, pp=1); while(1, forprime(p=2, ,if(n%p, if(2==p, return(m), break), n /= p; pp = p)); m *= pp); (m); };

A111701(n) = forprime(p=2, , if(n%p, return(n), n /= p));
A276084(n) = { for(i=1,oo,if(n%prime(i),return(i-1))); }
A329900(n) = if(n%2,1,prime(A276084(n))*A329900(A111701(n)));

For odd n, a(n) = 1, for even n, a(n) = A000040(A276084(n)) * a(A111701(n)).
For even n, a(n) = A000040(A276084(n)) * a(n/A002110(A276084(n))).
A108951(a(n)) = A328479(n), for n >= 1.
a(A108951(n)) = n, for n >= 1.
a(A328479(n)) = a(n),  for n >= 1.
a(A328478(n)) = 1, for n >= 1.
a(A002110(n)) = A000040(n), for n >= 1.
a(A000142(n)) = A307035(n), for n >= 0.
a(A283477(n)) = A019565(n), for n >= 0.
a(A329886(n)) = A005940(1+n), for n >= 0.
a(A329887(n)) = A163511(n), for n >= 0.
a(A329602(n)) = A329888(n), for n >= 1.
a(A025487(n)) = A181815(n), for n >= 1.
a(A124859(n)) = A181819(n), for n >= 1.
a(A181817(n)) = A025487(n), for n >= 1.
a(A181821(n)) = A122111(n), for n >= 1.
a(A002182(n)) = A329902(n), for n >= 1.
a(A260633(n)) = A329889(n), for n >= 1.
a(A033833(n)) = A330685(n), for n >= 1.
a(A307866(1+n)) = A330686(n), for n >= 1.
a(A330687(n)) = A330689(n), for n >= 1.

A319627 Primorial deflation of n (denominator): Let f be the completely multiplicative function over the positive rational numbers defined by f(p) = A034386(p) for any prime number p; f constitutes a permutation of the positive rational numbers; let g be the inverse of f; for any n > 0, a(n) is the denominator of g(n).

Original entry on oeis.org

1, 1, 2, 1, 3, 1, 5, 1, 4, 3, 7, 1, 11, 5, 2, 1, 13, 2, 17, 3, 10, 7, 19, 1, 9, 11, 8, 5, 23, 1, 29, 1, 14, 13, 3, 1, 31, 17, 22, 3, 37, 5, 41, 7, 4, 19, 43, 1, 25, 9, 26, 11, 47, 4, 21, 5, 34, 23, 53, 1, 59, 29, 20, 1, 33, 7, 61, 13, 38, 3, 67, 1, 71, 31, 6

Offset: 1

Rémy Sigrist, Sep 25 2018

See A319626 for the corresponding numerators and additional comments.

			f(21/5) = (2*3) * (2*3*5*7) / (2*3*5) = 42, hence g(42) = 21/5 and a(42) = 5.

Daniel Suteu, Table of n, a(n) for n = 1..10000

Cf. A025487 (positions of 1's), A064989, A329900, A358217 [= bigomega(a(n))].
Cf. A319626 (numerators, see comments there).
Cf. also A307035, A337377, A348990 [= a(A003961(n))], A349169 (odd numbers k such that A348993(k) = a(k)), A354365/A354366.

Array[#2/GCD[#1, #2] & @@ {#, Apply[Times, Map[If[#1 <= 2, 1, NextPrime[#1, -1]]^#2 & @@ # &, FactorInteger[#]]]} &, 120] (* Michael De Vlieger, Aug 27 2020 *)

a(n) = my (f=factor(n)); denominator(prod(i=1, #f~, my (p=f[i,1]); (p/if (p>2, precprime(p-1), 1))^f[i,2]))

a(n) = A064989(n) / gcd(n, A064989(n)).

a(n) = 1 iff n belongs to A025487.

"Primorial deflation" prefixed to the name by Antti Karttunen, Apr 29 2022

A325508 Product of primes indexed by the prime exponents of n!.

Original entry on oeis.org

1, 1, 2, 4, 10, 20, 42, 84, 204, 476, 798, 1596, 3828, 7656, 12276, 24180, 36660, 73320, 120840, 241680, 389424, 785680, 1294440, 2588880, 3848880, 7147920, 11264760, 15926040, 26057304, 52114608, 74421648, 148843296, 187159392, 340949280, 527531760, 926505360

Offset: 0

Gus Wiseman, May 08 2019

nonn

The prime indices of a(n) are the signature of n!, which is row n of A115627.

			We have 7! = 2^4 * 3^2 * 5^1 * 7^1, so a(7) = prime(4)*prime(2)*prime(1)*prime(1) = 84.
The sequence of terms together with their prime indices begins:
          1: {}
          1: {}
          2: {1}
          4: {1,1}
         10: {1,3}
         20: {1,1,3}
         42: {1,2,4}
         84: {1,1,2,4}
        204: {1,1,2,7}
        476: {1,1,4,7}
        798: {1,2,4,8}
       1596: {1,1,2,4,8}
       3828: {1,1,2,5,10}
       7656: {1,1,1,2,5,10}
      12276: {1,1,2,2,5,11}
      24180: {1,1,2,3,6,11}
      36660: {1,1,2,3,6,15}
      73320: {1,1,1,2,3,6,15}
     120840: {1,1,1,2,3,8,16}
     241680: {1,1,1,1,2,3,8,16}

Cf. A000142, A011371, A022559, A056171, A071626, A076934, A115627, A118914, A122111, A135291, A181819, A307035, A323014, A325272, A325273, A325276, A325509.

Table[Times@@Prime/@Last/@If[(n!)==1,{},FactorInteger[n!]],{n,0,30}]

a(n) = A181819(n!).
A001221(a(n)) = A071626(n).
A001222(a(n)) = A000720(n).
A056239(a(n)) = A022559(n).
A003963(a(n)) = A135291(n).
A061395(a(n)) = A011371(n).
A007814(a(n)) = A056171(n).
a(n) = A122111(A307035(n)). - Antti Karttunen, Nov 19 2019

cp's OEIS Frontend

A064989 Multiplicative with a(2^e) = 1 and a(p^e) = prevprime(p)^e for odd primes p.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

Python

Python

Formula

A108951 Primorial inflation of n: Fully multiplicative with a(p) = p# for prime p, where x# is the primorial A034386(x).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Sage

Formula

Extensions

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A329900 Primorial deflation of n: starting from x = n, repeatedly divide x by the largest primorial A002110(k) that divides it, until x is an odd number. Then a(n) = Product prime(k_i), for primorial indices k_1 >= k_2 >= ..., encountered in the process.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A325508 Product of primes indexed by the prime exponents of n!.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula