/************************************************************************************
1、从第三条边开始,判断编号小于这条边(不相邻)的所有边的两个端点与该边的位置关系
并判断该边的两个端点与编号小于该边的所有边的位置关系
2、如果两个端点分别在某边的上下区域,则把两端点坐标代入直线方程后相乘结果必小于零
3、面积s即任意两端点与原点的向量叉乘的行列式之和,化简后如程序中所示,若小于零取绝对值
**************************************************************************************/
#include
#include
using namespace std;
double judgeX2y2UpOrDown(double x1, double y1, double x2, double y2, double x3, double y3)
{
return y2 * x3 - y2 * x1 - y1 * x3 - x2 * y3 + x2 * y1 + x1 * y3;//判断(x2,y2)与通过(x1,y1),(x3,y3)的直线的位置关系
}
int main()
{
int i, j, n, t, cas = 1;
bool flag;
double s, v[1111][2];
while(scanf("%d", &n), n)
{
flag = 0;
s = 0;
for(i = 0; i < n; i++)
scanf("%lf%lf", &v[i][0], &v[i][1]);
printf("Figure %d: ", cas++);//冒号后面要有空格,否则WA
if(n < 3)
printf("Impossible\n\n");
else
{
for(i = 0; i < n; i++)
{
for(j = (i + 1) / n; j < i - 1; j++)
{
t = (i + 1) % n;
if(judgeX2y2UpOrDown(v[j][0], v[j][1], v[i][0], v[i][1], v[j + 1][0], v[j + 1][1]) *
judgeX2y2UpOrDown(v[j][0], v[j][1], v[t][0], v[t][1], v[j + 1][0], v[j + 1][1]) <= 0 &&
judgeX2y2UpOrDown(v[i][0], v[i][1], v[j][0], v[j][1], v[t][0], v[t][1]) *
judgeX2y2UpOrDown(v[i][0], v[i][1], v[j + 1][0], v[j + 1][1], v[t][0], v[t][1]) <= 0)
{
printf("Impossible\n");
flag = 1;
break;
}
}
if(flag)
break;
}
if(!flag)
{
for(i = 0; i < n; i++)
{
s += v[i][0] * v[(i + 1) % n][1];
s -= v[(i + 1) % n][0] * v[i][1];
}
if(s < 0)
s *= -1;
printf("%.2f\n\n", s / 2);
}
}
}
return 0;
}
zoj1010 Area
最新推荐文章于 2019-07-07 14:03:48 发布