%I #103 Nov 13 2024 13:27:52
%S 1,2,3,211,5,23,7,3331113965338635107,311,773,11,223,13,13367,1129,
%T 31636373,17,233,19,3318308475676071413,37,211,23,331319,773,3251,
%U 13367,227,29,547,31,241271,311,31397,1129,71129,37,373,313,3314192745739,41,379,43,22815088913,3411949,223,47,6161791591356884791277
%N Home primes: for n >= 2, a(n) = the prime that is finally reached when you start with n, concatenate its prime factors (A037276) and repeat until a prime is reached (a(n) = -1 if no prime is ever reached).
%C The initial 1 could have been omitted.
%C Probabilistic arguments give exactly zero for the chance that the sequence of integers starting at n contains no prime, the expected number of primes being given by a divergent sequence. - _J. H. Conway_
%C After over 100 iterations, a(49) is still composite - see A056938 for the latest information.
%C More terms:
%C a(50) to a(60) are 3517, 317, 2213, 53, 2333, 773, 37463, 1129, 229, 59, 35149;
%C a(61) to a(65) are 61, 31237, 337, 1272505013723, 1381321118321175157763339900357651;
%C a(66) to a(76) are 2311, 67, 3739, 33191, 257, 71, 1119179, 73, 379, 571, 333271.
%C This is different from A195264. Here 8 = 2^3 -> 222 -> ... -> 3331113965338635107 (a prime), whereas in A195264 8 = 2^3 -> 23 (a prime). - _N. J. A. Sloane_, Oct 12 2014
%D Jeffrey Heleen, Family Numbers: Mathemagical Black Holes, Recreational and Educational Computing, 5:5, pp. 6, 1990.
%D Jeffrey Heleen, Family numbers: Constructing Primes by Prime Factor Splicing, J. Recreational Math., Vol. 28 #2, 1996-97, pp. 116-119.
%H Christian N. K. Anderson, <a href="/A037274/a037274_1.txt">Table of known values of n, # of steps to reach a(n), and a(n)</a> or NA if a(n) has 30 digits or more. Also, the trajectory, with factors separated by a |, terminated by either "(end)" or "-> ?" if a(n) has 30 digits or more.
%H Patrick De Geest, <a href="http://www.worldofnumbers.com/topic1.htm">Home Primes < 100 and Beyond</a>
%H M. Herman and J. Schiffman, <a href="http://www.jstor.org/stable/10.5951/mathteacher.107.8.0606">Investigating home primes and their families</a>, Math. Teacher, 107 (No. 8, 2014), 606-614.
%H N. J. A. Sloane, <a href="/A195264/a195264.pdf">Confessions of a Sequence Addict (AofA2017)</a>, slides of invited talk given at AofA 2017, Jun 19 2017, Princeton. Mentions this sequence.
%H N. J. A. Sloane, Three (No, 8) Lovely Problems from the OEIS, Experimental Mathematics Seminar, Rutgers University, Oct 05 2017, <a href="https://vimeo.com/237029685">Part I</a>, <a href="https://vimeo.com/237030304">Part 2</a>, <a href="https://oeis.org/A290447/a290447_slides.pdf">Slides.</a> (Mentions this sequence)
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/HomePrime.html">Home Prime.</a>
%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Home_prime">Home prime</a>
%e 9 = 3*3 -> 33 = 3*11 -> 311, prime, so a(9) = 311.
%e The trajectory of 8 is more interesting:
%e 8 ->
%e 2 * 2 * 2 ->
%e 2 * 3 * 37 ->
%e 3 * 19 * 41 ->
%e 3 * 3 * 3 * 7 * 13 * 13 ->
%e 3 * 11123771 ->
%e 7 * 149 * 317 * 941 ->
%e 229 * 31219729 ->
%e 11 * 2084656339 ->
%e 3 * 347 * 911 * 118189 ->
%e 11 * 613 * 496501723 ->
%e 97 * 130517 * 917327 ->
%e 53 * 1832651281459 ->
%e 3 * 3 * 3 * 11 * 139 * 653 * 3863 * 5107
%e and 3331113965338635107 is prime, so a(8) = 3331113965338635107.
%p b:= n-> parse(cat(sort(map(i-> i[1]$i[2], ifactors(n)[2]))[])):
%p a:= n-> `if`(isprime(n) or n=1, n, a(b(n))):
%p seq(a(n), n=1..48); # _Alois P. Heinz_, Jan 09 2021
%t f[n_] := FromDigits@ Flatten[ IntegerDigits@ Table[ #[[1]], { #[[2]] }] & /@ FactorInteger@n, 2]; g[n_] := NestWhile[ f@# &, n, !PrimeQ@# &]; g[1] = 1; Array[g, 41] (* _Robert G. Wilson v_, Sep 22 2007 *)
%o (PARI) step(n)=my(f=factor(n),s="");for(i=1,#f~,for(j=1,f[i,2],s=Str(s,f[i,1]))); eval(s)
%o a(n)=if(n<4,return(n)); while(!isprime(n), n=step(n)); n \\ _Charles R Greathouse IV_, May 14 2015
%o (SageMath)
%o def digitLen(x,n):
%o r=0
%o while(x>0):
%o x//=n
%o r+=1
%o return r
%o def concatPf(x,n):
%o r=0
%o f=list(factor(x))
%o for c in range(len(f)):
%o for d in range(f[c][1]):
%o r*=(n**digitLen(f[c][0],n))
%o r+=f[c][0]
%o return r
%o def hp(x,n):
%o x1=concatPf(x,n)
%o while(x1!=x):
%o x=x1
%o x1=concatPf(x1,n)
%o return x
%o #example: prints the home prime of 8 in base 10
%o print(hp(8,10))
%o (Python)
%o from sympy import factorint, isprime
%o def f(n): return int("".join(str(p)*e for p, e in factorint(n).items()))
%o def a(n):
%o if n == 1: return 1
%o fn = n
%o while not isprime(fn): fn = f(fn)
%o return fn
%o print([a(n) for n in range(1, 40)]) # _Michael S. Branicky_, Jul 11 2022
%Y Cf. A006919, A037271, A037272, A037273, A037275, A037276, A037919-A037941, A048986, A056938.
%Y Cf. A195264 (use exponents instead of repeating primes).
%Y Cf. A084318 (use only one copy of each prime), A248713 (Fermi-Dirac analog: use unique representation of n>1 as a product of distinct terms of A050376).
%Y Cf. also A120716 and related sequences.
%K nonn,nice,base,changed
%O 1,2
%A _N. J. A. Sloane_, _Jeff Burch_
%E Corrected and extended by _Karl W. Heuer_, Sep 30 2003