• <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>

            poj 1269 Intersecting Lines

               題目意思是給出2條直線,然后判斷其是否相交,平行,還是重合。剛開始以為是判斷2條線段的關系,用了黑書的模板寫了,發(fā)現(xiàn)連樣例
            都過不了。后面改了很多才過了。先判斷2條直線所在的向量是否平行,這個可以判斷這2個向量的叉積是否為0,然后再判斷線段是否重合,
            可以選3點判斷叉積是否為0。如果向量不平行的話,直接求交點。求交點的公式是用了黑書里面的方法,先求出2個叉積代表2個三角形的
            有向面積,然后根據(jù)定比分點的關系(面積的比例等于交點分其中一條線段的比例)可以推出計算公式。
               有叉積和點積這2個工具確實能方便的解決很多問題。

               代碼如下:
            #include <stdio.h>
            #include <string.h>
            #include <math.h>
            struct Point
            {
                double fX;
                double fY;
            };
            Point beg[2], end[2];
            int nN;
            const double fPrecision = 1e-8;

            double Det(double fX1, double fY1, double fX2, double fY2)
            {
                return fX1 * fY2 - fX2 * fY1;
            }

            double Cross(Point a, Point b, Point c)
            {
                return Det(b.fX - a.fX, b.fY - a.fY, c.fX - a.fX, c.fY - a.fY);
            }

            int DblCmp(double fD)
            {
                if (fabs(fD) < fPrecision)
                {
                    return 0;
                }
                else
                {
                    return (fD > 0 ? 1 : -1);
                }
            }

            double DotDet(double fX1, double fY1, double fX2, double fY2)
            {
                return fX1 * fX2 + fY1 * fY2;
            }

            double Dot(Point a, Point b, Point c)
            {
                return DotDet(b.fX - a.fX, b.fY - a.fY, c.fX - a.fX, c.fY - a.fY);
            }

            int BetweenCmp(Point a, Point b, Point c)
            {
                return DblCmp(Dot(a, b, c));
            }

            int SegCross(Point a, Point b, Point c, Point d, Point& p)
            {
                double s1, s2, s3, s4;
                int d1, d2, d3, d4;
                d1 = DblCmp(s1 = Cross(a, b, c));
                d2 = DblCmp(s2 = Cross(a, b, d));
                d3 = DblCmp(s3 = Cross(c, d, a));
                d4 = DblCmp(s4 = Cross(c, d, b));
                
                Point e, f;
                e.fX = a.fX - b.fX;
                e.fY = a.fY - b.fY;
                f.fX = c.fX - d.fX;
                f.fY = c.fY - d.fY;
                if (DblCmp(Det(e.fX, e.fY, f.fX, f.fY)) == 0)//2個向量共線
                {
                    if (d1 * d2 > 0 && d3 * d4 > 0)//不在同一條直線上
                    {
                        return 0;
                    }
                    else
                    {
                        return 2;
                    }
                }
                
                //2條直線相交
                p.fX = (c.fX * s2 - d.fX * s1) / (s2 - s1);
                p.fY = (c.fY * s2 - d.fY * s1) / (s2 - s1);
                return 1;
            }

            int main()
            {
                //freopen("out.txt", "w", stdout);
                while (scanf("%d", &nN) == 1)
                {
                    printf("INTERSECTING LINES OUTPUT\n");
                    Point p;
                    for (int i = 0; i < nN; ++i)
                    {
                        scanf("%lf%lf%lf%lf", &beg[0].fX, &beg[0].fY, &end[0].fX, &end[0].fY);
                        scanf("%lf%lf%lf%lf", &beg[1].fX, &beg[1].fY, &end[1].fX, &end[1].fY);
                        int nRet = SegCross(beg[0], end[0], beg[1], end[1], p);
                        if (nRet == 0)
                        {
                            printf("NONE\n");
                        }
                        else if (nRet == 1)
                        {
                            printf("POINT %.2f %.2f\n", p.fX, p.fY);
                        }
                        else
                        {
                            printf("LINE\n");
                        }
                    }
                    printf("END OF OUTPUT\n");
                }
                
                return 0;
            }

            posted on 2012-07-17 15:20 yx 閱讀(1033) 評論(0)  編輯 收藏 引用 所屬分類: 計算幾何

            <2012年7月>
            24252627282930
            1234567
            891011121314
            15161718192021
            22232425262728
            2930311234

            導航

            統(tǒng)計

            公告

            常用鏈接

            留言簿(3)

            隨筆分類

            隨筆檔案

            me

            好友

            同學

            網友

            搜索

            最新評論

            閱讀排行榜

            評論排行榜

            精品人妻伦九区久久AAA片69| 亚洲国产精品一区二区久久hs| 亚洲欧美日韩久久精品第一区| 久久国产精品免费一区| 九九99精品久久久久久| 亚洲精品乱码久久久久久蜜桃不卡| 久久精品成人欧美大片| 亚洲精品国产成人99久久| 成人精品一区二区久久| 韩国三级中文字幕hd久久精品| 国内精品久久久久久久coent| 久久本道综合久久伊人| 中文精品99久久国产| 午夜精品久久久久久影视riav| 欧美伊人久久大香线蕉综合 | 国产精品久久久99| 久久se这里只有精品| 人妻无码αv中文字幕久久琪琪布 人妻无码精品久久亚瑟影视 | 日产精品99久久久久久| 久久精品国产亚洲精品2020| 国产亚洲美女精品久久久久狼| 久久er国产精品免费观看8| 亚洲国产成人久久一区WWW| 久久精品国产2020| 热久久这里只有精品| 欧美精品福利视频一区二区三区久久久精品 | 99久久精品国产一区二区蜜芽| 精品久久久久中文字| 久久人人添人人爽添人人片牛牛| 久久久久亚洲av无码专区导航| 99久久国产亚洲高清观看2024 | 国产精品久久久久久搜索| 色诱久久av| 国产亚洲精品美女久久久| 久久精品成人免费观看97| 色欲av伊人久久大香线蕉影院| 四虎国产永久免费久久| 人妻无码精品久久亚瑟影视| 亚洲午夜久久影院| 无码伊人66久久大杳蕉网站谷歌| 久久精品无码专区免费 |