2013芝加哥邀请赛H题-Overlapping Maps-计算几何

题意：一个高H，宽W的地图，还有一个同样地图，按比例缩小，右下角在地图中的位置是$(x,y)$,长宽缩小到原来的s%，逆时针旋转r度

如图：

$AM=\frac{hs}{100} AN=\frac{ws}{100}$

$AB=\frac{hs}{100}*\frac{b}{h}=\frac{sb}{100}$

同理$AC=\frac{sa}{100}$

$\overrightarrow{AB}=(-\frac{sb}{100}sin(r),\frac{sb}{100}cos(r))$

$\overrightarrow{AC}=(\frac{sa}{100}cos(r),\frac{sa}{100}sin(r))$

$\overrightarrow{AP}=(-\frac{sb}{100}sin(r)+\frac{sa}{100}cos(r),\frac{sb}{100}cos(r)+\frac{sa}{100}sin(r))$

列方程组

$$\begin{eqnarray} \begin{cases} x-\frac{sb}{100}sin(r)+\frac{sa}{100}cos(r)=a \cr y+\frac{sb}{100}cos(r)+\frac{sa}{100}sin(r)=b \end{cases} \end{eqnarray}$$

求解出a，b即可。

#include<bits/stdc++.h>
#define pi acos(-1)
using namespace std;
int main()
{
    double w,h,x,y,s,r;
    while(scanf("%lf%lf%lf%lf%lf%lf",&w,&h,&x,&y,&s,&r)&&(w+h+x+y+s+r)!=0) {
        r=r/180.0*pi;
       // cout << r << endl;
        double x1=1.0-s/100*cos(r);
        double x2=s/100*sin(r);
        double x4=-s/100*sin(r);
        double x5=1-s/100*cos(r);
        double px=(x*x5-x2*y)/(x1*x5-x2*x4);
        double py=(x1*y-x*x4)/(x1*x5-x2*x4);
        printf("%.2f %.2f\n",px,py);
    }
    return 0;
}