John D. Cook 2/15/2026

Wagon’s algorithm in Python

Read Original

This article provides a step-by-step implementation of Stan Wagon's algorithm in Python to find integer solutions for x² + y² = p, where p is a large odd prime. It covers finding a quadratic non-residue modulo p, computing a square root of -1 mod p, and applying a modified Euclidean algorithm with a precision-aware stopping condition using Python's isqrt. The tutorial includes complete code and demonstrates the solution for a 255-bit prime.

Wagon’s algorithm in Python

Comments

No comments yet

Be the first to share your thoughts!

Browser Extension

Get instant access to AllDevBlogs from your browser

Top of the Week

1
The Beautiful Web
Jens Oliver Meiert 2 votes
3
LLM Use in the Python Source Code
Miguel Grinberg 1 votes
4
Wagon’s algorithm in Python
John D. Cook 1 votes