poj2836 状压dp

    xiaoxiao2024-07-24  9

    题目链接:点击打开链接

    题意:

    给你 n 个点,然后用矩形去覆盖这些点;

    每个矩形最少覆盖两个点;

    点可以被重复覆盖;

    问最小的覆盖面积是多少;

    其中重复面积要重复计算;

    理解:

    状态压缩dp;

    1 超时的做法:

    递推式含义:dp[i] 表示 i 这种矩形状态下的最小值;

    递推式:dp[v[i] | v[k]] = max(dp[v[i] | v[k]], dp[v[i]] + dp[v[k]]);

    其中:

    v[i] 表示 i 符合条件;

    即:每个矩形最少覆盖两个点;

    即:i - (i & (-i)) > 0;

    初始值:就是计算构成最初符合条件的矩形状态的面积;

    但是毫无疑问会超时;

    其中 v[i] 的状态数太多;

    2 正确的做法:

    就是把 v[i] 重新定义;

    因为是两个点组成矩形;

    所以根据点枚举矩形就行;

    这样比上面符合条件的情况要少得多;

    递推式一样;

    就是把 v[i] 变一下;

    定义为:v[k] = (1 << i) | (1 << j);

    如果第 p 点包含在 v[k] 矩阵就 v[k] = v[k] | (1 << p);

    这样就不会超时了;

    代码如下:

    #include <cstdio> #include <cstdlib> #include <cstring> #include <cmath> #include <ctime> #include <iostream> #include <algorithm> #include <vector> #include <string> #include <map> #include <set> #include <queue> #include <stack> using namespace std; typedef long long LL; typedef pair<int, int> PII; const int MIN_INF = 1e-7; const int MAX_INF = (1e9) + 7; #define X first #define Y second int dp[(1 << 15) + 10]; int v[(1 << 15) + 10]; PII w[20]; int area(int n) { int x1 = 10000, y1 = 10000, x2 = -10000, y2 = -10000; for (int i = 0; i < 17; ++i) { if ((1 << i) & n) { x1 = min(x1, w[i].X); y1 = min(y1, w[i].Y); x2 = max(x2, w[i].X); y2 = max(y2, w[i].Y); } } int x = 1, y = 1; if (x1 != x2) { x = x2 - x1; } if (y1 != y2) { y = y2 - y1; }//cout << x * y << endl; return x * y; } int main() { int n; while (cin >> n && n) { memset(dp, 1, sizeof(dp)); for (int i = 0; i < n; ++i) { cin >> w[i].X >> w[i].Y; } int k = 0; for (int i = 0; i < n; ++i) { for (int j = i + 1; j < n; ++j) { v[k] = (1 << i) | (1 << j); for (int l = 0; l < n; ++l) { if (w[l].X <= max(w[i].X, w[j].X) && w[l].X >= min(w[i].X, w[j].X) && w[l].Y <= max(w[i].Y, w[j].Y) && w[l].Y >= min(w[i].Y, w[j].Y)) { v[k] |= (1 << l); } } dp[v[k]] = area(v[k]); ++k; } } for (int i = 0; i < (1 << n); ++i) { for (int j = 0; j < k; ++j) { dp[i | v[j]] = min(dp[i | v[j]], dp[i] + dp[v[j]]); } } cout << dp[(1 << n) - 1] << endl; } return 0; }

    转载请注明原文地址: https://ju.6miu.com/read-1290992.html
    最新回复(0)