cp's OEIS Frontend

A215029 To get a(n), start with m=n, let p = largest prime < m, set m = m-p if m>0, or m=m+p if m <= 0; repeat until p=2 has been processed; set a(n) = m.

%I A215029 #19 Dec 01 2018 09:28:13
%S A215029 0,1,2,1,-1,0,0,1,1,2,-1,0,0,1,1,2,-1,0,0,1,1,2,-1,0,0,1,0,1,0,1,1,2,
%T A215029 0,1,0,1,0,1,-1,0,1,2,0,1,-1,0,1,2,0,1,0,1,0,1,-1,0,1,2,-1,0,0,1,1,2,
%U A215029 -1,0,1,2,0,1,0,1,-1,0,0,1,0,1,0,1,1,2,-1,0,0,1,0,1,0,1,1,2,-1,0,1,2,-1,0,0,1,0,1,1,2,0,1,0,1,-1,0,0,1,0,1,1
%N A215029 To get a(n), start with m=n, let p = largest prime < m, set m = m-p if m>0, or m=m+p if m <= 0; repeat until p=2 has been processed; set a(n) = m.
%H A215029 Antti Karttunen, <a href="/A215029/b215029.txt">Table of n, a(n) for n = 0..20000</a>
%H A215029 Antti Karttunen, <a href="/A215029/a215029.txt">Data supplement: n, a(n) computed for n = 0..100003</a>
%H A215029 StackExchange, <a href="http://math.stackexchange.com/questions/176394/a-prime-number-pattern">A prime number pattern</a>, Jul 29 2012.
%p A215029 f:=proc(n) local m,a,i,p;
%p A215029 if n <= 2 then RETURN(n); fi;
%p A215029 m:=n; a:=n;
%p A215029 for i from 1 to n do
%p A215029 p:=prevprime(m);
%p A215029 if a>0 then a:=a-p else a:=a+p; fi;
%p A215029 m:=p;
%p A215029 if m <= 2 then RETURN(a); fi;
%p A215029 od;
%p A215029 # should never reach here
%p A215029 print("ERROR");
%p A215029 end;
%p A215029 [seq(f(i),i=0..120)];
%o A215029 (PARI) A215029(n) = if(n<=2,n,my(mp=precprime(n-1),d=n); while(mp>0, if(d>0, d -= mp, d += mp); mp = precprime(mp-1)); (d)); \\ _Antti Karttunen_, Nov 28 2018
%Y A215029 Cf. A215030-A215035, A214912.
%K A215029 sign
%O A215029 0,3
%A A215029 _N. J. A. Sloane_, Aug 05 2012