资讯

历史

科技

环境与自然

成长

游戏

财经

文学与艺术

美食

健康

家居

文化

情感

汽车

三农

军事

旅行

运动

教育

生活

星座命理

平面图最短路与对偶图网络流

创作时间:

作者:

@小白创作中心

平面图最短路与对偶图网络流

引用

来源

https://www.cnblogs.com/spdarkle/p/18464400

平面图的最小割与对偶图的最短路之间存在着密切的关系。这一理论在算法竞赛和计算机科学领域有着广泛的应用。本文将详细介绍平面图、对偶图的概念，并通过两个具体的OI题目来说明这一理论的应用。

平面图

平面图，指可以通过画图方式将使得边两两不相交的图。（无向图）

例如：

事实上是：

一些概念：

设(G)是平面图，由(G)的边将所在的平面划分成若干个区域，每个区域称为(G)的一个面，其中面积无限的面称为无界面或外部面，面积有限的称为内部面。包围每个面的所有边组成的回路称为该面的边界，边界的长度称为该面的次数（边数）。显然有次数之和为(2|E|)
阶数也就是(|V|)，边数也就是(|E|)。
若在简单平面图(G)的任意不相邻顶点间添加边，所得图为非平面图，称(G)为极大平面图。
一个(|V|=5)的完全图不是平面图，一个(|L|=|R|=3)的完全二分图不是平面图。一个图是平面图当且仅当不存在可以收缩到上述两种情况的子图
典型的平面图是网格图。
(|V|-|E|+C=2)，其中(C)是面数，欧拉公式，以下也采用(|V|,|E|,C)表示阶数边数面数。
S-T 平面图：图中的一个点为源点(s)，另外一个点为汇点(t)，且(s)和(t)都在图中的无界面的边界上

对偶图

平面图的对偶图也是一个平面图，具体构造方法如下：

将每个面视作一个点

按照如下方法构造：

(\forall (u,v,w)\in E)
该边同时属于面(f_1,f_2)，则连边((f_1,f_2,w))
该边仅属于面(f)，则连边((f,f,w))。也即自环。
注意都是无向边

性质：

对偶图的每一个环对应原图的一个割（一一对应）
对偶图与原图边数相同，对偶图点数等于原图面数，可得对偶图面数等于原图点数

建立对偶图与 S-T 平面图关系

容易想到一种构造方法：

连接(S-T)，边权为零，得到一个附加面，对这样的一个平面图建立对偶图

称原本附加面对应点为(s^)，原本无界面对应点为(t^)
我们让新图的边权等价为容量，则我们所求(s\to t)的最小割等价于该图包含(s^,t^)的一个边权和最小的环，由于(s^,t^)之间边权为零，因此这个边权和等价于删去((s^,t^))后(s^\to t^)的最短路。
我们使用 Dijkstra 算法求解最短路即可，复杂度得到大大降低。

实际应用

实际应用中这个平面图常常抽象为网格图，从((0,0)\to (n,m))。然后上右两个边界连向(S^)，下左两个边界连向(T^)，不过由于对称性反过来也没什么。

正向：CSP-S2021 交通规划

简要题意：给定一张网格图，每次给出(k)个边界上的交点作为关键点，每个点有颜色（黑白两色之一），你需要将非关键点染为黑白两色，使得(\sum_{(u,v,w),col_u\neq col_v}w)最小。求出这个最小值。

容易发现题目即为求解一个割掉边的方案使得所有黑色关键点与白色关键点不连通且割的边边权和最小。

很容易想到暴力网络流建图，可以获得 65 pts。

然后考虑(k=2)的部分分，此时(s,t)确定，可以利用转对偶图的方法通过。

两者综合可以获得相当可观的分数。

当一般(k)的时候，直接建立超级原点和超级汇点再转对偶图是错误的，它不一定是一个平面图（因为都在无界面上，而且可能黑白点有交叉，这就是非平面图了）。

但是我们受到了一些启发，我们将相邻的同色特殊点合并在一起，然后将合并后两个相邻非同色特殊点的部分依次划分为(s_i,t_i)（也即(0\to 1:s)，(1\to 0:t)），则显然最优方案一定是非同色块的两两配对（分类讨论即可证明），割掉这两个块之间的最短路即可。

我们可以(O(nmk\log nm))地跑最短路求出各个块两两间的最短路，也就是割掉的代价，那么问题转化为：有一个长度为(2n)的环，环上(i,j)连边的代价为(d_{i,j})，你需要对于每个奇数点与偶数点匹配（两两匹配），要求线段不交（相交定然不优），求最小匹配代价。

这个套路很多，可以将其视作括号序列等，利用破环为链后的区间DP可以轻松求解。

(f_{l,r}=\min(f_{l+1,r-1}+d_{l,r},f_{l,k}+f_{k+1,r}))

最终复杂度(O(nm(\sum k)\log nm+\sum k^3))

