%I #4 Dec 16 2024 16:02:36
%S 16141,36389,37907,69263,73189,86629,95651,103237
%N a(n) = If[Floor[2*Mod[(Prime[n + 1] - Prime[n])*PrimePi[n]/n, 8]] == 12, Prime[n]]
%t a = Flatten[Table[If[Floor[2*Mod[(Prime[n + 1] - Prime[n])*PrimePi[n]/n, 8]] == 12, Prime[n], {}], {n, 1, 10000}]]
%t Table[If[Floor[2*Mod[(Prime[n+1]-Prime[n])*PrimePi[n]/n,8]]==12,Prime[n],Nothing],{n,10000}] (* _Harvey P. Dale_, Dec 16 2024 *)
%Y Cf. A001359.
%K nonn
%O 0,1
%A _Roger L. Bagula_, Jun 27 2005