• <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>
            mmd 和 cz 的智慧結晶。某次比賽的時候寫的。


            /*
            1:兩個圓相離
            2:兩個圓外切
            3:兩個圓內切
            4:兩個圓相交
            5:如果兩個圓相含(包括兩圓坐標和半徑都相同)
            */


            int circle_intersect(circle_t A, circle_t B, point_t &ia, point_t &ib)
            {
                
            if (A.x == B.x && A.y == B.y)
                    
            return 5;
                
            double dd = dist (A.x, A.y, B.x, B.y);
                
            if (A.r + B.r + EP < dd)
                    
            return 1;
                
                
            double k, a, b, d, aa, bb, cc, c, drt;
                k 
            = A.r;
                a 
            = B.x - A.x;
                b 
            = B.y - A.y;
                c 
            = B.r;
                d 
            = sqr(c) - sqr(k) - sqr(a) - sqr(b);
                
                aa 
            = 4 * sqr(a) + 4 * sqr(b);
                bb 
            = 4 * b * d;
                cc 
            = sqr(d) - 4 * sqr(a) * sqr(k);
                
                drt 
            = sqr(bb) - 4 * aa * cc;
                
            if (drt < 0)
                    
            return 5;
                drt 
            = sqrt (drt);
                ia.y 
            = (-bb + drt) / 2 / aa;
                ib.y 
            = (-bb - drt) / 2 / aa;
                
            if (abs (a) < EP)
                
            {
                    ia.x 
            = sqrt (sqr (k) - sqr (ia.y));
                    ib.x 
            = -ia.x;
                }

                
            else
                
            {
                    ia.x 
            = (2 * b * ia.y + d) / -2 / a;
                    ib.x 
            = (2 * b * ib.y + d) / -2 / a;
                }

                ia.x 
            += A.x;
                ia.y 
            += A.y;
                ib.x 
            += A.x;
                ib.y 
            += A.y;
                
            if (abs (ia.y - ib.y) < EP)
                
            {
                    
            if (abs (A.r + B.r - dd) < EP)
                        
            return 2;
                    
            if (abs (dd - (max (A.r, B.r) - min (A.r, B.r))) < EP)
                        
            return 3;
                }

                
            return 4;
            }
            posted on 2007-09-20 18:34 Felicia 閱讀(5880) 評論(7)  編輯 收藏 引用 所屬分類: 計算幾何
            Comments
            • # re: [計算幾何]兩圓求交點
              Jie
              Posted @ 2007-12-03 18:47
              弱問EP是什么?  回復  更多評論   
            • # re: [計算幾何]兩圓求交點
              coder
              Posted @ 2008-04-22 11:46
              sqr,sqrt,乘法,等計算很慢的。能不能優化一下?  回復  更多評論   
            • # re: [計算幾何]兩圓求交點
              hehe
              Posted @ 2008-05-12 19:09
              EP是什么呢  回復  更多評論   
            • # re: [計算幾何]兩圓求交點
              lala
              Posted @ 2009-02-06 18:05
              我小白,請問sqr是什么?  回復  更多評論   
            • # re: [計算幾何]兩圓求交點[未登錄]
              Felicia
              Posted @ 2009-02-09 07:45
              平方  回復  更多評論   
            • # re: [計算幾何]兩圓求交點
              感謝
              Posted @ 2009-06-06 00:18
              那個方程真的難解,多謝你的代碼  回復  更多評論   
            • # re: [計算幾何]兩圓求交點[未登錄]
              路人甲
              Posted @ 2010-08-27 16:01
              @感謝
              感覺寫得有些復雜。令L為兩圓心間距,則可以構造三邊分別為R1 R2 L的三角形。進而解方程求得十字架各部分的長度(包括割線的長度),再通過向量運算就得到兩個交點了。  回復  更多評論   
             
            久久久久99精品成人片| 亚洲精品无码久久毛片| 久久99热只有频精品8| 大伊人青草狠狠久久| 香蕉久久一区二区不卡无毒影院| 91久久精品视频| 久久人人爽人人爽人人片AV麻烦| 亚洲中文久久精品无码| 大美女久久久久久j久久| 久久无码AV中文出轨人妻| 999久久久免费精品国产| 久久人人爽人人爽人人片AV东京热| 久久久久久午夜精品| 色综合久久最新中文字幕| 97久久国产综合精品女不卡| 国产精品成人久久久久三级午夜电影| 久久久久久免费视频| 精品久久久久久久久久中文字幕| 国内高清久久久久久| 久久久久国产精品麻豆AR影院| 亚洲精品美女久久久久99| 观看 国产综合久久久久鬼色 欧美 亚洲 一区二区 | 蜜桃麻豆www久久| 久久久久久国产精品无码下载| 久久99中文字幕久久| 久久成人国产精品免费软件| 久久99精品免费一区二区| 久久国产热精品波多野结衣AV| 久久天天躁狠狠躁夜夜avapp| 青青草原1769久久免费播放| 亚洲AV无码成人网站久久精品大| 欧美一区二区久久精品| 久久精品无码免费不卡| 久久精品男人影院| 伊人丁香狠狠色综合久久| 久久香蕉一级毛片| 中文精品久久久久国产网址| 91精品国产9l久久久久| 久久国产免费观看精品| 久久不射电影网| 婷婷综合久久中文字幕|