POJ 1279 半平面交

题意：
求一任意多边形的核的面积

题解：

半平面交。第一次搞。。。找了半天错，又和题解一样了。

#include <iostream>
#include <algorithm>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <cmath>

#define N 2222
#define EPS 1e-7
#define INF 1e9

using namespace std;

struct PO
{
    double x,y;
}p[N],tp[N],s[N],o;

int n;

inline void read()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%lf%lf",&p[i].x,&p[i].y);
    p[n+1]=p[1];
}

inline PO operator -(PO a,PO b)
{
    PO c;
    c.x=a.x-b.x;
    c.y=a.y-b.y;
    return c;
}

inline int dc(double x)
{
    if(x>EPS) return 1;
    else if(x<-EPS) return -1;
    return 0;
}

inline double cross(PO &a,PO &b,PO &c)
{
    return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
}

inline double getarea(PO *a,int n)//多边形面积，逆时针为正 
{
    double res=0.0;
    for(int i=1;i<=n;i++) res+=cross(o,a[i],a[i+1]);
    return res*0.5;
}

inline PO getpoint(PO &a,PO &b,PO &c,PO &d)//直线交点 
{
    PO ans,tp=b-a;
    double k1=cross(a,d,c);
    double k2=cross(b,c,d);
    ans.x=a.x+tp.x*k1/(k1+k2);
    ans.y=a.y+tp.y*k1/(k1+k2);
    return ans;
}

inline void change()//变顺时针为逆时针 
{
    for(int i=1;i<=(n>>1);i++) swap(p[i],p[n-i+1]);
    p[n+1]=p[1];
}

void prt(PO a)
{
    printf("%lf   %lf\n",a.x,a.y);
}

inline void getcut()
{
    tp[1].x=tp[5].x=-INF,tp[1].y=tp[5].y=-INF;
    tp[2].x=INF,tp[2].y=-INF;
    tp[3].x=INF,tp[3].y=INF;
    tp[4].x=-INF,tp[4].y=INF;
    
    int cp=4,tc;
    for(int i=1;i<=n;i++)
    {
        tc=0;
        for(int j=1;j<=cp;j++)
        {
            if(dc(cross(p[i],p[i+1],tp[j]))>=0) s[++tc]=tp[j];
            if(dc(cross(p[i],p[i+1],tp[j])*cross(p[i],p[i+1],tp[j+1]))<0)
                s[++tc]=getpoint(p[i],p[i+1],tp[j],tp[j+1]);
        }
        s[tc+1]=s[1];
        for(int j=1;j<=tc+1;j++) tp[j]=s[j];
        cp=tc;
    }
    n=cp;
}

inline void go()
{
    if(dc(getarea(p,n))<=0) change();
    getcut();
    printf("%.2lf\n",fabs(getarea(s,n)));
}

int main()
{
    int cas; scanf("%d",&cas);
    while(cas--) read(),go();
    return 0;
}