UVa 609 Metal Cutting

原创于 2026-06-27 08:19:49 发布 · 122 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

GEO检测

标签

#算法 #C++ #青少年编程 #深度优先

回溯法同时被 2 个专栏收录

119 篇文章

订阅专栏

计算几何

97 篇文章

订阅专栏

题目描述

给定一块 $\times m$ 的矩形金属板（四个顶点分别为 $(0, 0)$ 、 $(0, m)$ 、 $(n, m)$ 、 $(n, 0)$ ），需要从中切割出一个凸多边形（顶点数 $\le 8$ ）。切割机每次只能沿一条直线完全穿透当前金属板（不能中途转向），且所有切割必须沿着多边形的边进行。要求确定切割这些边的最优顺序，使得切割总长度最小。输出最小总长度，精确到小数点后 $3$ 位。

输入格式

第一行为整数 $N$ ，表示数据集个数。接下来每个数据集开头有一个空行（实际输入中可能存在，但用 cin 读取整数可忽略空白）。每个数据集第一行为两个整数 $n$ 和 $m$ （ $\le 500$ ）。第二行为整数 $p$ （ $\le p \le 8$ ），表示多边形顶点数。接下来 $p$ 行，每行两个整数 $x$ 和 $y$ （ $0 < x < n$ ， $0 < y < m$ ），按顺时针顺序给出多边形的顶点。多边形不自交，且不存在三个连续顶点共线。数据集之间有一个空行。

输出格式

对于每个数据集，输出一行：

Minimum total length = XXX

其中 XXX 为最小总长度，保留 $3$ 位小数。每个数据集输出之间用一个空行分隔。

样例

输入

输出

Minimum total length = 312.575

题目分析

切割操作要求每次沿多边形的一条边直线切穿当前金属板。由于多边形在矩形内部，且每次切割都沿着多边形的边进行，实际切割长度取决于该刀切割时，当前剩余金属区域（包含目标多边形）与切割线的交集长度。该交集是一条线段，其端点位于当前金属区域的边界（包括矩形边界和已经切过的多边形边）上。

由于切割顺序任意，我们需要尝试所有可能的边的切割顺序，计算每种顺序下的总切割长度，取最小值。因为 $\le 8$ ，枚举全排列 $p!$ 是可行的（ $8! = 40320$ ，每次模拟 $O(p^2)$ 级别的计算，完全可接受）。

解题思路

数据结构与几何计算

用 point 结构体表示二维点。
用 line 结构体表示一条有向线段，存储端点 p1、p2 以及直线方程系数 $a, b, c$ （ $a x + b y + c = 0$ ），便于求交点。
提供函数 getIntersection 计算两条直线的交点，若平行则返回 false。
提供函数 contains(a, b, c, d) 判断两个点 $a$ 和 $b$ 是否分别位于线段 $c d$ 的两侧（或端点处），用于筛选有效交点。

模拟一次切割

维护一个集合 metal，存储当前金属区域的边界线段。初始时包含矩形的四条边。对于当前要切割的多边形边 $e$ ，计算它与 metal 中所有线段的交点。所有交点中，位于 $e$ 上的两个端点之间的长度即为本次切割的有效长度。具体做法为：

对 metal 中每条线段，调用 getIntersection 求与 $e$ 所在直线的交点。
若存在交点，将其存入临时数组。
遍历交点对 $p_i, p_j)$ ，使用 contains 判断它们是否分别位于 $e$ 的两端，若是，则计算两点距离并取最小值（理论上只有一对有效，但为安全取最小）。
将该有效距离累加到总长度中，然后将 $e$ 加入 metal，供后续切割使用。

枚举所有顺序

将多边形的 $p$ 条边存入数组 sides，并按编号排序（便于 next_permutation 全排列）。
使用 do-while 循环枚举 sides 的每一种排列。
对每种排列，重置 metal 为矩形四边，然后按排列顺序依次执行切割模拟，累加长度。
若当前累加长度已经超过已知最优值，可提前剪枝（代码中已实现）。
记录所有排列中的最小总长度，输出。

正确性说明

由于每次切割必须沿多边形边进行，且不同顺序会影响每刀实际切割长度（因为已切割的边会改变剩余金属的边界），因此最优解必然对应某一种切割顺序。枚举所有顺序并模拟切割过程，即可找到全局最优。

复杂度分析

每条边切割时，与当前 metal 中最多 $(4 + p)$ 条线段求交点， $\le 8$ ，故每次切割复杂度 $O (p)$ 。
每种顺序有 $p$ 次切割，总复杂度 $\cdot p^2)$ 。当 $p = 8$ 时， $\times 64 \approx 2.6 \times 10^6$ ，运行时间可接受。
空间复杂度 $O (p)$ 。

