POJ 3358 凸包求面积
2013-07-26 16:16
246 查看
题目:题目链接
题目的意思就是给你N个点,然后让你求出包含所有树木在内的凸多边形,这就是凸包嘛。然后就是让你求出这个凸包的面积。50可以存放一只牛。问可以放多少牛。
分析:这中凸包的题目都大体不变的,直接上模板就OK了:
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <iostream>
using namespace std;
const int N = 1005;
const double PI = 3.1415927;
struct point
{
double x;
double y;
} p
, stack
;
double dis(point A, point B)
{
return sqrt((A.x-B.x)*(A.x-B.x)+(A.y-B.y)*(A.y-B.y));
}
double crossProd(point A, point B, point C)
{
return (B.x-A.x)*(C.y-A.y) - (B.y-A.y)*(C.x-A.x);
}
//以最左下的点为基准点,其他各点(逆时针方向)以极角从小到大的排序规则
int cmp(const void *a, const void *b)
{
point *c = (point *)a;
point *d = (point *)b;
double k = crossProd(p[0], *c, *d);//极角大小转化为求叉乘
if (k<0 || !k && dis(p[0], *c)>dis(p[0], *d)) return 1;
return -1;
}
double Graham(int n)
{
double x = p[0].x;
double y = p[0].y;
int mi = 0;
for (int i=1; i<n; ++i) //找到最左下的一个点
{
if (p[i].x<x || (p[i].x==x && p[i].y<y))
{
x = p[i].x;
y = p[i].y;
mi = i;
}
}
point tmp = p[mi];
p[mi] = p[0];
p[0] = tmp;
qsort(p+1, n-1, sizeof(point), cmp);
p
= p[0];
stack[0] = p[0];
stack[1] = p[1];
stack[2] = p[2];
int top = 2;
for (int i=3; i<=n; ++i) //加入一个点后,向右偏拐或共线,则上一个点不在凸包内,则--top,该过程直到不向右偏拐或没有三点共线的点
{
while (crossProd(stack[top-1], stack[top], p[i])<=0 && top>=2) --top;
stack[++top] = p[i];//在当前情况下符合凸包的点,入栈
}
double len = 0;
for(int i = 0; i <= top; ++i)
len +=(double)(stack[i].x*stack[(i+1)%(top+1)].y-stack[i].y*stack[(i+1)%(top+1)].x);
return len/2.0;
}
int main()
{
int n;
scanf ("%d", &n);
for (int i=0; i<n; ++i)
scanf ("%lf%lf", &p[i].x, &p[i].y);
double ans = Graham(n);
int ap = ans/50;
printf ("%d\n", ap);
return 0;
}
努力努力...
题目的意思就是给你N个点,然后让你求出包含所有树木在内的凸多边形,这就是凸包嘛。然后就是让你求出这个凸包的面积。50可以存放一只牛。问可以放多少牛。
分析:这中凸包的题目都大体不变的,直接上模板就OK了:
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <iostream>
using namespace std;
const int N = 1005;
const double PI = 3.1415927;
struct point
{
double x;
double y;
} p
, stack
;
double dis(point A, point B)
{
return sqrt((A.x-B.x)*(A.x-B.x)+(A.y-B.y)*(A.y-B.y));
}
double crossProd(point A, point B, point C)
{
return (B.x-A.x)*(C.y-A.y) - (B.y-A.y)*(C.x-A.x);
}
//以最左下的点为基准点,其他各点(逆时针方向)以极角从小到大的排序规则
int cmp(const void *a, const void *b)
{
point *c = (point *)a;
point *d = (point *)b;
double k = crossProd(p[0], *c, *d);//极角大小转化为求叉乘
if (k<0 || !k && dis(p[0], *c)>dis(p[0], *d)) return 1;
return -1;
}
double Graham(int n)
{
double x = p[0].x;
double y = p[0].y;
int mi = 0;
for (int i=1; i<n; ++i) //找到最左下的一个点
{
if (p[i].x<x || (p[i].x==x && p[i].y<y))
{
x = p[i].x;
y = p[i].y;
mi = i;
}
}
point tmp = p[mi];
p[mi] = p[0];
p[0] = tmp;
qsort(p+1, n-1, sizeof(point), cmp);
p
= p[0];
stack[0] = p[0];
stack[1] = p[1];
stack[2] = p[2];
int top = 2;
for (int i=3; i<=n; ++i) //加入一个点后,向右偏拐或共线,则上一个点不在凸包内,则--top,该过程直到不向右偏拐或没有三点共线的点
{
while (crossProd(stack[top-1], stack[top], p[i])<=0 && top>=2) --top;
stack[++top] = p[i];//在当前情况下符合凸包的点,入栈
}
double len = 0;
for(int i = 0; i <= top; ++i)
len +=(double)(stack[i].x*stack[(i+1)%(top+1)].y-stack[i].y*stack[(i+1)%(top+1)].x);
return len/2.0;
}
int main()
{
int n;
scanf ("%d", &n);
for (int i=0; i<n; ++i)
scanf ("%lf%lf", &p[i].x, &p[i].y);
double ans = Graham(n);
int ap = ans/50;
printf ("%d\n", ap);
return 0;
}
努力努力...
相关文章推荐
- poj2079凸包求最大三角形面积
- POJ 3348 Cows 寻找凸包 凸包面积
- POJ 3348 Cows (凸包面积)
- POJ 3348 Cows(凸包面积)
- Cows POJ - 3348 凸包面积
- poj 3348 Cows 求凸包以及凸包的面积
- POJ 3348 Cows(凸包面积)
- hdu 3934&&poj 2079 (凸包+旋转卡壳+求最大三角形面积)
- POJ 3348 求凸包面积
- POJ 3348 Cows(凸包+多边形面积)
- [poj] 3348 Cows || 求凸包面积
- Poj 3348 Cows (凸包面积)
- poj 3348 Cows(求凸包面积)
- POJ 3348 求凸包面积
- POJ 2079 求最大三角形面积 (凸包+旋转卡壳)
- poj 3348(凸包面积)
- POJ 3348 Cows 求凸包面积
- POJ 3348 || Cows(求凸包面积
- POJ 3348 Cows 凸包 + 多边形面积公式
- POJ 3348 Cows (计算凸包面积)