uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1058

　　半平面交求面积最值。直接枚举C(20,8)的所有情况即可。

代码如下：

 #include <cstdio>

 #include <cstring>

 #include <iostream>

 #include <algorithm>

 #include <cmath>

 #include <vector>

 using namespace std;

 const int N = ;

 struct Point {

     double x, y;

     Point() {}

     Point(double x, double y) : x(x), y(y) {}

 } pt[N], poly[N];

 template<class T> T sqr(T x) { return x * x;}

 typedef Point Vec;

 Vec operator + (Vec a, Vec b) { return Vec(a.x + b.x, a.y + b.y);}

 Vec operator - (Vec a, Vec b) { return Vec(a.x - b.x, a.y - b.y);}

 Vec operator * (Vec a, double p) { return Vec(a.x * p, a.y * p);}

 Vec operator / (Vec a, double p) { return Vec(a.x / p, a.y / p);}

 const double EPS = 1e-;

 inline int sgn(double x) { return (x > EPS) - (x < -EPS);}

 inline double dotDet(Vec a, Vec b) { return a.x * b.x + a.y * b.y;}

 inline double crossDet(Vec a, Vec b) { return a.x * b.y - a.y * b.x;}

 inline double dotDet(Point o, Point a, Point b) { return dotDet(a - o, b - o);}

 inline double crossDet(Point o, Point a, Point b) { return crossDet(a - o, b - o);}

 inline double vecLen(Vec x) { return sqrt(dotDet(x, x));}

 inline Vec normal(Vec x) { return Vec(-x.y, x.x) / vecLen(x);}

 struct DLine {

     Point p;

     Vec v;

     double ang;

     DLine() {}

     DLine(Point p, Vec v) : p(p), v(v) { ang = atan2(v.y, v.x);}

     bool operator < (const DLine &L) const { return ang < L.ang;}

     DLine move(double x) { return DLine(p + normal(v) * x, v);}

 } dl[N];

 inline bool onLeft(DLine L, Point p) { return crossDet(L.v, p - L.p) > ;}

 inline Point dLineIntersect(DLine a, DLine b) { return a.p + a.v * (crossDet(b.v, a.p - b.p) / crossDet(a.v, b.v));}

 struct Poly {

     vector<Point> pt;

     Poly() { pt.clear();}

     ~Poly() {}

     Poly(vector<Point> &pt) : pt(pt) {}

     Point operator [] (int x) const { return pt[x];}

     int size() { return pt.size();}

     double area() {

         if (pt.size() < ) return 0.0;

         double ret = 0.0;

         for (int i = , sz = pt.size(); i < sz; i++)

             ret += crossDet(pt[i], pt[(i + ) % sz]);

         return fabs(ret / 2.0);

     }

 } ;

 Point p[N];

 DLine q[N];

 int halfPlane(DLine *L, int n) {

     sort(L, L + n);

     int fi, la;

     q[fi = la = ] = L[];

     for (int i = ; i < n; i++) {

         while (fi < la && !onLeft(L[i], p[la - ])) la--;

         while (fi < la && !onLeft(L[i], p[fi])) fi++;

         q[++la] = L[i];

         if (fabs(crossDet(q[la].v, q[la - ].v)) < EPS) {

             la--;

             if (onLeft(q[la], L[i].p)) q[la] = L[i];

         }

         if (fi < la) p[la - ] = dLineIntersect(q[la - ], q[la]);

     }

     while (fi < la && !onLeft(q[fi], p[la - ])) la--;

     if (la < fi) return ;

     p[la] = dLineIntersect(q[la], q[fi]);

     int ret = ;

     for (int i = fi; i <= la; i++) poly[ret++] = p[i];

     return ret;

 }

 //Poly halfPlane(DLine *L, int n) {

 //    Poly ret = Poly();

 //    sort(L, L + n);

 //    int fi, la;

 //    Point *p = new Point[n];

 //    DLine *q = new DLine[n];

 //    q[fi = la = 0] = L[0];

 //    for (int i = 1; i < n; i++) {

 //        while (fi < la && !onLeft(L[i], p[la - 1])) la--;

 //        while (fi < la && !onLeft(L[i], p[fi])) fi++;

 //        q[++la] = L[i];

 //        if (fabs(crossDet(q[la].v, q[la - 1].v)) < EPS) {

 //            la--;

 //            if (onLeft(q[la], L[i].p)) q[la] = L[i];

 //        }

 //        if (fi < la) p[la - 1] = dLineIntersect(q[la - 1], q[la]);

 //    }

 //    while (fi < la && !onLeft(q[fi], p[la - 1])) la--;

 //    if (la < fi) return ret;

 //    p[la] = dLineIntersect(q[la], q[fi]);

 //    for (int i = fi; i <= la; i++) ret.pt.push_back(p[i]);

 //    return ret;

 //}

 double polyArea(Point *p, int n) {

     if (n < ) return 0.0;

     double sum = 0.0;

     p[n] = p[];

     Point O = Point(0.0, 0.0);

     for (int i = ; i < n; i++) sum += crossDet(O, p[i], p[i + ]);

     return fabs(sum / 2.0);

 }

 double work(int st, int n, double d) {

 //    cout << st << endl;

     for (int i = ; i < n; i++) dl[i] = DLine(pt[i], pt[i + ] - pt[i]).move((st & ( << i)) ==  ? 0.0 : d);

     int tmp = halfPlane(dl, n);

     return polyArea(poly, tmp);

 //    Poly tmp = halfPlane(dl, n);

 //    return tmp.area();

 }

 int cntBit(int x) {

     int cnt = ;

     while (x > ) {

         if (x & ) cnt++;

         x >>= ;

     }

     return cnt;

 }

 int main() {

 //    freopen("in", "r", stdin);

     int n, k;

     double d;

     while (~scanf("%d%d%lf", &n, &k, &d) && (n + k + d > EPS)) {

         for (int i = ; i < n; i++) scanf("%lf%lf", &pt[i].x, &pt[i].y);

         pt[n] = pt[];

         double ans = 0.0, area = polyArea(pt, n);

         for (int i = , end =  << n; i < end; i++) {

             if (cntBit(i) <= k) {

 //                cout << i << endl;

                 ans = max(ans, area - work(i, n, d));

             }

             if (sgn(ans - area) == ) break;

         }

         printf("%.2f\n", ans);

     }

     return ;

 }

