1.poj2318 TOYS

传送：http://poj.org/problem?id=2318

题意：有m个点落在n+1个区域内。问落在每个区域的个数。

分析：二分查找落在哪个区域内。叉积判断点与线段的位置。

  1.  #include<iostream>
  2.  #include<cstring>
  3.  #include<algorithm>
  4.  using namespace std;
  5.  const int maxn=;
  6.  struct point{
  7.      int x,y;
  8.      point(){ }
  9.      point(int _x,int _y):x(_x),y(_y){}
  10.      point operator +(const point &b) const{
  11.          return point(x+b.x,y+b.y);
  12.      }
  13.      point operator -(const point &b) const{
  14.          return point(x-b.x,y-b.y);
  15.      }
  16.      int operator *(const point &b) const{
  17.          return x*b.x,y*b.y;
  18.      }
  19.      int operator ^(const point &b) const{
  20.          return x*b.y-y*b.x;
  21.      }
  22.  };
  23.  struct line{
  24.      point s,e;
  25.      line(){}
  26.      line(point _s,point _e):s(_s),e(_e){}
  27.  }box[maxn];
  28.  int xmult(point p,point a,point b){
  29.      return (b-a)^(p-a);
  30.  }
  31.  int ans[maxn];
  32.  int main(){
  33.      ios::sync_with_stdio(false);
  34.      cin.tie();cout.tie();
  35.      int n,m,x1,y1,x2,y2,ui,li,x,y;
  36.      int _=;
  37.      while (cin >> n >> m >> x1 >> y1 >> x2 >> y2 && n){
  38.          if (_!=) cout << endl;
  39.          for (int i=;i<n;i++){
  40.              cin >> ui >> li;
  41.              box[i]=line(point(ui,y1),point(li,y2));
  42.          }
  43.          box[n]=line(point(x2,y1),point(x2,y2));
  44.          memset(ans,,sizeof(ans));
  45.          point p;
  46.          for (int i=;i<m;i++){
  47.              cin >> x >> y;
  48.              p=point(x,y);
  49.              int l=,r=n,mid,tmp;
  50.              while (l<=r){
  51.                  mid=(l+r)>>;
  52.                  if (xmult(p,box[mid].s,box[mid].e)<){
  53.                      tmp=mid;
  54.                      r=mid-;
  55.                  }
  56.                  else l=mid+;
  57.              }
  58.              ans[tmp]++;
  59.          }
  60.          for (int i=;i<=n;i++) cout << i << ": " << ans[i] << endl;
  61.          _++;
  62.      }
  63.      return ;
  64.  }

poj2318

2.poj2398 Toy Storage

传送：http://poj.org/problem?id=2398

题意：与poj2318相同，输出要求输出有t（t>0）个玩具的格子数。

分析：同2318

  1.  #include<iostream>
  2.  #include<cstring>
  3.  #include<algorithm>
  4.  using namespace std;
  5.  const int maxn=;
  6.  struct point{
  7.      int x,y;
  8.      point(){ }
  9.      point(int _x,int _y):x(_x),y(_y){}
  10.      point operator +(const point &b) const{
  11.          return point(x+b.x,y+b.y);
  12.      }
  13.      point operator -(const point &b) const{
  14.          return point(x-b.x,y-b.y);
  15.      }
  16.      int operator *(const point &b) const{
  17.          return x*b.x,y*b.y;
  18.      }
  19.      int operator ^(const point &b) const{
  20.          return x*b.y-y*b.x;
  21.      }
  22.  };
  23.  struct line{
  24.      point s,e;
  25.      line(){}
  26.      line(point _s,point _e):s(_s),e(_e){}
  27.      //
  28.      bool operator< (const line &b) const{
  29.          return s.x<b.s.x;
  30.      }
  31.  
  32.  }box[maxn];
  33.  int cmp(line p,line q){
  34.      return p.s.x<q.s.x;
  35.  }
  36.  int xmult(point p,point a,point b){
  37.      return (a-p)^(b-p);
  38.  }
  39.  int ans[maxn],vis[maxn];
  40.  int main(){
  41.      ios::sync_with_stdio(false);
  42.      cin.tie();cout.tie();
  43.      int n,m,x1,y1,x2,y2,ui,li,x,y;
  44.      while (cin >> n >> m >> x1 >> y1 >> x2 >> y2 && n){
  45.          for (int i=;i<n;i++){
  46.              cin >> ui >> li;
  47.              box[i]=line(point(ui,y1),point(li,y2));
  48.          }
  49.          box[n]=line(point(x2,y1),point(x2,y2));
  50.          sort(box,box++n);
  51.          memset(ans,,sizeof(ans));
  52.          point p;
  53.          for (int i=;i<m;i++){
  54.              cin >> x >> y;
  55.              p=point(x,y);
  56.              int l=,r=n,mid,tmp;
  57.              while (l<=r){
  58.                  mid=(l+r)>>;
  59.                  if (xmult(p,box[mid].s,box[mid].e)<){
  60.                      tmp=mid;
  61.                      r=mid-;
  62.                  }
  63.                  else l=mid+;
  64.              }
  65.              ans[tmp]++;
  66.          }
  67.          memset(vis,,sizeof(vis));
  68.          for (int i=;i<=n;i++) if (ans[i]) vis[ans[i]]++;
  69.          cout << "Box\n";
  70.          for (int i=;i<=n;i++){
  71.              if (!vis[i]) continue;
  72.              cout << i << ": " << vis[i] << endl;
  73.          }
  74.      }
  75.      return ;
  76.  }

