HDU 1665 or UVALive 3263 || That Nice Euler Circuit （欧拉定理求面数 F = E+2-V

N个端点的一笔画，第N个端点与点一个端点重合。

求这条闭合曲线，将平面分成多少部分。（包括封闭区域和无限大区域

欧拉定理：

设平面图的顶点数 V ，边数 E ，面数 F ，则 F = E+2-V

分析：

1：算出顶点数，包括交点，有可能三条线段交于一点，需要去重。

2：算出线段数，如果原来的线段新增加一个点，那么一条线段一为二，所以一条线段每新增一个点就多出一条线段。

代码细节见注释。

函数：

unique --》返回去重（需排序）后的数组最后一个元素的地址。（并不是删除元素而是将这些元素后移。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;

struct point
{
    double x,y,d;
    point(double a,double b):x(a),y(b){}
    point(){}
};
typedef point vec;
//向量+向量=向量
//点  +向量=点
vec operator +  ( vec a,   vec b)      {  return vec(a.x+b.x,a.y+b.y);}

vec operator -  ( point a, point b)    {  return vec(a.x-b.x,a.y-b.y);}

vec operator *  ( vec a,   double p)   {  return vec(a.x*p,a.y*p);}

vec operator /  ( vec a,   double p)   {  return vec(a.x/p,a.y/p);}

bool operator < ( const point &a,const point &b)
{
     return a.x<b.x || (a.x == b.x &&a.y<b.y);
}

const double eps =1e-6;
int dcmp( double x)
{
    if( fabs(x)<eps )return 0;

    return x < 0?-1:1;
}

bool operator == ( const point &a,const point &b)
{
    return dcmp(a.x-b.x) == 0 && dcmp(a.y-b.y)==0;
}
// 点积（向量 a，向量 b）
double dot( vec a,vec b)    {   return a.x*b.x + a.y*b.y;}


double cross( vec a,vec b ){    return a.x*b.y - a.y*b.x;}


//两条直线 p+tv 和 q + tw 有唯一交点，当且仅当 cross（v,w）！=0
point getlinetnersection(point p,vec v,point q,vec w)
{
    vec u = p-q;
    double t = cross(w,u) /cross(v,w);
    return p + v*t;
}

// 两条线段规范相交：每条线段的两个端点在另一条线段的两侧
bool segproperinter( point a1,point a2,point b1,point b2)
{
     double c1 = cross( a2-a1,  b1-a1 ),
            c2 = cross( a2-a1,  b2-a2 ),
            c3 = cross( b2-b1,  a1-b1 ),
            c4 = cross( b2-b1,  a2-b1 );

    return dcmp( c1 ) * dcmp( c2 ) < 0
            && dcmp( c3 )*dcmp( c4 )<0;
}
//一个点在另外一条线段上(不包括端点
bool onsegment( point p,point a1,point a2)
{
    return dcmp( cross(a1-p,a2-p)) == 0
            && dcmp( dot(a1-p,a2-p))<0;
}
const int MX = 300+5;
point p[MX],v[MX*MX/2+5];
int main()
{
    int n;
    int cas = 1;
    while( scanf("%d",&n) && n )
    {
        for( int i = 0 ; i < n ;++i)
        {
            scanf("%lf %lf",&p[i].x,&p[i].y);
            v[i] = p[i];
        }

        n--;// 第一个点与最后一个点一致
        int top = n,e = n;//顶点数 边数

        for( int i=0; i<n ;++i)//两两遍历是否有交点
            for( int j =1+i; j<n;++j)
        {
            if( segproperinter(p[i],p[i+1],p[j],p[j+1]) )//规范相交
                v[top++] = getlinetnersection(//覆盖最后重复的一点
                            p[i],p[i+1]-p[i],p[j],p[j+1]-p[j]);
        }

        sort(v,v+top);
        top = unique( v,v+top)-v;//去重后的最后数组长度

        //交点在线段上边数加一 ps:线段一定要使用原点集！！！
        for( int i = 0 ; i<top;++i)
            for( int j=0;j<n;++j)
            if( onsegment(v[i],p[j],p[j+1]))
                e++;
        printf("Case %d: There are %d pieces.\n",
               cas++,e+2-top);



    }
    return 0;
}