HDU4998 Rotate（旋转，线段相交模板）

Rotate

题目链接：

http://acm.hdu.edu.cn/showproblem.php?pid=4998

题目意思：

平面上有一个二维坐标轴上,进行n次操作，把坐标轴绕着(x,y) (这个坐标总是初始坐标轴的坐标) 逆时针转p弧度。最后的结果相当于进行一次操作，即绕着(X, Y) 逆时针旋转了P弧度。求 X,Y,P，题目保证总有解.

解题思路：

不难发现，最后的P是n次的p的和，因为这和绕什么点旋转无关。（在草稿纸上画一个图理解一下吧）

至于求最后C点的坐标，那就简单了。

此题相当于是一个模板题，模板：http://blog.csdn.net/piaocoder/article/details/39255451

对任意点x1,x2，绕一个坐标点v逆时针旋转p角度后的新的坐标设为t1,t2，x1=t1,x2=t2;

循环上述操作，得到最后的 x1,x2;令x3等于最开始的x1,x4等于最开始的x2;最后点C的坐标，即为直线x1x3的垂直平分线与直线x2x4的垂直平分线的交点。

AC代码：

#include <iostream>#include <cstdio>#include <cmath>using namespace std;////////////////////////////////////模板函数struct Point{    double x,y;    Point(double x=0,double y=0):x(x),y(y) {}};typedef Point Vector;const double pi=4*atan(1.0);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 p){    return Vector(A.x*p,A.y*p);}Vector operator / (Vector A,double p){    return Vector(A.x/p,A.y/p);}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){    return Vector(A.x*cos(rad)-A.y*sin(rad),A.x*sin(rad)+A.y*cos(rad));}Vector Normal(Vector A){    double L = Length(A);    return Vector(-A.y/L,A.x/L);}Point getlineintersection(Point P,Vector v,Point Q,Vector w){    Vector u = P-Q;    double t = Cross(w,u)/Cross(v,w);    return P+v*t;}////////////////////////////////////////int main(){    int T;    scanf("%d",&T);    while(T--){        int n,t;        Point v,t1,t2,t3,t4;        Vector a,b,tmp1,tmp2;        double p;        double x,y,p1=0;        t1.x=t1.y=t3.x=t3.y=0,t2.x=t2.y=t4.x=t4.y=1;        Point a1,a2,c;        scanf("%d",&n);        for(int i = 0; i < n; i++){            scanf("%lf%lf%lf",&v.x,&v.y,&p);            tmp1 = t1-v;            tmp2 = t2-v;            t1 = Rotate(tmp1,p)+v;            t2 = Rotate(tmp2,p)+v;            //printf("%.10lf %.10lf \n",t1.x,t1.y);            p1 += p;        }        t = p1/(2*pi);        p1 -= t*(2*pi);        a1 = (t1+t3)/2;a2=(t2+t4)/2;        a = t1-t3;b=t2-t4;        a = Normal(a);b = Normal(b);        c = getlineintersection(a1,a,a2,b);        printf("%.10lf %.10lf %.10lf\n",c.x,c.y,p1);    }    return 0;}