# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/

Search: id:a181103
Showing 1-1 of 1

%I A181103 #17 Sep 08 2022 08:45:54
%S A181103 2,3,7,13,17,19,23,29,41,47,59,83,89,101,107,137,173,179,239,251,257,
%T A181103 269,293,311,347,389,419,449,467,557,563,599,641,647,701,761,809,827,
%U A181103 839,857,887,947,953,971,977,1019,1049,1091,1109,1151,1163,1217,1229
%N A181103 a(n) is the least prime greater than a(n-1) such that abs(a(n-2) + a(n-1) - a(n)) is prime; a(1)=2, a(2)=3.
%C A181103 Partial products are in A181737.
%H A181103 Nathaniel Johnston, <a href="/A181103/b181103.txt">Table of n, a(n) for n = 1..10000</a>
%p A181103 A181103 := proc(n) option remember: local p: if(n<=2)then return n+1: fi: p:=procname(n-1): do p:=nextprime(p): if(isprime(abs(p-procname(n-2)-procname(n-1))))then return p: fi: od: end: seq(A181103(n),n=1..53); # _Nathaniel Johnston_, May 27 2011
%t A181103 nxt[{a_,b_}]:=Module[{lp=NextPrime[b]},While[!PrimeQ[Abs[a+b-lp]],lp= NextPrime[ lp]];{b,lp}]; Transpose[NestList[nxt,{2,3},60]][[1]] (* _Harvey P. Dale_, Feb 08 2015 *)
%o A181103 (Magma) A181103:=function(n); S:=[2, 3]; if n le 2 then return S[n]; end if; p:=5; for k in [3..n] do while not IsPrime(Abs(S[k-2]+S[k-1]-p)) do p:=NextPrime(p); end while; Append(~S, p); if k eq n then return p; end if; p:=NextPrime(p); end for; end function; [ A181103(n): n in [1..55] ]; // _Klaus Brockhaus_, Dec 19 2010
%Y A181103 Cf. A181737.
%K A181103 nonn,easy
%O A181103 1,1
%A A181103 _Giovanni Teofilatto_, Oct 03 2010
%E A181103 Definition refined and sequence extended by _R. J. Mathar_, Oct 12 2010
%E A181103 Index in definition corrected by _R. J. Mathar_, Oct 26 2010

# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE