BZOJ5251 八省联考2018劈配（网络流）

　　劈配，匹配，网络流。那么考虑怎么跑网络流。

　　先看第一问。首先套路的建出超源超汇。不用想也知道导师向汇连容量为战队人数上限的边。特别地，给出局也建一个点，向汇连容量inf的边（似乎没有必要）。对于一个新学员，假设我们已经知道了之前的学员的最优选择，可以把之前的每名学员和可以选择的导师连边，并由源向学员连容量为1的边。然后对于该名新学员，先只连最优选择的边，如果此时跑出的最大流不等于学员数，则表明这名学员无法选择最优，那么删掉最优边（此时这些边里一定没有流量，可以通过容量改为0实现）并连上次优边，次优边还是不行的话继续，一直这样下去直到最大流等于学员数。第一问就做完了。至于复杂度，O(能过)。

　　然后是第二问。对于每个人可以二分答案，然后跑最大流看其是否满足。似乎需要访问网络的历史状态？可持久化网络流！这玩意似乎没法可持久化啊……那就暴力记下来呗。由于有C的限制，这里面的边不会很多。于是就做完了，O(能过)。还有一种做法是先连上该学员可以选择的边，然后从第一名开始依次把最优选择加进去跑，直到无法满足，可能会快不少。

　

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
#define N 210
#define S 0
#define T 404
int test,c,n,m,b[N],s[N],p[N<<],cnt[N],t;
int cur[N<<],d[N<<],q[N<<],ans,cho[N];
vector<int> a[N][N];
struct data{int to,nxt,cap,flow;
}edge[N*N<<],hisedge[N][];
void addedge(int x,int y,int z)
{
    t++;edge[t].to=y,edge[t].nxt=p[x],edge[t].cap=z,edge[t].flow=,p[x]=t;
    t++;edge[t].to=x,edge[t].nxt=p[y],edge[t].cap=,edge[t].flow=,p[y]=t;
}
void addedge(int x,int y,int cap,int flow)
{
    t++;edge[t].to=y,edge[t].nxt=p[x],edge[t].cap=cap,edge[t].flow=flow,p[x]=t;
    t++;edge[t].to=x,edge[t].nxt=p[y],edge[t].cap=,edge[t].flow=-flow,p[y]=t;
}
void init()
{
    for (int i=;i<=m;i++) b[i]=read();
    for (int i=;i<=n;i++)
    {
        for (int j=;j<=m+;j++)
        a[i][j].clear();
        for (int j=;j<=m;j++)
        {
            int x=read();
            if (x) a[i][x].push_back(j);
        }
        a[i][m+].push_back(m+);
    }
    for (int i=;i<=n;i++) s[i]=read();
}
bool bfs()
{
    memset(d,,sizeof(d));d[S]=;
    int head=,tail=;q[]=S;
    do
    {
        int x=q[++head];
        for (int i=p[x];~i;i=edge[i].nxt)
        if (d[edge[i].to]==-&&edge[i].flow<edge[i].cap)
        {
            d[edge[i].to]=d[x]+;
            q[++tail]=edge[i].to;
        }
    }while (head<tail);
    return ~d[T];
}
int work(int k,int f)
{
    if (k==T) return f;
    int used=;
    for (int i=cur[k];~i;i=edge[i].nxt)
    if (d[k]+==d[edge[i].to])
    {
        int w=work(edge[i].to,min(f-used,edge[i].cap-edge[i].flow));
        edge[i].flow+=w,edge[i^].flow-=w;
        if (edge[i].flow<edge[i].cap) cur[k]=i;
        used+=w;if (used==f) return f;
    }
    if (used==) d[k]=-;
    return used;
}
void dinic()
{
    while (bfs())
    {
        memcpy(cur,p,sizeof(p));
        ans+=work(S,N);
    }
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj5251.in","r",stdin);
    freopen("bzoj5251.out","w",stdout);
    const char LL[]="%I64d";
#else
    const char LL[]="%lld";
#endif
    test=read(),c=read();
    while (test--)
    {
        n=read(),m=read();
        t=-;
        memset(p,,sizeof(p));
        memset(cnt,,sizeof(cnt));
        init();
        for (int i=;i<=m;i++) addedge(n+i,T,b[i]);
        addedge(n+m+,T,N);
        for (int j=;j<=t;j+=)
        if (edge[j].cap) hisedge[][++cnt[]]=edge[j],hisedge[][cnt[]].nxt=edge[j^].to;
        ans=;
        for (int i=;i<=n;i++)
        {
            addedge(S,i,);
            for (int j=;j<=m+;j++)
            {
                int tmp=t;
                for (int k=;k<a[i][j].size();k++)
                addedge(i,n+a[i][j][k],);
                dinic();
                if (ans==i) {cho[i]=j;break;}
                for (int k=tmp+;k<=t;k+=) edge[k].cap=;
            }
            for (int j=;j<=t;j+=)
            if (edge[j].cap) hisedge[i][++cnt[i]]=edge[j],hisedge[i][cnt[i]].nxt=edge[j^].to;
        }
        /*for (int i=0;i<=n;i++)
        {
            for (int j=1;j<=cnt[i];j++)
            cout<<hisedge[i][j].nxt<<' '<<hisedge[i][j].to<<' '<<hisedge[i][j].cap<<' '<<hisedge[i][j].flow<<endl;
            cout<<endl;
        }*/
        for (int i=;i<=n;i++) printf("%d ",cho[i]);
        cout<<endl;
        for (int i=;i<=n;i++)
        {
            int l=,r=i-,add=-;
            while (l<=r)
            {
                int mid=l+r>>;
                memset(p,,sizeof(p));
                t=-;
                addedge(S,i,);
                for (int j=;j<=s[i];j++)
                    for (int k=;k<a[i][j].size();k++)
                    addedge(i,n+a[i][j][k],);
                for (int j=;j<=cnt[mid];j++)
                addedge(hisedge[mid][j].nxt,hisedge[mid][j].to,hisedge[mid][j].cap,hisedge[mid][j].flow);
                ans=;
                dinic();
                if (!ans) r=mid-;
                else l=mid+,add=mid;
            }
            printf("%d ",i-add-);
        }
        cout<<endl;
    }
    return ;
}