A195264 Iterate x -> A080670(x) (replace x with the concatenation of the primes and exponents in its prime factorization) starting at n until reach 1 or a prime (which is then the value of a(n)); or a(n) = -1 if a prime is never reached.
1, 2, 3, 211, 5, 23, 7, 23, 2213, 2213, 11, 223, 13, 311, 1129, 233, 17, 17137, 19
Offset: 1
Examples
4 = 2^2 -> 22 =2*11 -> 211, prime, so a(4) = 211. 9 = 3^2 -> 32 = 2^5 -> 25 = 5^2 -> 52 = 2^2*13 -> 2213, prime, so a(9)=2213.
Links
- R. J. Mathar, Table of n, a(n) for n = 1..19
- J. H. Conway, Five $1000 Problems, Oct 10 2014 (recorded by Bill Cheswick).
- Alonso del Arte and Sean A. Irvine, Table of n, a(n) for n = 1..1000
- Hans Havermann, Table of n, a(n) for n = 1..10000 (includes links to lengthy (>40) and unknown-outcome evolutions, and a list of unfactored composites in the unknowns' last step)
- Hans Havermann, Table of n, a(n) for n = 1..10000 (includes links to lengthy (>40) and unknown-outcome evolutions, and a list of unfactored composites in the unknowns' last step) [Cached copy, pdf version as of May 08 2018, with permission]
- Hans Havermann, 13532385396179 precursors
- Hans Havermann, 13532385396179 precursors [Cached copy, pdf version as of May 08 2018, with permission]
- OEIS50 DIMACS Conference on Challenges of Identifying Integer Sequences, Problem Session 2, Oct 10 2014, J. H. Conway, Five $1000 Problems (starting at about 06.44). This sequence is mentioned in the fifth problem, starting at around 19:30.
- Tony Padilla and Brady Haran, 13532385396179, Numberphile video, 2017
- 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)
- Doron Zeilberger, Videos of Talks Delivered in SLOANE75/OEIS50 DIMACS Conference on Challenges of Identifying Integer Sequences (see Problem Session 2, Oct 10 2014)
Crossrefs
Programs
-
Mathematica
f[1] := 1; f[n_] := Block[{p = Flatten[FactorInteger[n]]}, k = Length[p]; While[k > 0, If[p[[k]] == 1, p = Delete[p, k]]; k--]; FromDigits[Flatten[IntegerDigits[p]]]]; Table[FixedPoint[f, n], {n, 19}] (* Alonso del Arte, based on the program for A080670, Sep 14 2011 *) fn[n_] := FromDigits[Flatten[IntegerDigits[DeleteCases[Flatten[ FactorInteger[n]], 1]]]]; Table[NestWhile[fn, n, # != 1 && ! PrimeQ[#] &], {n, 19}] (* Robert Price, Mar 15 2020 *)
-
PARI
a(n)={n>1 && while(!ispseudoprime(n), n=A080670(n));n} \\ M. F. Hasler, Oct 12 2014
Comments