BZOJ 3171 [Tjoi2013]循环格

Description
一个循环格就是一个矩阵，其中所有元素为箭头，指向相邻四个格子。每个元素有一个坐标（行，列），其中左上角元素坐标为（0,0）。给定一个起始位置（r，c），你可以沿着箭头防线在格子间行走。即如果（r,c）是一个左箭头，那么走到（r，c-1）;如果是右箭头那么走到（r，c+1）；如果是上箭头那么走到（r-1，c）；如果是下箭头那么走到（r+1，c）；每一行和每一列都是循环的，即如果走出边界，你会出现在另一侧。
一个完美的循环格是这样定义的：对于任意一个起始位置，你都可以i沿着箭头最终回到起始位置。如果一个循环格不满足完美，你可以随意修改任意一个元素的箭头直到完美。给定一个循环格，你需要计算最少需要修改多少个元素使其完美。

---

【题目分析】
题目大意就是要求每一个格子都有唯一的后继，互不相同，这样子的话就可以保证从每个格子出发都会有一个环。而且发现，两个格子连通时候的与不联通的情况下确定后继的代价是完全不同的，那么这就是一个二分图的最小权匹配，就可以转化成最小费用最大流的问题了。

---

【代码】

#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;

#define maxn 2500
#define maxp 20
#define inf (0x3f3f3f3f)
#define M(a) memset(a,-1,sizeof a)
#define M0(a) memset(a,0,sizeof a)
#define MB(a) memset(a,0x3f,sizeof a)

int h[maxn],to[maxn],ne[maxn],fl[maxn],cost[maxn],en=0;
int dis[maxn],a[maxp][maxp],minn[maxn],with[maxn],inq[maxn];
int mov[4][2]={{0,1},{0,-1},{-1,0},{1,0}};
int hasher[maxp][maxp][2],S,T,vis[maxn],n,m;
char s[maxp];

int read()
{
    int x=0,f=1; char ch=getchar();
    while (ch<'0'||ch>'9') { f=-1; ch=getchar(); }
    while (ch>='0'&&ch<='9') { x=x*10+ch-'0'; ch=getchar(); }
    return x*f;
}

void add(int a,int b,int r,int c)
{
//  printf("add %d to %d can %d is %d\n",a,b,r,c);
//  printf("add %d to %d can %d is %d\n",a,b,r,c);
    to[en]=b; ne[en]=h[a]; fl[en]=r; cost[en]=c; h[a]=en++;
    to[en]=a; ne[en]=h[b]; fl[en]=0; cost[en]=-c;h[b]=en++;
}

bool tell()
{
    queue <int> q;
    while (!q.empty()) q.pop();
    memset(dis,0x3f,sizeof dis);
    memset(minn,0x3f,sizeof minn);
    memset(inq,0,sizeof inq);
    memset(with,0,sizeof with);
    q.push(S);inq[S]=1;dis[S]=0;
    while (!q.empty())
    {
        int x=q.front(); q.pop(); inq[x]=0;
        for (int i=h[x];i>=0;i=ne[i])
        {
            if (fl[i]>0&&dis[to[i]]>dis[x]+cost[i])
            {
                dis[to[i]]=dis[x]+cost[i];
                minn[to[i]]=min(fl[i],minn[x]);
                with[to[i]]=i;
                if (!inq[to[i]]) q.push(to[i]),inq[to[i]]=1;
            }
        }
    }
//  for (int i=0;i<=n*m+1;++i) printf("%d ",dis[i]); printf("\n");

    if (dis[T]==inf) return false;
    else return true;
}

int zeng()
{
    for (int i=T;i!=S;i=to[with[i]^1])
    {
        fl[with[i]]-=minn[T];
        fl[with[i]^1]+=minn[T];
    }
    return dis[T]*minn[T];
}

int main()
{
    memset(h,-1,sizeof h);
    scanf("%d%d",&n,&m);
    int cnt=0;
    for (int i=0;i<n;++i)
    {
        scanf("%s",s);
        for (int j=0;j<m;++j)
        {
            if (s[j]=='D') a[i][j]=3;
            if (s[j]=='U') a[i][j]=2;
            if (s[j]=='L') a[i][j]=1;
            if (s[j]=='R') a[i][j]=0;
            hasher[i][j][0]=++cnt;
            hasher[i][j][1]=++cnt;
        }
    }
    S=0; T=++cnt;
    for (int i=0;i<n;++i)
        for (int j=0;j<m;++j)
        {
            add(S,hasher[i][j][0],1,0);
            add(hasher[i][j][1],T,1,0);
        }
    for (int i=0;i<n;++i)
        for (int j=0;j<m;++j)
        {
            for (int k=0;k<4;++k)
            {
                int tx=i+mov[k][0],ty=j+mov[k][1];
                tx=(tx+n)%n;
                ty=(ty+m)%m;
                if (k==a[i][j]) add(hasher[i][j][0],hasher[tx][ty][1],1,0);
                else add(hasher[i][j][0],hasher[tx][ty][1],1,1);
            }
        }
//  for (int i=1;i<cnt;++i) add(i*2,i*2+1,1,0);
    int ans=0;
    while (tell()) ans+=zeng();
    printf("%d\n",ans);
}