学习心得(code)

superlong@CoreCoder

  C++博客 :: 首页 :: 联系 :: 聚合  :: 管理
  74 Posts :: 0 Stories :: 5 Comments :: 0 Trackbacks

公告

文字可能放在http://blog.csdn.net/superlong100,此处存放代码

常用链接

留言簿(4)

我参与的团队

搜索

  •  

最新随笔

最新评论

  • 1. re: Poj 1279
  • 对于一个凹多边形用叉积计算面积 后能根据结果的正负来判断给的点集的时针方向?
  • --bsshanghai
  • 2. re: Poj 3691
  • 你写的这个get_fail() 好像并是真正的get_fail,也是说fail指向的串并不是当前结点的子串。为什么要这样弄呢?
  • --acmer1183
  • 3. re: HDU2295[未登录]
  • 这个是IDA* 也就是迭代加深@ylfdrib
  • --superlong
  • 4. re: HDU2295
  • 评论内容较长,点击标题查看
  • --ylfdrib
  • 5. re: HOJ 11482
  • 呵呵..把代码发在这里很不错..以后我也试试...百度的编辑器太烂了....
  • --csuft1

阅读排行榜

评论排行榜

#include <stdio.h>
#include 
<string.h>
#include 
<math.h>

#define eps 1e-8
#define zero(x) (((x)>0?(x):-(x))<eps)

struct point{
    
double x, y, z;
    
void read() {
        scanf(
"%lf %lf %lf"&x, &y, &z);
    }

    
void out() {
        printf(
"%lf %lf %lf\n", x, y, z);
    }

}
;

point unit(point u) 
{
    point ret;
    
double t = sqrt(u.x*u.x+u.y*u.y+u.z*u.z);
    ret.x 
= u.x / t;
    ret.y 
= u.y / t;
    ret.z 
= u.z / t;
    
return ret;
}


point xmult(point u, point v) 
{
    point ret;
    ret.x 
= u.y * v.z - v.y * u.z;
    ret.y 
= u.z * v.x - u.x * v.z;
    ret.z 
= u.x * v.y - u.y * v.x;
    
return ret;
}


double vlen(point p) {
    
return sqrt(p.x*p.x+p.y*p.y+p.z*p.z);
}


point subt(point u, point v) 
{
    point ret;
    ret.x 
= u.x - v.x;
    ret.y 
= u.y - v.y;
    ret.z 
= u.z - v.z;
    
return ret;
}


point pvec(point s1, point s2, point s3) 
{
    
return xmult(subt(s1, s2), subt(s2, s3));
}


double dmult(point u, point v) {
    
return u.x * v.x + u.y * v.y + u.z * v.z;
}


int dots_onplane(point a, point b, point c, point d) {
    
return zero( dmult( pvec(a, b, c), subt(d, a)));
}


int same_side(point p1, point p2, point s1, point s2, point s3) {
    
return dmult(pvec(s1, s2, s3), subt(p1, s1)) * dmult(pvec(s1, s2, s3), subt(p2, s1)) > eps;
}


int opposite_side(point p1, point p2, point s1, point s2, point s3) {
    
return dmult(pvec(s1, s2, s3), subt(p1, s1)) * dmult(pvec(s1, s2 ,s3), subt(p2, s1)) < -eps;
}


int intersect_ex(point l1, point l2, point s1, point s2, point s3) {
    
return opposite_side(l1, l2, s1, s2, s3) && opposite_side(s1, s2, l1, l2, s3) &&
           opposite_side(s2, s3, l1, l2, s1) 
&& opposite_side(s3, s1, l1, l2, s2);
}


int dot_inplane_in(point p, point s1, point s2, point s3) {
    
return zero(vlen(xmult(subt(s1, s2),subt(s1, s3)))-vlen(xmult(subt(p,s1),subt(p,s2)))-
                vlen(xmult(subt(p, s2), subt(p, s3)))
-vlen(xmult(subt(p,s3), subt(p,s1))));
}


point ini, dic, goal[
8], end, mid;

int main() {
    
int cases;
    scanf(
"%d"&cases);
    
for(int t = 1; t <= cases; t ++{
        ini.read();
        dic.read();
        
for(int i = 0; i < 8; i ++)
            goal[i].read();
        
bool flag = false;
        
if( opposite_side(ini, goal[7], goal[0], goal[1], goal[2]) > 0 ||
            dots_onplane(ini, goal[
0], goal[1], goal[2])  ) {
            
            end.x 
= ini.x + dic.x * 10000000;
            end.y 
= ini.y + dic.y * 10000000;
            end.z 
= ini.z + dic.z * 10000000;
            
            mid.x 
= (goal[0].x + goal[2].x) / 2;
            mid.y 
= (goal[0].y + goal[2].y) / 2;
            mid.z 
= (goal[0].z + goal[2].z) / 2;
            
            point tt, ts;
            tt.x 
= mid.x - ini.x; tt.y = mid.y - ini.y; tt.z = mid.z - ini.z;
            tt 
= unit(tt);
            ts 
= unit(dic);
            flag 
= false;
            
for(int i = 0; i < 4; i ++{
                
if( intersect_ex(ini, end, goal[i], goal[(i+1)%4],goal[(i+2)%4]) ) {
                        flag 
= true;
                }

            }

            
if!flag && !dots_onplane(ini, goal[0], goal[1], goal[2]) &&
                ts.x 
== tt.x && ts.y == tt.y && ts.z == tt.z )     flag = true;
            
            
if!flag && opposite_side(ini, end, goal[4], goal[5], goal[6]) ) {
                
for(int i = 0; i < 4; i ++{
                    
if( dot_inplane_in(ini, goal[i], goal[(i+1)%4],goal[(i+2)%4]) )
                        flag 
= true;
                }

            }

        }

        printf(
"Case %d: ", t);
        
if(flag) {
            puts(
"Stupid Larrionda!!!");
        }
 else {
            puts(
"Intelligent Larrionda!!!");
        }

    }
    
    
return 0;
}

posted on 2010-07-22 17:34 superlong 阅读(165) 评论(0)  编辑 收藏 引用

只有注册用户登录后才能发表评论。
网站导航: 博客园   IT新闻   BlogJava   知识库   博问   管理