能通过POJ的代码：

 #include <cstdio>

 #include <cstring>

 #include <iostream>

 #include <algorithm>

 #include <cmath>

 #include <vector>

 using namespace std;

 const int N = ;

 struct Point {

     double x, y;

     Point() {}

     Point(double x, double y) : x(x), y(y) {}

 } pt[N], poly[N];

 template<class T> T sqr(T x) { return x * x;}

 typedef Point Vec;

 Vec operator + (Vec a, Vec b) { return Vec(a.x + b.x, a.y + b.y);}

 Vec operator - (Vec a, Vec b) { return Vec(a.x - b.x, a.y - b.y);}

 Vec operator * (Vec a, double p) { return Vec(a.x * p, a.y * p);}

 Vec operator / (Vec a, double p) { return Vec(a.x / p, a.y / p);}

 const double EPS = 1e-;

 inline int sgn(double x) { return (x > EPS) - (x < -EPS);}

 inline double dotDet(Vec a, Vec b) { return a.x * b.x + a.y * b.y;}

 inline double crossDet(Vec a, Vec b) { return a.x * b.y - a.y * b.x;}

 inline double dotDet(Point o, Point a, Point b) { return dotDet(a - o, b - o);}

 inline double crossDet(Point o, Point a, Point b) { return crossDet(a - o, b - o);}

 inline double vecLen(Vec x) { return sqrt(dotDet(x, x));}

 inline Vec normal(Vec x) { return Vec(-x.y, x.x) / vecLen(x);}

 struct DLine {

     Point p;

     Vec v;

     double ang;

     DLine() {}

     DLine(Point p, Vec v) : p(p), v(v) { ang = atan2(v.y, v.x);}

     bool operator < (const DLine &L) const { return ang < L.ang;}

     DLine move(double x) { return DLine(p + normal(v) * x, v);}

 } dl[N];

 inline bool onLeft(DLine L, Point p) { return crossDet(L.v, p - L.p) > ;}

 inline Point dLineIntersect(DLine a, DLine b) { return a.p + a.v * (crossDet(b.v, a.p - b.p) / crossDet(a.v, b.v));}

 struct Poly {

     vector<Point> pt;

     Poly() { pt.clear();}

     ~Poly() {}

     Poly(vector<Point> &pt) : pt(pt) {}

     Point operator [] (int x) const { return pt[x];}

     int size() { return pt.size();}

     double area() {

         if (pt.size() < ) return 0.0;

         double ret = 0.0;

         for (int i = , sz = pt.size(); i < sz; i++)

             ret += crossDet(pt[i], pt[(i + ) % sz]);

         return fabs(ret / 2.0);

     }

 } ;

 Point p[N];

 DLine q[N], tmpDL[N];

 int halfPlane(DLine *L, int n) {

     for (int i = ; i < n; i++) tmpDL[i] = L[i];

     sort(L, L + n);

     int fi, la;

     q[fi = la = ] = L[];

     for (int i = ; i < n; i++) {

         while (fi < la && !onLeft(L[i], p[la - ])) la--;

         while (fi < la && !onLeft(L[i], p[fi])) fi++;

         q[++la] = L[i];

         if (fabs(crossDet(q[la].v, q[la - ].v)) < EPS) {

             la--;

             if (onLeft(q[la], L[i].p)) q[la] = L[i];

         }

         if (fi < la) p[la - ] = dLineIntersect(q[la - ], q[la]);

     }

     for (int i = ; i < n; i++) L[i] = tmpDL[i];

     while (fi < la && !onLeft(q[fi], p[la - ])) la--;

     if (la < fi) return ;

     p[la] = dLineIntersect(q[la], q[fi]);

     int ret = ;

     for (int i = fi; i <= la; i++) poly[ret++] = p[i];

     return ret;

 }

 double polyArea(Point *p, int n) {

     if (n < ) return 0.0;

     double sum = 0.0;

     p[n] = p[];

     Point O = Point(0.0, 0.0);

     for (int i = ; i < n; i++) sum += crossDet(O, p[i], p[i + ]);

     return fabs(sum / 2.0);

 }

 double work(int st, int n, double d) {

 //    cout << st << endl;

     for (int i = ; i < n; i++) dl[i] = DLine(pt[i], pt[i + ] - pt[i]).move((st & ( << i)) ==  ? 0.0 : d);

     int tmp = halfPlane(dl, n);

     return polyArea(poly, tmp);

 //    Poly tmp = halfPlane(dl, n);

 //    return tmp.area();

 }

 int cntBit(int x) {

     int cnt = ;

     while (x > ) {

         if (x & ) cnt++;

         x >>= ;

     }

     return cnt;

 }

 const double FINF = 1e100;

 double ans, d;

 void dfs(int id, int tt, int used, int k) {

     if (id >= tt) {

         int tmp = halfPlane(dl, tt);

         ans = min(ans, polyArea(poly, tmp));

         return ;

     }

     dl[id] = DLine(pt[id], pt[id + ] - pt[id]);

     dfs(id + , tt, used, k);

     if (used >= k) return ;

     dl[id] = DLine(pt[id], pt[id + ] - pt[id]).move(d);

     dfs(id + , tt, used + , k);

 }

 int main() {

 //    freopen("in", "r", stdin);

     int n, k;

     while (~scanf("%d%d%lf", &n, &k, &d) && (n + k + d > EPS)) {

         for (int i = ; i < n; i++) scanf("%lf%lf", &pt[i].x, &pt[i].y);

         pt[n] = pt[];

         double area = polyArea(pt, n);

         ans = FINF;

         dfs(, n, , k);

         printf("%.2f\n", area - ans);

     }

     return ;

 }

