亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频

? 歡迎來到蟲蟲下載站! | ?? 資源下載 ?? 資源專輯 ?? 關于我們
? 蟲蟲下載站

?? 1688 dolphin pool.cpp

?? heu暑期集訓中講到的一些ACM/ICPC培訓內容
?? CPP
字號:
#include<iostream>
#include<iomanip>
#include<string>
#include<algorithm>
#include<math.h>
#include<assert.h>
using namespace std;
#define NEG 1e-7  //浮點數比較中的允許誤差 

struct intersect{      //兩圓交點 
    int other;         //記錄同一交點在另一圓上的位置。other=另一圓序號*100+交點在圓上的序號。 
    char sign;         //相交弧與交點之間的相對位置(順時針或逆時針) 
    double x, y, x0, y0;  //(x,y)是交點坐標,(x0,y0)是圓心坐標 
}p[20][40];
struct circle{
    int x,y,r;
}cir[20];          //圓 
int n,nsec[20];    //圓的個數和每個圓上的交點數目 
bool v[20][40];    //交點是否訪問過 

/*下面三個函數是快速排序*/ 
bool operator<(const intersect& x, const intersect& y)
{
    return atan2(x.y-x.y0,x.x-x.x0)<atan2(y.y-y.y0,y.x-y.x0)-NEG;
}
int Partition(intersect b[],int l,int r,int a)
{
    intersect x=b[l],t;
    int i=l-1,j=r+1;
    while(1)
    {
        do{
            --j;
        }while(x<b[j]);
        do{
            ++i;
        }while(b[i]<x);
        if(i<j){
            p[b[i].other/100][b[i].other%100].other=a*100+j;  //排序中加入對other屬性的維護 
            p[b[j].other/100][b[j].other%100].other=a*100+i;
            t=b[i];b[i]=b[j];b[j]=t;
        }
        else return j;
    }
    return l;
}
void quicksort(intersect b[],int l,int r,int a)
{
    int n;
    if(l<r)
    {
        n=Partition(b,l,r,a);
        quicksort(b,l,n,a);
        quicksort(b,n+1,r,a);
    }
    return;
}

