题意:求n条线段的所有交点个数. 题目保证线段不会有重合部分.
直接求交点, 把所有的交点扔进一个集合统计一下就行了.
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <iostream>
#include <vector>
#include <set>
using namespace std;
const double eps = 1e-8;
const double INF = 1e20;
const double pi = acos (-1.0);
int dcmp (double x) {
if (fabs (x) < eps) return 0;
return (x < 0 ? -1 : 1);
}
inline double sqr (double x) {return x*x;}
struct Point {
double x, y;
Point (double _x = 0, double _y = 0):x(_x), y(_y) {}
void input () {scanf ("%lf%lf", &x, &y);}
void output () {printf ("%.2f %.2f\n", x, y);}
bool operator == (const Point &b) const {
return (dcmp (x-b.x) == 0 && dcmp (y-b.y) == 0);
}
bool operator < (const Point &b) const {
return (dcmp (x-b.x) == 0 ? dcmp (y-b.y) < 0 : x < b.x);
}
Point operator + (const Point &b) const {
return Point (x+b.x, y+b.y);
}
Point operator - (const Point &b) const {
return Point (x-b.x, y-b.y);
}
Point operator * (double a) {
return Point (x*a, y*a);
}
Point operator / (double a) {
return Point (x/a, y/a);
}
double len2 () {
return sqr (x) + sqr (y);
}
double len () {
return sqrt (len2 ());
}
Point change_len (double r) {
double l = len ();
if (dcmp (l) == 0) return *this;
r /= l;
return Point (x*r, y*r);
}
Point rotate_left () {
return Point (-y, x);
}
Point rotate_right () {
return Point (y, -x);
}
Point rotate (Point p, double ang) {
Point v = (*this)-p;
double c = cos (ang), s = sin (ang);
return Point (p.x + v.x*c - v.y*s, p.y + v.x*s + v.y*c);
}
Point normal () {
double l = len ();
return Point (-y/l, x/l);
}
};
double cross (Point a, Point b) {
return a.x*b.y-a.y*b.x;
}
double dot (Point a, Point b) {
return a.x*b.x + a.y*b.y;
}
double dis (Point a, Point b) {
Point p = b-a; return p.len ();
}
double rad_degree (double rad) {
return rad/pi*180;
}
double rad (Point a, Point b) {
return fabs (atan2 (fabs (cross (a, b)), dot (a, b)) );
}
bool parallel (Point a, Point b) {
double p = rad (a, b);
return dcmp (p) == 0 || dcmp (p-pi) == 0;
}
struct Line {
Point s, e;
Line () {}
Line (Point _s, Point _e) : s(_s), e(_e) {}
Line (Point p, double ang) {
s = p;
if (dcmp (ang-pi/2) == 0) {
e = s + Point (0, 1);
}
else
e = s + Point (1, tan (ang));
}
Line (double a, double b, double c) {
if (dcmp (a) == 0) {
s = Point (0, -c/b);
e = Point (1, -c/b);
}
else if (dcmp (b) == 0) {
s = Point (-c/a, 0);
e = Point (-c/a, 1);
}
else {
s = Point (0, -c/b);
e = Point (1, (-c-a)/b);
}
}
void input () {
s.input ();
e.input ();
}
void adjust () {
if (e < s) swap (e, s);
}
double length () {
return dis (s, e);
}
double angle () {
double k = atan2 (e.y-s.y, e.x-s.x);
if (dcmp (k) < 0) k += pi;
if (dcmp (k-pi) == 0) k -= pi;
return k;
}
};
int relation (Point p, Line l) {
int c = dcmp (cross (p-l.s, l.e-l.s));
if (c < 0) return 1;
else if (c > 0) return 2;
else return 3;
}
bool point_on_seg (Point p, Line l) {
return dcmp (cross (p-l.s, l.e-l.s)) == 0 &&
dcmp (dot (p-l.s, p-l.e) <= 0);
}
bool parallel (Line a, Line b) {
return parallel (a.e-a.s, b.e-b.s);
}
int seg_cross_seg (Line a, Line v) {
int d1 = dcmp (cross (a.e-a.s, v.s-a.s));
int d2 = dcmp (cross (a.e-a.s, v.e-a.s));
int d3 = dcmp (cross (v.e-v.s, a.s-v.s));
int d4 = dcmp (cross (v.e-v.s, a.e-v.s));
if ((d1^d2) == -2 && (d3^d4) == -2) return 2;
return (d1 == 0 && dcmp (dot (v.s-a.s, v.s-a.e)) <= 0) ||
(d2 == 0 && dcmp (dot (v.e-a.s, v.e-a.e)) <= 0) ||
(d3 == 0 && dcmp (dot (a.s-v.s, a.s-v.e)) <= 0) ||
(d4 == 0 && dcmp (dot (a.e-v.s, a.e-v.e)) <= 0);
}
int line_cross_seg (Line a, Line v) {
int d1 = dcmp (cross (a.e-a.s, v.s-a.s));
int d2 = dcmp (cross (a.e-a.s, v.e-a.s));
if (d1^d2 == -2) return 2;
return (d1 == 0 || d2 == 0);
}
int line_cross_line (Line a, Line v) {
if (parallel (a, v)) return relation (a.e, v) == 3;
return 2;
}
Point line_intersection (Line a, Line v) {
double a1 = cross (v.e-v.s, a.s-v.s);
double a2 = cross (v.e-v.s, a.e-v.s);
return Point ((a.s.x*a2-a.e.x*a1)/(a2-a1), (a.s.y*a2-a.e.y*a1)/(a2-a1));
}
int seg_intersectiong (Line a, Line b, Point &p) {
int rel = seg_cross_seg (a, b);
if (rel == 0) return 0;
int cnt = 0;
if (rel == 1) {
if (point_on_seg (a.e, b)) p = a.e, cnt++;
if (point_on_seg (a.s, b)) p = a.s, cnt++;
if (point_on_seg (b.e, a)) p = b.e, cnt++;
if (point_on_seg (b.s, a)) p = b.s, cnt++;
return 2;
}
p = line_intersection (a, b);
return 1;
}
double point_to_line (Point p, Line a) {
return fabs (cross (p-a.s, a.e-a.s) / a.length ());
}
double point_to_seg (Point p, Line a) {
if (dcmp (dot (p-a.s, a.e-a.s)) < 0 || dcmp (dot (p-a.e, a.s-a.e)) < 0)
return min (dis (p, a.e), dis (p, a.s));
return point_to_line (p, a);
}
Point projection (Point p, Line a) {
return a.s + (((a.e-a.s) * dot (a.e-a.s, p-a.s)) / (a.e-a.s).len2() );
}
Point symmetry (Point p, Line a) {
Point q = projection (p, a);
return Point (2*q.x-p.x, 2*q.y-p.y);
}
#define maxn 111
int n;
Line l[maxn];
set <Point> gg;
int main () {
while (cin >> n && n) {
gg.clear ();
for (int i = 0; i < n; i++) {
l[i].input ();
}
int ans = 0;
Point p;
for (int i = 0; i < n; i++) {
for (int j = i+1; j < n; j++) {
if (seg_intersectiong (l[i], l[j], p)) {
gg.insert (p);
}
}
}
cout << gg.size () << "\n";
}
return 0;
}