A276662 Home primes (A037274) for the odd composite numbers (A071904) (a(n) = -1 if no prime is ever reached).

311, 1129, 37, 773, 313, 311, 1129, 313, 3119014487, 31079, 317, 773, 1129, 3110647, 3103819425079, 310397, 5113, 31079, 3109, 3137, 310361, 31259, 331, 36389, 191176757654383, 31063, 337, 523, 324941, 31393, 127139, 33769, 31034567124791, 32369, 719, 5623, 347, 3371, 131777, 349, 31039, 34412909

Offset

1,1

Comments

a(n) = A037274(A071904).

Links

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

Example

The first entry is from 9 = 3*3. 33 = 3*11, and 311 is prime.
A longer 10 step progression is a(9) from 45. Specifically, 45=3*15 concatenating to 315=3*105 concatenating to 3105=3*1035 concatenating to 31035=3*10345 concatenating to 310345=5*62069 concatenating to 562069=41*13709 concatenating to 4113709=19*216511 concatenating to 19216511=17*1130383 concatenating to 171130383 = 3*57043461 concatenating to 357043461=3*119014487 concatenating to 3119014487 which is prime. a(9) then is 3119014487.

Mathematica

Map[NestWhile[Function[n, FromDigits@ Flatten@ IntegerDigits@ {#, n/#} &[FactorInteger[n][[1, 1]]]], #, ! PrimeQ@ # &] &, Select[Range[9, 157, 2], CompositeQ]] (* Michael De Vlieger, Sep 13 2016 *)

Prog

(PARI) genit(iend)={i5=9; while(i5<=iend, n=i5; while(isprime(n), n+=2); i5=n; endless=0; while(endless<99999, dun=0; z=divisors(n);
a=z[2]; b=n/a; k=length(digits(b)); q=a*10^k+b; if(isprime(q), dun=1; break); endless+=1; n=q); if(dun>0, print1(q, ", ")); i5+=2); }

Crossrefs

Cf. A037274, A071904.

Sequence in context: A067380 A046495 A061329 * A082862 A274236 A213814

Adjacent sequences: A276659 A276660 A276661 * A276663 A276664 A276665

Keyword

nonn,base

Author

Bill McEachen, Sep 11 2016

Extensions

Edited by N. J. A. Sloane, Oct 02 2016

Status

approved