/*求兩圓交點*/ 
void GetPoint(int a,int b,intersect& p1,intersect& p2)
{
    int A,B,C,D,E,F;
    double A1,B1,C1;
    A=-2*cir[a].x;B=-2*cir[a].y;
    C=cir[a].x*cir[a].x+cir[a].y*cir[a].y-cir[a].r*cir[a].r;
    D=-2*cir[b].x;E=-2*cir[b].y;
    F=cir[b].x*cir[b].x+cir[b].y*cir[b].y-cir[b].r*cir[b].r;
    if(A==D)                           //交線平行于x軸 
    {
        p1.y=p2.y=(F-C)/(double)(B-E);
        p1.x=(-A+sqrt((double)A*A-4*(p1.y*p1.y+B*p1.y+C)))/2.0;
        p2.x=(-A-sqrt((double)A*A-4*(p2.y*p2.y+B*p2.y+C)))/2.0;
    }
    else if(B==E)                      //交線平行于y軸 
    {
        p1.x=p2.x=(F-C)/(double)(A-D);
        p1.y=(-B+sqrt((double)B*B-4*(p1.x*p1.x+A*p1.x+C)))/2.0;
        p2.y=(-B-sqrt((double)B*B-4*(p2.x*p2.x+A*p2.x+C)))/2.0;
    }
    else                               //一般情形 
    {
        A1=(double)(E-B)*(E-B)+(double)(A-D)*(A-D);
        B1=2.0*(E-B)*(F-C)+((double)D)*(A-D)*(E-B)+(A-D)*(A-D)*(double)E;
        C1=(double)(F-C)*(F-C)+((double)D)*(A-D)*(F-C)+((double)F)*(A-D)*(A-D);
        p1.y=(-B1+sqrt(B1*B1-4*A1*C1))/(2*A1);
        p2.y=(-B1-sqrt(B1*B1-4*A1*C1))/(2*A1);
        p1.x=(p1.y*(E-B)+(F-C))/(A-D);
        p2.x=(p2.y*(E-B)+(F-C))/(A-D);
    }
    return;
}
/*判斷焦點p0是否被圓k覆蓋*/ 
bool Cover(int k,intersect& p0)
{
    if((p0.x-cir[k].x)*(p0.x-cir[k].x)+(p0.y-cir[k].y)*(p0.y-cir[k].y)>cir[k].r*cir[k].r+NEG)
        return false;
    else
        return true;
}
/*依序訪問交點,尋找小區域*/ 
int Visit(int i,int j)
{
    double lx,ly,s;
    s=0;lx=p[i][j].x;ly=p[i][j].y;
    while(!v[i][j])
    {
        v[i][j]=true;
        i=p[i][j].other;
        j=i%100;i/=100;           //找下一個圓 
        v[i][j]=true;
        if(j<nsec[i]-1)j++;       //找下一個交點 
        else j=0;
        s+=(lx*p[i][j].y-ly*p[i][j].x)/2;  //計算有向面積 
        lx=p[i][j].x;ly=p[i][j].y;
    }                             //循環y結束,找到一個由圓弧構成的環 
    if(s>-NEG)return 0;           //有向面積為正,該區域邊界是由多個圓組成的連通區域的外邊界 
    else return 1;                //有向區域為負,該區域是由多個圓圍成的小區域。 
}
int main()
{
    int t,index,i,j,k;
    intersect p1,p2;
    cin>>t;
    for(index=1;index<=t;index++)
    {
    /*讀入數據*/ 
        cin>>n;
        for(i=0;i<n;i++)
            cin>>cir[i].x>>cir[i].y>>cir[i].r;
        fill(nsec,nsec+n,0);
        
        /*求交點*/ 
        for(i=0;i<n;i++)
        {
            for(j=i+1;j<n;j++)
                if((cir[i].x-cir[j].x)*(cir[i].x-cir[j].x)+(cir[i].y-cir[j].y)*(cir[i].y-cir[j].y)<(cir[i].r+cir[j].r)*(cir[i].r+cir[j].r))
                {
                    GetPoint(i,j,p1,p2);            //求兩圓交點 
                    for(k=0;k<n;k++)
                        if(k!=i&&k!=j&&Cover(k,p1))break;     //判斷交點是否被覆蓋
                    if(k==n)
                    {
                        p[i][nsec[i]].x=p1.x;p[i][nsec[i]].y=p1.y;
                        p[j][nsec[j]].x=p1.x;p[j][nsec[j]].y=p1.y;
                        p[i][nsec[i]].x0=cir[i].x;p[i][nsec[i]].y0=cir[i].y;
                        p[j][nsec[j]].x0=cir[j].x;p[j][nsec[j]].y0=cir[j].y;
                        p[i][nsec[i]].other=j*100+nsec[j];
                        p[j][nsec[j]].other=i*100+nsec[i];
                        //根據兩交點對應半徑的叉積判斷相交弧在交點的順/逆時針方向 
                        if((p1.x-cir[i].x)*(p2.y-cir[i].y)-(p1.y-cir[i].y)*(p2.x-cir[i].x)<-NEG)  
                        {
                            p[i][nsec[i]].sign='+';p[j][nsec[j]].sign='-';
                        }
                        else
                        {
                            p[i][nsec[i]].sign='-';p[j][nsec[j]].sign='+';
                        }
                        nsec[i]++;nsec[j]++;
                    }
                    for(k=0;k<n;k++)      //對另一個交點操作,同上 
                        if(k!=i&&k!=j&&Cover(k,p2))break;
                    if(k==n)
                    {
                        p[i][nsec[i]].x=p2.x;p[i][nsec[i]].y=p2.y;
                        p[j][nsec[j]].x=p2.x;p[j][nsec[j]].y=p2.y;
                        p[i][nsec[i]].x0=cir[i].x;p[i][nsec[i]].y0=cir[i].y;
                        p[j][nsec[j]].x0=cir[j].x;p[j][nsec[j]].y0=cir[j].y;
                        p[i][nsec[i]].other=j*100+nsec[j];
                        p[j][nsec[j]].other=i*100+nsec[i];
                        if((p1.x-cir[i].x)*(p2.y-cir[i].y)-(p1.y-cir[i].y)*(p2.x-cir[i].x)<-NEG)
                        {
                            p[i][nsec[i]].sign='-';p[j][nsec[j]].sign='+';
                        }
                        else
                        {
                            p[i][nsec[i]].sign='+';p[j][nsec[j]].sign='-';
                        }
                        nsec[i]++;nsec[j]++;
                    }
                }
            quicksort(p[i],0,nsec[i]-1,i);
        }
        fill(&v[0][0],&v[n][0],false);
        int s=0;
        /*從每個未訪問的交點出發,尋找包含該點的區域*/
        for(i=0;i<n;i++)
            for(j=0;j<nsec[i];j++)
                if(!v[i][j])   
                    if(p[i][j].sign=='-')s+=Visit(i,j);
                    else s+=Visit(p[i][j].other/100,p[i][j].other%100);
        cout<<s<<endl;
    }
    return 0;
}

