题意
给出一个区间 ,允许进行如下操作:
- 将 转为 或者
- 将 转为 或者
且保证
但是给出了一系列的限制 ,表示当前区间为 时,限制当前的区间不能进行操作 (dir = L
)或者操作 (dir = R
),而启用这个限制则需要 的费用
你可以选择是否启用这个限制
询问最少需要花费多少来实现不能将区间 转变为 的区间。
分析
从 能否转变为 可以通过最短路来求算。
但是无法求知当最短路无法到达时(即题目要求的不能转变)最少需要多少的限制条件,而这些条件又是什么。
所以采用最大流来解决
最大流
画出网格图
将所有可以转换的两个状态之间用边连接,如果有提供限制的,将流量限制为费用,如果没有限制的,则设置为
对于整个矩阵而言,只需要一半的点用于建图,所以将汇点放在另外一半点中。所有 的点与汇点连接,而源点为
对于样例可以得到如下图
样例:
3 4
1 3 L 10
1 3 R 3
1 2 L 1
1 2 R 1
补充,图片漏画了的连线,其流量为
可以直接通过最大流求出答案
但是会TLE
对偶图
对偶图Wikipedia(https://en.wikipedia.org/wiki/Dual_graph)
通过对偶图,可以快速的将一张网格网络图求最大流转为求最短路
关于对偶图的解释请自行查阅资料
在原图上绘制对偶图得到
将对偶图中有用的元素将其分离出来得到
(图中未注明边权的边均为 )
可以通过最短路快速得到解
AC code
#include <bits/stdc++.h> using namespace std; #define ll long long const int maxn = 510; int n, m; ll dis[maxn * maxn]; char si; vector<pair<ll, int>> G[maxn * maxn]; void addedge(int u, int v, int cost) { G[u].push_back({cost, v}); } ll dijkstra(int s, int t) { memset(dis, 0x3f, sizeof(dis)); dis[s] = 0; priority_queue<pair<ll, int>, vector<pair<ll, int>>, greater<pair<ll, int>>> q; q.push({0ll, s}); while (!q.empty()) { ll u = q.top().second, c = q.top().first; q.pop(); if (dis[u] < c)continue; for (auto i : G[u]) { ll cc = i.first, v = i.second; if (dis[v] > dis[u] + cc) { dis[v] = dis[u] + cc; q.push({dis[v], v}); } } } return dis[t]; } inline int id(int x, int y) { return x * (n + 3) + y; } void solve() { cin >> n >> m; for (int i = 0; i < m; ++i) { int u, v, w; char c; cin >> u >> v >> c >> w; if (c == 'L') { addedge(id(u, v), id(u, v + 1), w); addedge(id(u, v + 1), id(u, v), w); } else { addedge(id(u, v), id(u - 1, v), w); addedge(id(u - 1, v), id(u, v), w); } } for (int i = 1; i <= n; ++i) { addedge(id(0, 0), id(0, i), 0); addedge(id(i, n + 1), id(n + 1, n + 1), 0); } dijkstra(id(0, 0), id(n + 1, n + 1)); if (dis[id(n + 1, n + 1)] >= 0x3f3f3f3f3f3f3f3f) cout << -1 << endl; else cout << dis[id(n + 1, n + 1)] << endl; } signed main() { ios_base::sync_with_stdio(false); cin.tie(nullptr); cout.tie(nullptr); #ifdef ACM_LOCAL freopen("in.txt", "r", stdin); freopen("out.txt", "w", stdout); int test_index_for_debug = 1; char acm_local_for_debug; while (cin >> acm_local_for_debug) { if (acm_local_for_debug == '$') exit(0); cin.putback(acm_local_for_debug); if (test_index_for_debug > 20) { throw runtime_error("Check the stdin!!!"); } auto start_clock_for_debug = clock(); solve(); auto end_clock_for_debug = clock(); cout << "Test " << test_index_for_debug << " successful" << endl; cerr << "Test " << test_index_for_debug++ << " Run Time: " << double(end_clock_for_debug - start_clock_for_debug) / CLOCKS_PER_SEC << "s" << endl; cout << "--------------------------------------------------" << endl; } #else solve(); #endif return 0; }