poj2398

3.poj3304 Segments

传送：http://poj.org/problem?id=3304

题意：给出n条线段两个端点的坐标，问所有线段投影到一条直线上，如果这些所有投影至少相交于一点就输出Yes！，否则输出No！。

分析：如果有存在这样的直线，过投影相交区域作直线的垂线，该垂线必定与每条线段相交，问题转化为问是否存在一条线和所有线段相交。

直线肯定经过两个端点。两两枚举端点（包括同一条线段），判断直线和线段是否相交。

  1.  #include<iostream>
  2.  #include<cstring>
  3.  #include<cmath>
  4.  #include<algorithm>
  5.  using namespace std;
  6.  const double eps=1e-;
  7.  const int maxn=;
  8.  struct point{
  9.      double x,y;
  10.      point(){}
  11.      point(double _x,double _y):x(_x),y(_y){}
  12.      point operator -(const point &b) const{
  13.          return point(x-b.x,y-b.y);
  14.      }
  15.      double operator ^(const point &b) const{
  16.          return x*b.y-y*b.x;
  17.      }
  18.      double operator *(const point &b) const{
  19.          return x*b.x+y*b.y;
  20.      }
  21.  };
  22.  struct line{
  23.      point s,t;
  24.      line(){}
  25.      line(point _s,point _t):s(_s),t(_t){}
  26.  }a[maxn];
  27.  int sgn(double x){
  28.      if (fabs(x)<eps) return ;
  29.      if (x<) return -;
  30.      return ;
  31.  }
  32.  double xmult(point p,point a,point b){  // PA X PB
  33.      return (a-p)^(b-p);
  34.  }
  35.  int seg_inter_line(line l1,line l2){
  36.      return sgn(xmult(l2.s,l1.s,l1.t))*sgn(xmult(l2.t,l1.s,l1.t))<=;
  37.  }
  38.  double dist(point a,point b){
  39.      return sqrt((b-a)*(b-a));
  40.  }
  41.  int check(line l1,int n){
  42.      if (sgn(dist(l1.s,l1.t))==) return ;
  43.      for (int i=;i<n;i++)
  44.          if (seg_inter_line(l1,a[i])==)    return ;
  45.      return ;
  46.  }
  47.  int main(){
  48.      int t,n; cin >> t;
  49.      double x1,x2,y1,y2;
  50.      while (t--){
  51.          cin >> n;
  52.          for (int i=;i<n;i++){
  53.              cin >> x1 >> y1 >> x2 >> y2;
  54.              a[i]=line(point(x1,y1),point(x2,y2));
  55.          }
  56.          int f=;
  57.          for (int i=;i<n;i++){
  58.              for (int j=;j<n;j++){
  59.                  if (check(line(a[i].s,a[j].s),n) || check(line(a[i].s,a[j].t),n)
  60.                      || check(line(a[i].t,a[j].s),n) || check(line(a[i].t,a[j].t),n)){
  61.                          f=; break;
  62.                      }
  63.              }
  64.          }
  65.          if (f) cout << "Yes!\n";
  66.          else cout << "No!\n";
  67.      }
  68.  }

poj3304

4.poj1269 Intersecting Lines

传送：http://poj.org/problem?id=1269

题意：给出两条直线，问两条直线的位置关系（平行，重合，相交）。相交的话输出交点。

  1.  #include<iostream>
  2.  #include<cstdio>
  3.  #include<algorithm>
  4.  #include<cmath>
  5.  #include<cstring>
  6.  using namespace std;
  7.  const double eps=1e-;
  8.  int sgn(double x){
  9.      if (fabs(x)<eps) return ;
  10.      if (x<) return -;
  11.      return ;
  12.  }
  13.  struct point{
  14.      double x,y;
  15.      point(){}
  16.      point(double _x,double _y):x(_x),y(_y){}
  17.      point operator -(const point &b) const{
  18.          return point(x-b.x,y-b.y);
  19.      }
  20.      double operator ^(const point &b) const{
  21.          return x*b.y-y*b.x;
  22.      }
  23.      double operator *(const point &b) const{
  24.          return x*b.x+y*b.y;
  25.      }
  26.  };
  27.  struct line{
  28.      point s,t;
  29.      line(){}
  30.      line(point _s,point _t):s(_s),t(_t){}
  31.      pair<point,int> operator &(const line &b) const{
  32.          point res=s;
  33.          if (sgn((s-t)^(b.s-b.t))==){
  34.              if (sgn((b.s-s)^(b.t-s))==) return make_pair(res,);  //两直线重合
  35.              else return make_pair(res,); //两直线平行
  36.          }
  37.          double k=((s-b.s)^(b.s-b.t))/((s-t)^(b.s-b.t));
  38.          res.x+=(t.x-s.x)*k;
  39.          res.y+=(t.y-s.y)*k;
  40.          return make_pair(res,); //两直线相交
  41.      }
  42.  };
  43.  int main(){
  44.      int n; cin >> n;
  45.      double x1,x2,y1,y2;
  46.      line l1,l2;
  47.      cout << "INTERSECTING LINES OUTPUT\n";
  48.      for (int i=;i<n;i++){
  49.          cin >> x1 >> y1 >> x2 >> y2;
  50.          l1=line(point(x1,y1),point(x2,y2));
  51.          cin >> x1 >> y1 >> x2 >> y2;
  52.          l2=line(point(x1,y1),point(x2,y2));
  53.          pair<point,int> ans=l1&l2;
  54.          if (ans.second==) printf("POINT %.2f %.2f\n",ans.first.x,ans.first.y);
  55.          else if (ans.second==) printf("NONE\n");
  56.          else printf("LINE\n");
  57.      }
  58.      cout << "END OF OUTPUT\n";
  59.      return ;
  60.  }

