平面图跑最大流 可以转换为其对偶图跑最短路 一个环对应一个割  找到最小环(即最短路)极为所求,注意辅助边的建立

加入读入优化  不过时间还是一般  估计是dij写的不好   大神勿喷~~~

/**************************************************************Problem: 1001User: 96655Language: C++Result: AcceptedTime:1724 msMemory:95120 kb
****************************************************************/#include<cstdio>
#include<cstring>
#include<algorithm>
#include<string>
#include<iostream>
#include<cstdlib>
#include<queue>
#include<map>
#include<set>
#include<vector>
#include<cmath>
#include<stack>
#include<utility>
using namespace std;
const int maxn=1001000;
const int INF=0x3f3f3f3f;
struct Edge
{int v,w,next;
} edge[maxn*6];
int head[maxn*2],p,n,m,y;
void addedge(int u,int v,int w)
{edge[p].w=w;edge[p].v=v;edge[p].next=head[u];head[u]=p++;
}
struct asd
{int x,d;asd (int a,int b):x(a),d(b) {}bool operator<(const asd &e)const{return d>e.d;}
};
priority_queue<asd>q;
int dis[maxn*2],vis[maxn*2];
void Dijkstra(int s)
{memset(vis,0,sizeof(vis));while(!q.empty())q.pop();for(int i=0; i<=y; i++)dis[i]=INF;dis[s]=0;q.push(asd(s,0));while(!q.empty()){asd e=q.top();q.pop();if(vis[e.x])continue;vis[e.x]=1;for(int i=head[e.x]; i!=-1; i=edge[i].next){int v=edge[i].v;if(dis[v]>dis[e.x]+edge[i].w){dis[v]=dis[e.x]+edge[i].w;q.push(asd(v,dis[v]));}}}
}
void read(int &x)
{char c;while((c=getchar())<'0' || c>'9');x=c-'0';while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
}
int main()
{read(n);read(m);if(n==1||m==1){if(n>m)swap(n,m);int ans=INF,a;for(int i=1; i<m; i++){scanf("%d",&a);ans=min(a,ans);}if(ans==INF)printf("0\n");else printf("%d\n",ans);return 0;}memset(head,-1,sizeof(head));p=0;y=(n-1)*(m-1)*2+1;int w,u,v;for(int i=1; i<=n; i++){for(int j=1; j<m; j++){read(w);if(i==1){u=((i-1)*(m-1)+j)*2;addedge(u,y,w);addedge(y,u,w);}else if(i==n){u=((i-2)*(m-1)+j)*2-1;addedge(0,u,w);}else{u=((i-1)*(m-1)+j)*2;v=((i-2)*(m-1)+j)*2-1;addedge(u,v,w);addedge(v,u,w);}}}for(int i=1; i<n; i++){for(int j=1; j<=m; j++){read(w);if(j==1){u=((i-1)*(m-1)+j)*2-1;addedge(0,u,w);addedge(u,0,w);}else if(j==m){u=((i-1)*(m-1)+j-1)*2;addedge(u,y,w);addedge(y,u,w);}else{u=((i-1)*(m-1)+j)*2-1;v=((i-1)*(m-1)+j-1)*2;addedge(u,v,w);addedge(v,u,w);}}}for(int i=1; i<n; i++){for(int j=1; j<m; j++){read(w);u=((i-1)*(m-1)+j)*2-1;v=((i-1)*(m-1)+j)*2;addedge(u,v,w);addedge(v,u,w);}}Dijkstra(0);printf("%d\n",dis[y]);return 0;
}

View Code

转载于:https://www.cnblogs.com/shuguangzw/p/4931470.html