代码实现

// Metal Cutting
// UVa ID: 609
// Verdict: Accepted
// Submission Date: 2019-01-24
// UVa Run Time: 0.800s
//
// 版权所有（C）2019，邱秋。metaphysis # yeah dot net

#include <bits/stdc++.h>
using namespace std;

const double EPSILON = 1e-7;

struct point {
    double x, y;
    point (double x = 0, double y = 0): x(x), y(y) {}
};

typedef vector<point> polygon;

struct line {
    int index;
    point p1, p2;
    double a, b, c;
    bool operator<(const line &u) const {
        return index < u.index;
    }
};

line getLine(int index, point p, point q) {
    line u;
    u.index = index;
    u.a = q.y - p.y, u.b = p.x - q.x, u.c = p.y * (q.x - p.x) - p.x * (q.y - p.y);
    u.p1 = p, u.p2 = q;
    return u;
}

bool getIntersection(line p, line q, point &pi) {
    if (fabs(p.a * q.b - q.a * p.b) < EPSILON) return false;  
    pi.x = (q.c * p.b - p.c * q.b) / (p.a * q.b - q.a * p.b);  
    pi.y = (q.c * p.a - p.c * q.a) / (p.b * q.a - q.b * p.a);  
    return true;
}

double getDist(point p1, point p2) {
    return sqrt(pow(p1.x - p2.x, 2) + pow(p1.y - p2.y, 2));
}

inline bool contains(point a, point b, point c, point d) {
    int xup, xdown, yup, ydown;
    xup = xdown = yup = ydown = 0;
    double maxx = max(c.x, d.x), minx = min(c.x, d.x);
    if (a.x - maxx > -EPSILON || b.x - maxx > -EPSILON) xup = 1;
    if (a.x - minx < EPSILON || b.x - minx < EPSILON) xdown = 1;
    double maxy = max(c.y, d.y), miny = min(c.y, d.y);
    if (a.y - maxy > -EPSILON || b.y - maxy > -EPSILON) yup = 1;
    if (a.y - miny < EPSILON || b.y - miny < EPSILON) ydown = 1;
    return xup + xdown + yup + ydown == 4;
}

vector<line> sides, metal;

double cut(line u) {
    point p;
    vector<point> ps;
    for (auto v : metal)
        if (getIntersection(u, v, p))
            ps.push_back(p);
    double cutted = 1e20;
    for (int i = 0; i < ps.size(); i++)
        for (int j = i + 1; j < ps.size(); j++)
            if (contains(ps[i], ps[j], u.p1, u.p2))
                cutted = min(cutted, getDist(ps[i], ps[j]));
    return cutted;
}

int main(int argc, char *argv[]) {
    cin.tie(0), cout.tie(0), ios::sync_with_stdio(false);
    int cases;
    cin >> cases;
    for (int cs = 1; cs <= cases; cs++) {
        int p;
        double n, m, x, y;
        polygon pg;
        cin >> n >> m >> p;
        for (int i = 0; i < p; i++) {
            cin >> x >> y;
            pg.push_back(point(x, y));
        }
        sides.clear();
        for (int i = 0; i < p; i++)
            sides.push_back(getLine(i, pg[i], pg[(i + 1) % p]));
        sort(sides.begin(), sides.end());
        double minL = 1e20;
        do {
            metal.clear();
            metal.push_back(getLine(0, point(0, 0), point(0, m)));
            metal.push_back(getLine(1, point(0, m), point(n, m)));
            metal.push_back(getLine(2, point(n, m), point(n, 0)));
            metal.push_back(getLine(3, point(n, 0), point(0, 0)));
            double L = 0;
            for (int i = 0; i < p; i++) {
                L += cut(sides[i]);
                metal.push_back(sides[i]);
                if (L > minL) break;
            }
            minL = min(minL, L);
        } while (next_permutation(sides.begin(), sides.end()));
        if (cs > 1) cout << '\n';
        cout << "Minimum total length = " << fixed << setprecision(3) << minL << '\n';
    }
    return 0;
}

总结

本题通过枚举切割顺序并模拟每刀的有效切割长度，求得最小总切割长度。核心是利用几何计算求直线与多边形边界的交点，并准确度量每次切割在剩余金属中的实际长度。由于 $p$ 很小，全排列枚举是可行且直接的解法。该题也启示我们，在几何切割问题中，顺序优化往往可以通过暴力枚举所有排列来解决，前提是问题规模足够小。