hdu4998:http://acm.hdu.edu.cn/showproblem.php?pid=4998
题意：给你n个点，以及绕每个点旋转的弧度。然后，问你经过这n次旋转，平面中的点总的效果是相当于哪个点旋转了多少弧度。
题解：我的第一道计算几何。可以选两个点，求出旋转之后的对应点，然后分别求出这两个点的中垂线，中垂线的交点就是要求的点，弧度就是所有弧度之和mod（2*pi）
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 #include<cmath> 6 const double pi=3.14159265; 7 using namespace std; 8 int n; 9 struct Point{ 10 double x; 11 double y; 12 double s; 13 }num[20]; 14 struct LINE{ 15 double a,b,c; 16 Point s; 17 }; 18 Point HHH(Point o,double alpha,Point p) {//点p绕着o旋转alpha弧度 19 Point tp; 20 p.x-=o.x; 21 p.y-=o.y; 22 tp.x=p.x*cos(alpha)-p.y*sin(alpha)+o.x; 23 tp.y=p.y*cos(alpha)+p.x*sin(alpha)+o.y; 24 return tp; 25 } 26 bool CCC(LINE l1,LINE l2,Point&p){//两条直线是否相交，如果相交，则交点为p 27 double d=l1.a*l2.b-l2.a*l1.b; 28 if(abs(d)<1e-6) return false; 29 p.x = (l2.c*l1.b-l1.c*l2.b)/d; 30 p.y = (l2.a*l1.c-l1.a*l2.c)/d; 31 return true; 32 } 33 LINE DDD(const Point &_a, const Point &_b){//求两点之间垂直平分线 34 LINE ret; 35 ret.s.x = (_a.x + _b.x)/2; 36 ret.s.y = (_a.y + _b.y)/2; 37 ret.a = _b.x - _a.x; 38 ret.b = _b.y - _a.y; 39 ret.c = (_a.y - _b.y) * ret.s.y + (_a.x - _b.x) * ret.s.x; 40 return ret; 41 } 42 43 LINE EEE(Point a,Point b){//求经过两点的中垂线 44 LINE ret; 45 if(abs(a.x-b.x)<1e-6){ 46 ret.a=1; 47 ret.b=0; 48 ret.c=-a.x; 49 return ret; 50 } 51 ret.a=a.y-b.y; 52 ret.b=b.x-a.x; 53 ret.c=a.x*b.y-a.y*b.x; 54 return ret; 55 } 56 int main(){ 57 int T; 58 scanf("%d",&T); 59 while(T--){ 60 scanf("%d",&n); 61 double sum=0; 62 for(int i=1;i<=n;i++){ 63 scanf("%lf%lf%lf",&num[i].x,&num[i].y,&num[i].s); 64 sum+=num[i].s; 65 } 66 Point s1,s2; 67 s1.x=2.345,s1.y=4.123; 68 s2.x=11.345,s2.y=12.123; 69 Point t1=s1,t2=s2; 70 for(int i=1;i<=n;i++){ 71 Point temp=HHH(num[i],num[i].s,t1); 72 t1=temp; 73 } 74 for(int i=1;i<=n;i++){ 75 Point temp=HHH(num[i],num[i].s,t2); 76 t2=temp; 77 } 78 LINE one1=DDD(s1,t1); 79 LINE one2=DDD(s2,t2); 80 Point a; 81 if(CCC(one1,one2,a)){ 82 if(sum>2*pi) 83 sum-=2*pi; 84 printf("%lf %lf %lf ",a.x,a.y,sum); 85 } 86 else{ 87 LINE ttt=EEE(s1,s2); 88 CCC(one1,ttt,a); 89 if(sum>2*pi) 90 sum-=2*pi; 91 printf("%lf %lf %lf ",a.x,a.y,sum); 92 } 93 } 94 }