HDU 3605 Escape
2012-08-17 22:19
169 查看
HDU_3605
之所以写这个题目是因为偶然看到了用匈牙利算法做二分图多重匹配这么个东西,看完之后发现只是在原有求二分图最大匹配的基础上对匈牙利算法做了一点拓展,代码还是很好写的,而且效率也很不错。
如果用网络流做的话,这么多点是承受不了的,但由于M很小,可以将N个点适应各个星球的M个数字综合起来看成一个二进制数,这样就能把N分成2^M类,这样图的顶点的数量就降下来了,于是再用网络流求解就OK了。不过即便这么费事地优化了原图,和裸着用匈牙利算法写的代码相比效率也才快100ms左右,可见匈牙利算法还是很高效的。
至于那些100ms左右的程序应该是用了读入优化以及各种其他的优化吧,这个题用读入优化效果还是很明显的。
View Code // 网络流算法
之所以写这个题目是因为偶然看到了用匈牙利算法做二分图多重匹配这么个东西,看完之后发现只是在原有求二分图最大匹配的基础上对匈牙利算法做了一点拓展,代码还是很好写的,而且效率也很不错。
如果用网络流做的话,这么多点是承受不了的,但由于M很小,可以将N个点适应各个星球的M个数字综合起来看成一个二进制数,这样就能把N分成2^M类,这样图的顶点的数量就降下来了,于是再用网络流求解就OK了。不过即便这么费事地优化了原图,和裸着用匈牙利算法写的代码相比效率也才快100ms左右,可见匈牙利算法还是很高效的。
至于那些100ms左右的程序应该是用了读入优化以及各种其他的优化吧,这个题用读入优化效果还是很明显的。
View Code // 网络流算法
#include<stdio.h> #include<string.h> #include<algorithm> #define MAXD 1044 #define MAXM 22558 #define INF 0x3f3f3f3f int N, M, first[MAXD], e, next[MAXM], v[MAXM], flow[MAXM]; int S, T, d[MAXD], q[MAXD], work[MAXD], num[MAXD]; void add(int x, int y, int z) { v[e] = y, flow[e] = z; next[e] = first[x], first[x] = e ++; } void init() { int i, j, k, x, cnt; memset(num, 0, sizeof(num)); cnt = 0; for(i = 1; i <= N; i ++) { x = 0; for(j = 0; j < M; j ++) scanf("%d", &k), x = x << 1 | k; if(!num[x]) ++ cnt; ++ num[x]; } S = 0, T = cnt + M + 1; memset(first, -1, sizeof(first[0]) * (T + 1)), e = 0; for(i = M; i >= 1; i --) { scanf("%d", &x); add(i, T, x), add(T, i, 0); } cnt = 0; for(i = 0; i < 1024; i ++) if(num[i]) { ++ cnt; add(S, M + cnt, num[i]), add(M + cnt, S, 0); for(j = 0; j < M; j ++) if(1 << j & i) add(M + cnt, j + 1, num[i]), add(j + 1, M + cnt, 0); } } int bfs() { int i, j, rear = 0; memset(d, -1, sizeof(d[0]) * (T + 1)); d[S] = 0, q[rear ++] = S; for(i = 0; i < rear; i ++) for(j = first[q[i]]; j != -1; j = next[j]) if(flow[j] && d[v[j]] == -1) { d[v[j]] = d[q[i]] + 1, q[rear ++] = v[j]; if(v[j] == T) return 1; } return 0; } int dfs(int cur, int a) { if(cur == T) return a; for(int &i = work[cur]; i != -1; i = next[i]) if(flow[i] && d[v[i]] == d[cur] + 1) if(int t = dfs(v[i], std::min(a, flow[i]))) { flow[i] -= t, flow[i ^ 1] += t; return t; } return 0; } int dinic() { int ans = 0, t; while(bfs()) { memcpy(work, first, sizeof(first[0]) * (T + 1)); while(t = dfs(S, INF)) ans += t; } return ans; } void solve() { printf("%s\n", dinic() == N ? "YES" : "NO"); } int main() { while(scanf("%d%d", &N, &M) == 2) { init(); solve(); } return 0; }
相关文章推荐
- 【网络流】 HDU 3605 Escape 缩图
- HDU 3605 Escape 建图状态压缩
- HDU 3605 Escape 最大流
- HDU 3605 Escape【最大流】
- HDU 3605 Escape
- HDU 3605 Escape (二分图的多重匹配)
- HDU 3605 Escape(最大流)
- HDU 3605 Escape 网络流最大流 合并点
- HDU 3605 Escape(最大流+状态压缩)
- HDU 3605 Escape(多重匹配||状压缩点的网络流)
- hdu 3605 Escape【图论-网络流-最大流-状态压缩】
- HDU 3605 Escape
- 二分图匹配 ( 多重匹配&&Hungarian)——Escape ( HDU 3605 )
- HDU_3605_Escape(多重匹配 / 网络流)
- HDU-3605 Escape(状态压缩+最大流求多重匹配、改版匈牙利算法)
- HDU 3605 Escape(多重KM,4级)
- HDU 3605 Escape(网络流 + 状压简化 )
- HDU 3605 Escape (缩点 + 最大流)
- hdu 3605 Escape【二分图多重匹配】
- hdu 3605 Escape