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\leq n \leq 10\,000$, $1 \leq q \leq 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|\leq 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
3 4
0 0
2 0
2 2
1 1 2
2 1
2 2
2 3

Output
3.4142135624 -1.4142135624
2.0000000000 0.0000000000
0.5857864376 -1.4142135624

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

Output
1.0000000000 -1.0000000000



-----Note-----

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

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.

Source: codeforces | Difficulty: VERY_HARD
Tags: Geometry


Solutions:
import math

def getCoordinate(gx, gy, alpha, x, y):
	x1 = gx + x * math.cos(alpha) - y * math.sin(alpha)
	y1 = gy + x * math.sin(alpha) + y * math.cos(alpha)
	return (x1, y1)

def getAngle(gx, gy, x, y):
	x = x - gx
	y = y - gy
	cos = x / math.sqrt(x ** 2 + y ** 2)
	alpha = math.acos(cos)
	if y < 0:
		alpha = -alpha
	return math.pi / 2 - alpha
(n, q) = list(map(int, input().split(' ')))
x = [0] * n
y = [0] * n
for i in range(n):
	(x[i], y[i]) = list(map(int, input().split(' ')))
r = [0] * q
f = [0] * q
t = [0] * q
v = [0] * q
for i in range(q):
	l = list(map(int, input().split(' ')))
	r[i] = l[0]
	if r[i] == 1:
		f[i] = l[1] - 1
		t[i] = l[2] - 1
	else:
		v[i] = l[1] - 1
gx = 0
gy = 0
s = 0
for i in range(n):
	ip = i + 1
	if ip == n:
		ip = 0
	ds = x[i] * y[ip] - x[ip] * y[i]
	s += ds
	gx += (x[i] + x[ip]) * ds
	gy += (y[i] + y[ip]) * ds
s /= 2
gx /= 6 * s
gy /= 6 * s
angles = [0] * n
for i in range(n):
	angles[i] = getAngle(gx, gy, x[i], y[i])
for i in range(n):
	x[i] -= gx
	y[i] -= gy
alpha = 0
fix = {0, 1}
for i in range(q):
	if r[i] == 2:
		(currX, currY) = getCoordinate(gx, gy, alpha, x[v[i]], y[v[i]])
		print('%.6f %.6f' % (currX, currY))
	else:
		if len(fix) == 2:
			fix.remove(f[i])
		for j in fix:
			(currX, currY) = getCoordinate(gx, gy, alpha, x[j], y[j])
			alpha = angles[j]
			(gx, gy) = (currX, currY - math.sqrt(x[j] ** 2 + y[j] ** 2))
		fix.add(t[i])


from sys import stdin, stdout
from math import atan2, sqrt, sin, cos, pi

def main():
	(n, q) = stdin.readline().split()
	n = int(n)
	q = int(q)
	lines = stdin.readlines()
	x = [0] * (n + 1)
	y = [0] * (n + 1)
	for i in range(n):
		(x[i], y[i]) = lines[i].split()
		x[i] = int(x[i])
		y[i] = int(y[i])
	x[n] = x[0]
	y[n] = y[0]
	a = 0
	cx = 0
	cy = 0
	for i in range(n):
		ii = i + 1
		a += x[i] * y[ii] - x[ii] * y[i]
		cx += (x[i] + x[ii]) * (x[i] * y[ii] - x[ii] * y[i])
		cy += (y[i] + y[ii]) * (x[i] * y[ii] - x[ii] * y[i])
	a *= 3
	r = [0.0] * n
	z = [0.0] * n
	for i in range(n):
		x[i] = x[i] * a - cx
		y[i] = y[i] * a - cy
		r[i] = sqrt(x[i] * x[i] + y[i] * y[i]) / a
		z[i] = atan2(y[i], x[i])
	cx /= a
	cy /= a
	p = [0, 1]
	mt = 0.0
	for i in range(q):
		tmp = lines[n + i].split()
		ta = int(tmp[0])
		tb = int(tmp[1])
		if ta == 1:
			ta = int(tmp[2])
			ta -= 1
			tb -= 1
			if p[0] == tb:
				(p[0], p[1]) = (p[1], p[0])
			p[1] = ta
			tc = z[p[0]] + mt
			cx += r[p[0]] * cos(tc)
			cy += r[p[0]] * (sin(tc) - 1)
			tc = pi / 2 - tc
			mt = pi / 2 - z[p[0]]
		else:
			tb -= 1
			tc = z[tb] + mt
			stdout.write('%.10lf %.10lf\n' % (cx + r[tb] * cos(tc), cy + r[tb] * sin(tc)))
