Gym - 247731E :room（最小费用流裸题）

题意：有N个宿舍(N<200)，给出第一年每个宿舍有哪4个同学。现在给出N个4元组y[][4]，表示这4个人想住一起，问最少多少人需要换宿舍。

思路：费用流，每个4元组y[]到每个宿舍连边，流量是1，费用是4-same。

#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;i++)
const int inf=<<;
const int maxn=;
using namespace std;
int To[maxn*],Laxt[maxn],Next[maxn*],cap[maxn*],cost[maxn*];
int S,T,cnt=,dis[maxn],ans;
bool inq[maxn],vis[maxn];
deque<int>q;
void add(int u,int v,int c,int cc)
{
    Next[++cnt]=Laxt[u];Laxt[u]=cnt;
    To[cnt]=v;cap[cnt]=c;cost[cnt]=cc;
}
bool spfa()
{
    for(int i=;i<=T;i++) inq[i]=;
    for(int i=;i<=T;i++) dis[i]=inf;
    inq[T]=; dis[T]=; q.push_back(T);
    while(!q.empty())
    {
        int u=q.front(); q.pop_front();
        inq[u]=;
        for(int i=Laxt[u];i;i=Next[i])
        {
            int v=To[i];
            if(cap[i^]&&dis[v]>dis[u]-cost[i])
            {
                dis[v]=dis[u]-cost[i];
                if(!inq[u]){
                    inq[v]=;
                    if(q.empty()||dis[v]>dis[q.front()]) q.push_back(v);
                    else q.push_front(v);
                }
            }
        }
    }
    return dis[S]<inf;
}
int dfs(int u,int flow)
{
    vis[u]=;
    if(u==T||flow==) return flow;
    int tmp,delta=;
    for(int i=Laxt[u];i;i=Next[i])
    {
        int v=To[i];
        if((!vis[v])&&cap[i]&&dis[v]==dis[u]-cost[i])
        {
            tmp=dfs(v,min(cap[i],flow-delta));
            delta+=tmp; cap[i]-=tmp; cap[i^]+=tmp;
        }
    }
    return delta;
}
int x[maxn][],y[maxn][];
int main()
{
    int N;
    scanf("%d",&N);
    rep(i,,N) rep(j,,) scanf("%d",&x[i][j]);
    rep(i,,N) rep(j,,) scanf("%d",&y[i][j]);
    rep(i,,N) sort(x[i],x[i]+);
    rep(i,,N) sort(y[i],y[i]+);
    S=,T=N+N+; cnt=;
    rep(i,,N) add(S,i,,),add(i,S,,);
    rep(i,,N) add(N+i,T,,),add(T,N+i,,);
    rep(i,,N) rep(j,,N) {
        int L1=,L2=,same=;
        while(L1<=&&L2<=) {
            if(x[i][L1]==y[j][L2]) same++,L1++,L2++;
            else if(x[i][L1]<y[j][L2]) L1++;
            else L2++;
        }
        add(i,N+j,,-same); add(N+j,i,,same-);
    }
    int ans=;
    while(spfa()){
        vis[T]=;
        while(vis[T]){
            for(int i=;i<=T;i++) vis[i]=;
            ans+=dis[S]*dfs(S,inf);
        }
    }
    printf("%d\n",ans);
    return ;
}