hdu 1086 判断交点

对于两条直线上的点， M(ax1,ay1) N (ax2,ay2) P(bx1,by1) Q(bx2,by2)

若两直线相交，则任意一条直线的两个端点，都在另一条直线的两侧。

我们可以使用向量外积 或者 线性规划 来判断。

此处使用外积。

MN 与 MP 的外积 MNP = (ax1 - ax2)*(ay1 - by1) - (ay1 - ay2)*(ax1 - bx1)

MN 与 MQ 的外积 MNQ = (ax1 - ax2)*(ay1 - by2) - (ay1 - ay2)*(ax1 - bx2)

若MNP 和 MNQ 异号，则点在线的两侧。

MNP 和 MNQ 的积 为0，则有点在线上，因此也算作有交点。

#include <stdio.h>
#include <stdlib.h>
struct node
{
double x,y;
}s1[105],s2[105];
int ans;
void wtf(int a,int b)
{
double abc = (s1[a].x - s2[a].x)*(s1[a].y - s1[b].y) - (s1[a].y - s2[a].y)*(s1[a].x - s1[b].x);
double abd = (s1[a].x - s2[a].x)*(s1[a].y - s2[b].y) - (s1[a].y - s2[a].y)*(s1[a].x - s2[b].x);
double cda = (s1[b].x - s2[b].x)*(s1[b].y - s2[a].y) - (s1[b].y - s2[b].y)*(s1[b].x - s2[a].x);
double cdb = (s1[b].x - s2[b].x)*(s1[b].y - s1[a].y) - (s1[b].y - s2[b].y)*(s1[b].x - s1[a].x);
if((abc * abd) <= 0 && (cda * cdb) <= 0) ans++;
}
int main()
{
//	freopen("t.txt","r",stdin);
int n,i,j;
while(~scanf("%d",&n) && n)
{
ans = 0;
for(i = 0; i < n; i++)
{
scanf("%lf%lf%lf%lf",&s1[i].x,&s1[i].y,&s2[i].x,&s2[i].y);
}
for(i = 1; i < n; i++)
{
for(j = 0; j < i; j++)
{
wtf(i,j);
}
}
printf("%d\n",ans);
}
return 0;
}