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
RRRD
URLL
LRRR
Sample Output
2
题解:
题目要求这个图要变成一个环,
环的话,出度入度都为1,
所以我们跑个费用流,
把每个点拆成两个,
st->x1 连一条流量为1,费用为0
x2->ed 连一条流量为1,费用为0
分别表示出入度。
每个点x1再向它周围的x2连边,如果方向相同,则建一条流量为1,费用为0
否则,建一条流量为1,费用为1的
代码;
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#define INF 1e9
using namespace std;
#define get(i,j) (i-1)*m+j
const int N=12000;
int n,m;
char s[20];
int map[20][20];
int fx[4]={-1,1,0,0},fy[4]={0,0,-1,1};
int st,ed;
struct node{
int x,y,z,cost,next,other;
}sa[N<<1];int len=0,first[N];
void ins(int x,int y,int z,int cost)
{
len++;
sa[len].x=x;
sa[len].y=y;
sa[len].z=z;
sa[len].cost=cost;
sa[len].next=first[x];
first[x]=len;
sa[len].other=len+1;
len++;
sa[len].x=y;
sa[len].y=x;
sa[len].z=0;
sa[len].cost=-cost;
sa[len].next=first[y];
first[y]=len;
sa[len].other=len-1;
}
int pre[N],frpe[N],dis[N];
bool vis[N];
queue<int>q;
bool spfa()
{
memset(pre,-1,sizeof(pre));pre[st]=0;
memset(vis,0,sizeof(vis));
for(int i=0;i<=ed;i++) dis[i]=INF;
dis[st]=0;
vis[st]=1;
q.push(st);
while(!q.empty())
{
int x=q.front();q.pop();
for(int i=first[x];i!=-1;i=sa[i].next)
{
int y=sa[i].y;
if(sa[i].z>0&&dis[y]>dis[x]+sa[i].cost)
{
dis[y]=dis[x]+sa[i].cost;
if(!vis[y])
{
q.push(y);
vis[y]=1;
}
pre[y]=x;
frpe[y]=i;
}
}
vis[x]=0;
}
if(dis[ed]>=INF) return false;
return true;
}
void MCMF()
{
int ans=0,liu=0;
while(spfa())
{
int minn=INF;
for(int i=ed;i!=st;i=pre[i])
{
if(minn>sa[frpe[i]].z)
minn=sa[frpe[i]].z;
}
ans+=minn*dis[ed];
liu+=minn;
for(int i=ed;i!=st;i=pre[i])
{
sa[frpe[i]].z-=minn;
sa[sa[frpe[i]].other].z+=minn;
}
}
printf("%d\n",ans);
}
int main()
{
scanf("%d%d",&n,&m);
st=0;ed=n*m*2+1;
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;
}
}
memset(first,-1,sizeof(first));
len=1;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{
ins(st,get(i,j),1,0);
ins(get(i,j)+m*n,ed,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(get(i,j),get(nx,ny)+m*n,1,0);
else ins(get(i,j),get(nx,ny)+m*n,1,1);
}
}
MCMF();
}