E. Hag's Khashba+计算几何+旋转+重心

E. Hag's Khashba

time limit per test

3 seconds

memory limit per test

256 megabytes

input

standard input

output

standard output

Hag is a very talented person. He has always had an artist inside him but his father forced him to study mechanical engineering.

Yesterday he spent all of his time cutting a giant piece of wood trying to make it look like a goose. Anyway, his dad found out that he was doing arts rather than studying mechanics and other boring subjects. He confronted Hag with the fact that he is a spoiled son that does not care about his future, and if he continues to do arts he will cut his 25 Lira monthly allowance.

Hag is trying to prove to his dad that the wooden piece is a project for mechanics subject. He also told his dad that the wooden piece is a strictly convex polygon with $n$

vertices.

Hag brought two pins and pinned the polygon with them in the $1$

-st and $2$-nd vertices to the wall. His dad has $q$

queries to Hag of two types.

• $1$

$f$ $t$: pull a pin from the vertex $f$, wait for the wooden polygon to rotate under the gravity force (if it will rotate) and stabilize. And then put the pin in vertex $t$

. $2$ $v$: answer what are the coordinates of the vertex $v$

• .

Please help Hag to answer his father's queries.

You can assume that the wood that forms the polygon has uniform density and the polygon has a positive thickness, same in all points. After every query of the 1-st type Hag's dad tries to move the polygon a bit and watches it stabilize again.

Input

The first line contains two integers $n$

and $q$ ( $3\le n\le 10000$, $1\le q\le 200000$

) — the number of vertices in the polygon and the number of queries.

The next $n$

lines describe the wooden polygon, the $i$-th line contains two integers $x_i$ and $y_i$ ( $|x_i|,|y_i|\le {10}^8$) — the coordinates of the $i$

-th vertex of the polygon. It is guaranteed that polygon is strictly convex and the vertices are given in the counter-clockwise order and all vertices are distinct.

The next $q$

lines describe the queries, one per line. Each query starts with its type $1$ or $2$. Each query of the first type continues with two integers $f$ and $t$ ( $1\le f,t\le n$) — the vertex the pin is taken from, and the vertex the pin is put to and the polygon finishes rotating. It is guaranteed that the vertex $f$ contains a pin. Each query of the second type continues with a single integer $v$ ( $1\le v\le n$

) — the vertex the coordinates of which Hag should tell his father.

It is guaranteed that there is at least one query of the second type.

Output

The output should contain the answer to each query of second type — two numbers in a separate line. Your answer is considered correct, if its absolute or relative error does not exceed ${10}^{-4}$

.

Formally, let your answer be $a$

, and the jury's answer be $b$. Your answer is considered correct if $\frac{|a-b|}{max(1,|b|)}\le {10}^{-4}$

Examples

Input

Copy

3 4
0 0
2 0
2 2
1 1 2
2 1
2 2
2 3

Output

Copy

3.4142135624 -1.4142135624
2.0000000000 0.0000000000
0.5857864376 -1.4142135624

Input

Copy

3 2
-1 1
0 0
1 1
1 1 2
2 1

Output

Copy

1.0000000000 -1.0000000000

Note

In the first test note the initial and the final state of the wooden polygon.

Red Triangle is the initial state and the green one is the triangle after rotation around $(2,0)$

.

In the second sample note that the polygon rotates $180$

degrees counter-clockwise or clockwise direction (it does not matter), because Hag's father makes sure that the polygon is stable and his son does not trick him.

这个图求重心你是绝对不能用坐标取中的

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define all(a) (a).begin(),(a).end()
#define pll pair<ll,ll>
#define vi vector<int>
#define pb push_back
const int inf=0x3f3f3f3f;
ll rd(){
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
const int N=1e4+10;
double x[N],y[N],r[N],a[N];
double D;
double sx,sy;
int main(){
#ifdef happy
    freopen("in.txt","r",stdin);
#endif
    int n=rd(),q=rd();
    rep(i,1,n)
    x[i]=rd(),y[i]=rd();
    double sa=0;
    rep(i,2,n-1){
        double px=(x[1]+x[i]+x[i+1])/3.0;
        double py=(y[1]+y[i]+y[i+1])/3.0;
        double a=abs((x[i]-x[1])*(y[i+1]-y[1])-(x[i+1]-x[1])*(y[i]-y[1]));
        sx+=a*px;
        sy+=a*py;
        sa+=a;
    }
    sx/=sa;
    sy/=sa;
    rep(i,1,n){
        r[i]=hypot(x[i]-sx,y[i]-sy);
        a[i]=atan2(y[i]-sy,x[i]-sx);
        while(a[i]<a[i-1])a[i]+=atan(1)*8;
    }
    int p[2]={1,2};
    double px[2]={x[1],x[2]};
    double py[2]={y[1],y[2]};

    while(q--){
        int k=rd();
        if(k==1){
            int i;
            int f=rd(),t=rd();
            if(p[0]==f)i=1,p[0]=t;
            else i=0,p[1]=t;
            sx=px[i];
            sy=py[i]-r[p[i]];
            D=atan(1)*2-a[p[i]];
            px[i^1]=sx+r[t]*cos(D+a[t]);
            py[i^1]=sy+r[t]*sin(D+a[t]);
        }else{
            int i=rd();
            double tx=sx+r[i]*cos(D+a[i]);
            double ty=sy+r[i]*sin(D+a[i]);
            printf("%.12lf %.12lf\n",tx,ty);
        }
    }
}