——written by Lyon

poj 1271 && uva 10117 Nice Milk （半平面交）的更多相关文章

  1. poj 2451 Uyuw's Concert （半平面交）

    2451 -- Uyuw's Concert 继续半平面交,这还是简单的半平面交求面积,不过输入用cin超时了一次. 代码如下: #include <cstdio> #include &l ...

  2. poj 2451 Uyuw's Concert（半平面交）

    Uyuw's Concert Time Limit: 6000MS   Memory Limit: 65536K Total Submissions: 8580   Accepted: 3227 De ...

  3. POJ - 1474 ：Video Surveillance （半平面交-求核）

    pro:顺时针给定多边形,问是否可以放一个监控,可以监控到所有地方,即问是否存在多边形的核. 此题如果两点在同一边界上(且没有被隔段),也可以相互看到. sol:求多边形是否有核.先给直线按角度排序, ...

  4. POJ 3384 Feng Shui（计算几何の半平面交+最远点对）

    Description Feng shui is the ancient Chinese practice of placement and arrangement of space to achie ...

  5. POJ 3525/UVA 1396 Most Distant Point from the Sea（二分+半平面交）

    Description The main land of Japan called Honshu is an island surrounded by the sea. In such an isla ...

  6. poj 3335（半平面交）

    链接:http://poj.org/problem?id=3335     //大牛们常说的测模板题 ------------------------------------------------- ...

  7. poj 1755 半平面交+不等式

    Triathlon Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 6461   Accepted: 1643 Descrip ...

  8. poj 1279 半平面交核面积

    Art Gallery Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 6668   Accepted: 2725 Descr ...

  9. poj 3335 Rotating Scoreboard（半平面交）

    Rotating Scoreboard Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 6420   Accepted: 25 ...