poj1269

5.poj1556 The Doors

传送：http://poj.org/problem?id=1556

题意：10×10的方阵内有若干个墙壁，墙壁之间有空隙（可通过），问从（0,5）到（10,5）的最短路。

分析：建图，线段两两判断是否相交，不相交就建边。然后跑最短路。（数据较小，随便最短路算法。

  1.  #include<iostream>
  2.  #include<cstdio>
  3.  #include<cmath>
  4.  #include<algorithm>
  5.  #include<cstring>
  6.  #include<algorithm>
  7.  using namespace std;
  8.  const double eps=1e-;
  9.  const double inf=1e20;
  10.  const int maxn=;
  11.  double mp[maxn][maxn];
  12.  int sgn(double x){
  13.      if (fabs(x)<eps) return ;
  14.      if (x<) return -;
  15.      return ;
  16.  }
  17.  struct point{
  18.      double x,y;
  19.      point(){}
  20.      point(double _x,double _y):x(_x),y(_y){}
  21.      point operator -(const point &b) const{
  22.          return point(x-b.x,y-b.y);
  23.      }
  24.      double operator ^(const point &b) const{
  25.          return x*b.y-y*b.x;
  26.      }
  27.      double operator *(const point &b) const{
  28.          return x*b.x+y*b.y;
  29.      }
  30.  };
  31.  struct line{
  32.      point s,e;
  33.      line(){}
  34.      line(point _s,point _e):s(_s),e(_e){}
  35.  }a[maxn];
  36.  double dist(point a,point b){
  37.      return sqrt((b-a)*(b-a));
  38.  }
  39.  int inter(line l1,line l2){
  40.      return
  41.          max(l1.s.x,l1.e.x) >= min(l2.s.x,l2.e.x)
  42.          && max(l2.s.x,l2.e.x) >= min(l1.s.x,l1.e.x)
  43.          && max(l1.s.y,l1.e.y) >= min(l2.s.y,l2.e.y)
  44.          && max(l2.s.y,l2.e.y) >= min(l1.s.y,l1.e.y)
  45.          && sgn((l2.s-l1.s)^(l1.e-l1.s))*sgn((l2.e-l1.s)^(l1.e-l1.s))<=
  46.          && sgn((l1.s-l2.s)^(l2.e-l2.s))*sgn((l1.e-l2.s)^(l2.e-l2.s))<=;
  47.  }
  48.  double floyd(int n){
  49.      for (int i=;i<=n;i++)
  50.          for (int j=;j<=n;j++)
  51.              for (int k=;k<=n;k++)
  52.                  if (mp[i][k]+mp[k][j]<mp[i][j]) mp[i][j]=mp[i][k]+mp[k][j];
  53.      return mp[][n];
  54.  }
  55.  int main(){
  56.      int n; double x,y1,y2,y3,y4;
  57.      while (cin >> n && n!=-){
  58.          for (int i=;i<=n;i++){
  59.              cin >> x >> y1 >> y2 >> y3 >> y4;
  60.              a[*i-]=line(point(x,y1),point(x,y2));
  61.              a[*i]=line(point(x,y3),point(x,y4));
  62.          }
  63.          for (int i=;i<=*n+;i++)
  64.              for (int j=;j<=*n+;j++)
  65.                  if (i==j) mp[i][j]=;else mp[i][j]=inf;
  66.          for (int i=;i<=*n;i++){
  67.              int l=(i+)/;
  68.              int f=;
  69.              point tmp;
  70.              if (i&) tmp=a[(i+)/].s; else tmp=a[(i+)/].e;
  71.              for (int j=;j<l;j++)
  72.                  if (inter(a[*j-],line(point(,),tmp))==
  73.                      && inter(a[*j],line(point(,),tmp))==) f=;
  74.              if (f) mp[][i]=mp[i][]=dist(point(,),tmp);
  75.              f=;
  76.              for (int j=l+;j<=n;j++)
  77.                  if (inter(a[*j-],line(point(,),tmp))==
  78.                      && inter(a[*j],line(point(,),tmp))==) f=;
  79.              if (f) mp[i][*n+]=mp[*n+][i]=dist(point(,),tmp);
  80.          }
  81.          for (int i=;i<=*n;i++){
  82.              for (int j=i+;j<=*n;j++){
  83.                  int l1=(i+)/,l2=(j+)/;
  84.                  int f=;
  85.                  point p1,p2;
  86.                  if (i&) p1=a[(i+)/].s; else p1=a[(i+)/].e;
  87.                  if (j&) p2=a[(j+)/].s; else p2=a[(j+)/].e;
  88.                  for (int k=l1+;k<l2;k++)
  89.                      if (inter(a[*k-],line(p1,p2))==
  90.                          && inter(a[*k],line(p1,p2))==) f=;
  91.                  if (f) mp[i][j]=mp[j][i]=dist(p1,p2);
  92.              }
  93.          }
  94.          int f=;
  95.          for (int i=;i<=n;i++)
  96.              if (inter(a[*i-],line(point(,),point(,)))==
  97.                  && inter(a[*i],line(point(,),point(,)))==) f=;
  98.          if (f) mp[][*n+]=mp[*n+][]=;
  99.          double ans=floyd(*n+);
  100.          printf("%.2f\n",ans);
  101.      }
  102.      return ;
  103.  }

