fzu acm 1015 土地划分

在Dukeswood这块土地上生活着一个富有的农庄主和他的几个孩子。在他临终时，他想把他的土地分给他的孩子。他有许多农场，每个农场都是一块矩形土地。他在农场地图上划上一些直线将矩形分成若干块。当他划直线时，他总是从矩形边界上的某一点划到另一个矩形边界上的点，这条线的结束点将成为下一条线的起始点。他划线时从不会让任三线共点。例如图1是某一种划分结果。

 

图1

 

划分的起始点和结束点均以五角星标记。当他完成划分后，他想要数一下划出的土地的块数以确保每个孩子都有一块地。例如，图1中土地被划分成18块。然而这个庄主由于年迈常会数错，因而他寻求你的帮助。

请写一个程序，输入原来的土地尺寸及线段的位置，输出划分出的土地块数。

输入输出格式

输入文件有多组数据组成。每组数据格式如下：
第一行输入地图的宽度w (1<=w<=1000)和高度 h (1<=h<=1000)，均为整数。
第二行输入线段数L (1<=L<=50)。
以下L+1行每行一个整数坐标(Xi,Yi)，庄主划的线段为(Xi,Yi)-(Xi+1,Yi+1)，i=1,2,…,L。当然(Xi,Yi)必定在矩形的边界上。
最后一组数据w=h=0，标志文件结束，不需要处理。
对于给定的输入，输出一行仅含一个数，即划分出的土地块数。

样例输入

18 12
8
2 0
6 12
10 0
18 9
15 12
0 6
14 0
10 12
0 9
7 6
6
2 0
5 6
7 3
0 3
3 0
3 6
0 5
0 0

样例输出

18
11

分析：此题一看就是计算几何题，要是不熟悉还真是无从下手，不知道怎么办。其实，这就是考我们判断线段焦点怎么判断

这就是用到数学上的向量叉乘。就是判断当前线段与之前的线段有几个交点，结果就是增加结点个数加一块。有了思路这题

接下来就是要认真的，耐心滴写下去……

#include<iostream>
#include<stdio.h>
using namespace std;

struct point //定义点结构体
{			//先接收所有的点
	double x;
	double y;
}Point[60];

struct vector
{//定义向量，存放边
	point a;
	point b;
}Vector[60];

bool VexJudge(vector Vec1,vector Vec2)
{//用于判断两条线段是否相交，当然这边只是判断点在两侧的情况实际是不太全面的，但这题来讲
	//是没问题的，因题目而定……
	int a,b;
	//行列式相乘
	a=((Vec2.a.x-Vec1.a.x)*(Vec1.b.y-Vec1.a.y)-(Vec1.b.x-Vec1.a.x)*(Vec2.a.y-Vec1.a.y))*
      ((Vec2.b.x-Vec1.a.x)*(Vec1.b.y-Vec1.a.y)-(Vec1.b.x-Vec1.a.x)*(Vec2.b.y-Vec1.a.y));
	b=((Vec1.b.x-Vec2.a.x)*(Vec2.b.y-Vec2.a.y)-(Vec2.b.x-Vec2.a.x)*(Vec1.b.y-Vec2.a.y))*
	  ((Vec1.a.x-Vec2.a.x)*(Vec2.b.y-Vec2.a.y)-(Vec2.b.x-Vec2.a.x)*(Vec1.a.y-Vec2.a.y));
	if(a<0&&b<0)
		return true;//相交返回真
	else 
		return false;
}

int function(vector Vec,int n)
{
	int i,number=0;
	for(i=0;i<n;i++)
	{//用这个循环，实现Vec向量与
	//之前的每条线段都判断一次
	//看有没有相交
		if(VexJudge(Vector[i],Vec))
			++number;	
	}
	return number;//返回相交的个数
}

int main()
{
	freopen("e://test.txt","r",stdin);
	int with,high,n,i,sum;
	while(scanf("%d%d",&with,&high)&&(with!=0||high!=0))
	{
		sum=0;
		scanf("%d",&n);
		for(i=0;i<=n;i++)//先用点结构体接收
			scanf("%lf%lf",&Point[i].x,&Point[i].y);
		if(n==1)//n为1，2时分开考虑
			sum=2;
		else if(n==2)
			sum=3;
		else
		{
			for(i=0;i<n;i++)//把点转换成向量
			{
				Vector[i].a.x=Point[i].x;
				Vector[i].a.y=Point[i].y;
				Vector[i].b.x=Point[i+1].x;
				Vector[i].b.y=Point[i+1].y;
			}
			sum=3;
			for(i=2;i<n;i++)//计算与之前的线段有没相交
			{//并计算个数
				sum+=function(Vector[i],i)+1;
			}
		}
		cout<<sum<<endl;
	}
	return 0;
}