【BZOJ 3171】 [Tjoi2013]循环格

Description

一个循环格就是一个矩阵，其中所有元素为箭头，指向相邻四个格子。每个元素有一个坐标（行，列），其中左上角元素坐标为（0,0）。给定一个起始位置（r，c）

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

 

Input

第一行两个整数R，C。表示行和列，接下来R行，每行C个字符LRUD，表示左右上下。

Output

一个整数，表示最少需要修改多少个元素使得给定的循环格完美

Sample Input

3 4
RRRD
URLL
LRRR

Sample Output

2

HINT

 

1<=R,L<=15

 

或许很明显，但我真真没看出这是一个网络流

hzwer题解

每个格子入度为一，出度为一

拆点，每个格子向它周围格子的右部点连边，费用为0或1，跑最小费用流

#include<cstdio>
#include<cstring>
using namespace std;
const int N=100,inf=1000000;
struct ee{int to,next,f,w;}e[N*N*200];
int S,T,cnt=1,n,k,ans,timer,m,u,v,w;
int head[N*N],dis[N*N],pre[N*N],q[N*N],map[N][N];
bool inq[N];
int fx[4]={-1,1,0,0},fy[4]={0,0,-1,1};
void ins(int u,int v,int f,int w){
    e[++cnt].to=v,e[cnt].next=head[u],e[cnt].f=f,e[cnt].w=w,head[u]=cnt;
    e[++cnt].to=u,e[cnt].next=head[v],e[cnt].f=0,e[cnt].w=-w,head[v]=cnt;
}
  
bool spfa(){
    for (int i=1;i<=T;i++) dis[i]=inf;
    int h=0,t=1;
    q[t]=S;dis[S]=0;inq[S]=1;
    while (h!=t){
        int now=q[++h];if(h==T) h=0;
        for (int i=head[now];i;i=e[i].next){
            int v=e[i].to;
            if (dis[v]>dis[now]+e[i].w&&e[i].f){
                dis[v]=dis[now]+e[i].w;
                pre[v]=i;
                if (!inq[v]){
                    q[++t]=v;if (t==T) t=0;
                    inq[v]=1;
                }
            }
        }
        inq[now]=0;
    }
    if (dis[T]==inf) return 0;
    return 1;
}
  
void updata(){
    int tmp=T;
    while (tmp!=S){
        int l=pre[tmp],v=e[l].to;
        e[l].f-=1;e[l^1].f+=1;
        tmp=e[l^1].to;
    }
    ans+=dis[T];
}
int main(){
    scanf("%d%d",&n,&m);
    T=n*m*2+1;
    char s[30];
    for(int i=1;i<=n;i++){
        scanf("%s",s+1);
        for (int j=1;j<=m;j++){
            if(s[j]=='U')map[i][j]=0;
            if(s[j]=='D')map[i][j]=1;
            if(s[j]=='L')map[i][j]=2;
            if(s[j]=='R')map[i][j]=3; 
        }
    }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            int u=(i-1)*m+j;
            ins(S,u,1,0);ins(u+n*m,T,1,0);
            for(int k=0;k<4;k++){
                int nx=i+fx[k],ny=j+fy[k];
                if(nx>n)nx=1;if(nx<1)nx=n;
                if(ny>m)ny=1;if(ny<1)ny=m;
                if(k==map[i][j]) ins(u,(nx-1)*m+ny+m*n,1,0);
                else ins(u,(nx-1)*m+ny+m*n,1,1);
            }
        }
    while (spfa()) 
    updata();
    printf("%d",ans);
}