BZOJ 3171 [Tjoi2013]循环格(费用流)
【题目链接】 http://www.lydsy.com/JudgeOnline/problem.php?id=3171
【题目大意】
一个循环格就是一个矩阵,其中所有元素为箭头,指向相邻四个格子。
每个元素有一个坐标(行,列),其中左上角元素坐标为(0,0)。给定一个起始位置(r,c)
你可以沿着箭头防线在格子间行走。即如果(r,c)是一个左箭头,那么走到(r,c-1);
如果是右箭头那么走到(r,c+1);如果是上箭头那么走到(r-1,c);
如果是下箭头那么走到(r+1,c);每一行和每一列都是循环的,
即如果走出边界,你会出现在另一侧。一个完美的循环格是这样定义的:
对于任意一个起始位置,你都可以i沿着箭头最终回到起始位置。
如果一个循环格不满足完美,你可以随意修改任意一个元素的箭头直到完美。
给定一个循环格,你需要计算最少需要修改多少个元素使其完美。
【题解】
我们发现一个完美的循环格每个元素都只有一个入度和一个出度,
因此我们将其拆分为入点和出现,入点连汇点,出点连源点,流量为1,费用为0,
对于每个箭头指向的反向,连相应的出点到入点,流量为1,费用为0,
对于相邻但是箭头不指向的地方,我们将其相互连接,流量为1,费用为1,
求最小费用最大流即答案。
【代码】
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int INF=0x7fffffff,N=1010;
int S,T,cnt,ans,d[N],q[N],from[N],g[N],flow,tot;
bool in[1010];
struct edge{int from,to,nxt,c,v;}e[100010];
void add(int u,int v,int w,int c){
e[++cnt].from=u;e[cnt].to=v;
e[cnt].nxt=g[u];g[u]=cnt;
e[cnt].c=c;e[cnt].v=w;
}void add_edge(int u,int v,int w,int c){add(u,v,w,c);add(v,u,0,-c);}
bool spfa(){
for(int i=S;i<=T;i++)d[i]=INF;
int t=0,w=1;d[S]=0;in[S]=1;q[0]=S;
while(t!=w){
int now=q[t];t++;if(t==T)t=0;
for(int i=g[now];i;i=e[i].nxt)
if(e[i].v&&d[e[i].to]>d[now]+e[i].c){
d[e[i].to]=d[now]+e[i].c;from[e[i].to]=i;
if(!in[e[i].to]){in[e[i].to]=1;q[w++]=e[i].to;if(w==T)w=0;}
}in[now]=0;
}return(d[T]!=INF);
}
void mcf(){
int x=INF;
for(int i=from[T];i;i=from[e[i].from])x=min(x,e[i].v);flow+=x;
for(int i=from[T];i;i=from[e[i].from]){e[i].v-=x;e[i^1].v+=x;ans+=e[i].c*x;}
}
const int dx[]={1,-1,0,0},dy[]={0,0,1,-1};
int R,C,mp[20][20];
char str[20][20];
int main(){
while(~scanf("%d%d",&R,&C)){
memset(g,0,sizeof(g));
cnt=1; ans=flow=0;
S=0,T=C*R*2+1;
for(int i=1;i<=R;i++)scanf("%s",str[i]+1);
for(int i=1;i<=R;i++)for(int j=1;j<=C;j++){
if(str[i][j]=='U')mp[i][j]=1;
if(str[i][j]=='D')mp[i][j]=0;
if(str[i][j]=='L')mp[i][j]=3;
if(str[i][j]=='R')mp[i][j]=2;
for(int k=0;k<4;k++){
int x=i+dx[k],y=j+dy[k];
if(x>R)x=1;if(x<1)x=R;
if(y>C)y=1;if(y<1)y=C;
if(k==mp[i][j])add_edge((i-1)*C+j,(x-1)*C+y+C*R,1,0);
else add_edge((i-1)*C+j,(x-1)*C+y+C*R,1,1);
}add_edge(S,(i-1)*C+j,1,0);
add_edge((i-1)*C+j+C*R,T,1,0);
}while(spfa())mcf();
printf("%d\n",ans);
}return 0;
}
BZOJ 3171 [Tjoi2013]循环格(费用流)的更多相关文章
- Bzoj 3171: [Tjoi2013]循环格 费用流
3171: [Tjoi2013]循环格 Time Limit: 1 Sec Memory Limit: 128 MBSubmit: 741 Solved: 463[Submit][Status][ ...
- bzoj 3171: [Tjoi2013]循环格
#include<cstdio> #include<iostream> #include<cstring> #define M 10000 #define inf ...
- [TJOI2013]循环格 费用流 BZOJ3171
题目背景 一个循环格就是一个矩阵,其中所有元素为箭头,指向相邻四个格子.每个元素有一个坐标(行,列),其中左上角元素坐标为(0,0).给定一个起始位(r,c),你可以沿着箭头方向在格子间行走.即:如果 ...
- bzoj 3171: [Tjoi2013]循环格 最小费用最大流
题目大意: http://www.lydsy.com/JudgeOnline/problem.php?id=3171 题解: 首先我们很容易发现一个结论: 出现完美循环当且仅当所有点的出入度均为1 所 ...
- bzoj 3171 [Tjoi2013]循环格(MCMF)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=3171 [题意] 给定一个方向矩阵,要求改变最少的格子,使得任意一个点都在一个环中. [ ...
- BZOJ 3171 循环格(费用流)
题意 一个循环格就是一个矩阵,其中所有元素为箭头,指向相邻四个格子.每个元素有一个坐标(行,列),其中左上角元素坐标为(0,0).给定一个起始位置(r,c),你可以沿着箭头防线在格子间行走.即如果(r ...
- 3171. [TJOI2013]循环格【费用流】
Description 一个循环格就是一个矩阵,其中所有元素为箭头,指向相邻四个格子.每个元素有一个坐标(行,列),其中左上角元素坐标为(0,0).给定一个起始位置(r,c) ,你可以沿着箭头防线在格 ...
- BZOJ_3171_[Tjoi2013]循环格_最小费用最大流
BZOJ_3171_[Tjoi2013]循环格_最小费用最大流 Description 一个循环格就是一个矩阵,其中所有元素为箭头,指向相邻四个格子.每个元素有一个坐标(行,列),其中左上角元素坐标为 ...
- [Tjoi2013]循环格
[Tjoi2013]循环格 2014年3月18日1,7500 Description Input 第一行两个整数R,C.表示行和列,接下来R行,每行C个字符LRUD,表示左右上下. Output 一个 ...
随机推荐
- C++之指针,引用与数组
引用只是对象的另一个名字,通过在变量名前面添加"&”符号来定义,而指针保存的是另一个对象的地址,它们两都提供了间接访问所服务变量的途径. 但是它们的差别还是挺大的: 先从它们的值说起 ...
- Lodash使用示例(比较全)
<html> <head> <meta name="viewport" content="width=device-width" ...
- c++设计模式系列----factory模式
问题: 假设我们要开发一个游戏--打怪物,首先,游戏有分等级,假设有初级,中级两个个等级(就不用flappy bird模式了,那个比较特殊,对一个玩家来说是难以具有持久吸引力的!),不同的等级怪物也是 ...
- Oracle 获取ddl语句
--得到所有表空间的ddl语句 SELECT DBMS_METADATA.GET_DDL('TABLESPACE', TS.tablespace_name)FROM DBA_TABLESPACES T ...
- expose a port on a living Docker container
if you have a container that with something running on its port 8000, you can run wget http://contai ...
- service XXX start启动报start: Rejected send message, 1 matche
转,原文地址:http://blog.sina.com.cn/s/blog_56d8ea9001018w1l.html [问题]start: Rejected send messag现象:crifan ...
- PHP的数据类型
原始类型共8种: 1, 4种标量类型:boolean(布尔型).integer(整形).float/double(浮点型).string(字符串型): 2, 2种复合型:array(数组).o ...
- mycncart自定义主题
本文是自己通过其他主题,自学的,如果有什么问题,可以提出建议? 参考资料:opencart官网 www.opencart.com 或 mycncart的官网上的一些教程 www.mycncart.c ...
- 取消myeclipse自动进入workspace
进入到C:\Program Files\MyEclipse 6.5\eclipse\configuration\.settings目录, 修改SHOW_WORKSPACE_SELECTION_DIAL ...
- Tomcat手动指定jdk路径
Windows系统: 修改bin下的setclasspath.bat,在开头处添加: rem set specific jdkset JAVA_HOME=%CATALINA_HOME%\jdk1.8. ...