Function mod_sqrt

pub fn mod_sqrt(a: u64, p: u64) -> Option<u64>

x² = a (mod p)を満たすxを一つ求める。