poj1556

6.poj2653 Pick-up sticks

传送：http://poj.org/problem?id=2653

题意：问有多少线段不与比它编号大的相交。

分析：线段交直接暴力判断。

  1.  #include<iostream>
  2.  #include<algorithm>
  3.  #include<cstring>
  4.  #include<cmath>
  5.  using namespace std;
  6.  const double eps=1e-;
  7.  const int maxn=;
  8.  int sgn(double x){
  9.      if (fabs(x)<eps) return ;
  10.      if (x<) return -;
  11.      return ;
  12.  }
  13.  struct point{
  14.      double x,y;
  15.      point(){}
  16.      point(double _x,double _y):x(_x),y(_y){}
  17.      point operator -(const point &b) const{
  18.          return point(x-b.x,y-b.y);
  19.      }
  20.      double operator ^(const point &b) const{
  21.          return x*b.y-y*b.x;
  22.      }
  23.      double operator *(const point &b) const{
  24.          return x*b.x+y*b.y;
  25.      }
  26.  };
  27.  struct line{
  28.      point s,e;
  29.      line(){}
  30.      line(point _s,point _e):s(_s),e(_e){}
  31.  }a[maxn];
  32.  double dist(point a,point b){
  33.      return sqrt((b-a)*(b-a));
  34.  }
  35.  int inter(line l1,line l2){
  36.      return
  37.          max(l1.s.x,l1.e.x) >= min(l2.s.x,l2.e.x)
  38.          && max(l2.s.x,l2.e.x) >= min(l1.s.x,l1.e.x)
  39.          && max(l1.s.y,l1.e.y) >= min(l2.s.y,l2.e.y)
  40.          && max(l2.s.y,l2.e.y) >= min(l1.s.y,l1.e.y)
  41.          && sgn((l2.s-l1.s)^(l1.e-l1.s))*sgn((l2.e-l1.s)^(l1.e-l1.s))<=
  42.          && sgn((l1.s-l2.s)^(l2.e-l2.s))*sgn((l1.e-l2.s)^(l2.e-l2.s))<=;
  43.  }
  44.  bool vis[maxn];
  45.  int main(){
  46.      ios::sync_with_stdio(false);
  47.      cin.tie();cout.tie();
  48.      int n; double x1,x2,y1,y2;
  49.      while (cin >> n && n){
  50.          for (int i=;i<n;i++){
  51.              cin >> x1 >> y1 >> x2 >> y2;
  52.              a[i]=line(point(x1,y1),point(x2,y2));
  53.              vis[i]=;
  54.          }
  55.          for (int i=;i<n;i++){
  56.              for (int j=i+;j<n;j++){
  57.                  if (inter(a[i],a[j])){vis[i]=;break;}
  58.              }
  59.          }
  60.          cout << "Top sticks: ";
  61.          int f=;
  62.          for (int i=;i<n;i++){
  63.              if (!vis[i]) continue;
  64.              if (f){cout << i+;f=;}
  65.              else cout << ", " << i+;
  66.          }
  67.          cout << ".\n";
  68.      }
  69.  }

poj2653

7.poj1066 Treasure Hunt

传送：http://poj.org/problem?id=1066

题意：在金字塔内有一个宝藏p(x,y);现在要取出这个宝藏在金字塔内有许多墙为了进入宝藏所在的房间必须把墙炸开但是炸墙只能炸每个房间墙的中点求将宝藏运出城堡所需要的最小炸墙数。

