1. 程式人生 > >[HAOI2011]向量

[HAOI2011]向量

input 復制 lag div nor 不一定 style des col

題目描述

給你一對數a,b,你可以任意使用(a,b), (a,-b), (-a,b), (-a,-b), (b,a), (b,-a), (-b,a), (-b,-a)這些向量,問你能不能拼出另一個向量(x,y)。

說明:這裏的拼就是使得你選出的向量之和為(x,y)

輸入輸出格式

輸入格式:

第一行數組組數t,(t<=50000)

接下來t行每行四個整數a,b,x,y (-2*10^9<=a,b,x,y<=2*10^9)

輸出格式:

t行每行為Y或者為N,分別表示可以拼出來,不能拼出來

輸入輸出樣例

輸入樣例#1: 復制
3
2 1 3 3
1 1 0 1
1 0 -2 3
輸出樣例#1: 復制
Y
N
Y

說明

樣例解釋:

第一組:(2,1)+(1,2)=(3,3)

第三組:(-1,0)+(-1,0)+(0,1)+(0,1)+(0,1)=(-2,3)

題解來自520dalao%%%

$k(a,b)+q(b,a)+w(a,-b)+c(b,-a)=(x,y)$

--> $(k+w)a+(q+c)b=x, (k-w)b+(q-c)a=y$ .

  由裴蜀定理可得:$ax+by=c$ ,x和y有整數解的充要條件是 $gcd(a,b)|c$ ,

證明:令$a=pgcd(a,b),\ b=qgcd(a,b)$ ,則原式=$(px+qy)gcd(a,b)=c$ ,顯然因為$gcd(a,b)$ 為整數,而要使x和y為整數,則$gcd(a,b)|c$ 。

  我們回到開始的方程組$(k+w)a+(q+c)b=x,(k-w)b+(q-c)a=y$ 。由裴蜀定理易得:

使$(k+w),(q+c),(k-w),(q-c)$ 均為整數的充要條件是:

$gcd(a,b)|x$ 且$gcd(a,b)|y$ 。但是註意到(k+w),(k-w)有整數解不一定k和w有整數解((q+c)和(q-c)是同理的)。此時不妨設$(k+w)=f,(k-w)=g$ ,則$k=(f+g)/2,w=(f-g)/2$ ,因為$2|(f+g)$ 且$2|(f-g)$ ,顯然要使k和w均為整數則f和g均為偶數或均為奇數((q+c)和(q-c)同理)

  於是我們考慮這四種情況:

  1、當(k+w)、(k-w)、(q+c)、(q-c)均為偶數時,$(k+w)a+(q+c)b=x$ 提公因數2結合$gcd(a,b)|x $ => $ 2gcd(a,b)|x$ 同理 $gcd(a,b)|y$

  2、當(k+w)和(k-w)為偶數,(q+c)和(q-c)為奇數時,$(k+w)a+(q+c)b=x$ 先左右兩邊同加b,再提公因數2 結合$gcd(a,b)|x$ --> $2gcd(a,b)|x+b$ 同理 $2gcd(a,b)|y+a$

  3、當(k+w)和(k-w)為奇數,(q+c)和(q-c)為偶數時,$(k+w)a+(q+c)b=x$ 先左右兩邊同加a,再提公因數2 結合$gcd(a,b)|x $ -->$2gcd(a,b)|x+a$ 同理 $2gcd(a,b)|y+b$

  4、當(k+w)、(k-w)、(q+c)、(q-c)均為奇數時,$(k+w)a+(q+c)b=x$ 先左右兩邊同加a+b,再提公因數2 結合$gcd(a,b)|x$ --> $2gcd(a,b)|x+a+b$ 同理 $2gcd(a,b)|y+a+b$

  只要滿足上述的任意一種情況,則說明本題k、w、q、c有整數解,說明可行,否則說明無解。

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5 #include<cmath>
 6 using namespace std;
 7 typedef long long lol;
 8 lol a,b,x,y;
 9 lol gcd(lol a,lol b)
10 {
11   if (!b) return a;
12   return gcd(b,a%b);
13 }
14 int main()
15 {int T,flag;
16   cin>>T;
17   while (T--)
18     {
19       scanf("%lld%lld%lld%lld",&a,&b,&x,&y);
20       lol d=gcd(2*a,2*b);
21       flag=0;
22       if (x%d==0&&y%d==0) flag=1;
23       if ((x+a)%d==0&&(y+b)%d==0) flag=1;
24       if ((x+b)%d==0&&(y+a)%d==0) flag=1;
25       if ((x+a+b)%d==0&&(y+a+b)%d==0) flag=1;
26       if (flag) printf("Y\n");
27       else printf("N\n");
28     }
29 }

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)

[HAOI2011]向量