我们可以知道每行最多可以有多少个格子不用建点,设为x[i],每列同理设为y[i],那么我们连接(source,i,x[i]),(i,sink,y[i])表示我们将一个格子不建点,那么(i,j,flag[i][j]),当i,j这个格子可以建点的时候连边表示我们不在这个格子建点,那么n*m-k-最大流就是答案。  

  因为我们考虑可以在哪一个位置不放点,使得整个矩阵仍然合法,这样我们就可以知道最多有多少个合法的不建点的合法格子。  

  备注:开始想写有下限的最小可行流的着。

/**************************************************************
Problem: 1458
User: BLADEVIL
Language: C++
Result: Accepted
Time:28 ms
Memory:1788 kb
****************************************************************/ //By BLADEVIL
#include <cstdio>
#include <cstring>
#include <algorithm>
#define maxn 400
#define maxm 30010
#define inf (~0U>>1) using namespace std; int n,m,k,source,sink,ans,l;
int need[maxn][];
int flag[maxn][maxn];
int que[maxn],dis[maxn],last[maxn],pre[maxm],other[maxm],len[maxm]; void connect(int x,int y,int z) {
pre[++l]=last[x];
last[x]=l;
other[l]=y;
len[l]=z;
//if (z) printf("%d %d %d\n",x,y,z);
} int bfs() {
memset(dis,,sizeof dis);
que[]=source; dis[source]=;
int h(),t();
while (h<t) {
int cur(que[++h]);
for (int p=last[cur];p;p=pre[p]) {
if (dis[other[p]]) continue;
if (!len[p]) continue;
dis[other[p]]=dis[cur]+;
que[++t]=other[p];
if (other[p]==sink) return ;
}
}
return ;
} int dinic(int x,int flow) {
if (x==sink) return flow;
int rest=flow;
for (int p=last[x];p;p=pre[p]) {
if (!len[p]) continue;
if (dis[other[p]]!=dis[x]+) continue;
if ((!rest)||(!len[p])) continue;
int tmp=dinic(other[p],min(rest,len[p]));
if (!tmp) dis[other[p]]=;
len[p]-=tmp; len[p^]+=tmp; rest-=tmp;
}
return flow-rest;
} int main() {
scanf("%d%d%d",&n,&m,&k); l=;
for (int i=;i<=n;i++) scanf("%d",&need[i][]),need[i][]=m-need[i][];
for (int i=;i<=m;i++) scanf("%d",&need[i][]),need[i][]=n-need[i][];
for (int i=;i<=k;i++) {
int x,y; scanf("%d%d",&x,&y);
flag[x][y]=;
need[x][]--; need[y][]--;
}
for (int i=;i<=n;i++) if (need[i][]<) {
printf("JIONG!\n");
return ;
}
for (int i=;i<=m;i++) if (need[i][]<) {
printf("JIONG!\n");
return ;
}
source=n+m+; sink=source+;
for (int i=;i<=n;i++) connect(source,i,need[i][]),connect(i,source,);
for (int i=;i<=m;i++) connect(i+n,sink,need[i][]),connect(sink,i+n,);
for (int i=;i<=n;i++)
for (int j=;j<=m;j++) if (!flag[i][j])
connect(i,j+n,),connect(j+n,i,);
ans=n*m-k;
while (bfs()) ans-=dinic(source,inf);//,printf("%d\n",ans);
printf("%d\n",ans);
return ;
}