bzoj 1001: [BeiJing2006]狼抓兔子 平面图最小割相关推荐

  1. BZOJ 1001: [BeiJing2006]狼抓兔子【最大流/SPFA+最小割,多解】

    1001: [BeiJing2006]狼抓兔子 Time Limit: 15 Sec  Memory Limit: 162 MB Submit: 23822  Solved: 6012 [Submit ...

  2. BZOJ 1001[BeiJing2006]狼抓兔子 最小割转最短路

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

  3. bzoj 1001: [BeiJing2006]狼抓兔子

    Time Limit: 15 Sec  Memory Limit: 162 MB Submit: 21007  Solved: 5251 [Submit][Status][Discuss] Descr ...

  4. bzoj1001 [BeiJing2006]狼抓兔子(最小割/平面图最大流转对偶图最短路)

    首先这题显然就是求一个最小割.然后跑dinic加点优化也是可以过的. 这题还有一种更妙的做法:st平面图最大流转对偶图最短路.大概做法就是首先在原图中加一条s->t的边,然后把原图的一个面当作一 ...

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

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

  6. 1001: [BeiJing2006]狼抓兔子(对偶图)

    1001: [BeiJing2006]狼抓兔子 Time Limit: 15 Sec  Memory Limit: 162 MB Submit: 23595  Solved: 5940 Descrip ...

  7. 1001: [BeiJing2006]狼抓兔子

    /**************************************************************Problem: 1001User: whymheLanguage: C+ ...

  8. 【BZOJ - 1001】狼抓兔子(无向图网络流,最小割,或平面图转对偶图求最短路SPFA)

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

  9. 对偶图 【BZOJ】1001: [BeiJing2006]狼抓兔子(对偶图+最短路)

    http://www.lydsy.com/JudgeOnline/problem.php?id=1001 可谓惨不忍睹,一下午就在调这题了. 很久以前看到这题是一眼最大流,看到n<=1000,我 ...

  10. 【bzoj 1001】狼抓兔子

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

最新文章

  1. 第十二届河南省ACM大学生程序设计竞赛赛后总结
  2. linux安装jdk、tomcat
  3. boost::locale::calendar用法的测试程序
  4. 微信公众开放平台开发04---百度BAE,java应用部署服务器,jetty了解
  5. zen3架构_AMD YES! Zen3架构亮相,消费者CPU市场的竞争热闹了
  6. linux php错误日志在哪里,宝塔的php错误日志在哪
  7. alarm/pause
  8. php将word转化成html,phpword转html
  9. \u5168\u56fd\u7f8e\u5bb9\u5927\u592b数据采集数据(\u82b1\u5bb9\u7f51 huaroo 公开数据),爬虫120例之26例
  10. 50个高端大气上档次的管理后台界面模板
  11. 简单好用又免费的百度翻译mac桌面端
  12. 机器学习实战之信用卡欺诈案列
  13. 1041. 考试座位号
  14. 用来在计算机各功能部件之间,广义相对论确认光在太阳附近会弯曲,引力透镜,水星进动等实验支持广义相对论。()...
  15. 斐波那契(兔子)数列python
  16. 软件设计的不同模型:瀑布式、快速原型法、迭代式(敏捷)
  17. Java培训完能不能就业?
  18. 博士申请 | 香港科技大学陈浩老师招收人工智能医疗方向全奖博士/博后
  19. 2019中国(深圳)IT领袖峰会大咖云集 引爆IT朋友圈...
  20. B860AV2.1-A/M/B_1g/2g-通刷_当贝纯净桌面-线刷固件包

热门文章

  1. 比特币 出块速度多少 为什么每秒7笔 以太坊15TPS
  2. thinkphp5 没有getField
  3. postges mysql_MySQL中on duplicate key update 在PostgesSQL中的实现方式
  4. redhat linux 设置ip,REDHAT LINUX企业版更改IP地址,网关,DNS和MAC地址----字符界面
  5. apo打印接口json参数_接口测试大全
  6. mysql表操作约束_MySQL操作表的约束
  7. 2.15.PHP7.1 狐教程-【PHP 抽象类、抽象方法】
  8. 2.4.PHP7.1 狐教程-【PHP常量】
  9. Kafka 安装部署、集群启动、命令行操作 与 可视化工具 Kafka Tool
  10. Java 23 种设计模式概述 与 简单工厂模式