usaco 3.4 Electric Fence 皮克定理

 

一条直线((0,0),(n,m))上的格点数(包含端点）等于n与m的最大公约数+1

皮克定理: 面积A和内部格点数目i、边上格点数目b的关系：A = i + b/2 - 1

 

 

1. #include <iostream>
2. #include <cmath>
3. #include <algorithm>
4. using namespace std;
6. /*
7. PROG: fence9
8. LANG: C++
9. ID: heben991
10. */
13. struct point
14. {
15.     long long x,y;
16. };
19. long long Abs(long long x)
20. {
21.     if(x<0)x=-x;
22.     return x;
23. }
25. long long cross(point a, point b)
26. {
27.     return a.x*b.y - a.y*b.x;
28. }
30. //多变形面积
31. long long area(point p[], int n = 3)
32. {
33.     int i, j, k;
34.     long long s, ss = 0;
36.     for(i = 0; i < n; i++)
37.     {
38.     j = (i+1)%n;
39.     s = cross(p[i],p[j]);
40.     ss += s;
41.     }
43.     return Abs(ss)/2;
44. }
46. // a,b整点之间的整点个数，包含a,b
47. long long gcd(long long a, long long b)
48. {
49.     while(a && b)
50.     {
51.         a%=b;
52.         if(a) b%=a;
53.     }
54.     return a+b;
55. }
58. //一条直线((0,0),(n,m))上的格点数(包含端点）等于n与m的最大公约数+1
59. long long countgrid(point a, point b)
60. {
61.     long long x = Abs(a.x-b.x), y = Abs(a.y-b.y);
62.     return gcd(x,y);
63. }
65. //多变形内整点个数，不包含边界pku 2954 1256
67. //皮克定理: 面积A和内部格点数目i、
68. //边上格点数目b的关系：A = i + b/2 - 1
69. long long internal_grid(point p[], int n = 3)
70. {
71.     long long a = area(p,n), in, b=0, t;
72.     int i;
74.     for(i = 0; i < n; i++)
75.     {
76.         t = countgrid(p[i], p[(i+1)%n]);
77.         b += t;
78.     }
80.     in = a - b/2 + 1;
81.     return in;
82. }
85. int main()
86. {
87.     freopen("fence9.in", "r", stdin);
88.     freopen("fence9.out","w",stdout);
90.     point t[3] = {0};
92.     cin >> t[1].x >> t[1].y >> t[2].x;
94.     cout << internal_grid(t) << endl;
95.     return 0;
96. }