UVALive 8275 M - Annual Congress of MUD 2017中国台湾花莲区域赛

网络流很显然，问题是直接暴力搞边数会爆炸
注意到本质不同的人（起终点不同）只有不到四百种，以此建点

#include<stdio.h>
#include<algorithm>
#include<cstring>
#define cint const int &
#define inf 2333333
#define N 505
using namespace std;

int tot,s[N],cur[N],tm[N],ts,lev[N],Q[N],S,T,n,c,ss[401],tt[21];

struct edge{int v,c,n;}e[N*N];
inline void push(cint u,cint v,cint c)
{
    e[++tot]=(edge){v,c,s[u]};s[u]=tot;
    e[++tot]=(edge){u,0,s[v]};s[v]=tot;
}

inline bool bfs()
{
    memcpy(cur,s,sizeof(cur));
    tm[Q[1]=S]=++ts;
    for (int l=1,r=1;l<=r;l++) for (int i=s[Q[l]];i;i=e[i].n) if (e[i].c && tm[e[i].v]<ts)
        tm[Q[++r]=e[i].v]=ts,lev[Q[r]]=lev[Q[l]]+1;
    return tm[T]==ts;
}

int dfs(cint u,int c)
{
    if (u==T) return c;
    int g=0,tmp;
    for (int i=cur[u];i && c;cur[u]=i,i=e[i].n) if (tm[e[i].v]==ts && lev[e[i].v]==lev[u]+1 && e[i].c && (tmp=dfs(e[i].v,min(c,e[i].c))))
        g+=tmp,c-=tmp,e[i].c-=tmp,e[i^1].c+=tmp;
    return g;
}

inline void prepare()
{
	tot=1;
	memset(s,0,sizeof(s));
	S=421;T=422;
	for (int i=1;i<=20;i++) tt[i]=tot+1,push(400+i,T,0);
	for (int i=1;i<=400;i++) ss[i]=tot+1,push(S,i,0);
	for (int i=1;i<=20;i++) for (int j=i;j<=20;j++) for (int k=i;k<=j;k++) push((i-1)*20+j,400+k,inf);
}

int main()
{
	while(scanf("%d%d",&n,&c)==2)
	{
		prepare();
		bool flag=0;
		for (int i=1,l,r;i<=n;i++) 
		{
			scanf("%d%d",&l,&r);
			e[ss[(l-1)*20+r]].c++;
			if (bfs()) dfs(S,inf);
			else
			{
				if (flag) printf(" ");
				printf("%d",i);
				for (int i=1;i<=20;i++) ++e[tt[i]].c;
				bfs();dfs(S,inf);
				flag=1;
			}
		}
		printf("\n");
	}
}