NBUT 1462 台球碰撞（物理运动）

台球碰撞

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

110 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

Hint

无

主要考虑物理运动，对于初速度进行x方向和 y 方向 的分速度，这样就可以对于整个斜方向转换为x和y两个方向的分方向，然后针对x和y两个分方向的进行分析，对于x和y分开算，然后算其在x和y方向的碰撞往返，然后对其全长求余，那么就能得到其最终的位置，另外，在这题中需要考虑球的半径，我们其实只用研究球心的运动就行了，我们可以推出，球心就在内部的距离边为R以内的范围内运动

附上代码：

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<string>
#define INF 0x3f3f3f3f
#define PI 3.14159265
using namespace std;

int main()
{
	double L,W,x,y,R,a,v,s;
	while(cin >> L >> W >> x >> y >> R >> a >> v >> s && L + W)
	{
		double Lr = L - 2 * R;        //  球心运动的范围
		double Wr = W - 2 * R;
		a *= acos(0.0)/90.0;         //   将角度转化成弧度
		double vx = v * cos(a);
		double vy = v * sin(a);
		double Lx = s * vx + x - R;
		double Wy = s * vy + y - R;
		int dx = int(Lx / Lr);
		int dy = int(Wy / Wr);
		Lx = fmod(Lx,Lr);
		Wy = fmod(Wy,Wr);
		if(Lx < 0)
		Lx = Lx * -1;
		if(Wy < 0)
		Wy = Wy * -1;
		if(dy % 2)
		{
			Wy = Wr - Wy;
		}
		if(dx % 2)
		{
			Lx = Lr - Lx;
		}
		printf("%.2lf %.2lf\n",Lx + R,Wy + R);
	} 
	return 0;
}