A330687 -id:A330687 - cp's OEIS frontend

A025487 Least integer of each prime signature A124832; also products of primorial numbers A002110.

1, 2, 4, 6, 8, 12, 16, 24, 30, 32, 36, 48, 60, 64, 72, 96, 120, 128, 144, 180, 192, 210, 216, 240, 256, 288, 360, 384, 420, 432, 480, 512, 576, 720, 768, 840, 864, 900, 960, 1024, 1080, 1152, 1260, 1296, 1440, 1536, 1680, 1728, 1800, 1920, 2048, 2160, 2304, 2310

Offset: 1

David W. Wilson

All numbers of the form 2^k1*3^k2*...*p_n^k_n, where k1 >= k2 >= ... >= k_n, sorted.
A111059 is a subsequence. - Reinhard Zumkeller, Jul 05 2010
Choie et al. (2007) call these "Hardy-Ramanujan integers". - Jean-François Alcover, Aug 14 2014
The exponents k1, k2, ... can be read off Abramowitz & Stegun p. 831, column labeled "pi".
For all such sequences b for which it holds that b(n) = b(A046523(n)), the sequence which gives the indices of records in b is a subsequence of this sequence. For example, A002182 which gives the indices of records for A000005, A002110 which gives them for A001221 and A000079 which gives them for A001222. - Antti Karttunen, Jan 18 2019
The prime signature corresponding to a(n) is given in row n of A124832. - M. F. Hasler, Jul 17 2019

			The first few terms are 1, 2, 2^2, 2*3, 2^3, 2^2*3, 2^4, 2^3*3, 2*3*5, ...

Franklin T. Adams-Watters, Table of n, a(n) for n = 1..10001 (first 291 terms from Will Nicholes)
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972.
Kevin Broughan, Equivalents of the Riemann Hypothesis, Vol. 1: Arithmetic Equivalents, Cambridge University Press, 2017. See section 8.2, "Hardy-Ramanujan Numbers".
YoungJu Choie, Nicolas Lichiardopol, Pieter Moree and Patrick Solé, On Robin's criterion for the Riemann hypothesis, Journal de théorie des nombres de Bordeaux, Vol. 19, No. 2 (2007), pp. 357-372. See section 5, p. 367.
Asaf Cohen Antonir and Asaf Shapira, An Elementary Proof of a Theorem of Hardy and Ramanujan (2022). arXiv:2207.09410 [math.NT]
Michael De Vlieger, Relations of A025487 to A002110, A002182, and A002201.
Steven R. Finch, Errata and Addenda to Mathematical Constants, arXiv:2001.00578 [math.HO], 2020, pp. 9-10.
G. H. Hardy and S. Ramanujan, Asymptotic formulae for the distribution of integers of various types, Proc. London Math. Soc, Ser. 2, Vol. 16 (1917), pp. 112-132. Also published in the book Collected Papers of Srinivasa Ramanujan, Chelsea, 1962, pages 245-261.
Jeffery Kline, On the eigenstructure of sparse matrices related to the prime number theorem, Linear Algebra and its Applications (2020) Vol. 584, 409-430.
L. B. Richmond, Asymptotic results for partitions (I) and the distribution of certain integers, Journal of Number Theory, Vol. 8, No. 4 (1976), pp. 372-389. See page 388.

Subsequence of A055932, A191743, and A324583.
Cf. A025488, A051282, A036041, A051466, A061394, A124832, A161360, A166469, A181815, A181817, A283980, A306802, A322584, A322585 (characteristic function), A329897, A329898, A329899, A329900, A329904, A330683.
Cf. A085089, A101296 (left inverses).
Equals range of values taken by A046523.
Cf. A178799 (first differences), A247451 (squarefree kernel), A146288 (number of divisors).
Subsequences of this sequence include: A000079, A000142, A000400, A001013, A001813, A002110, A002182, A005179, A006939, A025527, A056836, A061742, A064350, A066120, A087980, A097212, A097213, A111059, A119840, A119845, A126098, A129912, A140999, A166338, A166470, A166472, A166473, A166475, A167448, A168262, A168263, A168264, A179215, A181555, A181804, A181806, A181809, A181818, A181822, A181823, A181824, A181825, A181826, A181827, A182763, A182862, A182863, A212170, A220264, A220423, A250269, A250270, A260633, A266047, A284456, A300357, A304938, A329894, A330687; also A037019 and A330681 (when sorted), possibly also A289132.
Rearrangements of this sequence include A036035, A059901, A063008, A077569, A085988, A086141, A087443, A108951, A181821, A181822, A322827, A329886, A329887.
Cf. also array A124832 (row n = prime signature of a(n)) and A304886, A307056.

