cp's OEIS Frontend

A135025 Let b(1) = 2; and for n>= 2, if b(n-1) < prime(n) then b(n) = b(n-1) + prime(n) otherwise b(n) = b(n-1) - prime(n). The sequence gives the indices n where b(n-1) < b(n) < b(n+1).

%I A135025 #29 Aug 16 2022 07:39:47
%S A135025 4,9,22,57,146,367,946,2507,6634,17777,48522,133107,369020,1028405,
%T A135025 2880288,8100949,22877146,64823569,184274932,525282741,1501215194,
%U A135025 4299836187,12340952050,35486796313,102220582466,294917666855,852123981582,2465458792769
%N A135025 Let b(1) = 2; and for n>= 2, if b(n-1) < prime(n) then b(n) = b(n-1) + prime(n) otherwise b(n) = b(n-1) - prime(n). The sequence gives the indices n where b(n-1) < b(n) < b(n+1).
%C A135025 The b sequence, prefixed by 0, is A008348. The low points in b are 1 less than the terms of the present sequence, and are given in A309226. - _N. J. A. Sloane_, Aug 31 2019
%e A135025 b(1) = 2
%e A135025 b(2) = 5
%e A135025 b(3) = 0
%e A135025 b(4) = 7
%e A135025 b(5) = 18
%e A135025 b(3) < b(4) < b(5), so 4 is the first term of the sequence.
%p A135025 B := proc(n) option remember ; if n = 1 then 2; else if procname(n-1)-ithprime(n) < 0 then procname(n-1)+ithprime(n) ; else procname(n-1)-ithprime(n) ; fi; fi; end: A135025 := proc(n) option remember ; if n = 1 then 4; else for a from procname(n-1)+1 do if B(a-1) < B(a) and B(a) < B(a+1) then RETURN(a) ; fi; od: fi; end: for n from 1 do printf("%d,\n",A135025(n)) ; od: # _R. J. Mathar_, Feb 06 2009
%t A135025 B[n_] := B[n] = If[n == 1, 2, If[B[n-1] - Prime[n] < 0, B[n-1] + Prime[n], B[n-1] - Prime[n]]];
%t A135025 a[n_] := a[n] = If[n == 1, 4, For[k = a[n-1]+1, True, k++, If[B[k-1] < B[k] && B[k] < B[k+1], Return[k]]]];
%t A135025 Table[Print[n, " ", a[n]]; a[n], {n, 1, 16}] (* _Jean-François Alcover_, Aug 16 2022, after _R. J. Mathar_ *)
%Y A135025 Cf. A008348, A135026, A309226.
%K A135025 nonn,more
%O A135025 1,1
%A A135025 Lior Deutsch (liorde(AT)gmail.com), Feb 10 2008
%E A135025 New term added by Lior Deutsch (liorde(AT)gmail.com), Oct 17 2008
%E A135025 Definition corrected and entry revised by _Robert Israel_, _Michel Marcus_, and _N. J. A. Sloane_, Sep 29 2014
%E A135025 a(17)-a(28) from _Giovanni Resta_, Oct 02 2019