【对偶图转最短路】BZOJ2007 [Noi2010]海拔

题面在这里

首先注意到海拔只可能是0或1

首先海拔不可能大于1,然后考虑与h相连的路口均为0,1

则代价为\(h\cdot a+(1-h)\cdot b\ge min\{ a,b \}\),a,b分别为与0,1相连的总流量

把单点类比到联通块上,这样的话只有西北一片0和东南一片1了

代价只在分界线上产生,其实就是要找最小割

这个恰好是平面图,对偶图刷最短路即可

示例程序:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
#define cl(x,y) memset(x,y,sizeof(x))
using namespace std;
typedef long long ll;
inline char nc(){
static char buf[100000],*l=buf,*r=buf;
return l==r&&(r=(l=buf)+fread(buf,1,100000,stdin),l==r)?EOF:*l++;
}
inline int red(){
int res=0,f=1;char ch=nc();
while (ch<'0'||'9'<ch) {if (ch=='-') f=-f;ch=nc();}
while ('0'<=ch&&ch<='9') res=res*10+ch-48,ch=nc();
return res*f;
}

const int maxn=505,maxnn=250005,maxe=1000005,p[4][2]={{-1,0},{0,1},{1,0},{0,-1}};
int n,a[maxn][maxn][4],S,T;
int tot,son[maxe],nxt[maxe],lnk[maxnn],w[maxe];
inline void add(int x,int y,int z){
son[++tot]=y,nxt[tot]=lnk[x],lnk[x]=tot,w[tot]=z;
}
#define getid(x,y) (((x)-1)*n+(y))
ll dst[maxnn];
struct data{
int x;ll w;
data () {}
data (int _x,ll _w) {x=_x;w=_w;}
bool operator<(const data&b)const {return w>b.w;}
};
priority_queue<data> H;
void DIJ(){
cl(dst,63);
H.push(data(S,0));dst[S]=0;
while (!H.empty()){
if (H.top().w>dst[H.top().x]) {H.pop();continue;}
int k=H.top().x; H.pop();
for (int j=lnk[k];j;j=nxt[j])
if (dst[son[j]]>dst[k]+w[j])
dst[son[j]]=dst[k]+w[j],H.push(data(son[j],dst[son[j]]));
}
}
int main(){
n=red();S=n*n+1;T=S+1;
for (int i=0;i<=n;i++)
for (int j=1;j<=n;j++) a[i][j][0]=red();
for (int i=1;i<=n;i++)
for (int j=0;j<=n;j++) a[i][j][1]=red();
for (int i=0;i<=n;i++)
for (int j=1;j<=n;j++) a[i][j][2]=red();
for (int i=1;i<=n;i++)
for (int j=0;j<=n;j++) a[i][j][3]=red();
//make graph:
for (int i=1;i<=n;i++)
for (int j=1;j<=n;j++)
for (int t=0;t<4;t++){
int x=i+p[t][0],y=j+p[t][1],w;
if (y<1) {add(S,getid(i,j),a[x][y][(t+2)%4]);continue;}
if (x>n) {add(S,getid(i,j),a[i][j][(t+2)%4]);continue;}
if (t==1||t==2) w=a[i][j][t];else w=a[x][y][t];
if (x<1||y>n) add(getid(i,j),T,w);else
add(getid(i,j),getid(x,y),w);
}
DIJ();
printf("%lld",dst[T]);
return 0;
}