分析：枚举每个线段的端点与P点形成的直线，与多少线段相交。其次，可以从墙壁顶点进入，相交数+1与ans比较。
  1.  #include<iostream>
  2.  #include<cstring>
  3.  #include<algorithm>
  4.  #include<cmath>
  5.  #include<cstring>
  6.  using namespace std;
  7.  const double dx[]={,,,},dy[]={,,,};
  8.  const double eps=1e-;
  9.  const int maxn=;
  10.  int sgn(double x){
  11.      if (fabs(x)<eps) return ;
  12.      if (x<) return -;
  13.      return ;
  14.  }
  15.  struct point{
  16.      double x,y;
  17.      point(){}
  18.      point(double _x,double _y):x(_x),y(_y){}
  19.      point operator -(const point &b) const{
  20.          return point(x-b.x,y-b.y);
  21.      }
  22.      double operator ^(const point &b) const{
  23.          return x*b.y-y*b.x;
  24.      }
  25.      double operator *(const point &b) const{
  26.          return x*b.x+y*b.y;
  27.      }
  28.  }p[*maxn];
  29.  struct line{
  30.      point s,e;
  31.      line(){}
  32.      line(point _s,point _e):s(_s),e(_e){}
  33.  }a[maxn];
  34.  int inter(line l1,line l2){
  35.      return
  36.          max(l1.s.x,l1.e.x) >= min(l2.s.x,l2.e.x)
  37.          && max(l2.s.x,l2.e.x) >= min(l1.s.x,l1.e.x)
  38.          && max(l1.s.y,l1.e.y) >= min(l2.s.y,l2.e.y)
  39.          && max(l2.s.y,l2.e.y) >= min(l1.s.y,l1.e.y)
  40.          && sgn((l2.s-l1.s)^(l1.e-l1.s))*sgn((l2.e-l1.s)^(l1.e-l1.s))<=
  41.          && sgn((l1.s-l2.s)^(l2.e-l2.s))*sgn((l1.e-l2.s)^(l2.e-l2.s))<=;
  42.  }
  43.  int main(){
  44.      int n; double x1,y1,x2,y2;
  45.      while (cin >> n){
  46.          for (int i=;i<=n;i++){
  47.              cin >> x1 >> y1 >> x2 >> y2;
  48.              a[i]=line(point(x1,y1),point(x2,y2));
  49.              p[*i-]=point(x1,y1); p[*i]=point(x2,y2);
  50.          }
  51.          point s;
  52.          cin >> x1 >> y1;
  53.          s=point(x1,y1);
  54.          int ans=;
  55.          for (int i=;i<=*n;i++){
  56.              line l1=line(s,p[i]);
  57.              int res=;
  58.              for (int j=;j<=n;j++)
  59.                  if (inter(l1,a[j])) res++;
  60.              ans=min(res,ans);
  61.          }
  62.          for (int i=;i<;i++){
  63.              line l1=line(s,point(dx[i],dy[i]));
  64.              int res=;
  65.              for (int j=;j<=n;j++)
  66.                  if (inter(l1,a[j])) res++;
  67.              ans=min(ans,res+);
  68.          }
  69.          cout << "Number of doors = " << ans << endl;
  70.      }
  71.      return ;
  72.  }

poj1066

8.poj1410 Intersection

传送：http://poj.org/problem?id=1410

题意：给定一个矩阵和线段，若线段在矩阵内或与矩阵某条边相交，输出T，否则输出F。

分析：判断线段与四条边线段交。再判断是否在矩阵内。

  1.  #include<iostream>
  2.  #include<algorithm>
  3.  #include<cstring>
  4.  #include<cmath>
  5.  using namespace std;
  6.  const double eps=1e-;
  7.  const int maxn=;
  8.  int sgn(double x){
  9.      if (fabs(x)<eps) return ;
  10.      if (x<) return -;
  11.      return ;
  12.  }
  13.  struct point{
  14.      double x,y;
  15.      point(){}
  16.      point(double _x,double _y):x(_x),y(_y){}
  17.      point operator -(const point &b) const{
  18.          return point(x-b.x,y-b.y);
  19.      }
  20.      double operator ^(const point &b) const{
  21.          return x*b.y-y*b.x;
  22.      }
  23.      double operator *(const point &b) const{
  24.          return x*b.x+y*b.y;
  25.      }
  26.  }p[maxn];
  27.  struct line{
  28.      point s,e;
  29.      line(){}
  30.      line(point _s,point _e):s(_s),e(_e){}
  31.  };
  32.  int inter(line l1,line l2){
  33.      return
  34.          max(l1.s.x,l1.e.x) >= min(l2.s.x,l2.e.x)
  35.          && max(l2.s.x,l2.e.x) >= min(l1.s.x,l1.e.x)
  36.          && max(l1.s.y,l1.e.y) >= min(l2.s.y,l2.e.y)
  37.          && max(l2.s.y,l2.e.y) >= min(l1.s.y,l1.e.y)
  38.          && sgn((l2.s-l1.s)^(l1.e-l1.s))*sgn((l2.e-l1.s)^(l1.e-l1.s))<=
  39.          && sgn((l1.s-l2.s)^(l2.e-l2.s))*sgn((l1.e-l2.s)^(l2.e-l2.s))<=;
  40.  }
  41.  int onseg(point p,line l){
  42.      return
  43.          sgn((l.s-p)^(l.e-p))== &&
  44.          sgn((p.x-l.s.x)*(p.x-l.e.x))<= &&
  45.          sgn((p.y-l.s.y)*(p.y-l.e.y))<=;
  46.  }
  47.  //判断点在凸多边形内
  48.  //点形成一个凸包，而且按逆时针排序（如果是顺时针把里面的<0改为>0）
  49.  //点的编号:0~n-1
  50.  //返回值：
  51.  //-1:点在凸多边形外
  52.  //0:点在凸多边形边界上
  53.  //1:点在凸多边形内
  54.  int inConvexPoly(point a,point p[],int n){
  55.      for (int i=;i<n;i++){
  56.          if (sgn((p[i]-a)^(p[(i+)%n]-a))<) return -;
  57.          else if (onseg(a,line(p[i],p[(i+)%n]))) return ;
  58.      }
  59.      return ;
  60.  }
  61.  int main(){
  62.      int t; cin >> t;
  63.      double x1,y1,x2,y2;
  64.      while (t--){
  65.          cin >> x1 >> y1 >> x2 >> y2;
  66.          line l=line(point (x1,y1),point(x2,y2));
  67.          cin >> x1 >> y1 >> x2 >> y2;
  68.          if (x1>x2) swap(x1,x2); if (y1>y2) swap(y1,y2);
  69.          p[]=point(x1,y1);p[]=point(x2,y1);p[]=point(x2,y2);p[]=point(x1,y2);
  70.          if (inter(l,line(p[],p[])) || inter(l,line(p[],p[]))
  71.              || inter(l,line(p[],p[])) || inter(l,line(p[],p[]))){
  72.                  cout << "T\n";continue;
  73.              }
  74.          if (inConvexPoly(l.s,p,)>= || inConvexPoly(l.e,p,)>=){
  75.              cout << "T\n"; continue;
  76.          }
  77.          cout << "F\n";
  78.      }
  79.      return ;
  80.  }

