bzoj2007 NOI2010 网络流转对偶图

2007: [Noi2010]海拔

Time Limit: 20 Sec  Memory Limit: 552 MB
Submit: 2775  Solved: 1331
[Submit][Status][Discuss]

Description

YT市是一个规划良好的城市，城市被东西向和南北向的主干道划分为n×n个区域。简单起见，可以将YT市看作一个

正方形，每一个区域也可看作一个正方形。从而，YT城市中包括(n+1)×(n+1)个交叉路口和2n×(n+1)条双向道路

（简称道路），每条双向道路连接主干道上两个相邻的交叉路口。下图为一张YT市的地图(n = 2)，城市被划分为2

×2个区域，包括3×3个交叉路口和12条双向道路。 小Z作为该市的市长，他根据统计信息得到了每天上班高峰期

间YT市每条道路两个方向的人流量，即在高峰期间沿着该方向通过这条道路的人数。每一个交叉路口都有不同的海

拔高度值，YT市市民认为爬坡是一件非常累的事情，每向上爬h的高度，就需要消耗h的体力。如果是下坡的话，则

不需要耗费体力。因此如果一段道路的终点海拔减去起点海拔的值为h(注意h可能是负数)，那么一个人经过这段路

所消耗的体力是max{0, h}（这里max{a, b}表示取a, b两个值中的较大值）。 小Z还测量得到这个城市西北角的交

叉路口海拔为0，东南角的交叉路口海拔为1(如上图所示)，但其它交叉路口的海拔高度都无法得知。小Z想知道在

最理想的情况下（即你可以任意假设其他路口的海拔高度），每天上班高峰期间所有人爬坡所消耗的总体力和的最

小值。

Input

第一行包含一个整数n，含义如上文所示。接下来4n(n + 1)行，每行包含一个非负整数分别表示每一条道路每一个

方向的人流量信息。输入顺序：n(n + 1)个数表示所有从西到东方向的人流量，然后n(n + 1)个数表示所有从北到

南方向的人流量，n(n + 1)个数表示所有从东到西方向的人流量，最后是n(n + 1)个数表示所有从南到北方向的人

流量。对于每一个方向，输入顺序按照起点由北向南，若南北方向相同时由西到东的顺序给出(参见样例输入)。

Output

仅包含一个数，表示在最理想情况下每天上班高峰期间所有人爬坡所消耗的总体力和（即总体力和的最小值），结

果四舍五入到整数。

Sample Input

1
1
2
3
4
5
6
7
8

Sample Output

3
【样例说明】
样例数据见下图。
最理想情况下所有点的海拔如上图所示。
对于100%的数据：1 ≤ n ≤ 500，0 ≤ 流量 ≤ 1,000,000且所有流量均为整数。

 

80% n，m<=40

这个题的关键之处，就是得出一个结论：每个点要么为1要么为0，感性理解就可以了

然后就成了明显的最小割模型，割出每个点要么为1要么为0，然后就可以过80分了

 

对于全部数据，可以用对偶图做

考虑到这是一个emmmmmm 规则平面图（？应该是这个名字），可以构建其对偶图，用最短路来求最小割

 

对偶图

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
#define mp make_pair
#define fi first
#define se second
#define inf 0x3f3f3f3f
#define ll long long
#define N 260005
using namespace std;
int n,S,T,tot,hd[N],d[N],vis[N];
struct edge{int v,w,next;}e[N*10];
void adde(int u,int v,int w){
    e[++tot].v=v;
    e[tot].w=w;
    e[tot].next=hd[u];
    hd[u]=tot;
}
char gc(){
    static char s[1000000],*p1,*p2;
    if(p1==p2)p2=(p1=s)+fread(s,1,1000000,stdin);
    if(p1==p2)return EOF;
    return *p1++;
}

