hdu4998 Rotate（计算几何+模拟）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4998

题目大意：

n次旋转，每一次以（x,y）为中心点旋转弧度p。问最终旋转的效果相当于是围绕哪个点旋转弧度P？

输出这个点坐标以及旋转的弧度P

解：

这题就是模拟做的，因为数据量不大。模拟做就是写起来麻烦。

因为用的旋转模板是逆时针旋转，但最后求的P是在三角形里求解出来的，这样的P一定小于π，如果旋转角大于π需要判定下。

（和队友各种整才整过，菜~~）

代码：

#include <bits/stdc++.h>const double pi = acos (-1.0);const double eps=1e-6;using namespace std;struct Point{    double x,y,b;    Point(){}    Point(double _x,double _y)    {        x=_x;y=_y;    }    Point operator -(const Point &b)const    {        return Point(x-b.x,y-b.y);    }    double operator ^(const Point &b)const    {        return x*b.y-y*b.x;    }    void trans(double B)    {        double tx=x,ty=y;        x=tx*cos(B)-ty*sin(B);        y=tx*sin(B)+ty*cos(B);    }}a[15];void zhuan(Point &A,Point &B,double b){    A.x-=B.x;    A.y-=B.y;    A.trans(b);    A.x+=B.x;    A.y+=B.y;}Point AA(101,100),BB(101,1),CC(0,101);struct Line{    Point s,e;    Line(){}    Line(Point _s,Point _e)    {        s=_s;e=_e;    }    pair<int,Point> operator &(const Line &b)const{    Point res=s;    double t=((s-b.s)^(b.s-b.e))/((s-e)^(b.s-b.e));    res.x+=(e.x-s.x)*t;    res.y+=(e.y-s.y)*t;    return make_pair(2,res);    }};double dist(Point A,Point B){    return sqrt((A.x-B.x)*(A.x-B.x)+(A.y-B.y)*(A.y-B.y));}double qiuhudu(Point C,Point ans1,Point CC){    double a=dist(C,ans1);    double b=dist(CC,ans1);    double c=dist(C,CC);    double ans=(a*a+b*b-c*c)/2/a/b;    return acos(ans);}int main(){     int t,n;     scanf("%d",&t);     while(t--)     {         scanf("%d",&n);         Point A(101,100),B(101,1),C(0,101);         for(int i=0;i<n;i++)         {             scanf("%lf%lf%lf",&a[i].x,&a[i].y,&a[i].b);             zhuan(A,a[i],a[i].b);             zhuan(B,a[i],a[i].b);             zhuan(C,a[i],a[i].b);         }         Point mid1((A.x+AA.x)/2,(A.y+AA.y)/2);         Point mid2((B.x+BB.x)/2,(B.y+BB.y)/2);         zhuan(A,mid1,pi/2);         Line A1(A,mid1);         zhuan(B,mid2,pi/2);         Line A2(B,mid2);        pair<int,Point>s=A1&A2;        Point ans1(s.second.x,s.second.y);        double ans2=qiuhudu(C,ans1,CC);//cout<<"qiude ="<<ans2<<" ";        //cout<<ans2<<endl;        Point temp(CC.x,CC.y);        zhuan(temp,ans1,ans2);        if(fabs(temp.x-C.x)>eps||fabs(temp.y-C.y)>eps)            ans2=2*pi-ans2;            /*cout<<"CC"<<CC.x<<" "<<CC.y<<endl;            cout<<"C"<<C.x<<" "<<C.y<<endl;            cout<<"temp"<<temp.x<<" "<<temp.y<<endl;*/        printf("%.10f %.10f %.10f\n",s.second.x,s.second.y,ans2);     }    return 0;}