拆分 （网络流，构造）

可以构造使得答案就是下界值。调整
方法是，在加入航线的过程中，保证每个点都先填满k种不同颜色，然后清零从头开始。具体地，对于一条航线(u,v)，可以找到cx,cy是u,v的最小没有用过的颜色，假定cx<cy，那么从v开始一定有一条cx,cy,cx,cy…的交错链。将其整体取反就可以把航线染成cx了。
当一个点度数为k的时候，就可以清零（开新点）。

每次也可以用上下界网络流来完成这个过程。但是复杂度嘛…
具体地，对于每种颜色，跑一遍网络流，保证每个点流$di/k$的流量就行。（上下界分别是下整，上整）。每跑完一种颜色重构图，删掉已经染色的边。

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 5010;
int na,nb,m,k;
int f[N][310],kth[N][310],d[N],tot;
int co[N], ne[N], dd[N];

void kz(int x,int c1,int c2) {
	// printf("x=%d old:%d %d %d %d\n",x,f[x][c1],kth[x][c1],f[x][c2],kth[x][c2]);
	co[kth[x][c1]] = c2;
	co[kth[x][c2]] = c1;
	swap(f[x][c1], f[x][c2]);
	swap(kth[x][c1], kth[x][c2]);
	// printf("x=%d new:%d %d %d %d\n",x,f[x][c1],kth[x][c1],f[x][c2],kth[x][c2]);
	if (f[x][c2] == 0) return;
	kz(f[x][c2], c2, c1);
}

int ans;
int main() {
	freopen("b.in","r",stdin);
	freopen("b.out","w",stdout);
	cin>>na>>nb>>m>>k;
	for (int i = 1; i <= na + nb; i++) ne[i] = i;
	tot = na + nb;
	for (int i = 1; i <= m; i++) {
		for (int a = 1; a <= tot; a++) {
			for (int b = 1; b <= k; b++) {
				if (f[a][b] && f[f[a][b]][b] != a) {
					printf("WOCAO");
				}
			}
		}

		int u,v;
		scanf("%d %d",&u,&v);
		int uu=u,vv=v+na;
		v += na; dd[u]++,dd[v]++;
		u = ne[u], v = ne[v];

		if (d[u] == k) {
			u = ne[uu] = ++tot;
		}

		if (d[v] == k)  {
			v = ne[vv] = ++tot;
		}

		d[u]++,d[v]++;
		int cx,cy;
		for (int j = 1; j <= k; j++) if (f[u][j] == 0) {cx=j; break;}
		for (int j = 1; j <= k; j++) if (f[v][j] == 0) {cy=j; break;}

		if (cx > cy) swap(cx,cy), swap(u,v);
		kz(v, cx, cy);
		co[i] = cx;
		f[u][cx] = v, kth[u][cx] = i;
		f[v][cx] = u, kth[v][cx] = i;
	}
	for (int i = 1; i <= na + nb; i++) ans += dd[i] % k != 0;
	// cout<<ans<<endl;
	for (int i = 1; i <= m; i++) printf("%d ",co[i]); printf("\n");
}