Skip to content

Instantly share code, notes, and snippets.

Embed
What would you like to do?
(* On p=4n+1 := a^2+b^2 *)
sl = Table[n^2, {n, 1, 1500}];
t = Flatten@Table[a^2 + b^2, {a, 1, 1500}, {b, a + 1, 1501}];
Complement [Select[Prime[Range[100000]], Mod[#, 4] == 1 &], t]
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
You can’t perform that action at this time.