【SPOJ-BTCODE_A】Traversing Grid【GCD】

题意：

给出一个起点，一个终点。对于一个点(x, y)，有四种变换。

1 (y, x)

2 (x, -y)

3 (x + y, y)

4 (2 * x, y)

问能否从起点到达终点。

这题做了2个月...最后发现是少了一个特判。

我们发现，对于前3个操作，(x, y)可以变成(x, -y) (-x, y) (-x, -y) (x + y, y) (x+ y, -y) (-(x + y), y) (-(x + y), -y) ...

这些点的共同特性是绝对值的gcd一样。

对于第四个操作，gcd有可能变为2倍。

所以我们求出gcd(x1, y1)，gcd(x2, y2)，然后判断倍数是不是2的非负整数次幂就好了。

注意特判啊...

#include <cstdio>

typedef long long LL;

LL gcd(LL a, LL b) {
	return b ? gcd(b, a % b) : a;
}

LL _abs(LL a) {
	return a > 0 ? a : -a;
}

int main() {
	int T; scanf("%d", &T);
	while(T--) {
		LL x1, y1, x2, y2; scanf("%lld%lld%lld%lld", &x1, &y1, &x2, &y2);
		if(x1 == 0 && y1 == 0) {
			if(x2 == 0 && y2 == 0) printf("YES\n");
			else printf("NO\n");
			continue;
		}
		if(x2 == 0 && y2 == 0) {
			printf("NO\n");
			continue;
		}
		x1 = _abs(x1); y1 = _abs(y1); x2 = _abs(x2); y2 = _abs(y2);
		LL gcd1 = gcd(x1, y1), gcd2 = gcd(x2, y2);
		if(gcd2 % gcd1) printf("NO\n");
		else {
			LL d = gcd2 / gcd1;
			if(d == (d & -d)) printf("YES\n");
			else printf("NO\n");
		}
	}
	return 0;
}