ZYM/BH-int2048-2023
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

upd: bug located

Browse Source
This commit is contained in:
ZhuangYumin
2023-11-01 18:04:40 +08:00
parent 7c8904ff51
commit 5a67934594
2 changed files with 21 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

27
src/int2048.cpp
View File

@ -639,16 +639,27 @@ int2048 GetInv(const int2048 &B, int n) {
remain.val[2 * n / int2048::kNum] = kPow10[2 * n % int2048::kNum];
remain.num_length = (2 * n) + 1;
UnsignedMultiply(sub_soluton_copy_1, &current_B);
UnsignedMinus(remain, &sub_soluton_copy_1);
for (int i = 64; i > 0; i >>= 1) {
int2048 tmp_B(current_B);
tmp_B.UnsignedMultiplyByInt(i);
if (UnsignedCmp(remain, tmp_B) >= 0) {
res += i;
UnsignedMinus(remain, &tmp_B);
int cmp = UnsignedCmp(remain, sub_soluton_copy_1);
if (cmp == 0) return std::move(res);
if (cmp > 0) {
UnsignedMinus(remain, &sub_soluton_copy_1);
for (int i = 64; i > 0; i >>= 1) {
int2048 tmp_B(current_B);
tmp_B.UnsignedMultiplyByInt(i);
if (UnsignedCmp(remain, tmp_B) >= 0) {
res += i;
UnsignedMinus(remain, &tmp_B);
}
}
return std::move(res);
} else {
UnsignedMinus(sub_soluton_copy_1, &remain);
for (int i = 64; i > 0; i >>= 1) {
int2048 tmp_B(current_B);
tmp_B.UnsignedMultiplyByInt(i);
res -= i;
}
}
return std::move(res);
}
inline void UnsignedDivide(int2048 &A, const int2048 *pB) {
int2048 B(*pB);