bzoj 1458 网络流的更多相关文章

  1. BZOJ 1458: 士兵占领( 网络流 )

    先判无解 把整个棋盘都放上士兵, 只需求最多可以拿走多少个士兵即可.每一行看做一个点r(i), 每一列看做一个点c(i) S->r(i), c(i)->T 连边, 容量为可以拿走的最大士兵 ...

  2. BZOJ 1458 士兵占领

    http://www.lydsy.com/JudgeOnline/problem.php?id=1458 题意:n x m的棋盘,k个位置不能放,每行和每列都有要求至少的士兵,求能否有最少的满足条件的 ...

  3. bzoj 1458: 士兵占领 -- 最大流

    1458: 士兵占领 Time Limit: 10 Sec  Memory Limit: 64 MB Description 有一个M * N的棋盘,有的格子是障碍.现在你要选择一些格子来放置一些士兵 ...

  4. BZOJ 2929 网络流

    题意是啥--. 思路: 不是与1或n连起来的边 边权是1 否则是inf 跑网络流 //By SiriusRen #include <queue> #include <cstdio&g ...

  5. BZOJ 1797 网络流的可行边&必须边

    求完网络流以后 tarjan一发 判一判 //By SiriusRen #include <queue> #include <bitset> #include <cstd ...

  6. BZOJ 1458 / Luogu P4311 士兵占领 (上下界最小流 / 直接最大流)

    做法1:上下界最小流 先来一发上下界最小流,思路比较暴力,就是把行和列看作n+mn+mn+m个点,(i,j)(i,j)(i,j)如果能占领就从第iii行向第jjj列连一条边,上界为1下界为0;然后从s ...

  7. BZOJ 3931 (网络流+最短路)

    题面 传送门 分析 考虑网络流 注意到数据包走的是最短路,所以我们只需要考虑在最短路上的边 由于最短路可能有多条,我们先跑一遍Dijkstra,然后再\(O(m)\) 遍历每条边(u,v,w) 如果d ...

  8. bzoj 1458 士兵占领(最大流)

    [题意] n行m列,第i行必须放L[i],第j列必须放C[j],有障碍格,求满足条件至少需要放多少. [思路] 至少放多少等价于最多不放多少. 对行列分别建XY点,则连边(S,Xi,a)(Yi,T,b ...

  9. 【刷题】BZOJ 1458 士兵占领

    Description 有一个M * N的棋盘,有的格子是障碍.现在你要选择一些格子来放置一些士兵,一个格子里最多可以放置一个士兵,障碍格里不能放置士兵.我们称这些士兵占领了整个棋盘当满足第i行至少放 ...

随机推荐

  1. Unity4.6新UI系统初探(uGUI)

    一.引言 Unity终于在即将到来的4.6版本内集成了所见即所得的UI解决方案(视频).事实上从近几个版本开始,Unity就在为这套系统做技术扩展,以保证最终能实现较理想的UI系统.本文试图通过初步的 ...

  2. c#中结构体(struct)和类(class)的区别

    一.类与结构的示例比较: 结构示例: public struct Person { string Name; int height; int weight public bool overWeight ...

  3. javascript中的队列结构

    1.概念 队列和栈结构不同,栈是一种后进先出的结构,而队列是一种先进先出的结构.队列也是一种表结构,不同的是队列只能在队尾插入元素,在队首删除元素,可以将队列想象成一个在超时等待排队付钱的队伍,或者在 ...

  4. VisualStudio2013+EF6+MySql5.5环境下配置

    看院子里对EF框架和MySql的配置文章不少,但是几乎出自一篇文章的转载,而且这篇转载的文章的也比较坑爹,下面我将介绍一下我的配置过程: 第一步:安装mysql-connector-net-6.9.9 ...

  5. [tomcat7源码学习]初始化之catalina.home和catalina.base(转)

    我们在代码中为了获取某个配置文件路径下的文件经常会这么写 String tomcatPath = System.getProperty("catalina.home") + &qu ...

  6. 使用POI getCell 获取空的单元格之后在使用的时候报 NullPointerException

    解决办法,在得到cell之后先判断cell是否为空,然后再进行逻辑处理. 得到的cell建议使用去除策略(如左对齐,居中等)的cell,不然有可能受到策略影响而导致结果异常. org.apache.p ...

  7. Jdev Run Page 没有反应

    从旧电脑把原有的Jdeveloper完整的拷贝至新电脑,且已完整配置JDEV_USER_HOME,JAVA_HOME等环境变量, Run Page报以下错误. [Starting OC4J using ...

  8. 构建基于WCF Restful Service的服务

    前言 传统的Asmx服务,由于遵循SOAP协议,所以返回内容以xml方式组织.并且客户端需要添加服务端引用才能使用(虽然看到网络上已经提供了这方面的Dynamic Proxy,但是没有这种方式简便), ...

  9. 将matlab中数据输出保存为txt或dat格式

    :FID= FOPEN(filename,permission) 用指定的方式打开文件 FID=+N(N是正整数):表示文件打开成功,文件代号是N. FID=-1            : 表示文件打 ...

  10. [CareerCup] 10.5 Web Crawler 网络爬虫

    10.5 If you were designing a web crawler, how would you avoid getting into infinite loops? 这道题问如果让我们 ...