裸dinic就跑过去了,哪用得着平面图最小割=最短路……

#include <iostream>
#include <cstring>
#include <cstdio>
#include <queue>
using namespace std;
int n, m, w, ss, tt, hea[1000005], cnt, lev[1000005], cur[1000005], maxFlow;
const int oo=0x3f3f3f3f;
queue<int> d;
struct Edge{
int too, nxt, val;
}edge[8000005];
void add_edge(int fro, int too, int val){
edge[cnt].nxt = hea[fro];
edge[cnt].too = too;
edge[cnt].val = val;
hea[fro] = cnt++;
}
void addEdge(int fro, int too, int val){
add_edge(fro, too, val);
add_edge(too, fro, val);//无向边一开始就能增广
}
int f(int a, int b){
return (a-1)*m+b;
}
bool bfs(){
memset(lev, 0, sizeof(lev));
lev[ss] = 1;
d.push(ss);
while(!d.empty()){
int x=d.front();
d.pop();
for(int i=hea[x]; i!=-1; i=edge[i].nxt){
int t=edge[i].too;
if(!lev[t] && edge[i].val>0){
lev[t] = lev[x] + 1;
d.push(t);
}
}
}
return lev[tt]!=0;
}
int dfs(int x, int lim){
if(x==tt) return lim;
int addFlow=0;
for(int &i=cur[x]; i!=-1; i=edge[i].nxt){
int t=edge[i].too;
if(lev[t]==lev[x]+1 && edge[i].val>0){
int tmp=dfs(t, min(lim-addFlow, edge[i].val));
edge[i].val -= tmp;
edge[i^1].val += tmp;
addFlow += tmp;
if(addFlow==lim) break;
}
}
return addFlow;
}
void dinic(){
while(bfs()){
for(int i=ss; i<=tt; i++) cur[i] = hea[i];
maxFlow += dfs(ss, oo);
}
}
int main(){
memset(hea, -1, sizeof(hea));
cin>>n>>m;
for(int i=1; i<=n; i++)
for(int j=1; j<m; j++){
scanf("%d", &w);
addEdge(f(i,j), f(i,j+1), w);
}
for(int i=1; i<n; i++)
for(int j=1; j<=m; j++){
scanf("%d", &w);
addEdge(f(i,j), f(i+1,j), w);
}
for(int i=1; i<n; i++)
for(int j=1; j<m; j++){
scanf("%d", &w);
addEdge(f(i,j), f(i+1,j+1), w);
}
ss = f(1,1); tt = f(n,m);
dinic();
cout<<maxFlow<<endl;
return 0;
}

平面图最短路:

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
int n, m, ss, tt, uu, hea[2000005], cnt, dis[2000005], din;
bool vis[2000005];
const int oo=0x3f3f3f3f;
int p(int a, int b){
return (a-1)*(m-1)+b;
}
int q(int a, int b){
return p(a,b)+(n-1)*(m-1);
}
struct Edge{
int too, nxt, val;
}edge[8000005];
struct Node{
int idx, val;
}nd[9000005];
bool cmp(Node x, Node y){
return x.val>y.val;
}
void add_edge(int fro, int too, int val){
edge[++cnt].nxt = hea[fro];
edge[cnt].too = too;
edge[cnt].val = val;
hea[fro] = cnt;
}
void addEdge(int fro, int too, int val){
add_edge(fro, too, val);
add_edge(too, fro, val);
}
void dijk(){
memset(dis, 0x3f, sizeof(dis));
dis[ss] = 0;
nd[++din] = (Node){ss, 0};
while(din){
Node x=nd[1];
pop_heap(nd+1, nd+1+din, cmp);
din--;
if(vis[x.idx]) continue;
vis[x.idx] = true;
for(int i=hea[x.idx]; i; i=edge[i].nxt){
int t=edge[i].too;
if(!vis[t] && dis[t]>dis[x.idx]+edge[i].val){
dis[t] = dis[x.idx] + edge[i].val;
nd[++din] = (Node){t, dis[t]};
push_heap(nd+1, nd+1+din, cmp);
}
}
}
}
int main(){
cin>>n>>m;
if(n==1 || m==1){
int tmp=oo;
for(int i=1; i<=max(n,m); i++){
scanf("%d", &uu);
tmp = min(tmp, uu);
}
cout<<tmp<<endl;
return 0;
}
ss = 0; tt = q(n-1,m-1) + 1;
for(int i=1; i<=n; i++)
for(int j=1; j<m; j++){
scanf("%d", &uu);
if(i==1) addEdge(ss, p(i,j), uu);
else if(i==n) addEdge(q(i-1,j), tt, uu);
else addEdge(q(i-1,j), p(i,j), uu);
}
for(int i=1; i<n; i++)
for(int j=1; j<=m; j++){
scanf("%d", &uu);
if(j==1) addEdge(tt, q(i,j), uu);
else if(j==m) addEdge(p(i,j-1), ss, uu);
else addEdge(p(i,j-1), q(i,j), uu);
}
for(int i=1; i<n; i++)
for(int j=1; j<m; j++){
scanf("%d", &uu);
addEdge(p(i,j), q(i,j), uu);
}
dijk();
cout<<dis[tt]<<endl;
return 0;
}

