• <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>
            posts - 100,  comments - 15,  trackbacks - 0


            一.叉積

               設(shè) a(x1,y1), b(x2,y2)
            二維:a x b=x1*y2-x2*y1

               設(shè)p0(x0,y0), p1(x1,y1) ,p2(x2,y2), p3(x3,y3)
              < p0p1>= (p1-p0) =(x1-x0,y1-y0);
               <p2p3>= (p3-p2) =(x3-x2,y3-y2);
               <p0p1> x <p2p3> =(p1-p0) x (p3-p2) = (x1-x0)*(y3-y2)- (x3-x2)* (y1-y0); 

               (p1-p0) x (p3-p2)結(jié)果的意義:

               : <p0p1><p2,p3>順時(shí)針(0,pi)內(nèi)

               負(fù): <p0p1><p2,p3>逆時(shí)針(0,pi)內(nèi)

               0 : <p0p1>, <p2p3>共線,夾角為0pi

             

             

            double xmult(point p0,point p1,point p2){
                
            return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
            }

               可用于判p2與<p0 p1>的關(guān)系

             

            double xmult(point p0,point p1,point p2,point p3){
                
            return (p1.x-p0.x)*(p3.y-p2.y)-(p3.x-p2.x)*(p1.y-p0.y);
            }

               可用于判<p0p1>, <p2p3>位置關(guān)系

            double xmult(line u,line v){
                
            return (u.b.x-u.a.x)*(v.b.y-v.a.y)-(v.b.x-v.a.x)*(u.b.y-u.a.y);
            }


             

            //計(jì)算兩直線交點(diǎn),注意事先判斷直線是否平行!
            point intersection(line u,line v){
                point ret
            =u.a;
                
            double t=((u.a.x-v.a.x)*(v.a.y-v.b.y)-(u.a.y-v.a.y)*(v.a.x-v.b.x))
                        
            /((u.a.x-u.b.x)*(v.a.y-v.b.y)-(u.a.y-u.b.y)*(v.a.x-v.b.x));
                ret.x
            +=(u.b.x-u.a.x)*t;
                ret.y
            +=(u.b.y-u.a.y)*t;
                
            return ret;
            }



             

            int relat_and_intersection(line u,line v,point &ret)
            {
                
            /*
                LINE 重合
                NONE 平行
                INTERSECT 相交并返回交點(diǎn)
                
            */

                
            double a1,b1,c1,a2,b2,c2;
                
            //系數(shù)
                a1=u.a.y-u.b.y;
                b1
            =-(u.a.x-u.b.x);
                c1
            =u.a.x*u.b.y-u.a.y*u.b.x;

                a2
            =v.a.y-v.b.y;
                b2
            =-(v.a.x-v.b.x);
                c2
            =v.a.x*v.b.y-v.a.y*v.b.x;

                
            if(fabs(a1*b2-a2*b1)<eps)
                
            {
                    
            if( fabs(a1*c2-a2*c1)<eps && fabs(b1*c2-b2*c1)<eps )
                        
            return LINE;
                    
            else return NONE;
                }

                
            else 
                
            {
                    point ret;
                    ret.x
            =(b1*c2-b2*c1)/(a1*b2-a2*b1);
                    ret.y
            =(a2*c1-a1*c2)/(a1*b2-a2*b1);
                    
            return INTERSECT;
                }

            }
            posted on 2009-10-04 18:38 wyiu 閱讀(257) 評(píng)論(0)  編輯 收藏 引用 所屬分類(lèi): 常用模板和函數(shù)

            只有注冊(cè)用戶登錄后才能發(fā)表評(píng)論。
            相關(guān)文章:
            網(wǎng)站導(dǎo)航: 博客園   IT新聞   BlogJava   博問(wèn)   Chat2DB   管理


            久久国产精品久久久| 久久精品人妻中文系列| 国产精品久久自在自线观看| www.久久99| 无码精品久久一区二区三区| 精品无码久久久久国产动漫3d| 久久这里只有精品首页| 99久久精品国产麻豆| 日本精品一区二区久久久| 韩国免费A级毛片久久| 久久精品无码一区二区三区日韩| 久久香综合精品久久伊人| 九九99精品久久久久久| 久久人妻少妇嫩草AV蜜桃| 伊人久久大香线焦综合四虎| 伊人久久大香线蕉AV色婷婷色| 久久精品国产亚洲综合色| 久久综合久久自在自线精品自| 少妇被又大又粗又爽毛片久久黑人| 婷婷久久综合九色综合九七| 99久久精品费精品国产一区二区| 欧美精品九九99久久在观看| 亚洲精品99久久久久中文字幕 | 色婷婷久久综合中文久久蜜桃av| 亚洲嫩草影院久久精品| 久久天堂AV综合合色蜜桃网| 久久精品国产亚洲AV大全| 久久精品中文无码资源站| 久久精品夜夜夜夜夜久久| 精品久久久久中文字| 久久精品国产清高在天天线| 久久久高清免费视频| 99久久综合狠狠综合久久止| 国产精品久久久久一区二区三区| 性做久久久久久久久浪潮| 精品久久久久中文字幕一区| 国产精品女同久久久久电影院 | 久久亚洲国产精品一区二区| 久久久国产打桩机| 色天使久久综合网天天| 国产精品久久久久a影院|