hdu_1086 You can Solve a Geometry Problem too

题意：给定n条线段（n≤100)，求有多少个交点（若有多于两条线段交于一点，应重复计数（即可以视为任意两条线段不交于同一点））。

题解：判定线段相交的裸题。只需要枚举任意两条线段是否相交即可。（方法为判断线段AB的端点A,B是否在线段CD两侧，且C,D在线段AB两侧。注意特判一个端点在另一条线段上的情况。）

代码：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define eps 1e-15
using namespace std;
int n,t,ans;
struct vec{
    double x,y;
};
struct node{
    vec a,b;
}e[110];
double che(vec a,vec b){
    return a.x*b.y-a.y*b.x;
}
double ab(vec a){
    return a.x*a.x+a.y*a.y;
}
vec miu(vec a,vec b){
    return (vec){a.x-b.x,a.y-b.y};
}
bool on(node a,node b){
    double x=ab(miu(b.a,a.a)),y=ab(miu(a.b,b.a)),z=ab(miu(a.b,a.a));
    if(x+y==z)  return true;
    x=ab(miu(b.b,a.a));y=ab(miu(a.b,b.b));
    if(x+y==z)  return true;
    return false;
}
bool pan(node a,node b){
    vec aa=miu(a.b,a.a);
    double x=che(miu(a.b,a.a),miu(b.a,a.a))*che(miu(a.b,a.a),miu(b.b,a.a));
    double y=che(miu(b.b,b.a),miu(a.a,b.a))*che(miu(b.b,b.a),miu(a.b,b.a));
    if(x<0 && y<0)  return true;
    if(on(a,b))  return true;
    return false;
}
int main(){
    int i,j,k,l;
    scanf("%d",&n);
    while(n!=0){
        ans=0;memset(e,0,sizeof(e));
        for(i=1;i<=n;++i){
            scanf("%lf%lf%lf%lf",&e[i].a.x,&e[i].a.y,&e[i].b.x,&e[i].b.y);
            for(j=1;j<i;++j){
                if(pan(e[j],e[i]))  ans++;
            }
        }
        printf("%d\n",ans);
        scanf("%d",&n);
    }
    return 0;
}