Little Desprado2 is a student of Springfield Flowers Kindergarten. On this day, he had just learned how to draw triangles on grid coordinate paper. However, he soon found it very dull, so he came up with a more interesting question:
He had drawn two integral points of the triangle on the grid paper, and he denotes them (x1,y1)(x1,y1) and (x2,y2)(x2,y2). Now, he wanted to know the answer to the following question: where can he draw the third point (x3,y3)(x3,y3) so that the area of the triangle is positive but minimized?
Obviously, he can't solve this problem because he is too young and simple. Can you tell him the answer?
Please note that your answer's coordinates must consist of integers because he is drawing on grid paper, and the triangle shouldn't be a degenerated triangle to keep the area positive.
Input
The first line contains one integer TT (1≤T≤500001≤T≤50000), denoting the number of Little Desprado2's queries.
For each test case, there's a single line contains four integers x1, y1, x2, y2x1, y1, x2, y2 (−109≤x1, y1, x2, y2≤109−109≤x1, y1, x2, y2≤109) seperated by spaces, denoting two points are at (x1,y1)(x1,y1) and (x2,y2)(x2,y2), respectively.
It is guaranteed that the two points won't coincide.
Output
For each test case, print two integers x3, y3x3, y3 (−1018≤x3, y3≤1018−1018≤x3, y3≤1018) in a separated line, denoting your answer.
If there are multiple answers, you can print any one of them. It is guaranteed that there exists a solution in the above range.
Example
input
3
1 0 1 4
0 1 0 9
0 0 2 2
output
2 0 1 1 -1 0
题意: 给出两点坐标,需要输出第三个点坐标,使这三个点构成一个面积最小的三角形,要求点坐标必须为整数。
分析: 设已给出点A(x1, y1)和点B(x2, y2)的坐标,并且设向量AC为(x, y),则三角形面积为向量AB叉乘向量AC,将向量AB表示为(x2-x1, y2-y1),令a = x2-x1,b = y2-y1,则只需求出令a*y-b*x最小的整数解,而这个问题可以通过扩展欧几里得解决,根据扩展欧几里得可知a*y-b*x=c有整数解的条件是c = k*gcd(a, -b),也就是说该三角形面积只能是k*gcd(a, -b)中某一个取值,最小情况下自然是k = 1了,解出x和y后再加上A点坐标就是C点坐标了。
具体代码如下:
- #include
- #include
- #include
- #include
- #include
- #include
- #define int long long
- using namespace std;
-
- int exgcd(int a, int b, int &x, int &y) //扩展欧几里得
- {
- if(b == 0)
- {
- x = 1;
- y = 0;
- return a;
- }
- int eax = exgcd(b, a%b, x, y);
- int t = x;
- x = y;
- y = t - (a/b)*y;
- return eax;
- }
-
- signed main()
- {
- int T;
- cin >> T;
- while(T--){
- int x1, x2, y1, y2;
- scanf("%lld%lld%lld%lld", &x1, &y1, &x2, &y2);
- int a = x2-x1, b = y2-y1;
- if(a == 0){
- printf("%lld %lld\n", x1+1, y1);
- continue;
- }
- if(b == 0){
- printf("%lld %lld\n", x1, y1+1);
- continue;
- }
- int x, y;
- exgcd(-b, a, x, y);
- printf("%lld %lld\n", x+x1, y+y1);
- }
-
- return 0;
- }