A330683 -id:A330683 - cp's OEIS frontend

A329906 a(0) = 1; a(1) = 2; after which a(2n) = A329898(a(n)), a(2n+1) = A330683(a(n)).

1, 2, 3, 4, 5, 11, 6, 9, 7, 23, 15, 38, 8, 20, 13, 22, 10, 44, 30, 110, 19, 69, 49, 128, 12, 41, 27, 72, 17, 43, 29, 54, 14, 79, 56, 272, 37, 181, 136, 482, 26, 118, 86, 307, 61, 208, 156, 424, 16, 73, 52, 190, 34, 123, 89, 242, 24, 77, 55, 147, 36, 93, 66, 114, 18, 131, 97, 596, 68, 416, 323, 1448, 48, 286, 218, 990, 164, 711

Offset: 0

Antti Karttunen, Dec 24 2019

nonn

Note the indexing: domain begins from zero, but the range does not include it.

			This irregular table can be represented as a binary tree. Each child to the left is obtained by applying A329898 the parent, and each child to the right is obtained by applying A330683 to the parent:
                                      1
                                      |
                   ...................2...................
                  3                                       4
        5......../ \........11                  6......../ \........9
       / \                 / \                 / \                 / \
      /   \               /   \               /   \               /   \
     /     \             /     \             /     \             /     \
    7       23         15       38          8       20         13       22
  10 44   30  110    19  69    49 128     12 41   27  72     17  43   29  54
etc.

Index entries for sequences that are permutations of the natural numbers

Cf. A329905 (inverse permutation).

Cf. A163511, A181815, A329897, A329898, A329901, A330683.

A329906(n) = if(n<2,1+n,if(!(n%2),A329898(A329906(n/2)),A330683(A329906(n\2))));

a(0) = 1; a(1) = 2; after which a(2n) = A329898(a(n)), a(2n+1) = A330683(a(n)).

a(n) = A329901(A163511(n)).

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

Original entry on oeis.org

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

A181815 a(n) = largest integer such that, when any of its divisors divides A025487(n), the quotient is a member of A025487.

Original entry on oeis.org

1, 2, 4, 3, 8, 6, 16, 12, 5, 32, 9, 24, 10, 64, 18, 48, 20, 128, 36, 15, 96, 7, 27, 40, 256, 72, 30, 192, 14, 54, 80, 512, 144, 60, 384, 28, 108, 25, 160, 1024, 45, 288, 21, 81, 120, 768, 56, 216, 50, 320, 2048, 90, 576, 11, 42, 162, 240, 1536, 112, 432, 100, 640, 4096, 180, 1152

Offset: 1

Matthew Vandermast, Nov 30 2010

A permutation of the natural numbers.
The number of divisors of a(n) equals the number of ordered factorizations of A025487(n) as A025487(j)*A025487(k). Cf. A182762.
From Antti Karttunen, Dec 28 2019: (Start)
Rearranges terms of A108951 into ascending order, as A108951(a(n)) = A025487(n).
The scatter plot looks like a curtain of fractal spray, which is typical for many of the so-called "entanglement permutations". Indeed, according to the terminology I use in my 2016-2017 paper, this sequence is obtained by entangling the complementary pair (A329898, A330683) with the complementary pair (A005843, A003961), which gives the following implicit recurrence: a(A329898(n)) = 2*a(n) and a(A330683(n)) = A003961(a(n)). An explicit form is given in the formula section.
(End)

			For any divisor d of 9 (d = 1, 3, 9), 36/d (36, 12, 4) is a member of A025487. 9 is the largest number with this relationship to 36; therefore, since 36 = A025487(11), a(11) = 9.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Antti Karttunen, Entanglement Permutations, 2016-2017.
Index entries for sequences that are permutations of the natural numbers

