#include <iostream>
#include <cmath>
#include <complex>
using namespace std;

using vec = complex<double>;
using line = pair<vec,vec>;

double inner_product(vec u,vec v){
    return real(u*conj(v));
}
double cross_product(vec u,vec v){
    return imag(conj(u)*v);
}

double norm(vec& u){
    return abs(u);
}

vec projection(line l,vec p){//p???l???????°???±
    vec s=l.first, t=l.second;
    double k = inner_product(t-s,p-s)/inner_product(t-s,t-s);
    return (1.0-k)*s+k*t;
}

vec reflection(line l,vec p){
    return 2.0*projection(l,p)-p;
}


int main(void){
    int xp1,yp1,xp2,yp2;
    cin>>xp1>>yp1>>xp2>>yp2;
    line l(vec(xp1,yp1),vec(xp2,yp2));
    int q;
    cin>>q;
    cout<<fixed;
    cout.precision(10);
    for(int i=0;i<q;i++){
        int xp,yp;
        cin>>xp>>yp;
        vec X=reflection(l,vec(xp,yp));
        cout<<X.real()<<' '<<X.imag()<<endl;
    }

    return 0;
}
### Prompt

Please create a solution in cpp 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>
#include <complex>
using namespace std;

using vec = complex<double>;
using line = pair<vec,vec>;

double inner_product(vec u,vec v){
    return real(u*conj(v));
}
double cross_product(vec u,vec v){
    return imag(conj(u)*v);
}

double norm(vec& u){
    return abs(u);
}

vec projection(line l,vec p){//p???l???????°???±
    vec s=l.first, t=l.second;
    double k = inner_product(t-s,p-s)/inner_product(t-s,t-s);
    return (1.0-k)*s+k*t;
}

vec reflection(line l,vec p){
    return 2.0*projection(l,p)-p;
}


int main(void){
    int xp1,yp1,xp2,yp2;
    cin>>xp1>>yp1>>xp2>>yp2;
    line l(vec(xp1,yp1),vec(xp2,yp2));
    int q;
    cin>>q;
    cout<<fixed;
    cout.precision(10);
    for(int i=0;i<q;i++){
        int xp,yp;
        cin>>xp>>yp;
        vec X=reflection(l,vec(xp,yp));
        cout<<X.real()<<' '<<X.imag()<<endl;
    }

    return 0;
}
```