#include <iostream>
#include <vector>
#include <numeric>
#include <string>
#include <complex>
#include <iomanip>
using namespace std;
typedef complex<double> P;
const double EPS = 1e-10;

namespace std
{
bool operator<(P &a, P &b)
{
    return real(a) != real(b) ? real(a) < real(b) : imag(a) < imag(b);
}
} // namespace std

double dot(P a, P b)
{
    return real(a) * real(b) + imag(a) * imag(b);
}

double cross(P a, P b)
{
    return real(a) * imag(b) - imag(a) * real(b);
}

int main(void)
{

    double a, b, c, d;
    cin >> a >> b >> c >> d;

    P p1(a, b);
    P p2(c, d);
    int q;
    cin >> q;
    for (int i = 0; i < q; ++i)
    {
        double a, b;
        cin >> a >> b;
        P p(a, b);

        P x;
        x = p1 + dot(p - p1, p2 - p1) * (p2 - p1) / pow(abs(p2 - p1), 2.0);

        P reflection;
        reflection = x + x - p;

        cout << fixed << setprecision(15) << reflection.real() << " " << reflection.imag() << endl;
    }
}

### Prompt

Please formulate 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 <vector>
#include <numeric>
#include <string>
#include <complex>
#include <iomanip>
using namespace std;
typedef complex<double> P;
const double EPS = 1e-10;

namespace std
{
bool operator<(P &a, P &b)
{
    return real(a) != real(b) ? real(a) < real(b) : imag(a) < imag(b);
}
} // namespace std

double dot(P a, P b)
{
    return real(a) * real(b) + imag(a) * imag(b);
}

double cross(P a, P b)
{
    return real(a) * imag(b) - imag(a) * real(b);
}

int main(void)
{

    double a, b, c, d;
    cin >> a >> b >> c >> d;

    P p1(a, b);
    P p2(c, d);
    int q;
    cin >> q;
    for (int i = 0; i < q; ++i)
    {
        double a, b;
        cin >> a >> b;
        P p(a, b);

        P x;
        x = p1 + dot(p - p1, p2 - p1) * (p2 - p1) / pow(abs(p2 - p1), 2.0);

        P reflection;
        reflection = x + x - p;

        cout << fixed << setprecision(15) << reflection.real() << " " << reflection.imag() << endl;
    }
}

```