main()


import sys
read = lambda : map(int, sys.stdin.buffer.readline().split())
from math import *
(n, q) = read()
(x, y) = ([], [])
(cx, cy) = (0, 0)
deg = 0
(p1, p2) = (0, 1)

def pos(i):
	return (cx + x[i] * cos(deg) - y[i] * sin(deg), cy + y[i] * cos(deg) + x[i] * sin(deg))
for i in range(n):
	(_x, _y) = read()
	x.append(_x)
	y.append(_y)
s = 0
for i in range(2, n):
	_s = (x[i] - x[0]) * (y[i - 1] - y[0]) - (y[i] - y[0]) * (x[i - 1] - x[0])
	cx += (x[0] + x[i] + x[i - 1]) / 3 * _s
	cy += (y[0] + y[i] + y[i - 1]) / 3 * _s
	s += _s
cx /= s
cy /= s
for i in range(n):
	x[i] -= cx
	y[i] -= cy
for i in range(q):
	v = list(read())
	if v[0] == 1:
		if p1 == v[1] - 1:
			(p1, p2) = (p2, None)
		else:
			(p1, p2) = (p1, None)
		(cx, cy) = pos(p1)
		cy -= hypot(x[p1], y[p1])
		deg = pi / 2 - atan2(y[p1], x[p1])
		p2 = v[2] - 1
	else:
		print('%.15f %.15f' % pos(v[1] - 1))


from sys import stdin
from math import *
rl = lambda l: tuple(map(int, l.split()))
rd = lambda : rl(input())

class Point:

	def __init__(self, x, y):
		self.x = x
		self.y = y

	def __add__(self, pt):
		return Point(self.x + pt.x, self.y + pt.y)

	def __iadd__(self, pt):
		return self + pt

	def __mul__(self, n):
		return Point(self.x * n, self.y * n)

	def __imul__(self, n):
		return self * n

	def __truediv__(self, n):
		return Point(self.x / n, self.y / n)

	def __itruediv__(self, n):
		return self / n

	def __str__(self):
		return '%f %f' % (self.x, self.y)

	def distance(self, pt):
		return sqrt((self.x - pt.x) ** 2 + (self.y - pt.y) ** 2)

	def angle(self, pt):
		det_x = pt.x - self.x
		det_y = pt.y - self.y
		if det_x == 0:
			if det_y > 0:
				return pi / 2
			if det_y < 0:
				return pi / 2 * 3
			if det_y == 0:
				raise Exception('Fail to find angle between two identical points.')
		if det_x > 0:
			return atan(det_y / det_x)
		if det_x < 0:
			return atan(det_y / det_x) + pi

	def get_pt_by_dis_ang(self, dis, ang):
		return Point(self.x + cos(ang) * dis, self.y + sin(ang) * dis)

class Convex_polygon:

	def __init__(self):
		self.points = []
		self.core = None
		self.dis = None
		self.base_ang = 0
		self.ang = None
		self.pin = [0, 1]

	def add(self, pt):
		self.points.append(pt)

	def update(self):
		self.update_core()
		self.update_dis()
		self.update_ang()

	def cal_core_of_triangle(self, i, j, k):
		s = (self.points[j].x - self.points[i].x) * (self.points[k].y - self.points[i].y)
		s -= (self.points[j].y - self.points[i].y) * (self.points[k].x - self.points[i].x)
		return (s, (self.points[i] + self.points[j] + self.points[k]) / 3 * s)

	def update_core(self):
		self.core = Point(0, 0)
		s = 0
		for i in range(2, len(P.points)):
			(det_s, det_core) = self.cal_core_of_triangle(0, i, i - 1)
			self.core += det_core
			s += det_s
		self.core /= s

	def update_dis(self):
		self.dis = []
		for pt in self.points:
			self.dis.append(self.core.distance(pt))

	def update_ang(self):
		self.ang = []
		for pt in self.points:
			self.ang.append(self.core.angle(pt))

	def get_pt(self, i):
		return self.core.get_pt_by_dis_ang(self.dis[i], self.base_ang + self.ang[i])

	def __str__(self):
		return '\n'.join(('points: ' + '[' + ', '.join((str(pt) for pt in self.points)) + ']', 'core: ' + str(self.core), 'dis: ' + str(self.dis), 'base_ang ' + str(self.base_ang), 'ang ' + str(self.ang), 'pin ' + str(self.pin)))
