zcmu-1756 台球碰撞（正交分解）

Problem D: 台球碰撞

Time Limit: 1 Sec  Memory Limit: 128 MB
Submit: 91  Solved: 37
[Submit][Status][Web Board]

Description

在平面直角坐标系下，台球桌是一个左下角在(0,0)，右上角在(L,W)的矩形。有一个球心在(x,y)，半径为R的圆形母球放在台球桌上（整个球都在台球桌内）。受撞击后，球沿极角为a的射线（即：x正半轴逆时针旋转到此射线的角度为a）飞出，每次碰到球桌时均发生完全弹性碰撞（球的速率不变，反射角等于入射角）。 如果球的速率为v，s个时间单位之后球心在什么地方？

Input

输入文件最多包含25组测试数据，每个数据仅一行，包含8个正整数L,W,x,y,R,a,v,s（100<=L,W<=105, 1<=R<=5, R<=x<=L-R, R<=y<=W-R, 0<=a<360, 1<=v,s<=105），含义见题目描述。L=W=x=y=R=a=v=s=0表示输入结束，你的程序不应当处理这一行。

Output

对于每组数据，输出仅一行，包含两个实数x, y，表明球心坐标为(x,y)。x和y应四舍五入保留两位小数

Sample Input

100 100 80 10 5 90 2 23 1

10 100 70 10 5 180 1 9999

0 0 0 0 0 0 0 0

Sample Output

80.00 56.00

71.00 10.00

//将球的碰撞看做点的碰撞，则由相切得，矩形大小变为x的范围为[r,L-r]，y的范围为[r,W-r],然后考虑将速度正交分解，；

//%.2f会自动四舍五入

//100 100 80 10 5 270 2 23
#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cmath>
using namespace std;
const double PI=acos(-1.0);
int main()
{
    double l,w,x,y,r,a,v,s;
    while(scanf("%lf%lf%lf%lf%lf%lf%lf%lf",&l,&w,&x,&y,&r,&a,&v,&s))
    {
        if(!l&&!w&&!x&&!y&&!r&&!a&&!v&&!s)
            break;
        int flagx,flagy;//速度方向
        l=l-r;
        w=w-r;
        double vx,vy;
        if(a>=0&&a<=90)flagx=1,flagy=1;
        else if(a>=90&&a<=180)flagx=-1,flagy=1,a=180-a;
        else if(a>=180&&a<=270)flagx=-1,flagy=-1,a=a-180;
        else flagx=1,flagy=-1,a=360-a;
        if(a==90)vx=0,vy=v;
        else if(a==0)vx=v,vy=0;
        else
        {
            vx=v*cos(a/180*PI);
            vy=v*sin(a/180*PI);
        }
        double x1,y1;
        double time=0;
        double disx,disy;
        while(1)
        {
            if(flagx==1)disx=l-x;
            else disx=x-r;
            if(flagy==1)disy=w-y;
            else disy=y-r;
            if(disx/vx>disy/vy||vx==0)
            {
                double t=disy/vy;
                //cout<<disy<<"***"<<vy<<endl;
                if(time+t<s)
                {
                    x1=x+flagx*vx*t;
                    y1=y+flagy*vy*t;
                    time+=t;
                    x=x1;y=y1;flagy=-flagy;
                }else if(time+t==s)
                {
                    x1=x+flagx*vx*t;
                    y1=y+flagy*vy*t;
                    break;
                }else
                {
                    t=s-time;
                    x1=x+flagx*vx*t;
                    y1=y+flagy*vy*t;
                    break;
                }
            }else if(disx/vx<=disy/vy||vy==0)
            {
                double t=disx/vx;
                if(time+t<s)
                {
                    x1=x+flagx*vx*t;
                    y1=y+flagy*vy*t;
                    time+=t;
                    x=x1;y=y1;flagx=-flagx;
                }else if(time+t==s)
                {
                    x1=x+flagx*vx*t;
                    y1=y+flagy*vy*t;
                    break;
                }else
                {
                    t=s-time;
                    x1=x+flagx*vx*t;
                    y1=y+flagy*vy*t;
                    break;
                }
            }
            //cout<<x1<<" "<<y1<<time<<endl;
        }
        printf("%.2f %.2f\n",x1,y1);
    }
    return 0;
}