poj1265（Area）

题目地址：Area

 

题目大意：

     给你m个点x，y坐标的变化大小，形成一个封闭的图形，求出在图形上点的数目和在多边形内点的个数（每个点的距离间隔为“1”）和形成区域面积的大小。

 

解题思路：

    因为要求在多边形的边上和内部的点，想到pick定理：

pick定理：
                设F为平面上以格子点为定点的单纯多边形，则其面积为：S=on/2+in-1。 

                on为多边形边上点格点的个数，in为多边形内部格点的个数。 

                可用其计算多边形的面积，边界格点数或内部格点数。

　　首先可以通过叉积来求封闭多边形的面积，然后通过坐标变化的最大公约数可以确定移动到下一顶点时途经的顶点个数，最后通过pick定理求其内部的顶点。

 

代码：

#include <algorithm>
#include <iostream>
#include <sstream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <string>
#include <bitset>
#include <vector>
#include <queue>
#include <stack>
#include <cmath>
#include <list>
//#include <map>
#include <set>
using namespace std;
/***************************************/
#define ll long long
#define int64 __int64
/***************************************/
const int INF = 0x7f7f7f7f;
const double eps = 1e-8;
const double PIE=acos(-1.0);
const int d1x[]= {0,-1,0,1};
const int d1y[]= {-1,0,1,0};
const int d2x[]= {0,-1,0,1};
const int d2y[]= {1,0,-1,0};
const int fx[]= {-1,-1,-1,0,0,1,1,1};
const int fy[]= {-1,0,1,-1,1,-1,0,1};
/*vector <int>map[N];map[a].push_back(b);int len=map[v].size();*/
/***************************************/
void openfile()
{
    freopen("data.in","rb",stdin);
    freopen("data.out","wb",stdout);
}
priority_queue<int> qi1;
priority_queue<int, vector<int>, greater<int> >qi2;
/**********************华丽丽的分割线,以上为模板部分*****************/

struct node
{
    double x,y;
}p[200];
int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}
int EPS(double f)
{
    if (fabs(f)<eps)
        return 0;
    return f>0?1:-1;
}
double cha(node b1,node a1)
{
    return a1.x*b1.y-a1.y*b1.x;
}
int main()
{

    int cas,cnt=1;
    scanf("%d",&cas);
    while(cas--)
    {
        int n;
        memset(p,0,sizeof(p));
        scanf("%d",&n);
        double xx,yy;
        double area=0;
        int onp=0;
        int i;
        for(i=1;i<=n;i++)
        {
            scanf("%lf%lf",&xx,&yy);
            p[i].x=xx+p[i-1].x;
            p[i].y=yy+p[i-1].y;
            onp+=fabs(gcd(xx,yy));
            area+=cha(p[i],p[i-1]);
        }
        int inp=area/2+1-onp/2;
        printf("Scenario #%d:\n%d %d %.1f\n\n",cnt++,inp,onp,area/2.0);
    }
    return 0;
}