当在求 ( a / b ) % MOD 时由于有分数不能直接求
要转化成 a % MOD * ( x % MOD + MOD ) % MOD
x的求法如下:
void EGCD( LL a, LL b, LL &x, LL &y ){
if( b == 0 ){
x = 1;
y = 0;
return;
}
EGCD( b, a % b, x, y );
LL tp = x;
x = y;
y = tp - a / b * y;
}
LL &x, &y;
EGCD( b, MOD, x, y );