%I #21 Oct 01 2021 20:44:40
%S 2,0,23,53,409,293,211,1341,1343,2179,3967,15705,16033,19635,31425,
%T 24281,31429,31431,31433,155959,38501,58831,203713,268343,206699,
%U 370311,370313,370315,370317,1349591,1357261,1272749,1357265,1357267,2010801,2010803,2010805,2010807
%N Smallest number at distance 2n from nearest prime.
%C a(0) = 2. For n > 0, let f(m) = minimal distance from m to closest prime (excluding m itself). The a(n) = min { m : f(m) = 2n }.
%C f(m) is tabulated in A051700. - _R. J. Mathar_, Nov 18 2007
%F a(n) = A051652(2*n). - _Sean A. Irvine_, Oct 01 2021
%p A051700 := proc(m) if m <= 2 then op(m+1,[2,1,1]) ; else min(nextprime(m)-m,m-prevprime(m)) ; fi ; end: A051728 := proc(n) local m ; if n = 0 then RETURN(2); else for m from 0 do if A051700(m) = 2 * n then RETURN(m) ; fi ; od: fi ; end: seq(A051728(n),n=0..20) ; # _R. J. Mathar_, Nov 18 2007
%t a[n_] := Module[{m}, If[n == 0, Return[2], For[m = 0, True, m++, If[Min[NextPrime[m]-m, m-NextPrime[m, -1]] == 2*n, Return[m]]]]]; Table[Print[an = a[n]]; an, {n, 0, 33}] (* _Jean-François Alcover_, Feb 11 2014, after _R. J. Mathar_ *)
%t Join[{2},With[{t=Table[{n,Min[n-NextPrime[n,-1],NextPrime[n]-n]},{n,0,1358000}]},Table[SelectFirst[t,#[[2]]==2k&],{k,33}]][[All,1]]] (* _Harvey P. Dale_, Aug 13 2019 *)
%Y Related sequences: A023186-A023188, A046929-A046931, A051650, A051652, A051697-A051702, A051728-A051730.
%Y Cf. A132470.
%K nonn,easy,nice
%O 0,1
%A _N. J. A. Sloane_
%E More terms from _James A. Sellers_, Dec 07 1999
%E More terms from _Amiram Eldar_, Aug 28 2021