nyoj-1250-exgcd

时间:2022-06-24 17:22:58

机器人

时间限制:1000 ms  |  内存限制:65535 KB
难度:4
描述

Dr. Kong 设计的机器人卡尔非常活泼,既能原地蹦,又能跳远。由于受软硬件设计所限,机器人卡尔只能定点跳远。若机器人站在(X,Y)位置,它可以原地蹦,但只可以在(X,Y),(X,-Y),(-X,Y),(-X,-Y),(Y,X),(Y,-X),(-Y,X),(-Y,-X)八个点跳来跳去。

现在,Dr. Kong想在机器人卡尔身上设计一个计数器,记录它蹦蹦跳跳的数字变化(S,T),即,路过的位置坐标值之和。

你能帮助Dr. Kong判断机器人能否蹦蹦跳跳,拼出数字(S,T)吗?

假设机器人卡尔初始站在(0,0)位置上。

输入
第一行: K 表示有多少组测试数据。
接下来有K行,每行:X Y S T

1≤K≤10000 -2*109 <= X , Y, S, T <= 2*109
数据之间有一个空格。

输出
对于每组测试数据,输出一行:Y或者为N,分别表示可以拼出来,不能拼出来
样例输入
3
2 1 3 3
1 1 0 1
1 0 -2 3
样例输出
Y
N
Y
来源
第七届河南省程序设计大赛
  八种变换方式,有四对是呈相反状态的,例如(X,Y)和(-X,-Y)。所以只要对剩下的四个状态走若干次(可以是负数次表示走对立状态)
能达到(S,T)就好了。不妨令剩下的四种状态为(X,Y) (X,-Y) (Y,X) (Y,-X) ,对应的次数为a1,a2,a3,a4,我们有: S=(a1+a2)*X+(a3+a4)*Y
T=(a1-a2)*Y+(a3-a4)*X, 容易看出这两个线性方程可以用exgcd求解,如果S,T 不是gcd(X,Y)的整数倍显然不会成立。算出通解之后
注意到(a1+a2)+(a1-a2)=2*a1  (a3+a4)+(a3-a4)=2*a3 , 枚举一下系数的奇偶情况看是否对应的两项相加都可以是偶数即可。
  (不保证算法正确性,,但是AC了。
 #include<iostream>
#include<cstdio>
using namespace std;
#define LL long long
#define mp make_pair
#define pb push_back
#define inf 0x3f3f3f3f
void exgcd(LL a,LL b,LL &d,LL &x,LL &y){
if(!b){d=a;x=;y=;}
else{exgcd(b,a%b,d,y,x);y-=x*(a/b);}
}
int main(){
int t;
cin>>t;
while(t--){
LL X,Y,S,T,d;
LL A,B,C,D,x,y;
scanf("%lld%lld%lld%lld",&X,&Y,&S,&T);
exgcd(X,Y,d,A,B);
if(!(S%d==&&T%d==)){
puts("N");
}
else{
bool ok=;
LL d1=Y/d,d2=X/d;
for(int i=-;i<=;++i){
for(int j=-;j<=;++j){
LL _A=A*S/d+i*d1,_B=B*S/d-i*d2;
LL _C=A*T/d+j*d1,_D=B*T/d-j*d2; if((_A+_D)%==&&(_B+_C)%==)
ok=; }
}
ok?puts("Y"):puts("N");
}
}
return ;
}
/*
3
2 1 3 3
1 1 0 1
1 0 -2 3
*/

相关文章