poj1410

9.poj1696 Space Ant

传送：http://poj.org/problem?id=1696

题意：一只蚂蚁只会向左转，现在给出平面上很多个点，求解一种走法能使得蚂蚁能经过的点最多，每个顶点该蚂蚁只能经过一次，且所行走的路线不能发生交叉。

分析：对于题目所输入的点，先找出最左下方的顶点(即纵坐标最小的顶点)，然后对剩下的顶点按照对与左下点的极角排序，然后反复找最左下的点，反复进行极角排序，同时记录排序后左下的顶点。

极角排序方法：利用叉积，看向量p1和p2的叉积是否小于零，是则说明p1在p2逆时针方向，即p1的极角比p2的大，极角相同则按离p0距离降序排序。

  1.  #include<iostream>
  2.  #include<cstring>
  3.  #include<algorithm>
  4.  #include<cmath>
  5.  using namespace std;
  6.  const double eps=1e-;
  7.  const int maxn=;
  8.  int sgn(double x){
  9.      if (fabs(x)<eps) return ;
  10.      if (x<) return -;
  11.      return ;
  12.  }
  13.  struct point{
  14.      double x,y;
  15.      int id;
  16.      point(){}
  17.      point(double _x,double _y):x(_x),y(_y){}
  18.      point operator -(const point &b) const{
  19.          return point(x-b.x,y-b.y);
  20.      }
  21.      double operator ^(const point &b) const{
  22.          return x*b.y-y*b.x;
  23.      }
  24.      double operator *(const point &b) const{
  25.          return x*b.x+y*b.y;
  26.      }
  27.  }p[maxn];
  28.  double dist(point a,point b){
  29.      return sqrt((a-b)*(a-b));
  30.  }
  31.  int pos;
  32.  int cmp(point a,point b){
  33.      double tmp=(a-p[pos])^(b-p[pos]);
  34.      if (sgn(tmp)==) return dist(a,p[pos])<dist(b,p[pos]);
  35.      else if (sgn(tmp)<) return false;
  36.      else return true;
  37.  }
  38.  int main(){
  39.      int t,n; cin >> t;
  40.      double x,y;
  41.      while (t--){
  42.          cin >> n;
  43.          for (int i=;i<n;i++){
  44.              cin >> p[i].id >> p[i].x >> p[i].y;
  45.              if (p[i].y<p[].y || p[i].y==p[].y && p[i].x<p[].x) swap(p[i],p[]);
  46.          }
  47.          pos=;
  48.          for (int i=;i<n;i++){
  49.              sort(p+i,p+n,cmp);
  50.              pos++;
  51.          }
  52.          cout << n;
  53.          for (int i=;i<n;i++) cout << " " << p[i].id;
  54.          cout << endl;
  55.      }
  56.  }

poj1696

10.poj3347 Kadj Squares

传送：http://poj.org/problem?id=3347

题意：给出n个边长的正方形，要求尽可能贴着放。问最后从上往下看能看到哪几个正方形。

