Time Limit: 10 Sec Memory Limit: 256 MB
Description
给你一对数a,b,你可以任意使用(a,b), (a,-b), (-a,b), (-a,-b), (b,a), (b,-a), (-b,a), (-b,-a)这些向量,问你能不能拼出另一个向量(x,y)。
说明:这里的拼就是使得你选出的向量之和为(x,y)
Input
第一行数组组数t,(t<=50000)
接下来t行每行四个整数a,b,x,y (-2*109<=a,b,x,y<=2*109)
Output
t行每行为Y或者为N,分别表示可以拼出来,不能拼出来
Sample Input
3
2 1 3 3
1 1 0 1
1 0 -2 3
Sample Output
Y
N
Y
HINT
样例解释:
第一组:(2,1)+(1,2)=(3,3)
第三组:(-1,0)+(-1,0)+(0,1)+(0,1)+(0,1)=(-2,3)
这题我们可以发现向量
(−a,−b)=(a,b)∗−1
,
(−b,−a)=(b,a)∗−1
,
(b,−a)=−1∗(−b,a)
,
(−a,b)=−1∗(a,−b)
,所以我们可以发现我们只需要用4个向量
(a,b),(b,a),(−b,a),(a,−b)
就够了。
设
(x,y)=x1∗(a,b)+x2∗(b,a)+x3∗(−b,a)+x4∗(a,−b)
则我们可以得到方程
{(x1+x4)∗a+(x2−x3)∗b=x(x2+x3)∗a+(x1−x4)∗b=y
再令
f1=x1+x4,f2=x2−x3,f3=x2+x3,f4=x1−x4
则方程变为
{f1∗a+f2∗b=xf3∗a+f4∗b=y
我们可以发现,要使
x1,x2,x3,x4
有整数解,则
f1+f4
和
f2+f3
均应为偶数。
然后我们只需求出
f1,f2,f3,f4
即可,不难发现,我们可以用扩展欧几里得来解决。
则
⎧⎩⎨f1=X∗xgcd(a,b)+bgcd(a,b)∗tf2=Y∗xgcd(a,b)−agcd(a,b)∗t(t为整数)
其中的X,Y为扩欧(这里求的方程为aX+bY=gcd(a,b),所以需要将X,Y乘以
xgcd(a,b)
)得到的一组解,同时可以发现当
gcd(a,b)∤x
时,就没有整数解了。
同理可得第二个方程
⎧⎩⎨f3=X∗ygcd(a,b)+bgcd(a,b)∗pf4=Y∗ygcd(a,b)−agcd(a,b)∗p(p为整数)
然后我们就得到了
⎧⎩⎨f1+f4=X∗xgcd(a,b)+Y∗ygcd(a,b)+bgcd(a,b)∗t−agcd(a,b)∗pf2+f3=X∗ygcd(a,b)+Y∗xgcd(a,b)+bgcd(a,b)∗p−agcd(a,b)∗t
可以发现,
f1+f4和f2+f3
的奇偶性与
t和p
的值有关,然后我们又可以由奇偶性发现,我们只需要判断
t=0,p=0;t=1,p=0;t=0,p=1;t=1,p=1
这四种情况即可!
#include<iostream>
#include<cstdio>
using namespace std;
int a,b,c,d,X,Y;
int gcd(int x,int y)
{
if (y==0)
{
X=1;Y=0;
return x;
}
int z=gcd(y,x%y);
int t=X;
X=Y;
Y=t-Y*(x/y);
return z;
}
bool check(int t,int p)
{
if ((X*c+Y*d+b*t-a*p)%2==0&&(Y*c+X*d+a*t-b*p)%2==0) return true;
return false;
}
int main()
{
int T;
scanf("%d",&T);
while (T--)
{
scanf("%d%d%d%d",&a,&b,&c,&d);
if (a==0||b==0)
{
if (a==0&&b==0) printf("N\n");
if (a==0)
{
if (c%b==0&&d%b==0) printf("Y\n");else printf("N\n");
}
if (b==0)
{
if (c%a==0&&d%a==0) printf("Y\n");else printf("N\n");
}
continue;
}
int Gcd=gcd(a,b);
if (c%Gcd!=0||d%Gcd!=0) { printf("N\n");continue;}
a/=Gcd;b/=Gcd;c/=Gcd;d/=Gcd;
if (check(0,0)||check(1,0)||check(0,1)||check(1,1))
printf("Y\n");else printf("N\n");
}
return 0;
}