%I #34 Dec 10 2019 17:45:27
%S 0,1,1,2,2,3,0,0,3,0,0,1,0,0,4,0,1,1,0,0,2,0,1,0,0,0,0,1,0,0,0,0,0,0,
%T 0,0,1,0,2,0,1,0,0,1,0,0,0,0,1,0,0,0,0,0,0,1,0,0,1,0,0,0,0,0,2,0,0,0,
%U 1,0,0,2,0,0,0,1,0,1,0,0,0,1,0,0,0,0,0,1
%N Define a map from the primes to the primes by f(p) = (p-1)/2 if that is prime, or else (p+1)/2 if that is prime, and otherwise is undefined. Start with the n-th prime and iterate f until we cannot go any further; a(n) is the number of steps.
%C For each prime, the end of the trajectory is reached when one cannot generate another prime number from it.
%C For example, p(3) = 5 -> 2 (1 iteration), so a(3)=1. Also p(5) = 11 -> 5 -> 2 (2 iterations), 23 -> 11 -> 5 -> 2 (3 iterations) and 47 -> 23 -> 11 -> 5 -> 2 (4 iterations). Hence a(3) = 1, a(5) = 2, a(9) = 3 and a(15) = 4.
%C a(n) = 0 for n = 1, 7, 8, 10, 11, 13, 14, 16, 19, 20, 22, 24, 25, ... The corresponding primes are A176902(n) = 2, 17, 19, 29, 31, 41, 43, ... .
%C The sequence of the last terms of the trajectories begins with 2, 2, 2, 2, 2, 2, 17, 19, 2, 29, 31, 19, 41, 43, 2, 53, 29, 31, 67, ...
%C The following table gives the trajectories of the smallest prime requiring 0, 1, 2, 3, 4, 5, 6, iterations:
%C +------------+----------+------------------------------------------+
%C | Number of | smallest | trajectory |
%C | iterations | prime | |
%C +------------+----------+------------------------------------------+
%C | 0 | 2 | 2 |
%C | 1 | 3 | 3 -> 2 |
%C | 2 | 7 | 7 -> 3 -> 2 |
%C | 3 | 13 | 13 -> 7 -> 3 -> 2 |
%C | 4 | 47 | 47 -> 23 -> 11 -> 5 -> 2 |
%C | 5 | 2879 | 2879 -> 1439 -> 719 -> 359 -> 179 -> 89 |
%C | 6 | 1065601 | 1065601 -> 532801 -> 266401 -> 133201 -> |
%C | | | 66601 -> 33301 -> 16651 |
%C +------------+----------+------------------------------------------+
%e a(15) = 4 because prime(15) = 47 and 47 -> 23 -> 11 -> 5 -> 2 with 4 iterations.
%p for n from 1 to 100 do:
%p ii:=0:it:=0:p:=ithprime(n):
%p for i from 1 to 100 while(ii=0) :
%p p1:=(p-1)/2:p2:=(p+1)/2:
%p if type(p1,prime)=false and type(p2,prime)=false
%p then
%p ii:=1:printf(`%d, `,it):
%p else
%p it:=it+1:
%p if isprime(p1)
%p then
%p p:=p1:
%p else
%p p:=p2:
%p fi:
%p fi:
%p od:
%p od:
%t f[p_] := If[PrimeQ[(q = (p-1)/2)], q, If[PrimeQ[(r = (p+1)/2)], r, 0]]; g[n_] := -2 + Length @ NestWhileList[f, n, #>0 &]; g /@ Select[Range[457], PrimeQ] (* _Amiram Eldar_, Nov 16 2019 *)
%Y Cf. A000040, A005383, A005385, A176902.
%Y The underlying map is A330310.
%K nonn
%O 1,4
%A _Michel Lagneau_, Oct 25 2019
|