【模版】计算几何

#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
const double oo=1e19, pi=acos(-1), eps=1e-10;

struct dbl { double x,y; };
struct line { dbl P,v; };
struct circle { dbl O; double r; };

dbl operator + (dbl a, dbl b) { return (dbl){a.x+b.x, a.y+b.y}; }
dbl operator - (dbl a, dbl b) { return (dbl){a.x-b.x, a.y-b.y}; }
dbl operator * (dbl a, double b) { return (dbl){a.x*b, a.y*b}; }
dbl operator / (dbl a, double b) { return (dbl){a.x/b, a.y/b}; }

int dcmp(double d) { if(d<-eps) return -1; else return d>eps; }
double cj(dbl a, dbl b) { return a.x*b.y-a.y*b.x; }
double dj(dbl a, dbl b) { return a.x*b.x+a.y*b.y; }
double dis(dbl a, dbl b) { return sqrt( (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y) ); }
double jj(dbl v) { double rn=atan2(v.y, v.x); rn=(rn<0)?rn+2*pi:rn; return rn; }
void rotate(dbl& v, double g) { double C=cos(g),S=sin(g); v=(dbl){C*v.x-S*v.y, S*v.x+C*v.y}; }

int main()
{
  freopen("oj.in", "r", stdin);
  freopen("oj.out", "w", stdout);
  return 0;
}