_ := bigmod.NewModulus(addOne(decodeHex(t, B))) b, _ := bigmod.NewNat().SetBytes(decodeHex(t, B), q) gcd, err := bigmod.NewNat().GCDVarTime(a, b) // GCD doesn't work if a and b are both even, but LCM handles it. if err == nil { if got := strings.TrimLeft(hex.EncodeToString(gcd.Bytes(p)), "0"); got != GCD { t.Fatalf("unexpected GCD: got %s, want %s", got, GCD) } } lcm, err := totient(p, q) if oddDivisorLargerThan(decodeHex(t, GCD)) { if err != errDivisorTooLarge { t.Fatalf("expected divisor too large...