#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"); } |
English