POJ 4048 金华赛区E题

这题算是给自己敲了个警钟，模板必须高度可靠，原来我在各个OJ都能过题的模板，在现场塞被卡的很惨。。。

算法导论上的模板才是王道啊。。。（就差那么一点，就就，唉～～～～～～～～～～～～～～～～～）从此以后再不相信网上各种模板了。

还是贴代码吧。

#include<iostream>
#include<cstdio>
#include<algorithm>
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define MAXN 1504
typedef double datatype;
using namespace std;
//typedef complex<double> point;
struct point {                                             
	datatype x, y;
	point(){}
	point(datatype _x, datatype _y):x(_x),y(_y){}
	point(const point &temp)
	{
		x = temp.x;
		y = temp.y;
	}
}fir[MAXN],sec[MAXN];

datatype cmp(datatype x)
{
  if (x > 0)return 1;
  if (x < 0)return -1;
  return 0;  
}

point start,end;

datatype cross(const point &p1, const point &p2, const point &q1, const point &q2)
{
  return cmp((q2.y - q1.y)*(p2.x - p1.x) - (q2.x - q1.x)*(p2.y - p1.y));  
}

bool on_seg(const point &p1, const point &p2, const point &q)
{
  if (min(p1.x,p2.x) <= q.x && q.x <= max(p1.x,p2.y) &&
      min(p1.y,p2.y) <= q.y && q.y <= max(p1.y,p2.y) && cross(p1,p2,p1,q) == 0) {
    return true;      
  }
  return false;  
}

bool intersect(const point &p1, const point &p2, const point &q1, const point &q2)
{
  if (cross(p1,q1,p1,p2)*cross(p1,p2,p1,q2) > 0 && cross(q1,p1,q1,q2)*cross(q1,q2,q1,p2) > 0)
      return true;
  return on_seg(p1,p2,q1) || on_seg(p1,p2,q2) || on_seg(q1,q2,p1) || on_seg(q1,q2,p2);    
}

int get_cnt(int i, int n)
{
  int temp(1);
  for (int j(0); j<n; ++j) {
    if (i == j)continue;
    if (intersect(start,end,fir[j],sec[j]))++temp;  
  }
  return temp;      
}

void get_end(const point &p)
{
  end.x = start.x + (p.x - start.x)*100000.0;
  end.y = start.y + (p.y - start.y)*100000.0;  
}

int main()
{
  int t;
  cin>>t;
  while (t--) {
    int n;
    cin>>n;
    for (int i(0); i<n; ++i) {
      cin>>fir[i].x>>fir[i].y>>sec[i].x>>sec[i].y;  
    }
    cin>>start.x>>start.y;
    int cnt(1);
    for (int i(0); i<n; ++i) {
      get_end(fir[i]);
      cnt = max(cnt,get_cnt(i,n));
      get_end(sec[i]);
      cnt = max(cnt,get_cnt(i,n));
    }
    cout<<cnt<<endl;
  }
  return 0;
}