%I #6 Oct 15 2013 22:31:55
%S 13,89,509,1933,1637,1669,5953,27701,12853,19333,36389,28229,19609,
%T 82073,102701,89689,134513,31397,212701,815729,461717,1155733,927869,
%U 360653,396733,3422813,3279841,370261,5845193,6085441,6808273,9810653
%N a(n) is the smallest prime p of the form 4k+1 such that nextprime[p]-p=4n.
%e n=9: nextprime[a(9)]-a(9)=12889-12853=36=4n.
%t {m=4, r=1}; f[x_] := (Prime[x+1]-Prime[x])/m t=Table[0, {100}]; Do[s=f[n]; s1=Mod[Prime[n+1], m]; If[IntegerQ[s]&&Equal[s1, r]&&s<101&&t[[s]]==0, t[[s]]=Prime[n]], {n, 1, 1000000}]; t
%Y Cf. A002144, A002145, A082073-A082076, A082098.
%K nonn
%O 1,1
%A _Labos Elemer_, Apr 14 2003
|