#include <iostream>
#include <limits>
#include <iomanip>
using namespace std;

class Point{
public:
    long double x;
	long double y;
	Point(long double x=0,long double y=0) :x(x),y(y){ }
	Point operator +(Point &p) {
		return Point(x + p.x, y + p.y);
	}
	Point operator-(Point &p) {
		return Point(x - p.x, y - p.y);
	}
	Point operator *(long double a) {
		return Point(x*a, y*a);
	}
};

typedef struct {
	Point p1, p2;
}Line;

int main() {

    
	Line s;
	cin >> s.p1.x; cin >> s.p1.y; cin >> s.p2.x; cin >> s.p2.y;
	int q;
	cin >> q;
	for (int i = 0;i < q;i++) {
		Point p,q,r;
		cin >> p.x; cin >> p.y;
		 double prod = ((s.p2.x - s.p1.x)*(p.x - s.p1.x)) + ((s.p2.y - s.p1.y)*(p.y - s.p1.y));
		q.x = s.p1.x +(s.p2.x - s.p1.x)*(prod/((s.p1.x-s.p2.x)*(s.p1.x - s.p2.x)+ (s.p1.y - s.p2.y)*(s.p1.y - s.p2.y)));
		q.y = s.p1.y + (s.p2.y - s.p1.y)*(prod / ((s.p1.x - s.p2.x)*(s.p1.x - s.p2.x) + (s.p1.y - s.p2.y)*(s.p1.y - s.p2.y)));
		r.x = p.x + (q.x - p.x) * 2; r.y = p.y + (q.y - p.y) * 2;
		cout << fixed << setprecision(10)<< r.x << " " << r.y << endl;
	}

}
### 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 <iostream>
#include <limits>
#include <iomanip>
using namespace std;

class Point{
public:
    long double x;
	long double y;
	Point(long double x=0,long double y=0) :x(x),y(y){ }
	Point operator +(Point &p) {
		return Point(x + p.x, y + p.y);
	}
	Point operator-(Point &p) {
		return Point(x - p.x, y - p.y);
	}
	Point operator *(long double a) {
		return Point(x*a, y*a);
	}
};

typedef struct {
	Point p1, p2;
}Line;

int main() {

    
	Line s;
	cin >> s.p1.x; cin >> s.p1.y; cin >> s.p2.x; cin >> s.p2.y;
	int q;
	cin >> q;
	for (int i = 0;i < q;i++) {
		Point p,q,r;
		cin >> p.x; cin >> p.y;
		 double prod = ((s.p2.x - s.p1.x)*(p.x - s.p1.x)) + ((s.p2.y - s.p1.y)*(p.y - s.p1.y));
		q.x = s.p1.x +(s.p2.x - s.p1.x)*(prod/((s.p1.x-s.p2.x)*(s.p1.x - s.p2.x)+ (s.p1.y - s.p2.y)*(s.p1.y - s.p2.y)));
		q.y = s.p1.y + (s.p2.y - s.p1.y)*(prod / ((s.p1.x - s.p2.x)*(s.p1.x - s.p2.x) + (s.p1.y - s.p2.y)*(s.p1.y - s.p2.y)));
		r.x = p.x + (q.x - p.x) * 2; r.y = p.y + (q.y - p.y) * 2;
		cout << fixed << setprecision(10)<< r.x << " " << r.y << endl;
	}

}
```