• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>

            學(xué)習(xí)心得(code)

            superlong@CoreCoder

              C++博客 :: 首頁 :: 聯(lián)系 :: 聚合  :: 管理
              74 Posts :: 0 Stories :: 5 Comments :: 0 Trackbacks

            公告

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

            常用鏈接

            留言簿(4)

            我參與的團隊

            搜索

            •  

            最新隨筆

            最新評論

            • 1.?re: Poj 1279
            • 對于一個凹多邊形用叉積計算面積 后能根據(jù)結(jié)果的正負來判斷給的點集的時針方向?
            • --bsshanghai
            • 2.?re: Poj 3691
            • 你寫的這個get_fail() 好像并是真正的get_fail,也是說fail指向的串并不是當(dāng)前結(jié)點的子串。為什么要這樣弄呢?
            • --acmer1183
            • 3.?re: HDU2295[未登錄]
            • 這個是IDA* 也就是迭代加深@ylfdrib
            • --superlong
            • 4.?re: HDU2295
            • 評論內(nèi)容較長,點擊標題查看
            • --ylfdrib
            • 5.?re: HOJ 11482
            • 呵呵..把代碼發(fā)在這里很不錯..以后我也試試...百度的編輯器太爛了....
            • --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 閱讀(168) 評論(0)  編輯 收藏 引用
            韩国三级中文字幕hd久久精品| 四虎国产精品免费久久久| 一本伊大人香蕉久久网手机| 超级碰久久免费公开视频| 三级韩国一区久久二区综合| 久久人人爽人人爽人人av东京热| 7777久久久国产精品消防器材| 97久久精品国产精品青草| 成人亚洲欧美久久久久| 日韩久久久久中文字幕人妻| 久久久久久毛片免费播放| 国产成人综合久久综合| 色99久久久久高潮综合影院| 97热久久免费频精品99| 亚洲第一永久AV网站久久精品男人的天堂AV | www亚洲欲色成人久久精品| 国产激情久久久久影院| 亚洲国产精品无码久久久不卡| 久久超乳爆乳中文字幕| 色综合久久夜色精品国产| 久久免费国产精品一区二区| 久久99国产精品久久99小说| 青草影院天堂男人久久| 久久水蜜桃亚洲av无码精品麻豆| 久久精品国产精品亚洲下载| 91久久婷婷国产综合精品青草| 四虎国产精品成人免费久久| 久久久精品日本一区二区三区| 国产一级做a爰片久久毛片| 少妇久久久久久久久久| 亚洲欧美一区二区三区久久| 中文字幕一区二区三区久久网站| 久久久久人妻一区二区三区vr| 久久福利资源国产精品999| 久久午夜综合久久| 久久99精品久久久久久齐齐| 午夜精品久久久久9999高清| 亚洲AⅤ优女AV综合久久久| 婷婷久久综合九色综合绿巨人| 久久久久免费视频| 婷婷久久综合九色综合绿巨人|