POJ 2115

http://acm.pku.edu.cn/JudgeOnline/problem?id=2115

从题目可以推出 a+c*x = b (mod 2^k)    =>    c*x = (b-a) (mod 2^k)

这个一个非常典型的求解模线性方程

因为最后的解有多组 第一组是最小的 所以取第一个解

但第一解有时为负

因为所有解模2^k/d同于 所以最有解(x+2^k) % (2^k/d)  其中d=GCD(a, n)

 

#include <stdio.h> typedef long long llg; llg extgcd(llg a, llg b, llg &x, llg &y) { if(b == 0) { x = 1; y = 0; return a; } llg d = extgcd(b, a % b, x, y); llg t = x; x = y; y = t - a / b * y; return d; } llg mod_linear_solver(llg a, llg b, llg n) { llg x, y; llg d = extgcd(a, n, x, y); if(b % d == 0) { llg t = n / d; x = (x * (b / d)) % n; return (x + n) % t; } return -1; } int main(void) { llg a, b, c, k; while(scanf("%lld%lld%lld%lld", &a, &b, &c, &k) != EOF) { if(a == 0 && b == 0 && c == 0 && k == 0) break; llg m = 1LL << k; if(b >= m) { printf("FOREVER/n"); continue; } if(a-b == 0) { printf("0/n"); continue; } llg t = b - a; if(t < 0) t += m; llg n = mod_linear_solver(c, t, m); if(n < 0) printf("FOREVER/n"); else printf("%lld/n", n); } return 0; }