链接:

http://poj.org/problem?id=1265

题意:

给你一个多边形,求它的面积,内部格点数目,边上格点数目

题解:

pick公式:

给定顶点坐标均是整数点的简单多边形,有

          面积=内部格点数目+边上格点数目/2+1

边界上的格点数:

把每条边当做左开右闭的区间以避免重复,一条左开右闭的线段(x1,y1)->(x2,y2)上的格点数为:

          gcd(x2-x1,y2-y1)。

代码:

  1 #include <map>
  2 #include <set>
  3 #include <cmath>
  4 #include <queue>
  5 #include <stack>
  6 #include <cstdio>
  7 #include <string>
  8 #include <vector>
  9 #include <cstdlib>
 10 #include <cstring>
 11 #include <sstream>
 12 #include <iostream>
 13 #include <algorithm>
 14 #include <functional>
 15 using namespace std;
 16 #define rep(i,a,n) for (int i=a;i<n;i++)
 17 #define per(i,a,n) for (int i=n-1;i>=a;i--)
 18 #define all(x) (x).begin(),(x).end()
 19 #define pb push_back
 20 #define mp make_pair
 21 #define lson l,m,rt<<1
 22 #define rson m+1,r,rt<<1|1
 23 typedef long long ll;
 24 typedef vector<int> VI;
 25 typedef pair<int, int> PII;
 26 const ll MOD = 1e9 + 7;
 27 const int INF = 0x3f3f3f3f;
 28 const int MAXN = 2e4 + 7;
 29 // head
 30
 31 const double eps = 1e-8;
 32 int cmp(double x) {
 33     if (fabs(x) < eps) return 0;
 34     if (x > 0) return 1;
 35     return -1;
 36 }
 37
 38 const double pi = acos(-1);
 39 inline double sqr(double x) {
 40     return x*x;
 41 }
 42 struct point {
 43     double x, y;
 44     point() {}
 45     point(double a, double b) :x(a), y(b) {}
 46     void input() {
 47         scanf("%lf%lf", &x, &y);
 48     }
 49     friend point operator+(const point &a, const point &b) {
 50         return point(a.x + b.x, a.y + b.y);
 51     }
 52     friend point operator-(const point &a, const point &b) {
 53         return point(a.x - b.x, a.y - b.y);
 54     }
 55     friend point operator*(const double &a, const point &b) {
 56         return point(a*b.x, a*b.y);
 57     }
 58     friend point operator/(const point &a, const double &b) {
 59         return point(a.x / b, a.y / b);
 60     }
 61     double norm() {
 62         return sqrt(sqr(x) + sqr(y));
 63     }
 64 };
 65 double det(point a, point b) {
 66     return a.x*b.y - a.y*b.x;
 67 }
 68 double dot(point a, point b) {
 69     return a.x*b.x + a.y*b.y;
 70 }
 71 double dist(point a, point b) {
 72     return (a - b).norm();
 73 }
 74
 75 struct line {
 76     point a, b;
 77     line() {}
 78     line(point x, point y) :a(x), b(y) {}
 79 };
 80 double dis_point_segment(point p, point s, point t) {
 81     if (cmp(dot(p - s, t - s)) < 0) return (p - s).norm();
 82     if (cmp(dot(p - t, s - t)) < 0) return (p - t).norm();
 83     return fabs(det(s - p, t - p) / dist(s, t));
 84 }
 85 bool point_on_segment(point p, point s, point t) {
 86     return cmp(det(p - s, t - s)) == 0 && cmp(dot(p - s, p - t)) <= 0;
 87 }
 88 bool parallel(line a, line b) {
 89     return !cmp(det(a.a - a.b, b.a - b.b));
 90 }
 91 bool line_make_point(line a, line b,point &res) {
 92     if (parallel(a, b)) return false;
 93     double s1 = det(a.a - b.a, b.b - b.a);
 94     double s2 = det(a.b - b.a, b.b - b.a);
 95     res = (s1*a.b - s2*a.a) / (s1 - s2);
 96     return true;
 97 }
 98
 99 struct polygon {
100     int n;
101     point a[MAXN];
102     double area() {
103         double sum = 0;
104         a[n] = a[0];
105         rep(i, 0, n) sum += det(a[i], a[i + 1]);
106         return sum / 2;
107     }
108     point MassCenter() {
109         point ans = point(0, 0);
110         if (cmp(area()) == 0) return ans;
111         a[n] = a[0];
112         rep(i, 0, n) ans = ans + det(a[i + 1], a[i])*(a[i] + a[i + 1]);
113         return ans / area() / 6;
114     }
115     int gcd(int a, int b) {
116         return b == 0 ? a : gcd(b, a%b);
117     }
118     int Border_Int_Point_Num() {
119         int num = 0;
120         a[n] = a[0];
121         rep(i, 0, n) num += gcd(abs(int(a[i + 1].x - a[i].x)),
122             abs(int(a[i + 1].y - a[i].y)));
123         return num;
124     }
125     int Inside_Int_Point_Num() {
126         return int(area()) + 1 - Border_Int_Point_Num() / 2;
127     }
128 };
129
130 int n;
131 polygon p;
132
133 int main() {
134     int T;
135     cin >> T;
136     rep(cas, 1, T + 1) {
137         cin >> n;
138         p.n = n;
139         rep(i, 0, n) scanf("%lf%lf", &p.a[i].x, &p.a[i].y);
140         rep(i, 1, n) p.a[i] = p.a[i] + p.a[i - 1];
141         int a = p.Inside_Int_Point_Num();
142         int b = p.Border_Int_Point_Num();
143         double c = p.area();
144         printf("Scenario #%d:\n%d %d %.1f\n\n", cas, a, b, c);
145     }
146     return 0;
147 }