随机推荐

  1. Node.js调试技巧

    1. console.log 跟前端调试相同,通过一步步打印相关变量进行代码调试 2. 使用Node.js内置的调试器 通过node debug xxx.js来进行调试: [root@~/wade/n ...

  2. PHP配置环境中开启GD库

    下配置好的PHP环境中,GD库不像windows那样可以直接用,而是默认关闭,需要把它打开,去到php.ini文件中 找到php_gd2.dll把分号去掉即可.(注:GD库跟绘制二维码等有关)

  3. CSS Tools: Reset CSS

    样式初始化 /* http://meyerweb.com/eric/tools/css/reset/ v2.0 | 20110126 License: none (public domain) */ ...

  4. 【JZOJ3640】【COCI2014】utrka

    Mission 2<=N<=300,2<=M<=N∗(N−1) Solution SPFA. 由于只是二元关系,所以条件随便写. 具体来说,如果是u⇒v. 若v的最大领先时间还 ...

  5. Linux安装MariaDB（Mysql）和简单配置 mariadb

    Linux安装MariaDB(Mysql)和简单配置 1.安装MariaDB 安装命令 yum -y install mariadb mariadb-server 安装完成MariaDB,首先启动Ma ...

  6. OSGi教程：Resource API Specification

    此教程基于OSGi Core Release 7 OSGi Resource API规范 详细内容上面英文教程有详细解答 下面主要是一些个人见解,若有不当之处,欢迎指出: Resource:就是能够被 ...

  7. 下载安装APK(兼容Android7.0)

    我们使用手机的时候经常会看到应用程序提示升级,大部分应用内部都需要实现升级提醒和应用程序文件(APK文件)下载. 一般写法都差不多,比如在启动app的时候,通过api接口获得服务器最新的版本号,然后和 ...

  8. java8的stream系列教程之filter过滤集合的一些属性

    贴代码 List<Student> lists = new ArrayList<>(); Student student = new Student(); student.se ...

  9. oracle-ORA-01650错误

    Unable to extend rollback segment 原因:没有足够的撤销空间用来处理所有活动事务

  10. SDUT-3375_数据结构实验之查找三：树的种类统计

    数据结构实验之查找三:树的种类统计 Time Limit: 400 ms Memory Limit: 65536 KiB Problem Description 随着卫星成像技术的应用,自然资源研究机 ...