%I #2 Mar 30 2012 17:34:19
%S 19609,25471,35617,44293,58831,89689
%N a(n) = If[Floor[2*Mod[(Prime[n + 1] - Prime[n])*PrimePi[n]/n, 8]] == 15, Prime[n]]
%t a = Flatten[Table[If[Floor[2*Mod[(Prime[n + 1] - Prime[n])*PrimePi[n]/n, 8]] == 15, Prime[n], {}], {n, 1, 10000}]]
%Y Cf. A001359.
%K nonn
%O 0,1
%A _Roger L. Bagula_, Jun 27 2005