[TJOI2013] 循环格

题目描述：

qwq.

题目分析：

考虑到一个完美循环图每个点的入度与初度均为1
那么我们就可以做 最小权二分图匹配
把每个点拆成入点与出点
然后从每个格子向四周的格子连容量为1的边，与原方向相同则花费为0，否则为1
搞一下超级源汇，做最小权二分图匹配即可.

题目链接：

Luogu 3965
BZOJ 3171

Ac 代码：

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <queue>
const int inf=0x7fffffff;
const int maxm=2e6+100;
int head[maxm],to[maxm<<1],net[maxm<<1],cap[maxm<<1],cost[maxm<<1];
int cnt=1;
int n,m;
int id[200][200],tot;
inline void addedge(int u,int v,int c1,int c2)
{
    cnt++;
    to[cnt]=v,cap[cnt]=c1,cost[cnt]=c2,net[cnt]=head[u],head[u]=cnt;
    cnt++;
    to[cnt]=u,cap[cnt]=0,cost[cnt]=-c2,net[cnt]=head[v],head[v]=cnt;
}
namespace MCMF{
    int pre[maxm],id[maxm],dis[maxm],flow[maxm];
    int mincost,maxflow;
    bool vis[maxm];
    std::queue <int> dl;
    inline bool SPFA(int s,int t)
    {
        memset(pre,-1,sizeof(pre));
        memset(dis,127/3,sizeof(dis));
        dis[s]=0,pre[s]=0,flow[s]=inf,vis[s]=1;
        dl.push(s);
        while(!dl.empty())
        {
            int now=dl.front();
            dl.pop();
            vis[now]=0;
            for(int i=head[now];i;i=net[i])
            if(cap[i]&&dis[to[i]]>dis[now]+cost[i])
            {
                dis[to[i]]=dis[now]+cost[i];
                pre[to[i]]=now;
                id[to[i]]=i;
                flow[to[i]]=std::min(flow[now],cap[i]);
                if(!vis[to[i]]) vis[to[i]]=1,dl.push(to[i]);
            }
        }
        return pre[t]!=-1;
    }
    inline void change_cap(int s,int t,int x)
    {
        int now=t;
        while(now!=s)
        {
            cap[id[now]]-=x,cap[id[now]^1]+=x;
            now=pre[now];
        }
    }
    inline int mcmf(int s,int t)
    {
        maxflow=mincost=0;
        while(SPFA(s,t))
        {
            maxflow+=flow[t];
            mincost+=flow[t]*dis[t];
            change_cap(s,t,flow[t]);
        }
        return mincost;
    }
} 
int dx[]={0,0,1,0,-1},dy[]={0,1,0,-1,0};
int main()
{
    scanf("%d%d",&n,&m);
    int s=0,t=2*n*m+10; 
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++)
      {
        id[i][j]=++tot;
        addedge(s,id[i][j],1,0);
        addedge(id[i][j]+n*m,t,1,0);
      }
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++)
     {
        char c;
        std::cin>>c;
        int tmp;
        if(c=='R') tmp=1;
        if(c=='D') tmp=2;
        if(c=='L') tmp=3;
        if(c=='U') tmp=4;
        for(int k=1;k<=4;k++)
        {
            int x=i+dx[k],y=j+dy[k];
            if(x==0) x=n;
            if(x==n+1) x=1;
            if(y==0) y=m;
            if(y==m+1) y=1;
            addedge(id[i][j],id[x][y]+n*m,1,(tmp!=k));
        }
     }  
    printf("%d\n",MCMF::mcmf(s,t));
    return 0;
}