南邮 OJ 1180 H - 管道小球

H - 管道小球

时间限制(普通/Java) : 1000 MS/ 3000 MS          运行内存限制 : 65536 KByte
总提交 : 64            测试通过 : 18 

比赛描述

    有一个截面为长方形的管道，不放假设在笛卡尔坐标系中，该管道的中轴线恰好为x轴，管道的入口位于平面x=0，管道出口位于x=X。

    管道界面的长为a，宽为b，且长边平行于y轴，宽边平行于z轴。

    假设在管道的入口处（0，0，0）有一个半径可忽略并且不计重力的小球，现在给小球一个初速度v=（vx，vy，vz）使得小球运动。假设当小球与管道壁相碰，小球会完全的反弹。请问小球能否到达出口截面？截面坐标是什么？

输入

一行，包含6个整数：a，b，X，vx，vy，vz（0<b<a<10000,0<=x<=10000,-10000<=vx,vy,vz<=10000）

输出

若小球能到达管道出口截面，则输出一个实数：X，Y，Z，用单个空格隔开，表示到达是的坐标，否则，输出“impossible”，所有结果保留3位小数。

样例输入

10 10 10 10 0 0

样例输出

10.000 0.000 0.000

提示

undefined

题目来源

NJU 7th ACM contest

#include<iostream>
int main(){
// freopen("test.txt","r",stdin);
double a,b,X,vx,vy,vz,t,disy,disz;
scanf("%lf%lf%lf%lf%lf%lf",&a,&b,&X,&vx,&vy,&vz);
if(vx==0){
if(X==0){
printf("%.3lf %.3lf %.3lf\n",X,0,0);
}else{
printf("impossible\n");
}
return 0;
}
t = X/vx;
if(t<0){
printf("impossible\n");
return 0;
}
disy = vy*t;
while(disy<-a/2){
disy += 2*a;
}
while(disy>3*a/2){
disy -= 2*a;
}
if(disy>a/2){
disy = a-disy;
}

disz = vz*t;
while(disz<-b/2){
disz += 2*b;
}
while(disz>3*b/2){
disz -= 2*b;
}
if(disz>b/2){
disz = b-disz;
}
printf("%.3lf %.3lf %.3lf\n",X,disy,disz);
}