Jacobi symbol: Difference between revisions

<syntaxhighlight lang="ruby">func jacobi(na, kn) {

assert(kn > 0, "#{kn} must be positive")

assert(kn.is_odd, "#{kn} must be odd")

while (na %= kn) {

varif v = n(a.valuation(2is_even) {

t *= (-1)**v if (k%8var ~~v [3,5]= a.valuation(2)

n >> t *= (-1)**v if (n%8 ~~ [3,5])

(n,k) a >>= (k,n)v

t = -t if ([n%4, k%4] == [3,3])}

kn==1 ? t : 0

for na in (0..50), kn in (0..50) {

assert_eq(jacobi(na, 2*kn + 1), kronecker(na, 2*kn + 1))