login
a(n) = If[Floor[2*Mod[(Prime[n + 1] - Prime[n])*PrimePi[n]/n, 8]] == 15, Prime[n]]
2

%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