RosettaCodeData/Task/Curzon-numbers/Mathematica/curzon-numbers.math
2023-07-01 13:44:08 -04:00

21 lines
436 B
Text

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]]