分析：每新增一个正方形，就让他与左侧的每一个正方形贴紧，求的其左端坐标，最终结果一定是最大的那个。然后求的相应的最右端坐标。这样就转化为了线段，最后用朴素的n2算法求出每条线段没有被覆盖的长度，如果长度大于0，即可输出。

  1.  #include<iostream>
  2.  #include<algorithm>
  3.  #include<cmath>
  4.  using namespace std;
  5.  const int maxn=;
  6.  struct node{
  7.      int len,l,r;
  8.  }a[maxn];
  9.  int main(){
  10.      int n;
  11.      while (cin >> n && n){
  12.          for (int i=;i<n;i++){
  13.              cin >> a[i].len;  //将图形扩大根2倍，len为对角线长度的一半
  14.              a[i].l=;
  15.              for (int j=;j<i;j++){
  16.                  a[i].l=max(a[i].l,a[j].r-abs(a[i].len-a[j].len));
  17.              }
  18.              a[i].r=a[i].l+*a[i].len;
  19.          }
  20.          for (int i=;i<n;i++){
  21.              for (int j=;j<i;j++){
  22.                  if (a[i].l<a[i].r){
  23.                      if (a[i].len<a[j].len && a[i].l<a[j].r) a[i].l=a[j].r;
  24.                      else if (a[i].len>a[j].len && a[i].l<a[j].r) a[j].r=a[i].l;
  25.                  }
  26.              }
  27.          }
  28.          for (int i=;i<n;i++){
  29.              if (a[i].l<a[i].r) cout << i+ << " ";
  30.          }
  31.          cout << endl;
  32.      }
  33.      return ;
  34.  }

poj3347

11.poj1654 Area

传送：http://poj.org/problem?id=1654

题意：从原点开始，超八个方向走。问最后形成的多边形面积。

分析：rt。直接求多边形面积。

  1.  #include<iostream>
  2.  #include<algorithm>
  3.  #include<cstring>
  4.  #include<cmath>
  5.  using namespace std;
  6.  const int dx[]={,-,,,-,,,-,,},dy[]={,-,-,-,,,,,,};
  7.  const int maxn=;
  8.  typedef long long ll;
  9.  struct point{
  10.      int x,y;
  11.      point(){}
  12.      point(int _x,int _y):x(_x),y(_y){}
  13.      ll operator ^(const point &b)const{
  14.          return x*b.y-y*b.x;
  15.      }
  16.  };
  17.  struct polygon{
  18.      int n;
  19.      point p[maxn];
  20.      void add(point q){p[n++]=q;}
  21.      ll getarea(){
  22.          ll sum=;
  23.          for (int i=;i<n;i++) sum+=(p[i]^(p[(i+)%n]));
  24.          return sum<?-sum:sum;
  25.      }
  26.  };
  27.  polygon C;
  28.  int main(){
  29.      char ch;
  30.      int t,k; cin >> t;
  31.      while (t--){
  32.          C.n=;
  33.          ll x=,y=;
  34.          while (cin >> ch && ch!=''){
  35.              k=ch-'';
  36.              x+=dx[k]; y+=dy[k];
  37.              C.add(point(x,y));
  38.          }
  39.          ll ans=C.getarea();
  40.          if (ans&) cout << ans/ << ".5\n";
  41.          else cout << ans/ << endl;
  42.      }
  43.      return ;
  44.  }

poj1654

12.poj2954 Triangle

　　传送：http://poj.org/problem?id=2954

　　题意：求三角形内部有多少格点。

　　分析：pick定理。S=a+b/2-1。a代表多边形内部的格点，b代表多边形边上的格点。S代表多边形面积。

　　

  1.  #include<iostream>
  2.  #include<cstring>
  3.  #include<algorithm>
  4.  #include<cmath>
  5.  using namespace std;
  6.  typedef long long ll;
  7.  const int maxn=;
  8.  struct point{
  9.      int x,y;
  10.      point(){}
  11.      point(int _x,int _y):x(_x),y(_y){}
  12.      ll operator ^(const point &b)const{
  13.          return x*b.y-y*b.x;
  14.      }
  15.  };
  16.  struct polygon{
  17.      int n;
  18.      point p[maxn];
  19.      void add(point q){p[n++]=q;}
  20.      ll getarea(){
  21.          ll sum=;
  22.          for (int i=;i<n;i++) sum+=(p[i]^p[(i+)%n]);
  23.          return sum<?-sum:sum;
  24.      }
  25.  };
  26.  polygon C;
  27.  ll gcd(ll a,ll b){
  28.      if (b==) return a;
  29.      else return gcd(b,a%b);
  30.  }
  31.  int main(){
  32.      int x1,x2,x3,y1,y2,y3;
  33.      while (cin >> x1 >> y1 >> x2 >> y2 >> x3 >> y3){
  34.          if (!x1 && !x2 && !x3 && !y1 && !y2 && !y3) break;
  35.          C.n=;
  36.          C.add(point(x1,y1));C.add(point(x2,y2));C.add(point(x3,y3));
  37.          ll area=C.getarea();
  38.          int b=;
  39.          for (int i=;i<C.n;i++){
  40.              b+=gcd(abs(C.p[i].x-C.p[(i+)%C.n].x),abs(C.p[i].y-C.p[(i+)%C.n].y));
  41.          }
  42.          cout << (area+-b)/ << endl;    //S=a+b/2-1
  43.      }
  44.      return ;
  45.  }

poj2954

