#include<bits/stdc++.h> #define ff first #define ss second #define mp make_pair #define pb push_back using namespace std; typedef pair<int,int> pp; //skopiowane z https://cp-algorithms.com/graph/dinic.html struct FlowEdge { int v, u; long long cap, flow = 0; FlowEdge(int v, int u, long long cap) : v(v), u(u), cap(cap) {} }; struct Dinic { const long long flow_inf = 1e18; vector<FlowEdge> edges; vector<vector<int>> adj; int n, m = 0; int s, t; vector<int> level, ptr; queue<int> q; Dinic(int n, int s, int t) : n(n), s(s), t(t) { adj.resize(n); level.resize(n); ptr.resize(n); } void add_edge(int v, int u, long long cap) { edges.emplace_back(v, u, cap); edges.emplace_back(u, v, 0); adj[v].push_back(m); adj[u].push_back(m + 1); m += 2; } bool bfs() { while (!q.empty()) { int v = q.front(); q.pop(); for (int id : adj[v]) { if (edges[id].cap - edges[id].flow < 1) continue; if (level[edges[id].u] != -1) continue; level[edges[id].u] = level[v] + 1; q.push(edges[id].u); } } return level[t] != -1; } long long dfs(int v, long long pushed) { if (pushed == 0) return 0; if (v == t) return pushed; for (int& cid = ptr[v]; cid < (int)adj[v].size(); cid++) { int id = adj[v][cid]; int u = edges[id].u; if (level[v] + 1 != level[u] || edges[id].cap - edges[id].flow < 1) continue; long long tr = dfs(u, min(pushed, edges[id].cap - edges[id].flow)); if (tr == 0) continue; edges[id].flow += tr; edges[id ^ 1].flow -= tr; return tr; } return 0; } long long flow() { long long f = 0; while (true) { fill(level.begin(), level.end(), -1); level[s] = 0; q.push(s); if (!bfs()) break; fill(ptr.begin(), ptr.end(), 0); while (long long pushed = dfs(s, flow_inf)) { f += pushed; } } return f; } }; int main() { int n,m,k; scanf("%d%d%d",&n,&m,&k); vector<pp> ve; for(int i=0;i<m;++i) { int a,b; scanf("%d%d",&a,&b); ve.pb(mp(a,b)); } int cnt[k+1]; for(int i=0;i<=k;++i) cnt[i]=0; for(int b=k+1;b<=n;++b) { for(int e=b;e<=n;++e) { Dinic tmp(2*n+3,0,2*n+2); for(int i=1;i<=n;++i) tmp.add_edge(i,i+n,1); for(int i=1;i<=k;++i) tmp.add_edge(0,i,1); for(int i=b;i<=e;++i) tmp.add_edge(i+n,2*n+2,1); for(auto I : ve) tmp.add_edge(I.ff+n,I.ss,1); // printf("%d %d %lld\n",b,e,tmp.flow()); ++cnt[tmp.flow()]; } } for(int i=0;i<=k;++i) printf("%d\n",cnt[i]); // printf("\n"); }
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 | #include<bits/stdc++.h> #define ff first #define ss second #define mp make_pair #define pb push_back using namespace std; typedef pair<int,int> pp; //skopiowane z https://cp-algorithms.com/graph/dinic.html struct FlowEdge { int v, u; long long cap, flow = 0; FlowEdge(int v, int u, long long cap) : v(v), u(u), cap(cap) {} }; struct Dinic { const long long flow_inf = 1e18; vector<FlowEdge> edges; vector<vector<int>> adj; int n, m = 0; int s, t; vector<int> level, ptr; queue<int> q; Dinic(int n, int s, int t) : n(n), s(s), t(t) { adj.resize(n); level.resize(n); ptr.resize(n); } void add_edge(int v, int u, long long cap) { edges.emplace_back(v, u, cap); edges.emplace_back(u, v, 0); adj[v].push_back(m); adj[u].push_back(m + 1); m += 2; } bool bfs() { while (!q.empty()) { int v = q.front(); q.pop(); for (int id : adj[v]) { if (edges[id].cap - edges[id].flow < 1) continue; if (level[edges[id].u] != -1) continue; level[edges[id].u] = level[v] + 1; q.push(edges[id].u); } } return level[t] != -1; } long long dfs(int v, long long pushed) { if (pushed == 0) return 0; if (v == t) return pushed; for (int& cid = ptr[v]; cid < (int)adj[v].size(); cid++) { int id = adj[v][cid]; int u = edges[id].u; if (level[v] + 1 != level[u] || edges[id].cap - edges[id].flow < 1) continue; long long tr = dfs(u, min(pushed, edges[id].cap - edges[id].flow)); if (tr == 0) continue; edges[id].flow += tr; edges[id ^ 1].flow -= tr; return tr; } return 0; } long long flow() { long long f = 0; while (true) { fill(level.begin(), level.end(), -1); level[s] = 0; q.push(s); if (!bfs()) break; fill(ptr.begin(), ptr.end(), 0); while (long long pushed = dfs(s, flow_inf)) { f += pushed; } } return f; } }; int main() { int n,m,k; scanf("%d%d%d",&n,&m,&k); vector<pp> ve; for(int i=0;i<m;++i) { int a,b; scanf("%d%d",&a,&b); ve.pb(mp(a,b)); } int cnt[k+1]; for(int i=0;i<=k;++i) cnt[i]=0; for(int b=k+1;b<=n;++b) { for(int e=b;e<=n;++e) { Dinic tmp(2*n+3,0,2*n+2); for(int i=1;i<=n;++i) tmp.add_edge(i,i+n,1); for(int i=1;i<=k;++i) tmp.add_edge(0,i,1); for(int i=b;i<=e;++i) tmp.add_edge(i+n,2*n+2,1); for(auto I : ve) tmp.add_edge(I.ff+n,I.ss,1); // printf("%d %d %lld\n",b,e,tmp.flow()); ++cnt[tmp.flow()]; } } for(int i=0;i<=k;++i) printf("%d\n",cnt[i]); // printf("\n"); } |