#include <iostream>
#include <cmath>

using namespace std;
typedef double scalar;
class point{
	public:
	scalar x;
	scalar y;
	void print(){
		cout << x << " " << y; 
	}
	scalar norm2(void){
		return x * x + y * y;
	}
	scalar norm(void){
		return sqrt(norm2());
	}
	point operator =(point a){
		this->x = a.x;
		this->y = a.y;
		return *this;
	}
	point operator +(point a){
		point r;
		r.x = this->x + a.x;
		r.y = this->y + a.y;
		return r;
	}
	point operator -(point a){
		point r;
		r.x = this->x - a.x;
		r.y = this->y - a.y;
		return r;
	}
};

point operator *(scalar t, point a){
	point r;
	r.x = t * a.x; 
	r.y = t * a.y; 
	return r;
}

void projection(point p, point q, point r, point &x){
	point qp = q - p;
	point pr = p - r;
	scalar k = (qp.x * pr.x + qp.y * pr.y)/qp.norm2();
	x = p - (k * qp);
	
}

void reflection(point p, point q, point r, point &s){
	point x;
	projection(p,q,r,x);
	s = 2.0 * x - r;
}

int main(){
	point p1, p2, p, x;
	int q;
	cout.precision(20);
	cin >> p1.x >> p1.y >> p2.x >> p2.y;
	cin >> q;
	for(int i = 0;i < q;i++){
		cin >> p.x >> p.y;
		reflection(p1, p2, p, x);
		x.print();
		cout << endl;
	}
	return 0;
}
### Prompt

Your challenge is to write a Cpp solution to 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 <cmath>

using namespace std;
typedef double scalar;
class point{
	public:
	scalar x;
	scalar y;
	void print(){
		cout << x << " " << y; 
	}
	scalar norm2(void){
		return x * x + y * y;
	}
	scalar norm(void){
		return sqrt(norm2());
	}
	point operator =(point a){
		this->x = a.x;
		this->y = a.y;
		return *this;
	}
	point operator +(point a){
		point r;
		r.x = this->x + a.x;
		r.y = this->y + a.y;
		return r;
	}
	point operator -(point a){
		point r;
		r.x = this->x - a.x;
		r.y = this->y - a.y;
		return r;
	}
};

point operator *(scalar t, point a){
	point r;
	r.x = t * a.x; 
	r.y = t * a.y; 
	return r;
}

void projection(point p, point q, point r, point &x){
	point qp = q - p;
	point pr = p - r;
	scalar k = (qp.x * pr.x + qp.y * pr.y)/qp.norm2();
	x = p - (k * qp);
	
}

void reflection(point p, point q, point r, point &s){
	point x;
	projection(p,q,r,x);
	s = 2.0 * x - r;
}

int main(){
	point p1, p2, p, x;
	int q;
	cout.precision(20);
	cin >> p1.x >> p1.y >> p2.x >> p2.y;
	cin >> q;
	for(int i = 0;i < q;i++){
		cin >> p.x >> p.y;
		reflection(p1, p2, p, x);
		x.print();
		cout << endl;
	}
	return 0;
}
```