|
%I #16 Mar 10 2019 21:40:44
%S 2,3,5,7,11,13,17,19,23,25,29,31,35,37,41,43,47,49,53,55,59,61,67,71,
%T 73,77,79,83,85,89,91,95,97,101,103,107,109,113,115,119,121,127,131,
%U 133,137,139,143,145,149,151,155,157,161,163,167,169,173,175,179
%N Numbers m such that A099427(m) = 2.
%C A099427(a(n)) = 2;
%C primes and squares of primes greater than 9 are subsequences, cf. A000040, A001248, A000430;
%C GCD(A099427(a(n)-1), A099427(a(n))) = 1;
%C a(n) = A038179(n) for n <= 22.
%C The next term divisible by 3 is a(137)=429. - _Joe Slater_, Jan 10 2017
%C All terms after the first are odd, since A099427(n) == n+1 (mod 2) for n >= 3. - _Robert Israel_, Jan 10 2017
%H Reinhard Zumkeller, <a href="/A192489/b192489.txt">Table of n, a(n) for n = 1..10000</a>
%p A099427:= proc(n) option remember; 1 + igcd(n,procname(n-1)) end proc:
%p A099427(1):= 1:
%p select(A099427=2, [$1..1000]); # _Robert Israel_, Jan 10 2017
%t (* b = A099427 *) b[1] = 1; b[n_] := b[n] = GCD[n, b[n - 1]] + 1;
%t Select[Range[200], b[#] == 2&] (* _Jean-François Alcover_, Mar 10 2019 *)
%o (Haskell)
%o a192489 n = a192489_list !! (n-1)
%o a192489_list = f 2 1 where
%o f n x | x' == 2 = n : f (n+1) x'
%o | otherwise = f (n+1) x'
%o where x' = 1 + gcd n x
%K nonn
%O 1,1
%A _Reinhard Zumkeller_, Jul 02 2011
|
