Sampling a Uniform Random Solution of a Quadratic Equation Modulo pk
An n-ary integral quadratic form is a formal expression Q(x1,...,xn)=∑1≤i,j≤naijxixj in n-variables x1,...,xn, where aij=aji∈Z. We present a poly(n,k,logp,logt) randomized algorithm that given a quadratic form Q(x1,...,xn), a prime p, a positive integer k and an integer t, samples a uniform solution of Q(x1,...,xn)≡tmodpk.