Not Escaping（哈希+dp）

Not Escaping

[Link](Problem - E - Codeforces)

题意

给你一个$n\times m$的楼，一些梯子(只有梯子允许跨层操作)的信息例如$a,b,c,d,h$表示从$(a,b)$到$(c,d)$会获得$h$滴血，给你$n$个数$a_n$,$a_i:在第i层每走一步损失多少滴血$，例如在第$i$层从$x_i\to x_j$会损失$(x_j-x_i)\times a_i$滴血。

问你从$(1,1)$走到$(n,m)$最少损失多少滴血$(可以是负数)$，如果无法从$(1,1)$到$(n,m)$输出NO ESCAPE。

思路

首先只有起点、终点、有梯子的点有用，因为我们一定是从起点走到某个梯子往上爬并走到下一个梯子，可以将这些点映射成从1开始的点(用个$map$即可)。

我们设$f[i]:走到第i个点最少损失多少滴血$，发现第$i$个点的值可以由下面的梯子到当前这个点(如果有这样的梯子的话)和当前层的其他点通过同层走减血转移过来，也就是第$i$个点可以由前面的点推出最优解且不会被后继影响，因此考虑$dp$。

对于第$i$层的某个点$x_i$分类来看对于所有在它左边的点设$x_j$，从$j$到$i$的代价就是$f[j]+(x_i-x_j)\times a_i\to f[j]-a_i\times x_j+a_i\times x_i$,发现对于第$i$层点$x_i$来说$a_i\times x_i$为定值，因此我们只需要求一个$1\sim i-1$里的$f[j]-a_i\times x_j$的前缀最小值。

对于某个在它右边的点$x_j$,从$j$到$i$的代价就是$f[j]+(x_j-x_i)\times a_i\to f[j]+a_i\times x_j-a_i\times x_i$,同理我们只需要在$i+1\sim 总长$里求一个后缀最小值。

最后到了第$n$层枚举每个点向终点跑的代价即可。

对于每个用到的没有出现过的映射一个$idx$, 对于每一行出现的点用一个$vector$存对应的列和映射索引值，对于每个梯子来说用$vector$来存仿照邻接表，用起点$idx$表示地址，存终点的$idx$和这个梯子的权重。

Code

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <set>
#include <queue>
#include <vector>
#include <map>
#include <bitset>
#include <unordered_map>
#include <cmath> 
#include <stack>
#include <iomanip>
#include <deque> 
#include <sstream>
#define x first
#define y second
#define debug(x) cout<<#x<<":"<<x<<endl;
using namespace std;
typedef long double ld;
typedef long long LL;
typedef pair<int, int> PII;
typedef pair<double, double> PDD;
typedef unsigned long long ULL;
const int N = 2e5 + 10, M = 2 * N, INF = 0x3f3f3f3f, mod = 1e9 + 7;
const double eps = 1e-8, pi = acos(-1), inf = 1e20;
#define tpyeinput int
inline char nc() {static char buf[1000000],*p1=buf,*p2=buf;return p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++;}
inline void read(tpyeinput &sum) {char ch=nc();sum=0;while(!(ch>='0'&&ch<='9')) ch=nc();while(ch>='0'&&ch<='9') sum=(sum<<3)+(sum<<1)+(ch-48),ch=nc();}
int dx[] = {-1, 0, 1, 0}, dy[] = {0, 1, 0, -1};
int n, m, k, q;
map<PII, int> id;
vector<PII> r[N];
int a[N], cnt;
LL f[N];
vector<PII> w[N];
int main() {
    ios::sync_with_stdio(false), cin.tie(0);
    int T;
    cin >> T;
    while (T -- ) {
        cin >> n >> m >> q;
        for (int i = 1; i <= n; i ++) cin >> a[i], r[i].clear();
        cnt = 0;
        id.clear();

        while (q --) {
            int a, b, c, d, h; 
            cin >> a >> b >> c >> d >> h;
            PII sta = {a, b}, ed ={c, d};
            if (!id[sta]) 
                id[sta] = ++cnt, r[a].push_back({b, cnt});
            if (!id[ed])
                 id[ed] = ++cnt, r[c].push_back({d, cnt});
            w[id[sta]].push_back({id[ed], h});
        }
    
        for (int i = 1; i <= cnt; i ++)
            f[i] = 1ll << 60;
        for (auto i : r[1]) 
            f[i.y] = 1ll * (i.x - 1) * a[1];
        
        for (int i = 1; i <= n; i ++) 
            if (r[i].size()) {
               int m = r[i].size(); 
               sort(r[i].begin(), r[i].end());
               LL M = 1ll << 60;
               for (int j = 0; j < m; j ++) {
                   M = min(M, f[r[i][j].y] - 1ll * r[i][j].x * a[i]);
                   f[r[i][j].y] = min(f[r[i][j].y], M + 1ll * r[i][j].x * a[i]);
               }
               M = 1ll << 61; // 多一位，因为下面的转移式位+权重，60的话会导致有的值统计不到
               for (int j = m - 1; j >= 0; j --) {
                   M = min(M, f[r[i][j].y] + 1ll * r[i][j].x * a[i]);
                   f[r[i][j].y] = min(f[r[i][j].y], M - 1ll * r[i][j].x * a[i]);
               }
               for (int j = 0; j < m; j ++)
                    if (f[r[i][j].y] < 1ll << 60ll)
                        for (auto k : w[r[i][j].y])
                            f[k.x] = min(f[k.x], f[r[i][j].y] - k.y);
            }
        LL res = 1ll << 60;
        for (auto i : r[n])
            res = min(res, f[i.y] + 1ll * (m - i.x) * a[n]);
        if (res < 1ll << 60) cout << res << endl;
        else cout << "NO ESCAPE" << endl;
        for (int i = 1; i <= cnt; i ++) 
            w[i].clear();
    }
    return 0;
}