#include #include #include #include using namespace std; typedef long long lint; typedef vectorvi; typedef pairpii; #define rep(i,n)for(int i=0;i<(int)(n);++i) // https://github.com/koba-e964/contest/blob/master/comm/dinic.cpp /** * Dinic's algorithm for maximum flow problem. * Header requirement: vector, queue * Verified by: ABC010-D(http://abc010.contest.atcoder.jp/submissions/602810) * ARC031-D(http://arc031.contest.atcoder.jp/submissions/1050071) * POJ 3155(http://poj.org/problem?id=3155) */ template class Dinic { private: struct edge { int to; T cap; int rev; // rev is the position of reverse edge in graph[to] }; std::vector > graph; std::vector level; std::vector iter; /* Perform bfs and calculate distance from s */ void bfs(int s) { level.assign(level.size(), -1); std::queue que; level[s] = 0; que.push(s); while (! que.empty()) { int v = que.front(); que.pop(); for (int i = 0; i < graph[v].size(); ++i) { edge &e = graph[v][i]; if (e.cap > 0 && level[e.to] == -1) { level[e.to] = level[v] + 1; que.push(e.to); } } } } /* search augment path by dfs. if f == -1, f is treated as infinity. */ T dfs(int v, int t, T f) { if (v == t) { return f; } for (int &i = iter[v]; i < graph[v].size(); ++i) { edge &e = graph[v][i]; if (e.cap > 0 && level[v] < level[e.to]) { T newf = f == -1 ? e.cap : std::min(f, e.cap); T d = dfs(e.to, t, newf); if (d > 0) { e.cap -= d; graph[e.to][e.rev].cap += d; return d; } } } return 0; } public: /* v is the number of vertices (labeled from 0 .. v-1) */ Dinic(int v) : graph(v), level(v, -1), iter(v, 0) {} void add_edge(int from, int to, T cap) { graph[from].push_back((edge) {to, cap, graph[to].size()}); graph[to].push_back((edge) {from, 0, graph[from].size() - 1}); } T max_flow(int s, int t) { T flow = 0; while (1) { bfs(s); if (level[t] < 0) { return flow; } iter.assign(iter.size(), 0); T f; while ((f = dfs(s, t, -1)) > 0) { flow += f; } } } std::pair > max_flow_cut(int s, int t) { T flow = 0; while (1) { bfs(s); if (level[t] < 0) { std::vector ret; for (int i = 0; i < graph.size(); ++i) { if (level[i] < 0) { ret.push_back(i); } } return std::pair >(flow, ret); } iter.assign(iter.size(), 0); T f; while ((f = dfs(s, t, -1)) > 0) { flow += f; } } } }; // https://github.com/koba-e964/contest/blob/master/comm/MinCostFlow.cpp /* * Requirement of headers: vector, queue * Verified by: POJ 2135 (http://poj.org/problem?id=2135) */ class MinCostFlow { private: struct edge { int to, cap, cost, rev; // rev is the position of reverse edge in graph[to] }; typedef std::pair P; int v; // the number of vertices std::vector > graph; std::vector h; // potential std::vector dist; // minimum distance std::vector prevv, preve; // previous vertex and edge public: /* Initializes this solver. v is the number of vertices. */ MinCostFlow(int v) : v(v), graph(v), h(v), dist(v), prevv(v), preve(v) {} /* Initializes this solver with a existing instance. Only graph is copied. */ MinCostFlow(const MinCostFlow &ano) : v(ano.v), graph(), h(ano.v), dist(ano.v), prevv(ano.v), preve(ano.v) { for (int i = 0; i < ano.v; ++i) { std::vector tt; for (int j = 0; j < ano.graph[i].size(); ++j) { tt.push_back(ano.graph[i][j]); } graph.push_back(tt); } } /* Adds an edge. */ void add_edge(int from, int to, int cap, int cost) { graph[from].push_back((edge) {to, cap, cost, graph[to].size()}); graph[to].push_back((edge) {from, 0, -cost, graph[from].size() - 1}); } /* Calcucates the minimum cost flow whose source is s, sink is t, and flow is f. */ int min_cost_flow(int s, int t, int f) { const int inf = 0x3fffffff; int res = 0; std::fill(h.begin(), h.end(), 0); while (f > 0) { std::priority_queue, std::greater

> que; std::fill(dist.begin(), dist.end(), inf); dist[s] = 0; que.push(P(0, s)); while (! que.empty()) { P p(que.top()); que.pop(); int v = p.second; if (dist[v] < p.first) { continue; } for (int i = 0; i < graph[v].size(); ++i) { edge &e = graph[v][i]; if (e.cap > 0 && dist[e.to] > dist[v] + e.cost + h[v] - h[e.to]) { dist[e.to] = dist[v] + e.cost + h[v] - h[e.to]; prevv[e.to] = v; preve[e.to] = i; que.push(P(dist[e.to], e.to)); } } } if (dist[t] == inf) { return -1; // Cannot add flow anymore } for (int i = 0; i < v; ++i) { h[i] += dist[i]; } // Add flow fully int d = f; for (int i = t; i != s; i = prevv[i]) { d = std::min(d, graph[prevv[i]][preve[i]].cap); } f -= d; res += d * h[t]; for (int i = t; i != s; i = prevv[i]) { edge &e = graph[prevv[i]][preve[i]]; e.cap -= d; graph[i][e.rev].cap += d; } } // while (f > 0) return res; } }; int main(){ int n,a,c; cin>>n>>a; vi b(a); rep(i,a)cin>>b[i]; rep(i,a)b[i]--; cin>>c; vi d(c); rep(i,c)cin>>d[i]; rep(i,c)d[i]--; sort(b.rbegin(),b.rend()); sort(d.begin(),d.end()); vi tb(n),td(n); for(int i=0;itd[j]?0:1); } } } cout<