hdu 4216 Computational Geometry? · dp

题解

题意：公主在(0,0)，炸弹在(x₀,y₀)，给n个向量，可以选择向量之间的使用顺序，也可以选择取向量的反方向(-x,-y)，问一系列动作后，公主离炸弹最远多远

贪心失败：
原想法：既然向量之间的先后顺序不影响最终到达的点，那每次都根据条件使得公主离炸点距离最远，遂，wa
忘了距离和向量之间的夹角有关，噗，

后看别人题解：
正解dp
维护两个数组
mx[i][j]表示前i个向量里，横坐标为j的最大纵坐标
mi[i][j]表示前i个向量里，横坐标为j的最小纵坐标

并且由于向量坐标可能为负的，还需要加上一个基准值m，确保不会越界
转移的时候要看看上一个状态能否转移

错误总结：
吐血，
失策了，memset写不了-0x3f3f3f3f，
数组开小了，wa了一发，
试着能不能开一维的数组，wa了，后来想想，不对，必须二维，一维无法保证最远的时候是所有向量全部用上了

---

---

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
const int N = 5e4 + 10;
const int INF = 0x3f3f3f3f;
int n, m = 22500, k;
int mx[105][N], mi[105][N];//y轴正方向 y轴负方向

int dis(int x1, int y1, int x2, int y2) {
    return (x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2);
}

int main() {
    ios::sync_with_stdio(0);

    int T, x0, y0;
    cin >> T;
    for (int cs = 1; cs <= T; ++cs) {
        printf("Case %d: ", cs);
        cin >> n;
        for (int i = 0; i <= 100; ++i) {
            for (int j = 0; j <= 45000; ++j) {
                mx[i][j] = -INF;
                mi[i][j] = INF;
            }
        }
        cin >> x0 >> y0;
        
        mx[0][0 + m] = 0;
        mi[0][0 + m] = 0;

        for (int i = 1,x,y; i <= n; ++i) {
            cin >> x >> y;
            for (int j = -22000; j <= 22000; ++j) {
                if (mx[i - 1][j + m] != -INF)
                    mx[i][j + m + x] = max(mx[i][j + m + x], mx[i - 1][j + m] + y);

                if (mx[i - 1][j + m] != -INF)
                    mx[i][j + m - x] = max(mx[i][j + m - x], mx[i - 1][j + m] - y);

                if (mi[i - 1][j + m] != INF)
                    mi[i][j + m + x] = min(mi[i][j + m + x], mi[i - 1][j + m] + y);

                if (mi[i - 1][j + m] != INF)
                    mi[i][j + m - x] = min(mi[i][j + m - x], mi[i - 1][j + m] - y);

            }
        }

        int ans = 0;
        for (int i = -22000; i <= 22000; ++i) {
            if (mx[n][i + m] != -INF) {
                ans = max(ans, dis(i, mx[n][i + m], x0, y0));
            }
            if (mi[n][i + m] != INF) {
                ans = max(ans, dis(i, mi[n][i + m], x0, y0));
            }
        }
        printf("%.3f\n", sqrt(1.0 * ans));
    }
    return 0;
}