BZOJ 1779. [Usaco2010 Hol]Cowwar 奶牛战争

传送门

考虑构建网络流模型

把一个流量看成一只奶牛的攻击过程，那么答案就是最大流

因为每只奶牛只能操作一波，所以构造分层图，一层相当于一步

第一层就是初始状态，从 $S$ 向所有 $J$ 奶牛连一条流量为 $1$ 的边，表示只有一只 $J$

下一层，表示奶牛走一步后的状态，每只 $J$ 向下一层走一步可以到达的点连流量为 $1$ 的边，当然奶牛可以不走，所以也要向下一层原来的位置连流量为 $1$ 的边

再下一层就可以考虑奶牛攻击了，但是因为每个位置只能站一只奶牛，所以这一层还要再拆成两层，对于点之间连流量为 $1$ 的边，表示一个位置只能站一只奶牛，然后再把所有的 $T$牛 连一条流量为 $1$ 的边到 $T$（汇点），因为一只 $T$牛 只能被打翻一次

最后答案就是最大流

总结一下就是图分成四层，每一步慢慢考虑

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
using namespace std;
typedef long long ll;
inline int read()
{
    int x=,f=; char ch=getchar();
    while(ch<''||ch>'') { if(ch=='-') f=-; ch=getchar(); }
    while(ch>=''&&ch<='') { x=(x<<)+(x<<)+(ch^); ch=getchar(); }
    return x*f;
}
const int N=1e5+,M=4e6+,INF=1e9+;
int fir[N],from[M],to[M],val[M],cntt=;
inline void add(int a,int b,int c)
{
    from[++cntt]=fir[a]; fir[a]=cntt;
    to[cntt]=b; val[cntt]=c;
    from[++cntt]=fir[b]; fir[b]=cntt;
    to[cntt]=a; val[cntt]=;
}
int dep[N],Fir[N],S,T;
queue <int> q;
bool BFS()
{
    for(int i=S;i<=T;i++) Fir[i]=fir[i],dep[i]=;
    q.push(S); dep[S]=; int x;
    while(!q.empty())
    {
        x=q.front(); q.pop();
        for(int i=fir[x];i;i=from[i])
        {
            int &v=to[i]; if(dep[v]||!val[i]) continue;
            dep[v]=dep[x]+; q.push(v);
        }
    }
    return dep[T]>;
}
int DFS(int x,int mxf)
{
    if(x==T||!mxf) return mxf;
    int fl=,res;
    for(int &i=Fir[x];i;i=from[i])
    {
        int &v=to[i]; if(dep[v]!=dep[x]+||!val[i]) continue;
        if( res=DFS(v,min(mxf,val[i])) )
        {
            mxf-=res; fl+=res;
            val[i]-=res; val[i^]+=res;
            if(!mxf) break;
        }
    }
    return fl;
}
inline int Dinic() { int res=; while(BFS()) res+=DFS(S,INF); return res; }

int n,m;
char s[N];
int main()
{
    n=read(),m=read();
    scanf("%s",s+);
    S=,T=n*+; int a,b;
    for(int i=;i<=n;i++)
    {
        if(s[i]=='J') add(S,i,),add(i,n+i,);
        if(s[i]=='T') add(n*+i,T,);
        else add(n+i,n*+i,);
    }
    for(int i=;i<=m;i++)
    {
        a=read(),b=read();
        if(s[a]!='T'&&s[b]!='T') add(a,n+b,),add(b,n+a,);
        if(s[a]!='T'&&s[b]=='T') add(n*+a,n*+b,);
        if(s[b]!='T'&&s[a]=='T') add(n*+b,n*+a,);
    }
    printf("%d",Dinic());
    return ;
}