[BZOJ1779][Usaco2010 Hol]Cowwar 奶牛战争（网络流）

题目：

我是超链接

题解：

ZZ读题障碍患者：关于“不能先攻击再移动”这句话的理解大概是只能选择攻击一次/移动一次/移动+攻击一次
好久之前学长的胡策题（m=n-1），竟然良心的不用输出方案
网络流
建立超级源超级汇，建图见图

各边容量都是1，那么很奇怪呀，为什么n*2+i这一列不能省略呢？
其实我这个图就是栗子，显然这个图的意思是本来两头J都不能攻击到T，但是通过E的中转可以攻击到二者，但是E只能做一头牛的中转啊，然后那头牛就既不能攻击也不能动了，就卡在E这里，所以只能攻击到一头T
但如果我们省略掉那一列，E就只作为中转而并不能限制流量，会输出可以攻击两头T

代码：

#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#define INF 1e9
using namespace std;
const int N=4005;
int tot=-1,point[N],dis[N],cur[N],nxt[N*4+5],v[N*4+5],remind[N*4+5],lian[1005][1005];
void addline(int x,int y,int cap)
{
    ++tot; nxt[tot]=point[x]; point[x]=tot; v[tot]=y; remind[tot]=cap;
    ++tot; nxt[tot]=point[y]; point[y]=tot; v[tot]=x; remind[tot]=0;
}
int dfs(int now,int t,int limit)
{
    if (now==t || !limit) return limit;
    int flow=0,f;
    for (int i=cur[now];i!=-1;i=nxt[i])
    {
        cur[now]=i;
        if (dis[v[i]]==dis[now]+1 && (f=dfs(v[i],t,min(limit,remind[i]))))
        {
            flow+=f;
            limit-=f;
            remind[i]-=f;
            remind[i^1]+=f;
            if (!limit) break;
        }
    }
    return flow;
}
bool bfs(int s,int t)
{
    queue<int>q;
    for (int i=s;i<=t;i++) cur[i]=point[i];
    q.push(s);
    memset(dis,0x7f,sizeof(dis));
    dis[s]=0;
    while (!q.empty())
    {
        int x=q.front(); q.pop();
        for (int i=point[x];i!=-1;i=nxt[i])
          if (dis[v[i]]>INF && remind[i])
          {q.push(v[i]); dis[v[i]]=dis[x]+1;}
    }
    return dis[t]<INF;
}
int dinic(int s,int t)
{
    int ans=0;
    while (bfs(s,t)) ans+=dfs(s,t,INF);
    return ans;
}
int main()
{
    int n,m;char st[N];
    memset(nxt,-1,sizeof(nxt));
    memset(point,-1,sizeof(point));
    scanf("%d%d",&n,&m);
    scanf("%s",st+1);
    for (int i=1;i<=m;i++)
    {
        int x,y;scanf("%d%d",&x,&y);
        lian[x][++lian[x][0]]=y;
        lian[y][++lian[y][0]]=x;
    }
    for (int i=1;i<=n;i++)
    {
        if (st[i]=='J') 
        {
            addline(0,i,1);
            addline(i,i+n,1);
            for (int j=1;j<=lian[i][0];j++)
              if (st[lian[i][j]]!='T')
                addline(i,lian[i][j]+n,1);
        }
        if (st[i]=='T') addline(i+3*n,4*n+1,1);                                                 
    }     
    for (int i=1;i<=n;i++)
      if (st[i]!='T')
      {
        addline(i+n,i+2*n,1);
        for (int j=1;j<=lian[i][0];j++)
          if (st[lian[i][j]]=='T')
            addline(i+2*n,lian[i][j]+3*n,1);
      }
    printf("%d",dinic(0,4*n+1));
}