定義

１)

ある整数b以下の全ての素数p₁ , p₂ , ... , p_n を求めておく。 c=1 とする。(cはインクリメント用)

２)

とし

となる を計算する。

３)

がスムーズでないなら、cをインクリメントし 2) にリターン。

４)

スムーズな に対して 次のようなベクトル

がp_iを因子をして含む個数)

５)

2)〜4)を繰り返し スムーズなに対するベクトル

が 上一次従属になるまで集める。

６)

とすると

は 上一次従属なのでd_i=1 となるものが存在する。

d_i=1 となるような を集め、それらを とならべ直す。

７)

これはの中に各p_iが偶数個含まれているということである。

よって、

such that

８)

a₁²a₂²...a_k²=(a₁a₂...a_k)² より

such that (a₁a₂...a_k)²=r²

2)より

よって、

９)

として、0<x,y<m を満たすx,yを求める。