?? 快捷鍵說明

復制代碼 Ctrl + C
搜索代碼 Ctrl + F
全屏模式 F11
切換主題 Ctrl + Shift + D
顯示快捷鍵 ?
增大字號 Ctrl + =
減小字號 Ctrl + -
亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频
久久久91精品国产一区二区精品 | 亚洲欧美乱综合| 一区二区三区欧美| 久草在线在线精品观看| 欧美优质美女网站| 欧美国产一区在线| 久久99精品视频| 7777精品久久久大香线蕉| 亚洲欧美日韩一区| 成人av资源下载| 久久精品一区二区三区不卡| 美女一区二区三区在线观看| 欧美体内she精高潮| 亚洲欧美一区二区三区孕妇| 国产成人精品网址| 久久综合成人精品亚洲另类欧美| 五月天一区二区三区| 91国产成人在线| 一区二区三区日韩欧美| 成人avav在线| 国产精品国产三级国产专播品爱网| 看电影不卡的网站| 精品成人在线观看| 九一九一国产精品| 日韩精品一区二| 老鸭窝一区二区久久精品| 日韩女同互慰一区二区| 久久99久久久欧美国产| 日韩一区二区精品葵司在线 | 国产一区中文字幕| 日韩精品一区二区三区中文不卡| 日韩av电影天堂| 日韩视频一区二区三区| 美腿丝袜亚洲综合| 26uuu色噜噜精品一区| 国产91高潮流白浆在线麻豆| 国产精品视频你懂的| av一区二区三区| 亚洲免费伊人电影| 欧美日韩一级片网站| 日韩激情中文字幕| 精品久久久久久久久久久久包黑料| 久久aⅴ国产欧美74aaa| 久久精品亚洲麻豆av一区二区 | 欧美日韩国产美女| 热久久久久久久| 精品三级在线观看| 成人av资源站| 亚洲精品视频免费观看| 欧美日韩国产中文| 久久97超碰色| 中文字幕视频一区二区三区久| av亚洲精华国产精华精| 日韩伦理免费电影| 在线免费观看日韩欧美| 奇米色一区二区三区四区| 久久日一线二线三线suv| 波多野结衣91| 日日夜夜免费精品| 欧美激情一区在线| 欧美日韩一区精品| 国产成人自拍网| 亚洲色图清纯唯美| 欧美一区二区三区四区五区| 国产精品资源网站| 亚瑟在线精品视频| 久久嫩草精品久久久精品一| av成人免费在线观看| 午夜视频一区二区三区| 2020国产精品| 欧美亚洲精品一区| 国产精品一品二品| 亚洲一区av在线| 精品成人a区在线观看| 色综合久久久久网| 久久99精品国产麻豆不卡| 1000精品久久久久久久久| 日韩一区二区三区电影| 色激情天天射综合网| 激情综合网av| 亚洲一区二区三区在线播放| 国产丝袜在线精品| 日韩三级.com| 欧美在线视频不卡| av一本久道久久综合久久鬼色| 精品一区二区在线观看| 亚洲成av人片在www色猫咪| 国产精品欧美一区喷水| 日韩欧美一二区| 精品污污网站免费看| 国产精品一品二品| 国产综合色在线| 午夜精品123| 一区二区三区四区在线播放 | 精品国产a毛片| 欧美精品欧美精品系列| 91麻豆自制传媒国产之光| 国内精品嫩模私拍在线| 首页亚洲欧美制服丝腿| 亚洲激情综合网| 亚洲四区在线观看| 国产精品美女久久福利网站| 精品1区2区在线观看| 欧美一区二区免费观在线| 欧美日韩免费电影| 在线一区二区观看| 一本到三区不卡视频| 播五月开心婷婷综合| 成人性生交大合| 丁香激情综合五月| 国产成人av一区二区三区在线观看| 男男视频亚洲欧美| 蜜桃精品视频在线观看| 久久av老司机精品网站导航| 人人爽香蕉精品| 精品一区二区成人精品| 国精产品一区一区三区mba视频| 麻豆高清免费国产一区| 久久精品国产秦先生| 韩国女主播一区| 国产91精品免费| 91香蕉视频黄| 欧美日韩不卡视频| 欧美一二三区在线观看| wwwwww.欧美系列| 亚洲国产精品v| 亚洲人成影院在线观看| 亚洲在线成人精品| 日韩av一二三| 国产精品18久久久久久久久久久久 | 中文字幕免费不卡在线| 成人免费一区二区三区在线观看| 亚洲视频在线一区| 丝袜诱惑制服诱惑色一区在线观看| 亚洲a一区二区| 久久99精品久久只有精品| 国产91清纯白嫩初高中在线观看| www.久久久久久久久| 日本高清免费不卡视频| 日韩一区二区不卡| 久久精品视频一区二区三区| 中文字幕乱码日本亚洲一区二区 | 精品国产区一区| 国产精品色哟哟| 亚洲国产美女搞黄色| 韩日欧美一区二区三区| 日本黄色一区二区| 精品久久久久久亚洲综合网| 中文字幕一区二区三区在线播放 | xvideos.蜜桃一区二区| 亚洲人成网站精品片在线观看| 亚洲成av人片在www色猫咪| 国产精品综合视频| 欧美日韩视频第一区| 久久久亚洲精品石原莉奈| 一区二区三区不卡视频在线观看| 日韩福利视频导航| 99国产精品久| 精品国产人成亚洲区| 亚洲第一电影网| 成人手机在线视频| 欧美一区二区三区小说| 日韩毛片精品高清免费| 国产在线视频一区二区| 日本高清免费不卡视频| 国产亚洲成av人在线观看导航 | 亚洲一区二区偷拍精品| 激情文学综合插| 欧美日韩一区三区四区| 中文字幕亚洲成人| 国产在线精品一区二区| 欧美午夜不卡在线观看免费| 中文字幕免费在线观看视频一区| 日韩成人精品在线观看| 91官网在线观看| 综合久久国产九一剧情麻豆| 韩国成人精品a∨在线观看| 制服.丝袜.亚洲.中文.综合| 亚洲欧洲成人自拍| 国产成人h网站| 久久精品男人天堂av| 日韩电影在线免费观看| 欧美在线你懂的| 亚洲欧美日韩一区二区| 不卡欧美aaaaa| 国产精品蜜臀av| 国产成人日日夜夜| 久久婷婷国产综合精品青草| 久久精品国内一区二区三区| 欧美精品欧美精品系列| 亚洲成人动漫在线免费观看| 91国产视频在线观看| 亚洲精品日日夜夜| 在线视频综合导航| 亚洲欧美电影一区二区| 91蜜桃在线免费视频| 中文字幕中文在线不卡住| 国产成人亚洲综合a∨猫咪| 欧美激情一区不卡| 成人黄色av电影|