login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A137509 a(1)=2. For n >= 2, a(n) = the smallest integer > a(n-1) that has the same multiset of prime-factorization exponents as n has. 1
2, 3, 5, 9, 11, 14, 17, 27, 49, 51, 53, 63, 67, 69, 74, 81, 83, 92, 97, 98, 106, 111, 113, 135, 169, 177, 343, 356, 359, 366, 367, 3125, 3127, 3131, 3133, 3249, 3251, 3254, 3261, 3272, 3299, 3302, 3307, 3308, 3316, 3317, 3319, 3321, 3481, 3501, 3503, 3508 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

Starting the sequence at a(1)=1 instead leads to a(n) = n for every positive integer n.

LINKS

Table of n, a(n) for n=1..52.

EXAMPLE

12 = 2^2 * 3^1. So the multiset of exponents in the prime factorization of 12 is {1,2}. For a(12), we want the smallest integer > a(11)=53 of the form p^1 * q^2, where p and q are distinct primes. Checking: 54 = 2^1 *3^3, so 54 fails. 55 = 5^1*11^1. 56 = 2^3*7^1. 57 = 3^1*19^1. 58 = 2^1*29^1. 59=59^1. 60 = 2^2*3^1*5^1. 61 = 61^1. 62 = 2^1 *31^1. So 54 through 62 all fail. But 63 = 3^2 * 7^1, which has the same multiset of prime exponents, {1,2}, as 12 has. Therefore a(12) = 63.

MAPLE

pmset := proc(n) local e, a ; a := [] ; for e in ifactors(n)[2] do a := [op(a), e[2]] ; od: sort(a) ; end: A137509 := proc(n) option remember ; local nset, a ; if n = 1 then RETURN(2) ; fi ; nset := pmset(n) ; for a from A137509(n-1)+1 do if pmset(a) = nset then RETURN(a) ; fi ; od: end: seq(A137509(n), n=1..120) ; # R. J. Mathar, May 23 2008

CROSSREFS

Cf. A081761.

Sequence in context: A140329 A163292 A137518 * A014109 A102940 A157049

Adjacent sequences:  A137506 A137507 A137508 * A137510 A137511 A137512

KEYWORD

nonn

AUTHOR

Leroy Quet, Apr 23 2008

EXTENSIONS

More terms from R. J. Mathar, May 23 2008

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified November 17 03:06 EST 2019. Contains 329216 sequences. (Running on oeis4.)