In updating: "vx = (avx * 1ll * (vx + 10) + bvx) % 21 - 10" make sure you use real value of vx, not that divided by gcd(vx, vy);