• <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)的地方,就是測(cè)試數(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 閱讀(283) 評(píng)論(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

            最新隨筆

            搜索

            積分與排名

            最新評(píng)論

            閱讀排行榜

            好久久免费视频高清| 国产精品女同久久久久电影院| 精品久久久久久久国产潘金莲| 99久久精品免费看国产一区二区三区 | 狠狠色丁香婷婷久久综合不卡| 狠狠色噜噜狠狠狠狠狠色综合久久 | 久久婷婷五月综合色99啪ak| 国内精品伊人久久久久妇| 东京热TOKYO综合久久精品| 伊人色综合久久天天| 久久九九久精品国产免费直播| 久久91亚洲人成电影网站| 久久亚洲天堂| 日韩精品久久久久久| 国产亚洲精久久久久久无码77777| 久久久九九有精品国产| 亚洲AV无码久久精品色欲| 精品久久久久久99人妻| 久久99精品久久久久久hb无码| 久久精品极品盛宴观看| 99久久综合国产精品二区| 久久精品亚洲精品国产欧美| 99久久精品国产一区二区| 久久久噜噜噜久久| 精品999久久久久久中文字幕| 99久久99久久精品国产片果冻| 久久精品无码一区二区三区免费 | 国产精品亚洲美女久久久| 久久久久久夜精品精品免费啦| 久久99精品久久久大学生| 人妻系列无码专区久久五月天| 国产精品美女久久久久AV福利| 99久久精品国产高清一区二区 | 99久久精品免费看国产一区二区三区 | 亚洲色大成网站WWW久久九九| 久久久91人妻无码精品蜜桃HD | 日日狠狠久久偷偷色综合免费| 国产精品日韩欧美久久综合| 久久精品国产91久久麻豆自制| 国产精品99久久精品| 国内精品久久久久久99|