博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
解题报告 『[USACO07OPEN]Dining(网络最大流 + 拆点)』
阅读量:4609 次
发布时间:2019-06-09

本文共 2869 字,大约阅读时间需要 9 分钟。

网络流板子题 + 拆点,你谷评分一如既往的虚高。

 

代码实现如下:

#include 
using namespace std;#define rep(i, a, b) for (register int i = (a); i <= (b); i++)const int inf = 0x3f3f3f3f, maxn = 1e4 + 5;int n, f, d, S, T, tmp, ans = 0, f_sum, d_sum, num_edge = -1;int cur[maxn], dep[maxn], head[maxn];queue
q;struct node { int to, dis, nxt;}edge[maxn << 1];void origin(){memset(head, -1, sizeof(head));}int read() { int x = 0, flag = 0; char ch = ' '; while (ch != '-' && (ch < '0' || ch > '9')) ch = getchar(); if (ch == '-') { flag = 1; ch = getchar(); } while (ch >= '0' && ch <= '9') { x = (x << 1) + (x << 3) + ch - '0'; ch = getchar(); } return flag ? -x : x;}void addedge(int from, int to, int dis) { edge[++num_edge].nxt = head[from]; edge[num_edge].to = to; edge[num_edge].dis = dis; head[from] = num_edge;}int bfs(int S, int T) { memset(dep, 0, sizeof(dep)); while (!q.empty()) q.pop(); memcpy(cur, head, sizeof(cur)); dep[S] = 1; q.push(S); while (!q.empty()) { int u = q.front(); q.pop(); for (int i = head[u]; ~i; i = edge[i].nxt) { int v = edge[i].to; if (!dep[v] && edge[i].dis) { dep[v] = dep[u] + 1; q.push(v); } } } if (dep[T]) return 1; return 0;}int dfs(int u, int flow) { if (u == T || !flow) return flow; int d, used = 0; for (int i = cur[u]; ~i; i = edge[i].nxt) { cur[u] = i; int v = edge[i].to; if (dep[v] == dep[u] + 1 && (d = dfs(v, min(flow, edge[i].dis)))) { used += d; flow -= d; edge[i].dis -= d; edge[i ^ 1].dis += d; if (!flow) break; } } if (!used) dep[u] = -2; return used;}int dinic() { int ans = 0; while (bfs(S, T)) ans += dfs(S, inf); return ans;}void write(int x) { if (x < 0) { putchar('-'); x = -x; } if (x > 9) write(x / 10); putchar(x % 10 + '0');}int main() { origin(); n = read(), f = read(), d = read(); S = 0, T = f + (n << 1) + d + 1; rep(i, 1, f) { addedge(S, i, 1); addedge(i, S, 0); } rep(i, 1, n) { addedge(f + i, f + n + i, 1); addedge(f + n + i, f + i, 0); } rep(i, 1, d) { addedge(f + (n << 1) + i, T, 1); addedge(T, f + (n << 1) + i, 0); } rep(i, 1, n) { f_sum = read(), d_sum = read(); rep(j, 1, f_sum) { tmp = read(); addedge(tmp, f + i, 1); addedge(f + i, tmp, 0); } rep(j, 1, d_sum) { tmp = read(); addedge(f + n + i, f + (n << 1) + tmp, 1); addedge(f + (n << 1) + tmp, f + n + i, 0); } } ans = dinic(); write(ans); return 0;}
View Code

转载于:https://www.cnblogs.com/Kirisame-Marisa/p/10804429.html

你可能感兴趣的文章
How to Notify Command to evaluate in mvvmlight
查看>>
33. Search in Rotated Sorted Array
查看>>
461. Hamming Distance
查看>>
Python垃圾回收机制详解
查看>>
jquery 编程的最佳实践
查看>>
MeetMe
查看>>
IP报文格式及各字段意义
查看>>
(转载)rabbitmq与springboot的安装与集成
查看>>
C2. Power Transmission (Hard Edition)(线段相交)
查看>>
STM32F0使用LL库实现SHT70通讯
查看>>
Atitit. Xss 漏洞的原理and应用xss木马
查看>>
MySQL源码 数据结构array
查看>>
(文件过多时)删除目录下全部文件
查看>>
T-SQL函数总结
查看>>
python 序列:列表
查看>>
web移动端
查看>>
pythonchallenge闯关 第13题
查看>>
linux上很方便的上传下载文件工具rz和sz使用介绍
查看>>
React之特点及常见用法
查看>>
【WEB前端经验之谈】时间一年半,或沉淀、或从零开始。
查看>>