|
|
A237584
|
|
a(n) = number of iterations needed to stabilize the map m -> f(m) where f(m) is the concatenation in descending order of the prime distinct divisors of the numbers m = 2n+1, or = -1 if there is no such a(n).
|
|
1
|
|
|
0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 1, 0, 0, 1, 2, 0, 2, 0, 0, 1, 0, 1, 1, 0, 6, 1, 0, 0, 1, 2, 0, 1, 0, 0, 1, 3, 0, 1, 0, 3, 1, 0, 1, 1, 3, 0, 1, 0, 0, 7, 0, 0, 1, 0, 5, 2, 2, 1, 3, 1, 0, 1, 0, 1, 1, 0, 0, 3, 13, 9, 1, 0, 0, 1, 8, 0, 5, 3, 0, 1, 0, 1, 1, 0, 2
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,17
|
|
COMMENTS
|
Let f denote the map that replaces n by the concatenation of its prime divisors, written in decreasing order, each divisor being written in base 10 in the normal way. Then a(n) is the number of iterations needed to reach either a prime number or stabilize until a last composite number of a repetitive cycle (see the second example) when starting at 2n+1 and iterating f.
It seems that a(222) is unknown.
|
|
LINKS
|
|
|
EXAMPLE
|
a(57) = 5 because 2*57+1 = 115 -> 235 -> 475 -> 195 -> 1353 -> 41113 with 5 iterations, where the prime divisors of 115 in descending order are {23, 5} => 235, the prime divisors of 235 in descending order are {47, 5} => 475, the prime divisors of 475 in descending order are {19, 5} => 195, the prime divisors of 195 in descending order are {13,5,3} => 1353, and the prime divisors of 1353 are {41,11,3} => 41113 is prime and the last member of the trajectory.
a(1561)= 3 because 2*1561+1 = 3123 -> 3473 -> 15123 -> 713 with 3 iterations because 713 -> 3123 is already in the trajectory.
|
|
MAPLE
|
with(numtheory):for m from 1 to 500 do:n:=2*m+1:lst:={n}:nn:=n:for it from 1 to 10^4 while(type(nn, prime)=false) do:x:=factorset(nn):n1:=nops(x): s:=sum('x[n1-i+1]*10^(sum('length(x[n1-j+1])', 'j'=i+1..n1))', 'i'=1..n1):lst:=lst union {s}:nn:=s:od:n2:=nops(lst)-1: printf(`%d, `, n2):od:
|
|
PROG
|
(PARI) cmap(n) = {v = []; nb = 0; ok = 0; if (isprime(n), return (nb)); while (!ok, nb++; s = ""; f = factor(n); for (k=1, #f~, s = concat(Str(f[k, 1]), s)); new = eval(s); if (isprime(new), return (nb)); if (already(new, v), return(nb)); v = concat(v, new); n = new; ); }
|
|
CROSSREFS
|
|
|
KEYWORD
|
sign,base
|
|
AUTHOR
|
|
|
EXTENSIONS
|
Previous b-file transformed into an auxiliary file and new b-file restricted to known values by Michel Marcus, Sep 27 2015
|
|
STATUS
|
approved
|
|
|
|