If this sequence is considered the "primorial deflation" of A025487(n) (see first formula), the primorial inflation of n is A108951(n), and the primorial inflation of A025487(n) is A181817(n).
A181820(n) is another mapping from the members of A025487 to the positive integers.
Cf. A003961, A051282, A108951, A163511, A304886, A319626, A329897, A329898, A329900, A329901 (inverse), A329904, A329905, A329907, A330682 (reduced modulo 2), A330683.

(* First, load the program at A025487, then: *)
Map[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[g, ReplacePart[Table[0, {PrimePi[g[[-1, 1]]]}], #] &@ Map[PrimePi@ First@ # -> Last@ # &, g]]@ FactorInteger@ #], # > 0 &]}}, And[FreeQ[#[[-1, -1]], 0], Length[#[[-1, -1]] ] != 0] &][[All, 1]] ] &, Union@ Flatten@ f@ 6] (* Michael De Vlieger, Dec 28 2019 *)

A181815(n) = A329900(A025487(n)); \\ Antti Karttunen, Dec 24 2019

If A025487(n) is considered in its form as Product A002110(i)^e(i), then a(n) = Product p(i)^e(i). If A025487(n) is instead considered as Product p(i)^e(i), then a(n) = Product (p(i)/A008578(i))^e(i).
a(n) = A122111(A181820(n)). - Matthew Vandermast, May 21 2012
From Antti Karttunen, Dec 24-29 2019: (Start)
a(n) = Product_{i=1..A051282(n)} A000040(A304886(i)).
a(n) = A329900(A025487(n)) = A319626(A025487(n)).
a(n) = A163511(A329905(n)).
For n > 1, if A330682(n) = 1, then a(n) = A003961(a(A329904(n))), otherwise a(n) = 2*a(A329904(n)).
A252464(a(n)) = A329907(n).
A330690(a(n)) = A050378(n).
a(A306802(n)) = A329902(n).
(End)

A329898 a(n) is the position of 2*A025487(n) in A025487.

Original entry on oeis.org

2, 3, 5, 6, 7, 8, 10, 12, 13, 14, 15, 16, 17, 18, 19, 21, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 39, 40, 42, 45, 46, 47, 48, 49, 50, 51, 52, 53, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 70, 71, 74, 75, 76, 78, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 94, 95, 96, 97, 98, 99, 100

Offset: 1

Antti Karttunen, Dec 24 2019

nonn

Numbers k for which A007814(A025487(k)) > A007949(A025487(k)), i.e., numbers k for which the 2-adic valuation of A025487(k) is larger than its 3-adic valuation.

Numbers k for which A181815(k) is even.

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

Cf. A329897 (complement), A330683 (and its permutation).
Cf. A007814, A007949, A025487, A329904 (a left inverse), A329906.
Positions of even terms in A181815, zeros in A330682.

(* First, load the function f at A025487, then: *)
With[{s = Union@ Flatten@ f@ 6}, Map[If[2 # > Max@ s, Nothing, FirstPosition[s, 2 #][[1]] ] &, s]] (* Michael De Vlieger, Jan 11 2020 *)

upto_e = 64; \\ 64 -> 43608 terms.
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
A329898list(e) = { my(lista = List([1, 2]), i=2, u = 2^e, t, v025487); 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++); v025487 = vecsort(Vec(lista)); lista = List([]); for(i=1,oo,if(!(t=vecsearch(v025487,2*(v025487[i]))),return(Vec(lista)), listput(lista,t))); };
v329898 = A329898list(upto_e);
A329898(n) = v329898[n];

For all n >= 1, A329904(a(n)) = n.

A329897 Numbers k for which the 2-adic and 3-adic valuations of A025487(k) are equal, where A025487(k) is the k-th number which is a product of primorial numbers.

Original entry on oeis.org

1, 4, 9, 11, 20, 22, 23, 38, 41, 43, 44, 54, 69, 72, 73, 77, 79, 93, 110, 114, 118, 123, 124, 128, 129, 131, 147, 154, 181, 186, 190, 191, 199, 201, 208, 209, 212, 232, 242, 245, 246, 272, 279, 286, 294, 299, 300, 307, 310, 312, 321, 324, 327, 345, 359, 371, 374, 376, 416, 424, 425, 430, 434, 442, 446, 451, 454, 466, 469

Offset: 1

Antti Karttunen, Dec 24 2019

nonn

Numbers k for which A007814(A025487(k)) = A007949(A025487(k)).

Numbers k for which A181815(k) is odd.

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

Cf. A007814, A007949, A025487, A329898 (complement), A330682 (characteristic function).
Sequence A330683 sorted into ascending order.
Positions of odd terms in A181815.

s = {1}; k = 1; Do[If[GreaterEqual @@ (f = FactorInteger[n])[[;; , 2]] && PrimePi[f[[-1, 1]]] == Length[f], k++; If[Equal @@ IntegerExponent[n, {2, 3}], AppendTo[s, k]]], {n, 2, 10^5}]; s (* Amiram Eldar, Jul 28 2023 *)

A330681 a(n) = A283980(A025487(n)).

Original entry on oeis.org

1, 6, 36, 30, 216, 180, 1296, 1080, 210, 7776, 900, 6480, 1260, 46656, 5400, 38880, 7560, 279936, 32400, 6300, 233280, 2310, 27000, 45360, 1679616, 194400, 37800, 1399680, 13860, 162000, 272160, 10077696, 1166400, 226800, 8398080, 83160, 972000, 44100, 1632960, 60466176, 189000, 6998400, 69300, 810000, 1360800

Offset: 1

Antti Karttunen, Dec 26 2019

nonn

After 1, contains only the least representatives of such prime signatures where the maximal exponent occurs more than once. However, here the terms are not in ascending order.

			For example, 180 = 2^2 * 3^2 * 5^1 is present, because the maximal exponent in its prime factorization is 2, which occurs as an exponent of both 2 and 3, and because 180 is the minimal representative of the prime signature (2,2,1), as A046523(180) = 180.

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

Cf. A025487, A046523, A283980, A329899, A330683.

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(64);
A025487(n) = v025487[n];
A330681(n) = A283980(A025487(n));

a(n) = A283980(A025487(n)).
A046523(a(n)) = a(n).
A085089(a(n)) = A101296(a(n)) = A330683(n).

A341352 Inverse permutation to A341351.

Original entry on oeis.org

1, 2, 4, 9, 3, 22, 54, 6, 114, 246, 13, 488, 11, 5, 948, 1809, 29, 20, 3327, 66, 6020, 10624, 8, 18246, 38, 140, 30726, 43, 290, 51148, 84074, 17, 93, 135598, 570, 216398, 340886, 15, 72, 529051, 7, 814237, 186, 1090, 1240172, 147, 2057, 376, 1874464, 36, 2817289

Offset: 1

Antti Karttunen and Michael De Vlieger, Feb 14 2021

nonn

Index entries for sequences that are permutations of the natural numbers

Cf. A064216, A064989, A329898, A329901, A329906, A330683.

Cf. A341351 (inverse).

a(n) = A329901(A064216(n)).

cp's OEIS Frontend

A329906 a(0) = 1; a(1) = 2; after which a(2n) = A329898(a(n)), a(2n+1) = A330683(a(n)).

Views

Author

Keywords

Comments

Examples

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

A181815 a(n) = largest integer such that, when any of its divisors divides A025487(n), the quotient is a member of A025487.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A329898 a(n) is the position of 2*A025487(n) in A025487.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A329897 Numbers k for which the 2-adic and 3-adic valuations of A025487(k) are equal, where A025487(k) is the k-th number which is a product of primorial numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

A330681 a(n) = A283980(A025487(n)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Formula

A341352 Inverse permutation to A341351.

Views

Author

Keywords

Links