#include<bits/stdc++.h>
using namespace std;
#define pr pair<int,int>
#define mk make_pair
#define N 1005000
int id[505][505],num,n,m,dis[N],vis[N],f[505][505],wl[505][505],wh[505][505],ope[N],d[505][505],Id[N],col[N],cnt,nid[N];//操作涉及点编号
pr waiwei[N];//外围点对应格点
vector<pr >e[N];
void dijkstra(int s){
    for(int i=1;i<=num;++i)dis[i]=0x3f3f3f3f,vis[i]=0;
    priority_queue<pr >q;
    dis[s]=0,q.push(mk(0,s));
    while(!q.empty()){
        int u=q.top().second;q.pop();
        if(vis[u])continue;
        vis[u]=1;
        for(auto to:e[u]){
            int v=to.first,w=to.second;
            if(dis[v]>dis[u]+w){
                dis[v]=dis[u]+w;
                // cout<<u<<" "<<v<<' '<<w<<" "<<dis[v]<<"\n";
                q.push(mk(-dis[v],v));
            }
        }
    }
}
int dp(int len){
    for(int i=1;i<=len;++i)Id[i+len]=Id[i]=i;
    for(int i=0;i<=len<<1;++i)for(int j=0;j<=len<<1;++j)f[i][j]=0x3f3f3f3f;
    for(int i=1;i<=len<<1;++i)f[i+1][i]=0;
    for(int x=2;x<=len;x+=2)for(int l=1,r=x;r<=len<<1;++r,++l){
        f[l][r]=f[l+1][r-1]+d[Id[l]][Id[r]];
        for(int k=l+1;k<r-1;k+=2)f[l][r]=min(f[l][r],f[l][k]+f[k+1][r]);
    }
    int ans=0x3f3f3f3f;
    for(int l=1;l<len;++l)ans=min(ans,f[l][l+len-1]);
    return ans;
}
void add(int u,int v,int w){
    e[u].push_back(mk(v,w));e[v].push_back(mk(u,w));
}
int cost[N];
void build(){
    for(int i=1;i<=m;++i)waiwei[++num]=mk(1,i);
    for(int i=1;i<=n;++i)waiwei[++num]=mk(i,m);
    for(int i=m;i;--i)waiwei[++num]=mk(n,i);
    for(int i=n;i;--i)waiwei[++num]=mk(i,1);
    for(int i=1;i<=(m+n)*2;++i)cost[i]=0;
    for(int i=1;i<n;++i)for(int j=1;j<m;++j)id[i][j]=++num;
    for(int i=1;i<n;++i)for(int j=1;j<=m;++j){
        if(j==1)add(2*(m+n)-i+1,id[i][j],wl[i][j]);
        else if(j==m)add(m+i+1,id[i][j-1],wl[i][j]);
        else add(id[i][j],id[i][j-1],wl[i][j]);
    }
    for(int i=1;i<=n;++i)for(int j=1;j<m;++j){
        if(i==1)add(j+1,id[i][j],wh[i][j]);
        else if(i==n)add(2*m+n-j+1,id[i-1][j],wh[i][j]);
        else add(id[i][j],id[i-1][j],wh[i][j]);
    }
}
vector<int>S,T;
void sol(){
    cin>>cnt;int flg=1;
    for(int i=1;i<=cnt;++i){
        int w;cin>>w;cin>>ope[i];cin>>col[ope[i]];cost[ope[i]]=w;
    }
    for(int i=1;i<=(n+m)<<1;++i){
        int j=i+1;if(i==(n+m)<<1)j=1;
        add(i,j,cost[i]);cost[i]=0;
        // cout<<"Ad "<<i<<' '<<j<<' '<<cost[i]<<"\n";
    }
    int c=col[ope[1]];
    for(int i=1;i<=cnt;++i)flg&=(col[ope[i]]==c);
    if(flg){
        cout<<"0\n";
    }
    else {
        sort(ope+1,ope+cnt+1);
        for(int i=1;i<=cnt;++i)for(int j=1;j<=cnt;++j)d[i][j]=0x3f3f3f3f;
        for(int i=1;i<=cnt;++i)ope[i+cnt]=ope[i];int st=0;
        for(int i=1;i<=cnt;++i)if(col[ope[i]]!=col[ope[i+1]])st=i+1;
        int len=0;
        S.clear();T.clear();
        for(int i=st;i<st+cnt;++i){
            if(col[ope[i]]!=col[ope[i-1]]){
                if(col[ope[i]])S.push_back(ope[i]),++len,nid[ope[i]]=len;
                else T.push_back(ope[i]),++len,nid[ope[i]]=len;
            }
        }
        for(auto x:S){
            dijkstra(x);
            for(auto y:T)d[nid[x]][nid[y]]=d[nid[y]][nid[x]]=dis[y];
        }
        cout<<dp(len)<<"\n";
    }
    for(int i=1;i<=(n+m<<1);++i){
        e[i].pop_back();e[i].pop_back();
    }
}
int main(){
    freopen("traffic.in","r",stdin);
    freopen("traffic.out","w",stdout);
    ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
    cin>>n>>m;int T;cin>>T;
    for(int i=1;i<n;++i)for(int j=1;j<=m;++j)cin>>wl[i][j];
    for(int i=1;i<=n;++i)for(int j=1;j<m;++j)cin>>wh[i][j];
    build();
    while(T--){
        sol();
    }
}