lightoj 1320

输入：

给组数；

给边数N，给X最大，给Y最大，

给N条边，X1Y1，X2Y2;

问区域内有几个块。

思路清晰，对于每条边，有自己的影响，影响就是和之前的边有多少交点。有一个交点加一块，自身是一块。

但是有个技巧，就是浮点数的点，怎么存；

发现直线相交的交点是

X=（l2.b*l1.c - l1.b*l2.c）/det；

Y=（l1.a*l2.c - l2.a*l1.c）/det；

所以存的时候直接存个PAIR表示点对，注意GCD，否则会重点；

由此而来的就是对于边界的比较。很机智。

#include<bits/stdc++.h>
using namespace std;

int n,L,W;

struct line{
int a,b,c;
line(){}
line(int x1,int y1,int x2, int y2){
a = y1 - y2;
b = x2 - x1;
c = a*x1 + b*y1;
}
};

struct point{
int x,y;
point(){}
point(int a): x(a),y(1){}
point(int a, int b){
x = a/__gcd(a,b);
y = b/__gcd(a,b);
}
bool operator < (const point &o) const{
return x*o.y < y*o.x;
}

};

bool intersection(line l1,line l2, point &x,point &y){
int det = l1.a*l2.b - l2.a*l1.b;
if(det == 0) // lines are parallel
return false;
x = point(l2.b*l1.c - l1.b*l2.c,det);
y = point(l1.a*l2.c - l2.a*l1.c,det);
return true;
}

vector<line> lines;

int solve(){
point minx(0);
point miny(0);
point maxx(L);
point maxy(W);
int ans = 1;
for(int i = 0; i< n; ++i){
set<pair<point,point> > points;
for(int j = 0; j< i; ++j){
point x,y;
if(intersection(lines[i],lines[j],x,y)){
if(minx < x and x < maxx and miny < y and y < maxy)
points.insert(make_pair(x,y));
}
}
ans += points.size() + 1;
}
return ans;
}

int main(){
int tc;cin>>tc;
int a,b,c,d;
for(int tid = 1; tid <= tc; ++tid){
cin>>n>>L>>W;
lines.clear();
for(int i = 0 ; i< n; ++i){
cin>>a>>b>>c>>d;
lines.push_back(line(a,b,c,d));
}
printf("Case %d: %d\n",tid,solve());

}

return 0;
}