次のコードを検討してください。
mpz_t x, n, out;
mpz_init_set_ui(x, 2UL);
mpz_init_set_ui(n, 7UL);
mpz_init(out);
mpz_invert(out, x, n);
gmp_printf ("%Zd\n", out);//prints 4. 2 * 4 (mod 7) = 1. OK
mpz_powm_ui(out, x, -1, n);//prints 1. 2 * 1 (mod 7) = 2. How come?
gmp_printf ("%Zd\n", out);
mpz_clear(x);
mpz_clear(n);
mpz_clear(out);
ドキュメントによると、負の指数をサポートするはずですが、mpz_powm 関数が負の指数をどのように処理するのか理解できません。n を法として数値を -1 に上げるのは、n を法としてそれを反転することと同じだと思います。ここで何が欠けていますか?