(n, q) = rd()
P = Convex_polygon()
for _ in range(n):
	P.add(Point(*rd()))
P.update()
for _ in stdin.readlines():
	s = rl(_)
	if s[0] == 1:
		(__, f, t) = s
		f -= 1
		t -= 1
		P.pin.remove(f)
		i = P.pin[0]
		top_pt = P.get_pt(i)
		P.core = Point(top_pt.x, top_pt.y - P.dis[i])
		P.base_ang = pi / 2 - P.ang[i]
		P.pin.append(t)
	else:
		(__, v) = s
		v -= 1
		print(P.get_pt(v))


from math import hypot
[n, q] = map(int, input().strip().split())
xys = [tuple(map(int, input().strip().split())) for _ in range(n)]
qis = [tuple(map(int, input().strip().split())) for _ in range(q)]
dxys = [(xys[(i + 1) % n][0] - xys[i][0], xys[(i + 1) % n][1] - xys[i][1]) for i in range(n)]
S = 3 * sum((x * dy - y * dx for ((x, y), (dx, dy)) in zip(xys, dxys)))
Sx = sum(((dx + 2 * x) * (x * dy - y * dx) for ((x, y), (dx, dy)) in zip(xys, dxys)))
Sy = sum(((dy + 2 * y) * (x * dy - y * dx) for ((x, y), (dx, dy)) in zip(xys, dxys)))
for p in [2, 3]:
	while S % p == Sx % p == Sy % p == 0:
		S //= p
		Sx //= p
		Sy //= p
xyms = [(S * x - Sx, S * y - Sy) for (x, y) in xys]
hs = [hypot(x, y) for (x, y) in xyms]

def to_coord(x, y):
	return ((x + Sx) / S, (y + Sy) / S)
hangs = (0, 1)
hang_on = None
(cx, cy) = (0.0, 0.0)

def get_v(v):
	if hang_on is None:
		return xyms[v]
	else:
		(ux, uy) = xyms[hang_on]
		(vx, vy) = xyms[v]
		h = hs[hang_on]
		return ((uy * vx - ux * vy) / h, (ux * vx + uy * vy) / h)
for qi in qis:
	if qi[0] == 1:
		(_, f, t) = qi
		s = hangs[1 - hangs.index(f - 1)]
		(dx, dy) = get_v(s)
		cx += dx
		cy += dy - hs[s]
		hang_on = s
		hangs = (s, t - 1)
	else:
		(_, v) = qi
		(dx, dy) = get_v(v - 1)
		print(*to_coord(cx + dx, cy + dy))


from math import hypot, atan2, cos, sin, sqrt, pi
import sys
mapInt = lambda : map(int, sys.stdin.buffer.readline().split())
(N, Q) = mapInt()
points = {}
for n in range(N):
	(x, y) = mapInt()
	points[n + 1] = [x, y]
centroid = [0, 0]
A = 0
for n in range(1, N + 1):
	w = points[n][0] * points[n % N + 1][1] - points[n % N + 1][0] * points[n][1]
	centroid[0] += (points[n][0] + points[n % N + 1][0]) * w
	centroid[1] += (points[n][1] + points[n % N + 1][1]) * w
	A += w
centroid[0] /= 3 * A
centroid[1] /= 3 * A
for i in points.keys():
	points[i][0] -= centroid[0]
	points[i][1] -= centroid[1]
totAng = 0
pins = [1, 2]
for q in range(Q):
	inp = list(mapInt())
	if inp[0] == 1:
		(f, t) = (inp[1], inp[2])
		if pins[0] == pins[1] == f:
			pins[0] = t
			continue
		elif pins[0] == f:
			pivot = pins[1]
			pins[0] = t
		else:
			pivot = pins[0]
			pins[1] = t
		(x, y) = points[pivot]
		centroid[0] += sin(totAng + atan2(x, -y)) * hypot(x, y)
		centroid[1] += (-cos(totAng + atan2(x, -y)) - 1) * hypot(x, y)
		totAng = pi - atan2(x, -y)
	elif inp[0] == 2:
		(x, y) = points[inp[1]]
		pointX = sin(totAng + atan2(x, -y)) * hypot(x, y) + centroid[0]
		pointY = -cos(totAng + atan2(x, -y)) * hypot(x, y) + centroid[1]
		print('%0.20f %0.20f' % (pointX, pointY))