#include <bits/stdc++.h>
using namespace std;
typedef complex<double> Point;
typedef pair<Point, Point> Line;
const double EPS = 1e-9;

// 内積 |a||b|cosθ
double dot(Point a, Point b){
	// return a.real() * b.real() + a.imag() * b.imag();
	return (conj(a) * b).real();
}
// 射影　直線に点pを投影　点pから直線に下ろした垂線の足
Point projection(Line l, Point p){
    Point a = l.first, b = l.second, c = p;
    double t = dot(p - a, b - a) / norm(b - a); 
    return a + t * (b - a);
}

int main(){
    double x, y;
    Point p1, p2;
    cin >> x >> y;
    p1 = {x, y};
    cin >> x >> y;
    p2 = {x, y};
    Line l = {p1, p2};
    int q;
    cin >> q;
    for(int i = 0; i < q; i++){
        cin >> x >> y;
        Point c = {x, y};
        Point res = projection(l, c);
        res = c + 2.0 * (res - c);
        printf("%.10f %.10f\n", res.real(), res.imag());
    }
}

### Prompt

Create a solution in CPP for the following problem:

For given three points p1, p2, p, find the reflection point x of p onto p1p2.

<image>

Constraints

* 1 ≤ q ≤ 1000
* -10000 ≤ xi, yi ≤ 10000
* p1 and p2 are not identical.

Input


xp1 yp1 xp2 yp2
q
xp0 yp0
xp1 yp1
...
xpq−1 ypq−1


In the first line, integer coordinates of p1 and p2 are given. Then, q queries are given for integer coordinates of p.

Output

For each query, print the coordinate of the reflection point x. The output values should be in a decimal fraction with an error less than 0.00000001.

Examples

Input

0 0 2 0
3
-1 1
0 1
1 1


Output

-1.0000000000 -1.0000000000
0.0000000000 -1.0000000000
1.0000000000 -1.0000000000


Input

0 0 3 4
3
2 5
1 4
0 3


Output

4.2400000000 3.3200000000
3.5600000000 2.0800000000
2.8800000000 0.8400000000

### Response

```cpp
#include <bits/stdc++.h>
using namespace std;
typedef complex<double> Point;
typedef pair<Point, Point> Line;
const double EPS = 1e-9;

// 内積 |a||b|cosθ
double dot(Point a, Point b){
	// return a.real() * b.real() + a.imag() * b.imag();
	return (conj(a) * b).real();
}
// 射影　直線に点pを投影　点pから直線に下ろした垂線の足
Point projection(Line l, Point p){
    Point a = l.first, b = l.second, c = p;
    double t = dot(p - a, b - a) / norm(b - a); 
    return a + t * (b - a);
}

int main(){
    double x, y;
    Point p1, p2;
    cin >> x >> y;
    p1 = {x, y};
    cin >> x >> y;
    p2 = {x, y};
    Line l = {p1, p2};
    int q;
    cin >> q;
    for(int i = 0; i < q; i++){
        cin >> x >> y;
        Point c = {x, y};
        Point res = projection(l, c);
        res = c + 2.0 * (res - c);
        printf("%.10f %.10f\n", res.real(), res.imag());
    }
}

```