//投影
#include<bits/stdc++.h>
using namespace std;
const double PI=acos(-1);
const double EPS=1e-6;
bool is_Zero(double &a){
	return a<EPS&&a>-EPS;
}
const int INFTY=1<<29;
struct Point{
	double x,y;
	Point(double x,double y):x(x),y(y){}
	Point(){}
	Point(const Point&a){x=a.x;y=a.y;}
	Point operator+(const Point&a){
		return Point(x+a.x,y+a.y);
	}
	Point operator-(const Point&a){
		return Point(x-a.x,y-a.y);
	}
	Point operator*(double k){
		return Point(k*x,k*y);
	}
	Point operator/(double k){
		return Point(x/k,y/k);
	} 
};
typedef Point Vector;
Point operator*(double k,const Point&a){
	return Point(k*a.x,k*a.y);
}
double norm(const Point&a){
	return a.x*a.x+a.y*a.y;
}
double abs(const Point&a){
	return sqrt(norm(a));
}
double dot_product(const Point&a,const Point&b){
	return a.x*b.x+a.y*b.y;
}
double cross_product(const Point&a,const Point&b){
	return a.x*b.y-a.y*b.x;
}
Point projection_Point(Point&p1,Point&p2,Point&p){
	Vector base=p2-p1;
	double r=dot_product(base,p-p1)/norm(base);
	return p1+base*r;
}

Point reflection_Point(Point&p1,Point&p2,Point&p){
	Point temp=projection_Point(p1,p2,p);
	return 2*temp-p;
}

int main(){
	Point p1,p2;
	Point temp;
	scanf("%lf%lf%lf%lf",&p1.x,&p1.y,&p2.x,&p2.y);
	int q;
	scanf("%d",&q);
	while(q--){
		scanf("%lf%lf",&temp.x,&temp.y);
		Point a=reflection_Point(p1,p2,temp);
		printf("%.10lf %.10lf\n",a.x,a.y);
	}
	return 0;
}

### 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;
const double PI=acos(-1);
const double EPS=1e-6;
bool is_Zero(double &a){
	return a<EPS&&a>-EPS;
}
const int INFTY=1<<29;
struct Point{
	double x,y;
	Point(double x,double y):x(x),y(y){}
	Point(){}
	Point(const Point&a){x=a.x;y=a.y;}
	Point operator+(const Point&a){
		return Point(x+a.x,y+a.y);
	}
	Point operator-(const Point&a){
		return Point(x-a.x,y-a.y);
	}
	Point operator*(double k){
		return Point(k*x,k*y);
	}
	Point operator/(double k){
		return Point(x/k,y/k);
	} 
};
typedef Point Vector;
Point operator*(double k,const Point&a){
	return Point(k*a.x,k*a.y);
}
double norm(const Point&a){
	return a.x*a.x+a.y*a.y;
}
double abs(const Point&a){
	return sqrt(norm(a));
}
double dot_product(const Point&a,const Point&b){
	return a.x*b.x+a.y*b.y;
}
double cross_product(const Point&a,const Point&b){
	return a.x*b.y-a.y*b.x;
}
Point projection_Point(Point&p1,Point&p2,Point&p){
	Vector base=p2-p1;
	double r=dot_product(base,p-p1)/norm(base);
	return p1+base*r;
}

Point reflection_Point(Point&p1,Point&p2,Point&p){
	Point temp=projection_Point(p1,p2,p);
	return 2*temp-p;
}

int main(){
	Point p1,p2;
	Point temp;
	scanf("%lf%lf%lf%lf",&p1.x,&p1.y,&p2.x,&p2.y);
	int q;
	scanf("%d",&q);
	while(q--){
		scanf("%lf%lf",&temp.x,&temp.y);
		Point a=reflection_Point(p1,p2,temp);
		printf("%.10lf %.10lf\n",a.x,a.y);
	}
	return 0;
}

```