cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

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

Views

Author

Antti Karttunen, Dec 22 2019

Keywords

Comments

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)

Examples

			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

Links

Crossrefs

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.

Programs

  • Mathematica
    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));
  • PARI
    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

Formula

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)

Extensions

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

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

Views

Author

Paul Boddington, Jul 21 2005

Keywords

Comments

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

Examples

			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).

Links

Crossrefs

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.

Programs

  • Mathematica
    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 *)
  • PARI
    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
  • Python
    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
  • Sage
    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

Formula

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

Extensions

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

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

Views

Author

Antti Karttunen, Dec 22 2019

Keywords

Comments

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

Links

Crossrefs

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.

Programs

  • Mathematica
    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 *)
  • PARI
    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); };
  • PARI
    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)));

Formula

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.
Showing 1-3 of 3 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.