【kuangbin专题】计算几何基础的更多相关文章

  1. [kuangbin]专题六 最小生成树 题解+总结

    kuangbin专题链接:https://vjudge.net/article/752 kuangbin专题十二 基础DP1 题解+总结:https://www.cnblogs.com/RioTian ...

  2. nyis oj 68 三点顺序 （计算几何基础）

    三点顺序 时间限制:1000 ms  |  内存限制:65535 KB 难度:3 描写叙述 如今给你不共线的三个点A,B,C的坐标,它们一定能组成一个三角形,如今让你推断A,B,C是顺时针给出的还是逆 ...

  3. kuangbin专题 数论基础 part1？

    线段树专题太难了,那我来做数学吧! 但数学太难了,我......(扯 这两天想了做了查了整理了几道数学. 除了一些进阶的知识,像莫比乌斯反演,杜教筛,min25学不会我跳了,一些基础的思维还是可以记录 ...

  4. 【BZOJ】1043: [HAOI2008]下落的圆盘（计算几何基础+贪心）

    http://www.lydsy.com/JudgeOnline/problem.php?id=1043 唯一让我不会的就是怎么求圆的周长并QAAQ... 然后发现好神!我们可以将圆弧变成$[0, 2 ...

  5. 计算几何基础——矢量和叉积 && 叉积、线段相交判断、凸包（转载）

    转载自 http://blog.csdn.net/william001zs/article/details/6213485 矢量 如果一条线段的端点是有次序之分的话,那么这种线段就称为 有向线段,如果 ...

  6. BZOJ_1610_[Usaco2008_Feb]_Line连线游戏_(计算几何基础+暴力)

    描述 http://www.lydsy.com/JudgeOnline/problem.php?id=1610 给出n个点,问两两确定的直线中,斜率不同的共有多少条. 分析 暴力枚举直线,算出来斜率放 ...

  7. 二维计算几何基础题目泛做（SYX第一轮）

    题目1: POJ 2318 TOYS 题目大意: 给一个有n个挡板的盒子,从左到右空格编号为0...n.有好多玩具,问每个玩具在哪个空格里面. 算法讨论: 直接叉积判断就可以.注意在盒子的边界上面也算 ...

  8. 最小生成树 kuangbin专题最后一个题

    题目链接:https://cn.vjudge.net/contest/66965#problem/N 注释:这道题需要用krustra,用prim的话可能会超时.并且在计算距离的时候要尽量减少步骤,具 ...

  9. kuangbin专题十六 KMP&&扩展KMP HDU2087 剪花布条

    一块花布条,里面有些图案,另有一块直接可用的小饰条,里面也有一些图案.对于给定的花布条和小饰条,计算一下能从花布条中尽可能剪出几块小饰条来呢? Input输入中含有一些数据,分别是成对出现的花布条和小 ...

随机推荐

  1. 如何安全管理windows系统日志，windows系统日志的报表和告警

    如何安全管理windows系统日志,windows系统日志的报表和告警 无论大小,每个拥有IT基础设施的组织都容易发生内部安全攻击.您的损失等同于黑客的收益:访问机密数据.滥用检索到的信息.系统崩溃, ...

  2. vue2.0 移动端，下拉刷新，上拉加载更多插件 转：

    转自:http://www.cnblogs.com/gmajip/p/7204072.html <template lang="html"> <div class ...

  3. Python sys.argv[] 的用法

    sys.argv变量是一个list, 执行 python abc.py a b c 时, sys.argv[0]为 abc.py sys.argv[1]为 a sys.argv[2]为 b sys.a ...

  4. set集合的排序

    在hibernate的OneToMany的实体关联的时候,one端的set是无序的,可是需要按照顺序来搞的话就比较麻烦了. 下面给出一个例子. Set<DiaryPicture> diar ...

  5. tomcat7 安装 windows 服务

    tomcat 可以安装成windows 服务,这样 每次启动就不需要启动tomcat了. 具体配置: 1.修改 service.bat 在行首添加 set "JAVA_HOME=E:\jdk ...

  6. java常用设计模式八：代理模式

    一.概述 代理模式是指客户端并不直接调用实际的对象,而是通过调用代理,来间接的调用实际的对象. 其特征是代理类与委托类有同样的接口,真正的核心业务逻辑还是在实际对象里面. 二.为什么要使用代理模式 当 ...

  7. openstack之镜像管理

    概览 [root@cc07 fast-pulsar]# glance help | grep image [--os-image-url OS_IMAGE_URL] [--os-image-api-v ...

  8. XML xmlns

    xmlns xml namespaces 参考 http://www.w3school.com.cn/tags/tag_prop_xmlns.asp http://www.w3school.com.c ...

  9. 用Html创建简历

    用最基本的html实现简历的创建 代码 <!doctype html> <html> <head><meta charset="utf-8" ...

  10. HDU 3078 LCA转RMQ

    题意: n个点 m个询问 下面n个数字表示点权值 n-1行给定一棵树 m个询问 k u v k为0时把u点权值改为v 或者问 u-v的路径上 第k大的数 思路: LCA转RMQ求出 LCA(u,v) ...