Rotate (hdu4998) 几何

来源:互联网 发布:ubuntu命令行解压rar 编辑:程序博客网 时间:2024/06/03 21:54

题目:点击打开链接

题意:一共T组数据,给出n,代表再输入n行,每行x,y代表坐标,P代表这个坐标围绕这个点逆时针旋转P弧度,最后求出一个点和弧度能完成以上效果;

解题思路:其实旋转角度就是总旋转角度对2π取余,然后随意找两个启始点,求出按照题目意思旋转后的终点,这两个终点分别都是到两个启始点距离相等,然后解出方程,我就是以(0.0)(1.1)为启始点解的方程

#include <bits/stdc++.h>using namespace std;const double PI=acos(-1);int main(){    int T,n;    scanf("%d",&T);    while(T--)    {       scanf("%d",&n);       double r=0;       double tx1,ty1,tx2,ty2,x1,x2,y1,y2,a,rx0,ry0;       x1=0;y1=0;       x2=1;y2=1;       for(int i=1;i<=n;i++)       {           scanf("%lf%lf%lf",&rx0,&ry0,&a);           r+=a;           if(r>PI*2)            r-=PI*2;            tx1 = (x1 - rx0) * cos(a) - (y1 - ry0) * sin(a) + rx0;            ty1 = (x1 - rx0) * sin(a) + (y1 - ry0) * cos(a) + ry0;            tx2 = (x2 - rx0) * cos(a) - (y2 - ry0) * sin(a) + rx0;            ty2 = (x2 - rx0) * sin(a) + (y2 - ry0) * cos(a) + ry0;            x1=tx1;            x2=tx2;            y1=ty1;            y2=ty2;       }       double tmp1=(2*x2-2)*(-2*y1)/(2*x1);       tmp1=tmp1+(2*y2-2);       double tmp2=(x2*x2+y2*y2-2);       double tmp3=(2*x2-2)*(x1*x1+y1*y1)/(2*x1);       double y=(tmp2-tmp3)/tmp1;       double tmp4=x1*x1+y1*y1-2*y*y1;       double x=tmp4/(2*x1);       printf("%.6lf %.6lf %.6lf\n",x,y,r);    }    return 0;}


原创粉丝点击