zjuoj 3608 Signal Detection

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3608

Signal Detection

---

Time Limit: 2 Seconds      Memory Limit: 65536 KB

---

Parallelepiped
Type	Prism
Faces	6 parallelograms
Edges	12
Vertices	8

Dr. Gale is testing his laser system. He uses a detector to collect the signals from a laser generator. He also puts a special prism in the system so that he can filter the noise. But he is not very sure where to put the detector to collect the signals. Can you help him with this problem?

In order to simplify the problem, here we assume the prism is a parallelepiped, which is a three-dimensional figure formed by six parallelograms. The laser goes in one direction and the detector can receive signals from any direction. The detector is placed on the ground where the z-coordinate is zero. There is no energy lost in the refraction. That is to say, there is no reflection in the signal transmission. You don't need to consider the situation of total reflection or the degenerate situation.

Input

There are multiple test cases. The first line of input contains an integer T (T ≤ 50) indicating the number of test cases. Then T test cases follow.

The first line of each test case contains 3 integers, indicating the coordinates of the laser generator. The second line contains 3 integers describing a point the laser will go through without the prism. The third line contains 3 integers describing a vertex A of the prism. The fourth to the sixth lines contain 3 integers each, describing an adjacent vertex of A. The seventh line contains a real number u, the refractive index of the prism.(1 < u ≤ 10) The absolute value of the coordinates will not exceed 1000. The z coordinates are all nonnegative. The prism and the laser generator are strictly above the ground and the laser generator will not be inside the prism.

Output

If there is no place in the ground that can receive the signals output "Error". Otherwise, output the x and y coordinates the place accurate to 0.001.

Sample Input

2
0 0 10
0 0 0
-5 -5 1
5 -5 11
-5 5 1
-6 -5 2
1.4142136
0 0 10
0 0 11
-5 -5 1
5 -5 1
-5 5 1
-5 -5 2
2

Sample Output

0.423 0.000
Error

References

• http://en.wikipedia.org/wiki/Parallelepiped

---

Author: GUAN, Yao
Contest: The 9th Zhejiang Provincial Collegiate Programming Contest

AC代码：

 #include <cstdlib>
 #include <cstdio>
 #include <cstring>
 #include <cmath>
 #include <iostream>
 #include <algorithm>
 #include <string>
 #include <vector>
 #include <set>
 #include <map>
 #include <queue>
 #include <time.h>
 using namespace std;
 typedef long long llong;

 inline int bit(int x, int i){ return (x>>i) & ;}
 inline int setb(int x, int i){ return x | ( << i);}
 inline int clrb(int x, int i){ return x & (~( << i));}
 inline int lowb(int x){return x & (-x);}
 const double eps = 1e-;
 struct Point{
     double x, y, z;
     Point(){}
     Point(double x, double y, double z):x(x), y(y), z(z){}
     double len(){
         return sqrt(x * x + y * y + z * z);
     }
     Point shrink(double l = 1.0){
         double s = l / len();
         return Point(x * s, y * s, z * s);
     }
     void input(){
         scanf("%lf %lf %lf", &x, &y, &z);
     }
 };

 Point operator +(const Point &p, const Point &q){
     return Point(p.x + q.x, p.y + q.y, p.z + q.z);
 }
 Point operator -(const Point &p, const Point &q){
     return Point(p.x - q.x, p.y - q.y, p.z - q.z);
 }
 Point operator *(const Point &p, const double &s){
     return Point(p.x * s, p.y * s, p.z * s);
 }
 Point operator *(const Point &p, const Point &q){
     return Point(p.y * q.z - p.z * q.y,
                  p.z * q.x - p.x * q.z,
                  p.x * q.y - p.y * q.x);
 }
 double operator &(const Point &p, const Point &q){
     return p.x * q.x + p.y * q.y + p.z * q.z;
 }

 Point p, q, v[];
 double r;
 int face[][] ={
 {, , , },
 {, , , },
 {, , , },
 {, , , },
 {, , , },
 {, , , },
 };
 Point fn[];
 double fb[];

 double inter(const Point &p, const Point &d, const Point &n, double b){
     return (b - (n & p)) / (n & d);
 }
 bool collide(double R){
     int hf = -;
     double ht;
     Point d = (q - p);
     for(int i = ;i < ; ++i){
         if(fabs(fn[i] & d) < eps) continue;
         double t = inter(p, d, fn[i], fb[i]);
         if(t <= eps) continue;
         Point hit = p + d * t;
         bool ok = true;
         for(int j = ;j < ; ++j){
             if((((v[face[i][j]]-hit) * (v[face[i][(j + )%]] - hit)) & fn[i]) <= ) ok = false;
         }
         if(ok && (hf <  || t < ht)){
             hf = i;
             ht = t;
         }
     }
     if(hf < ) return false;

     Point hit = p + d * ht;
     Point vx = fn[hf].shrink();
     if(fabs((vx * d).len()) < eps){
         p = hit;
         q = hit + d;
         return true;
     }
     double dx = vx & d;
     if(dx < ) vx = vx * -, dx = -dx;
     Point vy = ((vx * d) * vx).shrink();
     double dy = vy & d;
     if(dy < ) vy = vy * -, dy = -dy;
     double theta0 = atan2(dy, dx);
     double theta = asin(sin(theta0) / R);
     d = vx * cos(theta) + vy * sin(theta);
     p = hit;
     q = hit + d;
     return true;
 }
 int main(){
     int TT;
     scanf("%d", &TT);
     for(int cas = ; cas <= TT; ++cas){
         p.input();
         q.input();

         v[].input();
         v[].input();
         v[].input();
         v[] = v[] + v[] - v[];

         v[].input();
         v[] = v[] + v[] - v[];
         v[] = v[] + v[] - v[];
         v[] = v[] + v[] - v[];

         scanf("%lf", &r);
         for(int i = ;i < ; ++i){
             fn[i] = (v[face[i][]] - v[face[i][]]) *
                     (v[face[i][]] - v[face[i][]]);
             fb[i] = fn[i] & v[face[i][]];
         }
         if(collide(r)){
             collide(./r);
         }
         Point norm(, , );
         Point dir = q - p;
         if(fabs(norm & dir) < eps) puts("Error");
         else{
             double t = inter(p, dir, norm, );
             if(t < ) puts("Error");
             else{
                 Point hit = p + dir * t;
                 printf("%.3f %.3f\n", hit.x, hit.y);
             }
         }
     }
     return ;
 }