【BZOJ4439】[Swerc2015]Landscaping【最小割】
2016-05-13 16:58
453 查看
【题目链接】
题解:
第一眼觉得是裸最小割模型,写了之后发现样例跑出个5000...
问了神犇,才知道:
1 相同地之间也得加容量为A的边,因为一开始的高低地不能决定最后策略的高低地。
2 要加双向边,因为一开始并不知道某个点应该属于S集还是T集。
复杂度:
时间复杂度:O(maxflow(n, m)),空间复杂度O(n + m)
1WA:
相同地之间也得加容量为A的边,而且得加双向边。
GET:
原来我一直不会最小割呀。
/* Telekinetic Forest Guard */
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long LL;
const int maxg = 55, maxn = 2505, maxm = 50005, maxq = 100000, inf = 0x3f3f3f3f;
int n, m, head[maxn], cur[maxn], cnt, bg, ed, A, B, depth[maxn], q[maxq];
char str[maxg];
struct _edge {
int v, w, next;
} g[maxm << 1];
inline void add(int u, int v, int w) {
g[cnt] = (_edge){v, w, head[u]};
head[u] = cnt++;
}
inline void insert(int u, int v, int w) {
add(u, v, w); add(v, u, 0);
}
inline bool bfs() {
for(int i = 0; i <= ed; i++) depth[i] = -1;
int h = 0, t = 0, u, i; depth[q[t++] = bg] = 1;
while(h != t) for(i = head[u = q[h++]]; ~i; i = g[i].next) if(g[i].w && !~depth[g[i].v]) {
depth[g[i].v] = depth[u] + 1;
if(g[i].v == ed) return 1;
q[t++] = g[i].v;
}
return 0;
}
inline int dfs(int x, int flow) {
if(x == ed) return flow;
int left = flow;
for(int i = cur[x]; ~i; i = g[i].next) if(g[i].w && depth[g[i].v] == depth[x] + 1) {
int tmp = dfs(g[i].v, min(g[i].w, left));
left -= tmp; g[i].w -= tmp; g[i ^ 1].w += tmp;
if(g[i].w) cur[x] = i;
if(!left) return flow;
}
if(left == flow) depth[x] = -1;
return flow - left;
}
int main() {
scanf("%d%d%d%d", &n, &m, &A, &B); bg = 0; ed = n * m + 1;
for(int i = 0; i <= ed; i++) head[i] = -1; cnt = 0;
for(int i = 1; i <= n; i++) {
scanf("%s", str + 1);
for(int j = 1; j <= m; j++) {
if(str[j] == '#') insert(bg, (i - 1) * m + j, B);
else insert((i - 1) * m + j, ed, B);
if(i != n) insert((i - 1) * m + j, i * m + j, A), insert(i * m + j, (i - 1) * m + j, A);
if(j != m) insert((i - 1) * m + j, (i - 1) * m + j + 1, A), insert((i - 1) * m + j + 1, (i - 1) * m + j, A);
}
}
int ans = 0;
while(bfs()) {
for(int i = 0; i <= ed; i++) cur[i] = head[i];
ans += dfs(bg, inf);
}
printf("%d\n", ans);
return 0;
}
题解:
第一眼觉得是裸最小割模型,写了之后发现样例跑出个5000...
问了神犇,才知道:
1 相同地之间也得加容量为A的边,因为一开始的高低地不能决定最后策略的高低地。
2 要加双向边,因为一开始并不知道某个点应该属于S集还是T集。
复杂度:
时间复杂度:O(maxflow(n, m)),空间复杂度O(n + m)
1WA:
相同地之间也得加容量为A的边,而且得加双向边。
GET:
原来我一直不会最小割呀。
/* Telekinetic Forest Guard */
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long LL;
const int maxg = 55, maxn = 2505, maxm = 50005, maxq = 100000, inf = 0x3f3f3f3f;
int n, m, head[maxn], cur[maxn], cnt, bg, ed, A, B, depth[maxn], q[maxq];
char str[maxg];
struct _edge {
int v, w, next;
} g[maxm << 1];
inline void add(int u, int v, int w) {
g[cnt] = (_edge){v, w, head[u]};
head[u] = cnt++;
}
inline void insert(int u, int v, int w) {
add(u, v, w); add(v, u, 0);
}
inline bool bfs() {
for(int i = 0; i <= ed; i++) depth[i] = -1;
int h = 0, t = 0, u, i; depth[q[t++] = bg] = 1;
while(h != t) for(i = head[u = q[h++]]; ~i; i = g[i].next) if(g[i].w && !~depth[g[i].v]) {
depth[g[i].v] = depth[u] + 1;
if(g[i].v == ed) return 1;
q[t++] = g[i].v;
}
return 0;
}
inline int dfs(int x, int flow) {
if(x == ed) return flow;
int left = flow;
for(int i = cur[x]; ~i; i = g[i].next) if(g[i].w && depth[g[i].v] == depth[x] + 1) {
int tmp = dfs(g[i].v, min(g[i].w, left));
left -= tmp; g[i].w -= tmp; g[i ^ 1].w += tmp;
if(g[i].w) cur[x] = i;
if(!left) return flow;
}
if(left == flow) depth[x] = -1;
return flow - left;
}
int main() {
scanf("%d%d%d%d", &n, &m, &A, &B); bg = 0; ed = n * m + 1;
for(int i = 0; i <= ed; i++) head[i] = -1; cnt = 0;
for(int i = 1; i <= n; i++) {
scanf("%s", str + 1);
for(int j = 1; j <= m; j++) {
if(str[j] == '#') insert(bg, (i - 1) * m + j, B);
else insert((i - 1) * m + j, ed, B);
if(i != n) insert((i - 1) * m + j, i * m + j, A), insert(i * m + j, (i - 1) * m + j, A);
if(j != m) insert((i - 1) * m + j, (i - 1) * m + j + 1, A), insert((i - 1) * m + j + 1, (i - 1) * m + j, A);
}
}
int ans = 0;
while(bfs()) {
for(int i = 0; i <= ed; i++) cur[i] = head[i];
ans += dfs(bg, inf);
}
printf("%d\n", ans);
return 0;
}
相关文章推荐
- CodeForces - 96B Lucky Numbers (easy)
- POJ 1696 Space Ant Graham-scan
- 励志学习文章精选
- hdu1495非常可乐(BFS)
- 机器学习开源代码大集合
- 深度学习BP算法的推导附加RNN-LSTM的推导
- myeclipse中打开java文件中文乱码
- 微信WeixinJSBridge API
- linux tomcat 配置
- 测试用例与软件测试基础
- java并发编程实战-对象的共享
- Problem K
- 表连接 join和(+)、union和uion all
- easyUi --修改
- http-server 超轻量级web服务器
- 阻塞、非阻塞、异步、同步以及select/poll和epoll
- 在线画时序图的工具
- Ubuntu SSH 私匙和公匙的产生原理阐述
- mysql tcp 4层负载
- windows异常处理__try __except