19/11/26 22:19:47.56 sE/nea3J.net
>>256
では、BIT SHIFT ではなく、浮動小数点演算にして、以下の様にすれば速くなるかもしれません。
(i) 初期条件
η = 1/a; // 多倍長の浮動小数点
x_0 = 1;
y_0 = a * x_0 - b;
(ii) 漸化式
x_{k+1} = x_k - (int_N)(y_k * η);
y_{k+1} = a * x_{k+1} - b;
但し、int_N は、多倍長の浮動小数点を多倍長整数に直す cast。