一、题目描述
二、算法分析说明与代码编写指导
设要经过 x 次循环才可以跳出循环体,则有 a + cx ≡ b (% 2^k)
设 m = 2^k,则上式可转化为:cx - 2^k · y = b - a。
套用扩展欧几里得算法并将 x 调整为最小正整数解即为答案。如果 gcd(c, -2^k) 不能整除 (b - a) 则无解。
三、AC 代码(0 ms)
尝试提交时发现,无论是 cmath 还是 algorithm 的 abs() 函数,POJ 的古董编译器均无法正确重载。
#include<cstdio>
#pragma warning(disable:4996)
//Solve ax + by = GCD(a, b) for x, y.
template<class _Ty> _Ty extgcd(const _Ty& a, const _Ty&