[loj2477]劈配

考虑依次选择每一位考生，设当前选到第$i+1$位，前i个分别为$p1,p2,……pi$（注意：这里只确定了导师的志愿编号），然后枚举第$p_{i+1}$，通过网络流建图+判定，复杂度为$o(nm*f(n+m,cn))$（$f(n,m)$表示点数为n、边数为m的网络流），无法通过

考虑优化，由于网络流可以先任意流，那么直接在之前的残余网络上加上第k份志愿并判断能否联通即可，复杂度$o(cn^{2}m)$

第二个问题类似，存储下每一次跑完后的残余网络，然后从下往上枚举其排名（志愿为si及其以上的）并建图，复杂度同样为$o(cn^{2}m)$

 

#include<bits/stdc++.h>
using namespace std;
#define N 205
struct ji{
    int nex,to,len;
}edge[30*N],e[N][30*N];
vector<int>v[N][N];
int E,t,n,m,x,EE[N],head[N<<1],vis[N<<1],h[N][N<<1];
void add(int x,int y,int z){
    edge[E].nex=head[x];
    edge[E].to=y;
    edge[E].len=z;
    head[x]=E++;
}
bool dfs(int k){
    if (k==n+m+1)return 1;
    if (vis[k])return 0;
    vis[k]=1;
    for(int i=head[k];i!=-1;i=edge[i].nex)
        if ((edge[i].len)&&(dfs(edge[i].to))){
            edge[i].len--;
            edge[i^1].len++;
            return 1;
        }
    return 0;
}
int main(){
    scanf("%d%*d",&t);
    while (t--){
        scanf("%d%d",&n,&m);
        E=0;
        memset(head,-1,sizeof(head));
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)v[i][j].clear();
        for(int i=1;i<=m;i++){
            scanf("%d",&x);
            add(i+n,n+m+1,x);
            add(n+m+1,i+n,0);
        }
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++){
                scanf("%d",&x);
                if (x)v[i][x].push_back(j);
            }
        EE[0]=E;
        memcpy(h[0],head,sizeof(head));
        memcpy(e[0],edge,sizeof(edge));
        for(int i=1;i<=n;i++){
            bool flag=0;
            EE[i]=E;
            memcpy(h[i],head,sizeof(head));
            memcpy(e[i],edge,sizeof(edge));
            for(int j=1;j<=m;j++){
                E=EE[i];
                memset(vis,0,sizeof(vis)); 
                memcpy(head,h[i],sizeof(head));
                memcpy(edge,e[i],sizeof(edge));
                for(int k=0;k<v[i][j].size();k++){
                    add(i,v[i][j][k]+n,1);
                    add(v[i][j][k]+n,i,0);
                }
                if (dfs(i)){
                    flag=1;
                    printf("%d ",j);
                    break;
                }
            }
            if (!flag)printf("%d ",m+1);
            EE[i]=E;
            memcpy(h[i],head,sizeof(head));
            memcpy(e[i],edge,sizeof(edge));
        }
        printf("
");
        for(int i=1;i<=n;i++){
            scanf("%d",&x);
            bool flag=0;
            for(int j=i;j;j--){
                E=EE[j-1];
                memset(vis,0,sizeof(vis)); 
                memcpy(head,h[j-1],sizeof(head));
                memcpy(edge,e[j-1],sizeof(edge));
                for(int k=1;k<=x;k++)
                    for(int l=0;l<v[i][k].size();l++){
                        add(i,v[i][k][l]+n,1);
                        add(v[i][k][l]+n,i,0);
                    }
                if (dfs(i)){
                    flag=1;
                    printf("%d ",i-j);
                    break;
                }
            }
            if (!flag)printf("%d ",i);
        }
        printf("
");
    }
}