login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A037276 Start with 1; for n>1, replace n with the concatenation of its prime factors in increasing order. 35
1, 2, 3, 22, 5, 23, 7, 222, 33, 25, 11, 223, 13, 27, 35, 2222, 17, 233, 19, 225, 37, 211, 23, 2223, 55, 213, 333, 227, 29, 235, 31, 22222, 311, 217, 57, 2233, 37, 219, 313, 2225, 41, 237, 43, 2211, 335, 223, 47, 22223, 77, 255, 317, 2213, 53, 2333 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
LINKS
N. J. A. Sloane, Table of n, a(n) for n = 1..20000 [First 10000 terms from T. D. Noe]
Patrick De Geest, Home Primes
N. J. A. Sloane, Confessions of a Sequence Addict (AofA2017), slides of invited talk given at AofA 2017, Jun 19 2017, Princeton. Mentions this sequence.
N. J. A. Sloane, Three (No, 8) Lovely Problems from the OEIS, Experimental Mathematics Seminar, Rutgers University, Oct 05 2017, Part I, Part 2, Slides. (Mentions this sequence)
Eric Weisstein's World of Mathematics, Home Prime
EXAMPLE
If n = 2^3*5^5*11^2 = 3025000, a(n) = 222555551111 (n=2*2*2*5*5*5*5*5*11*11, then remove the multiplication signs).
MAPLE
# This is for n>1
read("transforms") ;
A037276 := proc(n)
local L, p ;
L := [] ;
for p in ifactors(n)[2] do
L := [op(L), seq(op(1, p), i=1..op(2, p))] ;
end do:
digcatL(L) ;
end proc: # R. J. Mathar, Oct 29 2012
MATHEMATICA
co[n_, k_] := Nest[Flatten[IntegerDigits[{#, n}]] &, n, k - 1]; Table[FromDigits[Flatten[IntegerDigits[co @@@ FactorInteger[n]]]], {n, 54}] (* Jayanta Basu, Jul 04 2013 *)
FromDigits@ Flatten@ IntegerDigits[Table[#1, {#2}] & @@@ FactorInteger@ #] & /@ Range@ 54 (* Michael De Vlieger, Jul 14 2015 *)
PROG
(PARI) a(n)={ n<4 & return(n); for(i=1, #n=factor(n)~, n[1, i]=concat(vector(n[2, i], j, Str(n[1, i])))); eval(concat(n[1, ]))} \\ M. F. Hasler, Jun 19 2011
(Haskell)
a037276 = read . concatMap show . a027746_row
-- Reinhard Zumkeller, Apr 03 2012
(Python)
from sympy import factorint
def a(n):
f=factorint(n)
l=sorted(f)
return 1 if n==1 else int("".join(str(i)*f[i] for i in l))
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jun 23 2017
CROSSREFS
Cf. A037274, A048985, A067599, A080670, A084796. Different from A073646.
Cf. also A027746, A289660 (a(n)-n).
Sequence in context: A288532 A073647 A073646 * A068190 A084796 A084797
KEYWORD
nonn,easy,base
AUTHOR
STATUS
approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 26 21:53 EDT 2024. Contains 372004 sequences. (Running on oeis4.)