#include #include #include using namespace atcoder; using mint = modint1000000007; using namespace std; #define rep(i,n) for(int i=0;i<(n);i++) #define Inf 1000000001 const double eps = 1e-10; template struct point{ T x,y,rad,dir; point(T a=0.0,T b=0.0){ x = a; y = b; fix_rd(); } void update_x(T a,T b){ x = a*x + b; fix_rd(); } void update_x(T a){ update_x(0.0,a); } void update_y(T a,T b){ y = a*y + b; fix_rd(); } void update_y(T a){ update_y(0.0,a); } void update_rad(T a,T b){ rad = a*rad + b; fix_xy(); } void update_rad(T a){ update_rad(0.0,a); } void update_dir(T a,T b){ dir = a*dir + b; fix_xy(); } void update_dir(T a){ update_dir(0.0,a); } void fix_xy(){ x = rad * cos(dir); y = rad * sin(dir); fix_rd(); } void fix_rd(){ rad = hypot(x,y); if(rad==0.0)dir=0.0; else dir = atan2(y,x); fix_zero(); } void fix_zero(){ if(abs(x) V){ return hypot(x-V.x,y-V.y); } T size(){ return get_dis(point()); } T angle_difference(point V){ double ret = dir - V.dir; if(ret<-acos(-1.0))ret = acos(-1.0)*2.0+ret; if(ret>acos(-1.0))ret=-acos(-1.0)*2.0+ret; return ret; } //中点 point get_midpoint(point V){ V.update_x(0.5,x/2.0); V.update_y(0.5,y/2.0); return V; } T get_inner_product(point V){ return x*V.x+y*V.y; } T get_cross_product(point V){ return x*V.y-y*V.x; } point &operator+=(const point &another){ update_x(1,another.x); update_y(1,another.y); return (*this); } point &operator-=(const point &another){ update_x(1,-another.x); update_y(1,-another.y); return (*this); } point operator+(const point &another)const{ return (point(*this)+=another); } point operator-(const point &another)const{ return (point(*this)-=another); } void show(){ cout< struct line{ point a,t; line(point V1,point V2){ a=V1; t=V2-V1; } T get_signed_dis(point V){ point PA = a-V; return PA.get_cross_product(t)/t.size(); } T get_dis(point V){ return abs(get_signed_dis(V)); } point get_projection(point P){ T r = t.get_inner_product(P-a)/t.size(); point temp = t; temp.update_rad(0.0,r); return a+temp; } point get_cross_point(line L){ point ret(1e20,1e20); if(abs(t.get_cross_product(L.t)) temp = t; temp.x *= d0/(d1-d0); temp.y *= d0/(d1-d0); ret = a - temp; return ret; } }; template struct segment{ point V1,V2; segment(point a=point(),point b=point()){ V1=a; V2=b; } T get_dis(point P){ T ret = 1e20; line L(V1,V2); point Q = L.get_projection(P); if(Q.x+eps>min(V1.x,V2.x)&&Q.y+eps>min(V1.y,V2.y) &&Q.x l){ if(get_cross_point(l).x<1e20)return 0.0; return min({get_dis(l.V1),get_dis(l.V2),l.get_dis(V1),l.get_dis(V2)}); } point get_cross_point(segment l){ line L1(V1,V2),L2(l.V1,l.V2); point P = L1.get_cross_point(L2); if(get_dis(P) (1e20,1e20); } }; int main(){ int n; cin>>n; vector> S; vector> t; rep(i,n){ double x,y;cin>>x>>y; point a(x,y); cin>>x>>y; point b(x,y); S.push_back(segment(a,b)); t.push_back(a); t.push_back(b); } int ans = 0; rep(i,t.size()){ for(int j=i+1;j l(t[i],t[j]); int tt = 0; rep(k,n){ if(l.get_signed_dis(S[k].V1) * l.get_signed_dis(S[k].V2)