|
%I #31 Jan 01 2022 18:06:16
%S 9,10,15,16,21,22,25,26,28,33,34,36,39,40,45,46,49,50,52,55,56,58,63,
%T 64,66,69,70,75,76,78,81,82,85,86,88,91,92,94,96,99,100,105,106,111,
%U 112,115,116,118,120,122,123,124,126,129,130,133,134,136,141,142
%N Composites c with result 0 under iterations of {r mod (max prime p <= r)} starting at r = c.
%C Intersection of A002808 and A175089.
%C Composites c such that A121559(c) = 0. - _Michel Marcus_, Aug 22 2014
%H Robert Israel, <a href="/A175090/b175090.txt">Table of n, a(n) for n = 1..10000</a>
%e Iteration procedure for a(3) = 15: 15 mod 13 = 2, 2 mod 2 = 0.
%p g:= proc(x) if isprime(x) then 0 else x mod prevprime(x) fi end proc:
%p f:= proc(x) local y; y:= x; while y > 1 do y:= g(y) od; y = 0 end proc:
%p select(not(isprime) and f, [$4..200]); # _Robert Israel_, Feb 09 2015
%t Composites := Select[Range[2, 200], ! PrimeQ[#] &]; Select[Composites, PrimeQ[# - NextPrime[#, -1]] &] (* _Carlos Eduardo Olivieri_, Feb 09 2015 *)
%Y Cf. A007917 and A064722 (both for the iterations).
%K nonn
%O 1,1
%A _Jaroslav Krizek_, Jan 28 2010
%E Missing term 55 inserted, more terms added, _Michel Marcus_, Aug 22 2014
|
