Description
有一个w*h的矩形,它两边与坐标轴平行,且矩形中心与坐标原点重合。
现将它沿原点逆时针旋转a度产生一个新的矩形,求这个新矩形与原矩形相交的面积。
详见对应样例2的图片。
Input
第一行三个整数w,h,a。
Output
一个实数表示两个矩形相交的面积。
Solutions
暴力出旋转之后的矩形与原矩形所有交点,
分两种情况:8个交点和4个交点,
特判一下90°和180°就OK了。
代码
1 uses math; 2 var 3 w,h,aa:int64; 4 x1,y1,x2,y2,k,b,op,po:array [1..4] of real; 5 function min(o,p:longint):longint; 6 begin 7 if o<p then exit(o); 8 exit(p); 9 end; 10 11 procedure main; 12 var 13 i:longint; 14 num,top1,top2:int64; 15 a,ty,s,xx,yy,t,tk,s1,s2,xo,yo:real; 16 begin 17 ty:=aa; 18 a:=aa*3.1415926535/180; num:=0; top1:=0; top2:=0; 19 x1[1]:=-w/2; x1[2]:=-w/2; x1[3]:=w/2; x1[4]:=w/2; 20 y1[1]:=h/2; y1[2]:=-h/2; y1[3]:=-h/2; y1[4]:=h/2; 21 xx:=-w/2; yy:=h/2; 22 if a>90 then a:=(180-ty)*3.1415926535/180; 23 k[1]:=tan(a); k[2]:=-1/tan(a); 24 k[3]:=tan(a); k[4]:=-1/tan(a); 25 ty:=ty*3.1415926535/180; xo:=0; yo:=0; 26 for i:=1 to 4 do 27 begin 28 x2[i]:=x1[i]*cos(ty)-y1[i]*sin(ty); 29 y2[i]:=x1[i]*sin(ty)+y1[i]*cos(ty); 30 b[i]:=y2[i]-k[i]*x2[i]; 31 if x2[i]<=x1[1] then xo:=x1[1]-x2[i]; 32 if y2[i]>=y1[1] then yo:=y2[i]-y1[1]; 33 if ((yy-b[i])/k[i]>=xx) and ((yy-b[i])/k[i]<=xx+w) then 34 begin 35 inc(num); 36 inc(top2); po[top2]:=(yy-b[i])/k[i]; 37 end; 38 if ((yy-b[i]-h)/k[i]>=xx) and ((yy-b[i]-h)/k[i]<=xx+w) then inc(num); 39 if (k[i]*xx+b[i]<=yy) and (k[i]*xx+b[i]>=yy-h) then 40 begin 41 inc(num); 42 inc(top1); op[top1]:=(k[i]*xx+b[i]); 43 end; 44 if (k[i]*(xx+w)+b[i]<=yy) and (k[i]*(xx+w)+b[i]>=yy-h) then inc(num); 45 end; 46 if num=8 then 47 begin 48 if op[1]>op[2] then t:=op[1]-op[2] 49 else t:=op[2]-op[1]; 50 if po[1]>po[2] then tk:=po[1]-po[2] 51 else tk:=po[2]-po[1]; 52 s1:=t*xo; s2:=tk*yo; 53 s:=w*h-s1-s2; 54 end else 55 if top1>0 then 56 begin 57 if op[1]>op[2] then t:=op[1]-op[2] 58 else t:=op[2]-op[1]; 59 s:=t*w; 60 end else 61 begin 62 if po[1]>po[2] then t:=po[1]-po[2] 63 else t:=po[2]-po[1]; 64 s:=t*h; 65 end; 66 writeln(s:0:9); 67 end; 68 69 begin 70 readln(w,h,aa); 71 if aa mod 90=0 then 72 begin 73 if aa=90 then writeln(min(w,h)*min(w,h)) 74 else writeln(w*h); 75 end else 76 main; 77 end.