Codeforce 975E

E. Hag's Khashba
time limit per test
3 seconds
memory limit per test
256 megabytes
standard input
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 nn vertices.

Hag brought two pins and pinned the polygon with them in the 11-st and 22-nd vertices to the wall. His dad has qq queries to Hag of two types.

  • 11 ff tt: pull a pin from the vertex ff, wait for the wooden polygon to rotate under the gravity force (if it will rotate) and stabilize. And then put the pin in vertex tt.
  • 22 vv: answer what are the coordinates of the vertex vv.

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.


The first line contains two integers nn and qq (3n100003≤n≤100001q2000001≤q≤200000) — the number of vertices in the polygon and the number of queries.

The next nn lines describe the wooden polygon, the ii-th line contains two integers xixi and yiyi (|xi|,|yi|108|xi|,|yi|≤108) — the coordinates of the ii-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 qq lines describe the queries, one per line. Each query starts with its type 11 or 22. Each query of the first type continues with two integers ff and tt (1f,tn1≤f,t≤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 ff contains a pin. Each query of the second type continues with a single integer vv (1vn1≤v≤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.


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 10410−4.

Formally, let your answer be aa, and the jury's answer be bb. Your answer is considered correct if |ab|max(1,|b|)104|a−b|max(1,|b|)≤10−4

3 4
0 0
2 0
2 2
1 1 2
2 1
2 2
2 3
3.4142135624 -1.4142135624
2.0000000000 0.0000000000
0.5857864376 -1.4142135624
3 2
-1 1
0 0
1 1
1 1 2
2 1
1.0000000000 -1.0000000000

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)(2,0).

In the second sample note that the polygon rotates 180180 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.





#define LL long long
#define N 10010
const double eps=1e-3;
using namespace std;
int dcmp(double x){if (fabs(x)<eps)return 0;else return x<0?-1:1;}
struct Point
    double x,y;
    Point(double x=0,double y=0):x(x),y(y){ }
    void read(){scanf("%lf%lf",&x,&y);}
typedef Point Vector;
Point a[N];
int g[N];

Vector operator + (Vector a,Vector b){return Vector(a.x+b.x,a.y+b.y);}
Vector operator - (Vector a,Vector b){return Vector(a.x-b.x,a.y-b.y);}
Vector operator * (Vector a,double b){return Vector(a.x*b,a.y*b);}
Vector operator / (Vector a,double b){return Vector(a.x/b,a.y/b);}
double Dot(Vector a,Vector b){return a.x*b.x+a.y*b.y;}  //点积
double Length(Vector a){return sqrt(Dot(a,a));}
double Cross(Vector a,Vector b){return a.x*b.y-a.y*b.x;} //叉积
Vector Rotate(Vector a,double rad)// 向量 a 逆时针旋转 rad
{return Vector(a.x*cos(rad)-a.y*sin(rad),a.x*sin(rad)+a.y*cos(rad));}
Point Zhongxin(Point *a,int n)
    Point ans=Point(0,0);double ss=0;
    for (int i=1;i<n-1;i++)
        double s=fabs(Cross(a[i]-a[0],a[i+1]-a[0]));
    if (dcmp(ss)==0) return Point(0,0);
    return ans/ss;

double GetAngle(Vector a)
    return atan2(a.x,a.y);

int main()
    int n,m;double ang=0; g[0]=g[1]=1; int a1=0,a2=1;
    for (int i=0;i<n;i++) a[i].read();
    Point cc=Zhongxin(a,n);
    for (int i=0;i<n;i++) a[i]=a[i]-cc;
    for (int i=1;i<=m;i++)
        int fg,j,k;
        if (fg&1)
            if (!g[j])
                int k=j==a1?a2:a1;
                Point t=cc+Rotate(a[k],ang);
                double r=GetAngle(t-cc);
            if (j==a1) a1=k;else a2=k;
            Point t=cc+Rotate(a[j-1],ang);
            double x=t.x,y=t.y;
            printf("%.8f %.8f\n",x,y);

文章标签: Codeforces 计算几何
个人分类: 计算几何 Codeforces
上一篇HDU 6230 CCPC 2017 哈尔滨
下一篇CSU 2212 Particle Collision
想对作者说点什么? 我来说一句