#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; } |