转载于:https://www.cnblogs.com/baocong/p/6731151.html

POJ 1265 计算几何 多边形面积 内部格点数 边上格点数相关推荐

  1. poj 1654 Area 多边形面积

    /* poj 1654 Area 多边形面积题目意思很简单,但是1000000的point开不了 */ #include<stdio.h> #include<math.h> # ...

  2. POJ 3855 计算几何·多边形重心

    思路: 多边形面积->任选一个点,把多边形拆成三角,叉积一下 三角形重心->(x1+x2+x3)/3,(y1+y2+y3)/3 多边形重心公式题目中有,套一下就好了 计算多边形重心方法: ...

  3. POJ 1265-Area(计算几何+皮克定理+多边形面积公式)

    题目地址:POJ 1265 题意:给定一个格点多边形,求出内部点数in,边上点数on,和面积S. 思路:运用的定理很多. 1.皮克定理:S=in+on/2-1,即in=(2*S+2-on)/2. 2. ...

  4. pku 2954 Triangle pku 1265 Area Pick定理的应用 + 叉积求多边形面积

    Pick定理证明:http://translate.google.com/translate?u=http://episte.math.ntu.edu.tw/articles/sm/sm_25_10_ ...

  5. 【SSL.1213】多边形面积(difficult)(计算几何)

    多边形面积(difficult) Description 中学生杰利陶醉于数学研究.他思考的问题对于专家而言可能太多简单,但作为一个15岁的业余爱好者,他做得非常棒.他太热衷思考数学问题,以至于轻易就 ...

  6. 多边形面积_ssl1213_计算几何

    Description 中学生杰利陶醉于数学研究.他思考的问题对于专家而言可能太多简单,但作为一个15岁的业余爱好者,他做得非常棒.他太热衷思考数学问题,以至于轻易就学会了试用数学方法来解决碰到的问题 ...

  7. 观星(计算几何/凸包/多边形面积)

    观星 对于平面上有n个点分为三类,要求寻找一个三角形,三个顶点分别属于这三类,求解最大面积. N<=3000 首先考虑到O(n2)O(n^2)O(n2)的枚举,然后对于另外一种考虑这个点的位置, ...

  8. POJ 1265 Area(Pick定理)

    POJ 1265 Area(Pick定理) 根据pick定理来做 皮克定理是指一个计算点阵中顶点在格点上的多边形面积公式该公式可以表示为 2S=2a+b−2 其中a表示多边形内部的点数,b表示多边形边 ...

  9. poj 1265 Area(pick 定理)

    链接:poj 1265 题意:从原点出发.给出一些dx,dy移动增量,终于形成一个多边形, 求多边形内部的格点数目,边上的格点数目 .以及面积. 补充知识: 1.以格子点为顶点的线段.覆盖的点的个数为 ...

  10. POJ 1265 pick 定理

    首先来脑补一下几个几何知识 非常奇妙的证明:图形必在格点之外 从cut-the-knot上看到的. 问题:     设想一个平面上布满间距为1的横纵直线,形成由一个个1×1正方形组成的网格.任意给一个 ...

最新文章

  1. python和linux哪个难学-“慢”下来的 Python 要怎么竞争?
  2. PAT甲级1012 The Best Rank :[C++题解]4个成绩取排名最低:排序、二分(好题)
  3. SQLite 创建表(http://www.w3cschool.cc/sqlite/sqlite-create-table.html)
  4. 文档丨Oracle 20c 概念手册
  5. HC大会,华为联合合作伙伴发布一站式物联网IoT开发工具小熊派BearPi
  6. 达内python第二阶段月考_oracle练习题 达内第二次月考题
  7. 一、初识函数定义与调用
  8. Pandas Series interpolate
  9. 析砂性土层php泥浆护壁,砂卵石层钻探护壁工艺分析
  10. 初中计算机ps教程,初中信息技术《认识Photoshop CS2的工作界面》教案
  11. 使用npm发布自己的npm组件包
  12. 安装AdventureWorks2008R2
  13. drupal 电商网站_为什么小型企业应该为电子商务网站选择Drupal
  14. python自动化(六)持续集成:2.Jenkins技术讲解
  15. 使用 html 写一个表达爱意的网页
  16. 服务端开发基础知识点
  17. 内事不决问张昭,外事不决问周瑜,“ 排序 ”不决问威少
  18. IE浏览器快速切换各个版本
  19. 逻辑代数的基本公式和常用公式基本定理
  20. 大数据应用及未来展望

热门文章

  1. error: x264_bit_depth undeclared (first use in this function) did you mean x264_picture_t
  2. 泰山游记:道阻且长,活着走出
  3. Linux系统剪贴板getContents()非常之慢
  4. MySql解决办法:2004:Can't create TCP/IP socket (24)
  5. AlphaGo程序出BUG了——第四局观感
  6. 公交车宜配备逃生绳索
  7. ios不能保存png_ios转一加8t使用体验随用随更新
  8. android 高德amap开发一(地图创建)
  9. matlab 输入时间,求助:在MATLAB里如何输入时间序列中的时间
  10. data1 op data2 c语言,编写一个简单计算器程序,输入格式为:data1 op data2。其中data1和data2是参加...