BZOJ2299 [HAOI2011]向量 【裴蜀定理】

时间:2022-03-02 21:40:13

题目链接

BZOJ2299

题解

题意就是给我们四个方向的向量\((a,b),(b,a),(-a,b),(b,-a)\),求能否凑出\((x,y)\)

显然我们就可以得到一对四元方程组,用裴蜀定理判断一下方程有没有解即可

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<map>
#define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt)
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define mp(a,b) make_pair<int,int>(a,b)
#define cls(s) memset(s,0,sizeof(s))
#define cp pair<int,int>
#define LL long long int
using namespace std;
const int maxn = 100005,maxm = 100005,INF = 1000000000;
inline int read(){
int out = 0,flag = 1; char c = getchar();
while (c < 48 || c > 57){if (c == '-') flag = -1; c = getchar();}
while (c >= 48 && c <= 57){out = (out << 3) + (out << 1) + c - 48; c = getchar();}
return out * flag;
}
LL a,b,x,y,tmp;
LL gcd(LL a,LL b){return b ? gcd(b,a % b) : a;}
int main(){
int T = read();
while (T--){
a = read(); b = read(); x = read(); y = read();
if (a < b) swap(a,b);
tmp = gcd(gcd(2 * a * b,a * a + b * b),a * a - b * b);
if ((b * x - a * y) % tmp || (b * x + a * y) % tmp || (a * x - b * y) % tmp || (a * x + b * y) % tmp)
puts("N");
else puts("Y");
}
return 0;
}