OFFSET
1,1
COMMENTS
The cycle that gets entered consists of the primes in A287865. It appears that the mapping p -> gpf(2*p+1) produces no other cycles.
Conjecture: under repeated mapping all primes ultimately enter the same cycle.
EXAMPLE
Prime 2 is in the sequence as it maps to 5. And so is 17 as it maps to 7. The primes 3, 5, 7, 11, 13, 19, 23 and 47 are not included, as they are part of the cycle itself (and hence considered zero iterations away from the cycle).
MAPLE
gpf:= n -> max(numtheory:-factorset(n)):
filter:= proc(n) local S, t, x;
t:= gpf(2*n+1);
if t = n then return false fi;
S:= {n, t};
x:= t;
do
x:= gpf(2*x+1);
if member(x, S) then return (x = t) fi;
S:= S union {x};
od;
end proc:
select(filter, [seq(ithprime(i), i=1..1000)]); # Robert Israel, Feb 03 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Johannes M.V.A. Koelman, Jan 28 2025
STATUS
approved