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

            pku 1228 Grandpa's Estate 凸多邊形的唯一性(凸包)

            題意:
            一個(gè)凸多邊形的邊界上有若干木樁,現(xiàn)丟失部分木樁,問由剩下的木樁能否唯一確定這個(gè)多邊形

            解法:
            首先能夠唯一確定的條件是由剩下的木樁確定的凸包的每條邊上至少包含3個(gè)木樁,這個(gè)自己畫圖比劃下就知道了- -
            然后就是求一個(gè)凸包了。在這種坐標(biāo)都是整數(shù)的情況下,凸包最好不要用atan2函數(shù),而是用叉積來比較。我特地用純C寫了個(gè),有要的童鞋可以拿去當(dāng)模板
            有個(gè)陰險(xiǎn)的地方,就是測試數(shù)據(jù)只有3個(gè)點(diǎn),而且3點(diǎn)一線。。。你懂的

            代碼
             1#  include <stdio.h>
             2#  include <stdlib.h>
             3# define N 1200
             4# define cross(x1,y1,x2,y2) ((x1)*(y2)-(x2)*(y1))
             5# define min(a,b) ((a)<(b)?(a):(b))
             6# define max(a,b) ((a)>(b)?(a):(b))
             7typedef struct
             8{
             9    int x,y;
            10}
            point;
            11int n,c;
            12point data[N],ans[N],std;
            13int dis(point *pos)
            14{
            15    return (pos->x-std.x)*(pos->x-std.x)+(pos->y-std.y)*(pos->y-std.y);
            16}

            17int isin(point *a,point *b,point *pos)
            18{
            19    if(pos->x>max(a->x,b->x)||pos->x<min(a->x,b->x)||pos->y>max(a->y,b->y)||pos->y<min(a->y,b->y)) return 0;
            20    else if(cross(pos->x-a->x,pos->y-a->y,b->x-a->x,b->y-a->y)!=0return 0;
            21    else return 1;
            22}

            23int cmp(const void *a,const void *b)
            24{
            25    point *aa=(point *)a,*bb=(point *)b;
            26    if(cross(bb->x-std.x,bb->y-std.y,aa->x-std.x,aa->y-std.y))
            27            return cross(bb->x-std.x,bb->y-std.y,aa->x-std.x,aa->y-std.y);
            28    else 
            29            return dis(aa)-dis(bb);
            30}

            31void sort()
            32{
            33    int i;
            34    int x=0xfffffff,y=0xfffffff;
            35    for(i=0;i<n;i++)
            36        if(data[i].y<y||data[i].y==y&&data[i].x<x)
            37            y=data[i].y,x=data[i].x;
            38    std.x=x;
            39    std.y=y;
            40    qsort(data,n,sizeof(point),cmp);
            41}

            42void build()
            43{
            44    int i;
            45    c=0;
            46    sort();
            47    for(i=0;i<n;i++)
            48    {
            49        while(c>=2&&cross(data[i].x-ans[c-1].x,data[i].y-ans[c-1].y,ans[c-1].x-ans[c-2].x,ans[c-1].y-ans[c-2].y)>=0) c--;
            50        ans[c++]=data[i];
            51    }

            52    if(c>0) ans[c++]=ans[0];
            53}

            54int chk()
            55{
            56    int i;
            57    for(i=0;i<c-1;i++)
            58    {
            59        int count=0,j;
            60        for(j=0;j<n;j++)
            61            if(isin(&ans[i],&ans[i+1],&data[j]))
            62                count++;
            63        if(count<3return 0;
            64    }

            65    return 1;
            66}

            67int main()
            68{
            69    int test;
            70    scanf("%d",&test);
            71    while(test--)
            72    {
            73        int i;
            74        scanf("%d",&n);
            75        for(i=0;i<n;i++)
            76            scanf("%d %d",&data[i].x,&data[i].y);
            77        build();
            78        if(c>3&&chk()) printf("YES\n");
            79        else printf("NO\n");
            80    }

            81    return 0;
            82}

            83

            posted on 2011-01-15 02:27 yzhw 閱讀(296) 評論(0)  編輯 收藏 引用 所屬分類: geometry&phycise

            <2011年1月>
            2627282930311
            2345678
            9101112131415
            16171819202122
            23242526272829
            303112345

            導(dǎo)航

            統(tǒng)計(jì)

            公告

            統(tǒng)計(jì)系統(tǒng)

            留言簿(1)

            隨筆分類(227)

            文章分類(2)

            OJ

            最新隨筆

            搜索

            積分與排名

            最新評論

            閱讀排行榜

            亚洲AV成人无码久久精品老人 | 久久久精品免费国产四虎| 久久午夜综合久久| 久久久高清免费视频| 久久香蕉超碰97国产精品| 精品久久人人爽天天玩人人妻| 日日狠狠久久偷偷色综合0| 伊人久久综合成人网| 狠狠精品久久久无码中文字幕| 中文字幕无码久久人妻| 久久99精品国产麻豆宅宅| 一级女性全黄久久生活片免费| 国产一区二区三区久久| 亚洲国产精品久久久天堂| 久久精品国产精品亚洲| 九九99精品久久久久久| 2021最新久久久视精品爱| 久久精品国产亚洲5555| 国产精品99久久免费观看| 久久精品国产亚洲av麻豆图片| 久久人人超碰精品CAOPOREN| 日本三级久久网| 久久免费精品视频| 99久久人妻无码精品系列| 国内精品九九久久精品| 免费精品久久久久久中文字幕| 国产高清国内精品福利99久久| 99久久超碰中文字幕伊人 | 国产精品久久久福利| 久久国产精品无码一区二区三区 | 久久人妻AV中文字幕| 久久精品国产亚洲Aⅴ香蕉 | 久久精品国产99国产精品澳门 | 国产欧美久久一区二区| 中文字幕久久久久人妻| 亚洲综合伊人久久综合| 亚洲人成伊人成综合网久久久| 亚洲精品蜜桃久久久久久| 亚洲va久久久噜噜噜久久男同| 少妇内射兰兰久久| 久久国产欧美日韩精品|