#include <bits/stdc++.h>

using namespace std;
int q;
struct point
{
    double x;
    double y;
};
point vt(point a,point b)
{
    point vt;
    vt.x=b.x-a.x;
    vt.y=b.y-a.y;
    return vt;
}
double scalar(point a,point b)
{
    double scalar=0;
    scalar=sqrt((b.x-a.x)*(b.x-a.x)+(b.y-a.y)*(b.y-a.y));
    return scalar;
}
double mul1(point a,point b)
{
    double mul=a.x*b.x+a.y*b.y;
    return mul;
}
point mul2(point a,double b)
{
    point c;
    c.x=a.x*b;
    c.y=a.y*b;
    if (c.x==-0) c.x=0;
    if (c.y==-0) c.y=0;
    return c;
}
int main()
{
    point p1,p2;
    cin>>p1.x>>p1.y;
    cin>>p2.x>>p2.y;
    cin>>q;
    while (q--)
    {
        point p;
        cin>>p.x>>p.y;
        point v=vt(p1,p);
        point u=vt(p1,p2);
        //cout<<v.x<< " "<<v.y<<endl;
        //cout<<u.x<< " "<<u.y<<endl;
        double su=scalar(p1,p2);
        //cout<<su<<endl;
        point ans;
        ans=mul2(u,mul1(u,v)/(su*su));
        point x1;
        x1.x=ans.x+p1.x;
        x1.y=ans.y+p1.y;
        point px1=vt(p,x1);
        px1=mul2(px1,2);
        cout<<setprecision(10)<<fixed<<px1.x+p.x<<" "<<px1.y+p.y<<endl;
    }
    return 0;
}


### Prompt

Please provide a cpp coded solution to the problem described below:

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;
int q;
struct point
{
    double x;
    double y;
};
point vt(point a,point b)
{
    point vt;
    vt.x=b.x-a.x;
    vt.y=b.y-a.y;
    return vt;
}
double scalar(point a,point b)
{
    double scalar=0;
    scalar=sqrt((b.x-a.x)*(b.x-a.x)+(b.y-a.y)*(b.y-a.y));
    return scalar;
}
double mul1(point a,point b)
{
    double mul=a.x*b.x+a.y*b.y;
    return mul;
}
point mul2(point a,double b)
{
    point c;
    c.x=a.x*b;
    c.y=a.y*b;
    if (c.x==-0) c.x=0;
    if (c.y==-0) c.y=0;
    return c;
}
int main()
{
    point p1,p2;
    cin>>p1.x>>p1.y;
    cin>>p2.x>>p2.y;
    cin>>q;
    while (q--)
    {
        point p;
        cin>>p.x>>p.y;
        point v=vt(p1,p);
        point u=vt(p1,p2);
        //cout<<v.x<< " "<<v.y<<endl;
        //cout<<u.x<< " "<<u.y<<endl;
        double su=scalar(p1,p2);
        //cout<<su<<endl;
        point ans;
        ans=mul2(u,mul1(u,v)/(su*su));
        point x1;
        x1.x=ans.x+p1.x;
        x1.y=ans.y+p1.y;
        point px1=vt(p,x1);
        px1=mul2(px1,2);
        cout<<setprecision(10)<<fixed<<px1.x+p.x<<" "<<px1.y+p.y<<endl;
    }
    return 0;
}


```