import Data.Set (singleton, fromList, deleteFindMin, union)
a025487 n = a025487_list !! (n-1)
a025487_list = 1 : h [b] (singleton b) bs where
   (_ : b : bs) = a002110_list
   h cs s xs'@(x:xs)
     | m <= x    = m : h (m:cs) (s' `union` fromList (map (* m) cs)) xs'
     | otherwise = x : h (x:cs) (s  `union` fromList (map (* x) (x:cs))) xs
     where (m, s') = deleteFindMin s
-- Reinhard Zumkeller, Apr 06 2013

isA025487 := proc(n)
    local pset,omega ;
    pset := sort(convert(numtheory[factorset](n),list)) ;
    omega := nops(pset) ;
    if op(-1,pset) <> ithprime(omega) then
        return false;
    end if;
    for i from 1 to omega-1 do
        if padic[ordp](n,ithprime(i)) < padic[ordp](n,ithprime(i+1)) then
            return false;
        end if;
    end do:
    true ;
end proc:
A025487 := proc(n)
    option remember ;
    local a;
    if n = 1 then
        1 ;
    else
        for a from procname(n-1)+1 do
            if isA025487(a) then
                return a;
            end if;
        end do:
    end if;
end proc:
seq(A025487(n),n=1..100) ; # R. J. Mathar, May 25 2017

PrimeExponents[n_] := Last /@ FactorInteger[n]; lpe = {}; ln = {1}; Do[pe = Sort@PrimeExponents@n; If[ FreeQ[lpe, pe], AppendTo[lpe, pe]; AppendTo[ln, n]], {n, 2, 2350}]; ln (* Robert G. Wilson v, Aug 14 2004 *)
(* Second program: generate all terms m <= A002110(n): *)
f[n_] := {{1}}~Join~
  Block[{lim = Product[Prime@ i, {i, n}],
   ww = NestList[Append[#, 1] &, {1}, n - 1], dec},
   dec[x_] := Apply[Times, MapIndexed[Prime[First@ #2]^#1 &, x]];
   Map[Block[{w = #, k = 1},
      Sort@ Prepend[If[Length@ # == 0, #, #[[1]]],
        Product[Prime@ i, {i, Length@ w}] ] &@ Reap[
         Do[
          If[# < lim,
             Sow[#]; k = 1,
             If[k >= Length@ w, Break[], k++]] &@ dec@ Set[w,
             If[k == 1,
               MapAt[# + 1 &, w, k],
               PadLeft[#, Length@ w, First@ #] &@
                 Drop[MapAt[# + Boole[i > 1] &, w, k], k - 1] ]],
           {i, Infinity}] ][[-1]]
] &, ww]]; Sort[Join @@ f@ 13] (* Michael De Vlieger, May 19 2018 *)

isA025487(n)=my(k=valuation(n,2),t);n>>=k;forprime(p=3,default(primelimit),t=valuation(n,p);if(t>k,return(0),k=t);if(k,n/=p^k,return(n==1))) \\ Charles R Greathouse IV, Jun 10 2011

factfollow(n)={local(fm, np, n2);
  fm=factor(n); np=matsize(fm)[1];
  if(np==0,return([2]));
  n2=n*nextprime(fm[np,1]+1);
  if(np==1||fm[np,2]Franklin T. Adams-Watters, Dec 01 2011 */

is(n) = {if(n==1, return(1)); my(f = factor(n));  f[#f~, 1] == prime(#f~) && vecsort(f[, 2],,4) == f[, 2]} \\ David A. Corneth, Feb 14 2019

upto(Nmax)=vecsort(concat(vector(logint(Nmax,2),n,select(t->t<=Nmax,if(n>1,[factorback(primes(#p),Vecrev(p)) || p<-partitions(n)],[1,2]))))) \\ M. F. Hasler, Jul 17 2019

\\ For fast generation of large number of terms, use this program:
A283980(n) = {my(f=factor(n)); prod(i=1, #f~, my(p=f[i, 1], e=f[i, 2]); if(p==2, 6, nextprime(p+1))^e)}; \\ From A283980
A025487list(e) = { my(lista = List([1, 2]), i=2, u = 2^e, t); while(lista[i] != u, if(2*lista[i] <= u, listput(lista,2*lista[i]); t = A283980(lista[i]); if(t <= u, listput(lista,t))); i++); vecsort(Vec(lista)); }; \\ Returns a list of terms up to the term 2^e.
v025487 = A025487list(101);
A025487(n) = v025487[n];
for(n=1,#v025487,print1(A025487(n), ", ")); \\ Antti Karttunen, Dec 24 2019

def sharp_primorial(n): return sloane.A002110(prime_pi(n))
N = 2310
nmax = 2^floor(log(N,2))
sorted([j for j in (prod(sharp_primorial(t[0])^t[1] for k, t in enumerate(factor(n))) for n in (1..nmax)) if j <= N])
# Giuseppe Coppoletta, Jan 26 2015

What can be said about the asymptotic behavior of this sequence? - Franklin T. Adams-Watters, Jan 06 2010
Hardy & Ramanujan prove that there are exp((2 Pi + o(1))/sqrt(3) * sqrt(log x/log log x)) members of this sequence up to x. - Charles R Greathouse IV, Dec 05 2012
From Antti Karttunen, Jan 18 & Dec 24 2019: (Start)
A085089(a(n)) = n.
A101296(a(n)) = n [which is the first occurrence of n in A101296, and thus also a record.]
A001221(a(n)) = A061395(a(n)) = A061394(n).
A007814(a(n)) = A051903(a(n)) = A051282(n).
a(A101296(n)) = A046523(n).
a(A306802(n)) = A002182(n).
a(n) = A108951(A181815(n)) = A329900(A181817(n)).
If A181815(n) is odd, a(n) = A283980(a(A329904(n))), otherwise a(n) = 2*a(A329904(n)).
(End)
Sum_{n>=1} 1/a(n) = Product_{n>=1} 1/(1 - 1/A002110(n)) = A161360. - Amiram Eldar, Oct 20 2020

Offset corrected by Matthew Vandermast, Oct 19 2008

Minor correction by Charles R Greathouse IV, Sep 03 2010

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.

A050377 Number of ways to factor n into "Fermi-Dirac primes" (members of A050376).

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 1, 1, 4, 1, 2, 1, 2, 1, 1, 1, 2, 2, 1, 2, 2, 1, 1, 1, 4, 1, 1, 1, 4, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 4, 2, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1, 2, 6, 1, 1, 1, 2, 1, 1, 1, 4, 1, 1, 2, 2, 1, 1, 1, 4, 4, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 4, 1, 2, 2, 4, 1, 1, 1, 2, 1

Offset: 1

Christian G. Bower, Nov 15 1999

a(n) depends only on prime signature of n (cf. A025487). So a(24) = a(375) since 24 = 2^3 * 3 and 375 = 3 * 5^3 both have prime signature (3,1).

Antti Karttunen, Table of n, a(n) for n = 1..100000 (the first 10000 terms from Reinhard Zumkeller)
Index entries for sequences computed from exponents in factorization of n.

Cf. A000123, A001055, A002110, A005117, A018819, A050376-A050380, A025487.

Cf. A108951, A330687 (positions of records), A330688 (record values), A330689, A330690, A382295.

A018819:= proc(n) option remember;
  if n::odd then procname(n-1)
  else procname(n-1) + procname(n/2)
  fi
end proc:
A018819(0):= 1:
f:= n -> mul(A018819(s[2]),s=ifactors(n)[2]):
seq(f(n),n=1..100); # Robert Israel, Jan 14 2016

b[0] = 1; b[n_] := b[n] = b[n - 1] + If[EvenQ[n], b[n/2], 0];
a[n_] := Times @@ (b /@ FactorInteger[n][[All, 2]]);
Array[a, 102] (* Jean-François Alcover, Jan 27 2018 *)

A018819(n) = if( n<1, n==0, if( n%2, A018819(n-1), A018819(n/2)+A018819(n-1))); \\ From A018819
A050377(n) = factorback(apply(e -> A018819(e), factor(n)[, 2])); \\ Antti Karttunen, Dec 28 2019

Dirichlet g.f.: Product_{n in A050376} (1/(1-1/n^s)).
a(p^k) = A000123([k/2]) for all primes p.
a(A002110(n)) = 1.
Multiplicative with a(p^e) = A018819(e). - Christian G. Bower and David W. Wilson, May 22 2005
a(n) = Sum{a(d): d^2 divides n}, a(1) = 1. - Reinhard Zumkeller, Jul 12 2007
a(A108951(n)) = A330690(n). - Antti Karttunen, Dec 28 2019
a(n) = 1 for all squarefree values of n (A005117). - Eric Fox, Feb 03 2020
G.f.: Sum_{k>=1} a(k) * x^(k^2) / (1 - x^(k^2)). - Ilya Gutkovskiy, Nov 25 2020
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p prime} f(1/p) = 1.7876368001694456669... (A382295), where f(x) = (1-x) / Product_{k>=0} (1 - x^(2^k)). - Amiram Eldar, Oct 03 2023

More terms from Antti Karttunen, Dec 28 2019

cp's OEIS Frontend

A330689 Primorial deflation of A330687 (record positions in A050377): a(n) is the unique integer x such that A108951(x) = A330687(n).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A025487 Least integer of each prime signature A124832; also products of primorial numbers A002110.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

PARI

PARI

Sage

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

A050377 Number of ways to factor n into "Fermi-Dirac primes" (members of A050376).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A330688 Record values in A050377, number of ways to factor n into "Fermi-Dirac primes" (A050376).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A330684 Square roots of the positions of records in A050377, number of ways to factor n into "Fermi-Dirac primes" (A050376).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

A382295 Decimal expansion of the asymptotic mean of the number of ways to factor k into "Fermi-Dirac primes" when k runs over the positive integers.

Views

Author

Keywords

Examples

Crossrefs