pub fn eval(params: &VdfParams, x: &BigUint) -> VdfOutput
Evaluate the VDF: compute y = x^(2^T) mod N and proof. This is the slow part — T sequential squarings.