luogu4001 [BJOI2006]狼抓兔子的更多相关文章

  1. P4001 [BJOI2006]狼抓兔子(对偶图)

    P4001 [BJOI2006]狼抓兔子 最短路+对偶图 看这题最容易想到的就是网络流.Dinic可以过,据说还跑得比正解快. 如果不写网络流,那么需要知道2个前置知识:平面图和对偶图(右转baidu ...

  2. BJOI2006狼抓兔子

    1001: [BeiJing2006]狼抓兔子 Time Limit: 15 Sec  Memory Limit: 162 MBSubmit: 9967  Solved: 2267[Submit][S ...

  3. 752. [BJOI2006] 狼抓兔子

    ★★★   输入文件:bjrabbit.in   输出文件:bjrabbit.out   简单对比时间限制:1 s   内存限制:162 MB Description   Source: Beijin ...

  4. [BJOI2006]狼抓兔子

    题目描述 现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到,但抓兔子还是比较在行的,而且现在的兔子还比较笨,它们只有两个窝,现在你做为狼王,面对下面这样一个网格的地形: ...

  5. 1001. [BJOI2006]狼抓兔子【最小割】

    Description 现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到,但抓兔子还是比较在行的, 而且现在的兔子还比较笨,它们只有两个窝,现在你做为狼王,面对下面这样一 ...

  6. BZOJ1001 BJOI2006 狼抓兔子

    Description 现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到,但抓兔子还是比较在行的,而且现在的兔子还比较笨,它们只有两个窝,现在你做为狼王,面对下面这样一个 ...

  7. [BJOI2006]狼抓兔子(网络流)

    题目描述 现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到,但抓兔子还是比较在行的,而且现在的兔子还比较笨,它们只有两个窝,现在你做为狼王,面对下面这样一个网格的地形: ...

  8. BZOJ1001 洛谷4001 [BJOI2006]狼抓兔子 题解

    题目 这个题目有多种解法,这个题也是一个比较经典的题了,正是因为他的多样的做法,这个题主要难在建图和优化,因为这是一个网格图,所以spfa肯定过不去,所以用最短路解法的话,只能用dij,而网络流也是要 ...

  9. BZOJ1001 BJOI2006狼抓兔子(最小割+最短路)

    显然答案就是最小割.直接跑dinic也能过,不过显得不太靠谱. 考虑更正确的做法.作为一个平面图,如果要把他割成两半,那么显然可以用一条曲线覆盖且仅覆盖所有割边.于是我们把空白区域看成点,隔开他们的边 ...

随机推荐

  1. 轻量级的绘制图表js库--Morris.js

    Morris.js 是一个轻量级的 JS 库,使用 jQuery 和 Raphaël 来生成各种时序图. 虽说现在移动手机网络已经到了4G,但是在移动web客户端开发过中,为了达到良好的体验效果,需要 ...

  2. VC 对话框设置背景颜色和图片

    改变背景颜色,有两种方法: 1.在app的初始化函数中调用:void SetDialogBkColor( COLORREF clrCtlBk = RGB(192, 192, 192), COLORRE ...

  3. FMDB浅析(思想)

    http://www.cnblogs.com/OTgiraffe/p/5931800.html 一.FMDB介绍 FMDB是一种第三方的开源库,FMDB就是对SQLite的API进行了封装,加上了面向 ...

  4. mybatis(一):思维导图

  5. linux readahead

    blockdev --getra /dev/sda blockdev --setra 2048 /dev/sda 必须将其写入配置文件/etc/rc.local,否则重启就会失效.[root@loca ...

  6. iOS--UIScrollView基本用法和代理方法

    主要是为了记录下UIScrollView的代理方法吧 在帮信息学院的学长做东西的时候需要大量用到分块浏览,所以就涉及到很多的关于scrollview,所以也就有了这篇文章   - (void)view ...

  7. C++ 学习笔记 开篇

    从大一开始学习C语言,大学期间做了许多嵌入式的开发项目,毕业后从事嵌入式开发工作主要的开发语言也是C语言.虽然期间断断续续的学习过C++,做过QT.C#上位机但也只是在其他语言的外壳下使用C在开发,始 ...

  8. vmware虚拟机安装Windows 7后虚拟机自动挂起

    vmware虚拟机安装windows7后在一段时间中没有操作,虚拟机会自动挂起,是因为windows7中的设置的自动睡眠,打开[控制面板]=>[电源选项]=>[选择关闭显示器时间]将下面两 ...

  9. ubuntu14.04搭建LAMP环境(nginx,php,mysql,linux)详解

    最近更换开发环境至ubuntu,整理开发环境和常用软件的安装配置(更新排版) 以下安装过程经过多次操作得出,参照步骤进行操作即可 一.LAMP基本环境搭建 1 切换root账号 sudo su 2,安 ...

  10. Golang tcp转发 remoteAddr错误

    Golang tcp 转发 第一版本 accept获取的Conn里的localAddr做为源地址,remoteAddr来做为目的地址 // tcpForward package main import ...