• <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>
            隨筆-19  評論-1  文章-0  trackbacks-0
             1.  一種是用矢量叉乘法:
                    由三個頂點向所求的點引出矢量(注意方向),然后任意用其中兩個矢量形成平面,再用這個平面和剩下的矢量叉乘,得出一個新矢量,方向向里,則在三角形外,反之在里面。
            2.用面積方法

            #include<stdio.h>
            #include
            <math.h>
            struct TPoint {
                
            float x;
                
            float y;
            }
            ;

            //求叉積
            float mul(struct TPoint p1, struct TPoint p2, struct TPoint p0) {
                
            return ((p1.x - p0.x)*(p2.y - p0.y)-(p2.x - p0.x)*(p1.y - p0.y));
            }

            /*由三個頂點向所求的點引出矢量(注意方向),然后任意用其中兩個矢量形成平面,
             * 再用這個平面和剩下的矢量叉乘,得出一個新矢量,方向向里,則在三角形外,反之在里面。
             
            */

            int inside(struct TPoint tr[], struct TPoint p) {
                
            int i;
                
            for (i = 0; i < 3; i++)
                    
            if (mul(p, tr[i], tr[(i + 1% 3]) * mul(p, tr[(i + 2% 3], tr[(i + 1% 3]) > 0)
                        
            return 0;
                
            return 1;
            }


            float area(struct TPoint p1, struct TPoint p2, struct TPoint p3) {
                
            return fabs((p1.x - p3.x)*(p2.y - p3.y)-(p2.x - p3.x)*(p1.y - p3.y));
            }

            //用面積判斷p是否在三角形內
            int inside2(struct TPoint tr[], struct TPoint p) {
                
            if (fabs(area(tr[0], tr[1], tr[2]) -
                        area(p, tr[
            1], tr[2]) -
                        area(tr[
            0], p, tr[2]) -
                        area(tr[
            0], tr[1], p)) < 1.0e-20)
                    
            return 1;
                
            else
                    
            return 0;
            }


            int main() {
                
            struct TPoint tr[3= {{-11},{10},{30}},  p = {12};

                
            //方法一
                printf("algorithm   1:");
                
            if (inside(tr, p))
                    printf(
            "In\n");
                
            else
                    printf(
            "Out\n");

                
            //方法一
                printf("algorithm   2:");
                
            if (inside2(tr, p))
                    printf(
            "In\n");
                
            else
                    printf(
            "Out\n");
            }
            posted on 2010-10-12 09:40 孟起 閱讀(1883) 評論(0)  編輯 收藏 引用 所屬分類: 計算幾何
            精品综合久久久久久97| 94久久国产乱子伦精品免费| 色综合久久天天综合| 狠狠人妻久久久久久综合蜜桃| 久久久久99精品成人片| 亚洲精品高清国产一线久久| 久久99国内精品自在现线| 色偷偷888欧美精品久久久| 久久久久婷婷| 精品久久久久久无码专区| 久久夜色精品国产亚洲| 欧美日韩精品久久免费| 91精品日韩人妻无码久久不卡 | 中文字幕日本人妻久久久免费| 久久久久亚洲av无码专区导航 | 欧美精品一区二区久久| 成人免费网站久久久| 狠狠色丁香婷婷久久综合 | 久久er热视频在这里精品| 久久久久久久综合狠狠综合| 久久久中文字幕| 精品蜜臀久久久久99网站| 国产精品乱码久久久久久软件| 久久亚洲综合色一区二区三区| 久久久久人妻精品一区二区三区| 污污内射久久一区二区欧美日韩| 久久久精品一区二区三区| 97久久精品人妻人人搡人人玩| 97久久婷婷五月综合色d啪蜜芽 | 久久精品桃花综合| 久久青青草原精品国产不卡| 久久精品中文字幕久久| 国产一区二区三区久久| 久久亚洲私人国产精品vA| 精品久久久久久国产| 18禁黄久久久AAA片| 99久久免费国产精品特黄| 久久亚洲AV无码精品色午夜麻豆 | 久久精品国产亚洲一区二区| 狠狠狠色丁香婷婷综合久久五月| 国产美女久久精品香蕉69|