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

            最新隨筆

            搜索

            積分與排名

            最新評論

            閱讀排行榜

            99久久人妻无码精品系列蜜桃| 日韩精品久久久久久久电影蜜臀| 99久久国产亚洲高清观看2024| 91亚洲国产成人久久精品网址| 久久精品国产99久久丝袜| 久久久久久久免费视频| 久久精品无码午夜福利理论片 | 久久久久九国产精品| 一本一本久久aa综合精品| 青青青国产成人久久111网站| 亚洲国产婷婷香蕉久久久久久| 狠狠色丁香婷婷久久综合不卡 | 国产成人综合久久久久久| 狠狠色丁香久久婷婷综合_中| 91精品国产高清久久久久久io| 人妻无码精品久久亚瑟影视| 久久亚洲精品中文字幕三区| 色播久久人人爽人人爽人人片AV| 99久久精品这里只有精品| 人妻精品久久久久中文字幕69| 久久青青草原亚洲av无码| 精品久久久久久久久久中文字幕| 国产精品岛国久久久久| 久久人做人爽一区二区三区 | 久久伊人精品青青草原高清| 久久免费的精品国产V∧| 国产精品99久久久久久宅男小说| 国产亚州精品女人久久久久久 | 伊人久久综合精品无码AV专区| 久久强奷乱码老熟女| 久久久精品视频免费观看| 93精91精品国产综合久久香蕉| 国产午夜免费高清久久影院| 久久久久久久97| 国内精品伊人久久久久AV影院| 久久精品国产亚洲av影院| 大伊人青草狠狠久久| 91精品国产高清久久久久久91| 日本精品久久久中文字幕| 久久精品国产清自在天天线| 久久亚洲电影|