#include <bits/stdc++.h>
using namespace std;
const int MX=200200,ME=2200200;
int n,m,k,num,i,x,y,cnt,nle,le,ri,ed,it,fi,fr,q[MX],u[MX],pre[MX],se[2][MX],f[2][ME],t[ME];
vector<int> g[MX];
long long ans[55];
long long sm(long long x) {
return (x*(x+1LL))/2;
}
void addedge(int fr, int tt) {
t[ed]=fr; f[0][ed]=0; g[tt].push_back(ed++);
t[ed]=tt; f[0][ed]=1; g[fr].push_back(ed++);
}
int bfs(int* f) {
int fr=1;
q[0]=0;
u[0]=++it;
int mn=num+1;
for (int fi=0; fi<fr; fi++) {
int i=q[fi];
if (i>=nle) mn=min(mn,i);
for (int ed: g[i]) if (f[ed]>0) {
int nxt=t[ed];
if (u[nxt]!=it) {
u[nxt]=it;
pre[nxt]=ed;
q[fr++]=nxt;
}
}
}
return mn;
}
void addpath(int cur, int i) {
while (i) {
int ed=pre[i];
if (ed&1) se[cur][i]=ed;
--f[cur][ed];
ed^=1;
++f[cur][ed];
i=t[ed];
}
}
void rmpath(int cur, int i) {
while (i) {
int ed=se[cur][i];
se[cur][i]=-1;
++f[cur][ed];
ed^=1;
--f[cur][ed];
i=t[ed];
}
}
int main() {
scanf("%d%d%d",&n,&m,&k);
num=n+n-k;
for (i=1; i<=k; i++) addedge(0,i);
for (i=k+1; i<=n; i++) addedge(i,n+i-k);
for (i=0; i<m; i++) {
scanf("%d%d",&x,&y);
addedge((x<=k)?x:(x+n-k),y);
}
for (i=0; i<=num; i++) se[0][i]=-1;
for (cnt=0; cnt<k; cnt++) {
le=k+1;
nle=n+1;
int cur=(cnt&1);
int oth=(cur^1);
int mn=bfs(f[cur]);
if (mn>num) break;
ri=mn-n+k;
addpath(cur,mn);
if (cnt<k) {
for (i=0; i<ed; i++) f[oth][i]=f[cur][i];
for (i=0; i<=num; i++) se[oth][i]=se[cur][i];
}
ans[cnt+1]+=n-ri+1;
for (++le, ++nle; le<=n; ++le, ++nle) {
if (se[cur][n+le-k-1]>=0) {
rmpath(cur,n+le-k-1);
mn=bfs(f[cur]);
if (mn>num) break;
addpath(cur,mn);
ri=max(ri,mn-n+k);
}
ans[cnt+1]+=n-ri+1;
}
}
ans[0]=sm(n-k);
for (i=0; i<=k; i++) printf("%lld\n",ans[i]-ans[i+1]);
return 0;
}
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 | #include <bits/stdc++.h> using namespace std; const int MX=200200,ME=2200200; int n,m,k,num,i,x,y,cnt,nle,le,ri,ed,it,fi,fr,q[MX],u[MX],pre[MX],se[2][MX],f[2][ME],t[ME]; vector<int> g[MX]; long long ans[55]; long long sm(long long x) { return (x*(x+1LL))/2; } void addedge(int fr, int tt) { t[ed]=fr; f[0][ed]=0; g[tt].push_back(ed++); t[ed]=tt; f[0][ed]=1; g[fr].push_back(ed++); } int bfs(int* f) { int fr=1; q[0]=0; u[0]=++it; int mn=num+1; for (int fi=0; fi<fr; fi++) { int i=q[fi]; if (i>=nle) mn=min(mn,i); for (int ed: g[i]) if (f[ed]>0) { int nxt=t[ed]; if (u[nxt]!=it) { u[nxt]=it; pre[nxt]=ed; q[fr++]=nxt; } } } return mn; } void addpath(int cur, int i) { while (i) { int ed=pre[i]; if (ed&1) se[cur][i]=ed; --f[cur][ed]; ed^=1; ++f[cur][ed]; i=t[ed]; } } void rmpath(int cur, int i) { while (i) { int ed=se[cur][i]; se[cur][i]=-1; ++f[cur][ed]; ed^=1; --f[cur][ed]; i=t[ed]; } } int main() { scanf("%d%d%d",&n,&m,&k); num=n+n-k; for (i=1; i<=k; i++) addedge(0,i); for (i=k+1; i<=n; i++) addedge(i,n+i-k); for (i=0; i<m; i++) { scanf("%d%d",&x,&y); addedge((x<=k)?x:(x+n-k),y); } for (i=0; i<=num; i++) se[0][i]=-1; for (cnt=0; cnt<k; cnt++) { le=k+1; nle=n+1; int cur=(cnt&1); int oth=(cur^1); int mn=bfs(f[cur]); if (mn>num) break; ri=mn-n+k; addpath(cur,mn); if (cnt<k) { for (i=0; i<ed; i++) f[oth][i]=f[cur][i]; for (i=0; i<=num; i++) se[oth][i]=se[cur][i]; } ans[cnt+1]+=n-ri+1; for (++le, ++nle; le<=n; ++le, ++nle) { if (se[cur][n+le-k-1]>=0) { rmpath(cur,n+le-k-1); mn=bfs(f[cur]); if (mn>num) break; addpath(cur,mn); ri=max(ri,mn-n+k); } ans[cnt+1]+=n-ri+1; } } ans[0]=sm(n-k); for (i=0; i<=k; i++) printf("%lld\n",ans[i]-ans[i+1]); return 0; } |
English