• <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 凸多邊形的唯一性(凸包)

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

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

            代碼
             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 閱讀(289) 評論(0)  編輯 收藏 引用 所屬分類: geometry&phycise

            <2010年10月>
            262728293012
            3456789
            10111213141516
            17181920212223
            24252627282930
            31123456

            導航

            統計

            公告

            統計系統

            留言簿(1)

            隨筆分類(227)

            文章分類(2)

            OJ

            最新隨筆

            搜索

            積分與排名

            最新評論

            閱讀排行榜

            久久久久久久波多野结衣高潮| 亚洲伊人久久精品影院| 理论片午午伦夜理片久久 | 久久无码高潮喷水| 久久夜色精品国产噜噜麻豆 | 午夜精品久久久内射近拍高清 | 久久久久亚洲精品男人的天堂| 日本WV一本一道久久香蕉| 国产综合久久久久| 一本综合久久国产二区| 久久精品国产半推半就| 久久九九兔免费精品6| 久久99久久成人免费播放| 久久亚洲春色中文字幕久久久| 久久久久97国产精华液好用吗| 国产精品久久久久久搜索 | 国产午夜精品久久久久九九| 久久亚洲AV成人无码| 看全色黄大色大片免费久久久| 久久亚洲国产午夜精品理论片| 久久精品久久久久观看99水蜜桃| 精品久久久久久国产免费了| 天天爽天天狠久久久综合麻豆| 女人高潮久久久叫人喷水| 久久免费99精品国产自在现线| 香蕉久久一区二区不卡无毒影院 | 国产福利电影一区二区三区久久久久成人精品综合 | 无码伊人66久久大杳蕉网站谷歌| 日韩十八禁一区二区久久| 国产—久久香蕉国产线看观看| 97久久精品无码一区二区天美| 亚洲中文字幕无码一久久区| 久久久精品久久久久影院| 伊人久久一区二区三区无码| 久久久久国产一区二区| 精品免费久久久久国产一区 | 久久综合给合久久狠狠狠97色| 精品国产乱码久久久久久人妻| 久久久久久国产精品无码下载| 久久无码专区国产精品发布| 亚洲日韩中文无码久久|