void r(int &x){
    x=0;int f=1;char ch=gc();
    while(ch>'9'||ch<'0'){
        if(ch=='-')f=-1;
        ch=gc();
    }
    while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=gc();
    x*=f;
}
typedef pair<int,int>pii;
priority_queue<pii,vector<pii>,greater<pii> >q;
void dijkstra(){
    for(int i=S;i<=T;i++)d[i]=inf;
    d[S]=0;q.push(mp(d[S],S));
    while(!q.empty()){
        pii u=q.top();q.pop();
        if(vis[u.se])continue;
        vis[u.se]=1;
        for(int i=hd[u.se];i;i=e[i].next){
            int v=e[i].v;
            if(d[v]<=d[u.se]+e[i].w)continue;
            d[v]=d[u.se]+e[i].w;
            q.push(mp(d[v],v));
        }
    }
    printf("%d",d[T]);
}
int main(){
#ifdef wsy
    freopen("data.in","r",stdin);
#else
    freopen("altitude.in","r",stdin);
    freopen("altitude.out","w",stdout);
#endif
    r(n);
    S=0;T=n*n+1;
    int x,s,t;
    
    for(int i=1;i<=n;i++)
    r(x),s=i,t=T,
    adde(i,T,x);
    for(int i=2;i<=n;i++)
    for(int j=1;j<=n;j++)
    r(x),
    s=(i-1)*n+j,t=(i-2)*n+j,
    adde(s,t,x);
    for(int i=1;i<=n;i++)
    r(x),
    s=S,t=(n-1)*n+i,
    adde(S,t,x);
    
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n+1;j++){
        r(x);
        if(j==1){
            s=S;t=(i-1)*n+j;
            adde(S,t,x);
        }
        else if(j==n+1){
            s=i*n;t=T;
            adde(i*n,T,x);
        }
        else{
            s=j+(i-1)*n-1;t=j+(i-1)*n;
            adde(s,t,x);
        } 
    }
    
    for(int i=1;i<=n;i++)
    r(x),s=T,t=i,
    adde(s,t,x);
    for(int i=2;i<=n;i++)
    for(int j=1;j<=n;j++)
    r(x),
    s=(i-2)*n+j,t=(i-1)*n+j,
    adde(s,t,x);
    for(int i=1;i<=n;i++)
    r(x),
    s=(n-1)*n+i,t=S,
    adde(S,t,x);
    
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n+1;j++){
        r(x);
        if(j==1){
            s=(i-1)*n+j;t=S;
            adde(S,t,x);
        }
        else if(j==n+1){
            s=T;t=i*n;
            adde(s,t,x);
        }
        else{
            s=j+(i-1)*n;t=j+(i-1)*n-1;
            adde(s,t,x);
        }
    }
    dijkstra();
    return 0;
}

 

 

网络流

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
#define inf 0x3f3f3f3f
#define ll long long
#define N 1705
using namespace std;
int n,S,T,tot,mp[N][N],cur[N],vis[N],d[N],hd[N];
int dx[]={0,0,1,-1};
int dy[]={1,-1,0,0};
bool check(int x,int y){return x<=n&&y<=n&&x>=1&&y>=1;}
struct edge{int v,cap,next;}e[N*N*4];
void adde(int u,int v,int c){
    e[tot].v=v;
    e[tot].next=hd[u];
    e[tot].cap=c;
    hd[u]=tot++;
}
bool bfs(){
    queue<int>q;
    memset(vis,0,sizeof(vis));
    q.push(S);d[S]=0;vis[S]=1;
    while(!q.empty()){
        int u=q.front();q.pop();
        for(int i=hd[u];~i;i=e[i].next){
            int v=e[i].v;
            if(e[i].cap&&!vis[v]){
                vis[v]=1;
                d[v]=d[u]+1;
                q.push(v);
            }
        }
    }
    return vis[T];
}

int dfs(int u,int a){
    if(u==T||!a)return a;
    int fl=0,f;
    for(int &i=cur[u];~i;i=e[i].next){
        int v=e[i].v;
        if(e[i].cap&&d[v]==d[u]+1&&(f=dfs(v,min(a,e[i].cap)))){
            fl+=f;a-=f;
            e[i].cap-=f;
            e[i^1].cap+=f;
            if(!a)break;
        }
    }
    return fl;
}
int main(){
#ifdef wsy
    freopen("data.in","r",stdin);
#else
    freopen("altitude.in","r",stdin);
    freopen("altitude.out","w",stdout);
#endif
    scanf("%d",&n);
    memset(hd,-1,sizeof(hd));
    for(int i=1;i<=n+1;i++){
        for(int j=1;j<=n;j++){
            int cnt=(i-1)*(n+1)+j; 
            scanf("%d",&mp[cnt][cnt+1]);
        }
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n+1;j++){
            int cnt=(i-1)*(n+1)+j;
            scanf("%d",&mp[cnt][cnt+n+1]);
        }
    }
    for(int i=1;i<=n+1;i++){
        for(int j=1;j<=n;j++){
            int cnt=(i-1)*(n+1)+j; 
            scanf("%d",&mp[cnt+1][cnt]);
        }
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n+1;j++){
            int cnt=(i-1)*(n+1)+j;
            scanf("%d",&mp[cnt+n+1][cnt]);
        }
    }
    n++;S=0;T=n*n+1;
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++){
        int id=(i-1)*n+j;
        if(id==n*n)continue;
        for(int o=0;o<4;o++){
            int nx=i+dx[o];
            int ny=j+dy[o];
            if(!check(nx,ny))continue;
            int nd=(nx-1)*n+ny;
            //if(nd==1)continue;
            adde(id,nd,mp[id][nd]);
            adde(nd,id,0);
        //    printf("%d %d %d\n",id,nd,mp[id][nd]);
        }
    }
    adde(S,1,inf);adde(1,S,0);
//    printf("%d %d %d\n",S,1,inf);
    adde(n*n,T,inf);adde(T,n*n,0);
//    printf("%d %d %d\n",n*n,T,inf);
    int flow=0;
    while(bfs()){
        for(int i=S;i<=T;i++)cur[i]=hd[i];
        flow+=dfs(S,inf);
    }
    printf("%d",flow);
    return 0;
}