首页 > 代码库 > uva 1331 - Minimax Triangulation(dp)
uva 1331 - Minimax Triangulation(dp)
题目链接:uva 1331 - Minimax Triangulation
题目大意:按照顺时针或者逆时针的顺序给出多边的点,要将这个多边形分解成n-2个三角形,要求使得这些三角行中面积最大的三角形面积尽量小,求最小值。
解题思路:状态很好想,dp[i][j]表示从第i个点到第j个点,划分成j-i-1个三角形的最优解,然后每次转移时,枚举长度和左边界始点,那么根据长度和左边界点就可以知道右边界点,然后枚举左边界和右边界中间的点k,dp[i][j] = min(dp[i][j], max(max(dp[i][k], dp[k][j]), Area(i, k, j)).但是有一个问题,即i,k,j三点围成的三角形是否符合要求,判断的条件即为是否存在除i,k,j三点外的一点位于三角形中,有面积法判断。
然后其实我还一直纠结说凹多边形的时候怎么处理掉得,如图
BCD的组成的三角形式合法的,但是后来想了想,虽然BCD组成的三角形是合法的,但是不会影响到最后的答案。
。
#include <cstdio> #include <cstring> #include <cmath> #include <algorithm> using namespace std; const int N = 100; const double INF = 0x3f3f3f3f3f3f; const double eps = 1e-9; struct point { double x, y; void get() { scanf("%lf%lf", &x, &y); } }p[N]; int n; double dp[N][N]; double area (point a, point b, point c) { return fabs((b.x-a.x)*(c.y-a.y) - (c.x-a.x)*(b.y-a.y))/2; } bool judge (int a, int b, int c) { double cur = area(p[a], p[b], p[c]); for (int i = 0; i < n; i++) { if (i == a || i == b || i == c) continue; double tmp = area(p[a], p[b], p[i]) + area(p[b], p[c], p[i]) + area(p[c], p[a], p[i]); if (fabs(tmp - cur) < eps) return false; } return true; } double solve () { for (int i = 0; i < 2; i++) { for (int j = 0; j < n; j++) dp[j][(j+i)%n] = 0; } for (int i = 0; i < n; i++) dp[i][(i+2)%n] = area(p[i], p[(i+1)%n], p[(i+2)%n]); for (int k = 3; k < n; k++) { for (int i = 0; i < n; i++) { int t = (i + k) % n; dp[i][t] = INF; for (int j = (i+1)%n; j != t; j = (j+1)%n) { if (judge(i, t, j)) dp[i][t] = min(dp[i][t], max(max(dp[i][j], dp[j][t]), area(p[i], p[j], p[t]))); } } } double ans = INF; for (int i = 0; i < n; i++) ans = min (ans, dp[i][(i+n-1)%n]); return ans; } int main () { int cas; scanf("%d", &cas); while (cas--) { scanf("%d", &n); for (int i = 0; i < n; i++) p[i].get(); printf("%.1lf\n", solve()); } return 0; }
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。