1001: [BeiJing2006]狼抓兔子
Time Limit: 15 Sec Memory Limit: 162 MBSubmit: 25653 Solved: 6512
[Submit][Status][Discuss]
Description
现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到，但抓兔子还是比较在行的，
而且现在的兔子还比较笨，它们只有两个窝，现在你做为狼王，面对下面这样一个网格的地形：
左上角点为(1,1),右下角点为(N,M)(上图中N=4,M=5).有以下三种类型的道路
1:(x,y)<==>(x+1,y)
2:(x,y)<==>(x,y+1)
3:(x,y)<==>(x+1,y+1)
道路上的权值表示这条路上最多能够通过的兔子数，道路是无向的. 左上角和右下角为兔子的两个窝，
开始时所有的兔子都聚集在左上角(1,1)的窝里，现在它们要跑到右下解(N,M)的窝中去，狼王开始伏击
这些兔子.当然为了保险起见，如果一条道路上最多通过的兔子数为K，狼王需要安排同样数量的K只狼，
才能完全封锁这条道路，你需要帮助狼王安排一个伏击方案，使得在将兔子一网打尽的前提下，参与的
狼的数量要最小。因为狼还要去找喜羊羊麻烦.
Input
第一行为N,M.表示网格的大小，N,M均小于等于1000.
接下来分三部分
第一部分共N行，每行M-1个数，表示横向道路的权值.
第二部分共N-1行，每行M个数，表示纵向道路的权值.
第三部分共N-1行，每行M-1个数，表示斜向道路的权值.
输入文件保证不超过10M
Output
输出一个整数，表示参与伏击的狼的最小数量.
Sample Input
3 4
5 6 4
4 3 1
7 5 3
5 6 7 8
8 7 6 5
5 5 5
6 6 6
5 6 4
4 3 1
7 5 3
5 6 7 8
8 7 6 5
5 5 5
6 6 6
Sample Output
14
HINT
2015.4.16新加数据一组，可能会卡掉从前可以过的程序。
Source
这道题所有人看的第一眼都会认为这不是一道网络流的水题么
可是数据范围告诉我们会超时的
这个图是对称的，兔子过不去就相当于从左下角到右上角把路封死，就像这样
我们可以从左下角往右上角一次连边
最后跑最短路即可
具体实现看代码
#include <bits/stdc++.h>
#define ll long long
using namespace std;
inline int read(){
int x=0;int f=1;char ch=getchar();
while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
while(isdigit(ch)) {x=x*10+ch-'0';ch=getchar();}
return x*f;
}
const int MAXN=2e6+10;
namespace zhangenming{
struct node{
int y,next;
}e[MAXN<<2];
int linkk[MAXN<<2],len=0,n,m,sum1,sum2,a[4000010];
inline void insert(int xx,int yy){
e[++len].y=yy;e[len].next=linkk[xx];linkk[xx]=len;
}
void init(){
n=read();m=read();
}
int q[MAXN],vis[MAXN<<1]={},dis[MAXN<<1],head=0,tail=0;
void spfa(){
memset(dis,10,sizeof(dis));
for(int i=1;i<m;i++){
int tn=(n-1)*(m-1)+i;
q[++tail]=tn;dis[tn]=0;vis[tn]=1;
}
for(int i=1;i<n;i++){
int tn=(m-1)*n+(i-1)*m+1;
q[++tail]=tn;dis[tn]=0;vis[tn]=1;
}
int maxn=2000000;
while(head!=tail){
if(head==maxn) head=0;
int tn=q[++head];
for(int i=linkk[tn];i;i=e[i].next){
if(dis[e[i].y]>dis[tn]+a[tn]){
dis[e[i].y]=dis[tn]+a[tn];
if(!vis[e[i].y]){
q[++tail]=e[i].y;
if(tail==maxn) tail=0;
vis[e[i].y]=1;
}
}
}
vis[tn]=0;
}
}
inline int heng(int xx,int yy){
return (xx-1)*(m-1)+yy;
}
inline int shu(int xx,int yy){
return (sum1+(xx-1)*m+yy);
}
inline int xie(int xx,int yy){
return sum2+(xx-1)*(m-1)+yy;
}
void solve(){
sum1=n*(m-1);
sum2=sum1+(n-1)*m;
for(int i=1;i<=n;i++){
for(int j=1;j<m;j++){
int tn=heng(i,j);
a[tn]=read();
}
}
for(int i=1;i<n;i++){
for(int j=1;j<=m;j++){
int tn=shu(i,j);
a[tn]=read();
}
}
for(int i=1;i<n;i++){
for(int j=1;j<m;j++){
int tn=xie(i,j);
a[tn]=read();
}
}
for(int i=1;i<n;i++){
for(int j=1;j<m;j++){
int tn=heng(i,j);int k1=shu(i,j+1);int k2=xie(i,j);
insert(k1,tn);insert(k2,tn);
}
}
for(int i=n;i>1;i--){
for(int j=1;j<m;j++){
int tn=heng(i,j);int k1=shu(i-1,j);int k2=xie(i-1,j);
insert(tn,k1);insert(tn,k2);
}
}
for(int i=1;i<n;i++){
for(int j=1;j<m;j++){
int tn=xie(i,j);
int k1=shu(i,j);int k2=shu(i,j+1);
insert(k1,tn);insert(tn,k2);
}
}
spfa();
int minn=10000000;
for(int i=1;i<n;i++){
int tn=shu(i,m);
minn=min(minn,dis[tn]+a[tn]);
}
for(int i=1;i<m;i++){
int tn=heng(1,i);
minn=min(minn,dis[tn]+a[tn]);
}
cout<<minn<<endl;
}
}
int main(){
//freopen("All.in","r",stdin);
//freopen("aa.out","w",stdout);
using namespace zhangenming;
init();
solve();
return 0;
}