ClearAll[CurzonNumberQ]
CurzonNumberQ[b_Integer][n_Integer]:=PowerMod[b,n,b n+1]==b n
val=Select[Range[100000],CurzonNumberQ[2]];
Take[val,50]
val[[1000]]

val=Select[Range[100000],CurzonNumberQ[4]];
Take[val,50]
val[[1000]]

val=Select[Range[100000],CurzonNumberQ[6]];
Take[val,50]
val[[1000]]

val=Select[Range[100000],CurzonNumberQ[8]];
Take[val,50]
val[[1000]]

val=Select[Range[100000],CurzonNumberQ[10]];
Take[val,50]
val[[1000]]