POJ1279-求解多边形内核面积

题目：题目链接

 

题意：题目意思就是给你一个多边形，然后让你求出内核的面积。

 

分析：还是求内核的问题，先按照求内核的方法求出点，然后又是一个多边形求面积。对于顺时针还是逆时针，最后的

面积判个绝对值就可以了：

 

代码：

 

#include <iostream>
#include <cstdio>
#include <string>
#include <string.h>
#include <map>
#include <vector>
#include <cstdlib>
#include <algorithm>
#include <cmath>
#include <queue>
#include <set>
#include <stack>
#include <functional>
#include <fstream>
#include <sstream>
#include <iomanip>
#include <numeric>
#include <cassert>
#include <bitset>
#include <stack>
#include <ctime>
#include <list>
#define INF 0x7fffffff
#define max3(a,b,c) (max(a,b)>c?max(a,b):c)
#define min3(a,b,c) (min(a,b)<c?min(a,b):c)
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;

#define exp 1e-10

struct node
{
    double x;
    double y;
};

node point[1505];//记录最开始的多边形
node q[1505]; //临时保存新切割的多边形
node p[1505]; //保存新切割出的多边形
int n,m;//n的原先的点数,m是新切割出的多边形的点数
double a,b,c;

void getline(node x,node y)  //获取直线ax+by+c==0
{
    a=y.y-x.y;
    b=x.x-y.x;
    c=y.x*x.y-x.x*y.y;
}

node intersect(node x,node y) //获取直线ax+by+c==0  和点x和y所连直线的交点
{
    double u=fabs(a*x.x+b*x.y+c);
    double v=fabs(a*y.x+b*y.y+c);
    node ans;
    ans.x=(x.x*v+y.x*u)/(u+v);
    ans.y=(x.y*v+y.y*u)/(u+v);
    return ans;
}

void cut()                            //用直线ax+by+c==0切割多边形
{
    int cutm=0,i;
    for(i=1; i<=m; i++)
    {
        if(a*p[i].x + b*p[i].y+c>=0)  //暂时定位顺时钟给出点的
        {
            //所以一个点在直线右边的话，那么带入值就会大于等于0
            q[++cutm]=p[i];           //说明这个点还在切割后的多边形内，将其保留
        }
        else
        {
            if(a*p[i-1].x+b*p[i-1].y+c>0) //该点不在多边形内，但是它和它相邻的点构成直线与
            {
                //ax+by+c==0所构成的交点可能在新切割出的多边形内，
                q[++cutm]=intersect(p[i-1],p[i]); //所以保留交点
            }
            if(a*p[i+1].x+b*p[i+1].y+c>0)
            {
                q[++cutm]=intersect(p[i+1],p[i]);
            }
        }
    }
    for(i = 1; i <= cutm; i++)
        p[i] = q[i];
    p[cutm+1] = q[1];
    p[0] = q[cutm];
    m = cutm;
}
//求核的面积
double multi(node p0, node p1, node p2)
{
    return (p1.x-p0.x)*(p2.y-p0.y)-(p1.y-p0.y)*(p2.x-p0.x);
}

double solve()
{
    int i;
    for(i = 1; i <= n; i++)
        p[i] = point[i];
    point[n+1] = point[1];
    p[n+1]=p[1];
    p[0]=p[n];
    m = n;
    for(i=1; i<=n; i++)
    {
        getline(point[i], point[i+1]); //根据point[i]和point[i+1]确定直线ax+by+c==0
        cut();                         //用直线ax+by+c==0切割多边形
    }
    if (m < 3)
        return 0;
    double area = 0;
    for (int i = 1; i < m-1; i++)
        area += multi(p[0], p[i], p[i+1]);
    if (area < 0) area = -area;
    return area/2;
}

int main()
{
    int cas,i;
    scanf("%d",&cas);
    while(cas--)
    {
        scanf("%d",&n);
        for(i=1; i<=n; i++)
            scanf("%lf%lf",&point[i].x,&point[i].y);
        double ans = solve();
        printf("%.2lf\n", ans);
    }
    return 0;
}