一.缩点的概念
缩点,也称为点缩法(Vertex Contraction),是图论中的一种操作,通常用于缩小图的规模,同时保持了图的某些性质。这个操作的目标是将图中的一些节点合并为一个超级节点,同时调整相关边,以便保持图的连通性和其他性质。
具体步骤如下:
-
选择一个要缩点的节点:选择图中的一个节点,将它合并到另一个节点上。
-
合并节点:将选定的节点合并到另一个节点上,形成一个新的超级节点。通常情况下,选择入度或出度较小的节点进行合并,以减小新图的规模。
-
调整边:将与被合并节点相邻的边重新连接到新的超级节点上。注意要避免重复边和自环。
-
重复步骤1~3:继续选择节点进行缩点,直到不满足合并条件为止。
缩点操作通常用于算法设计和图分析中,有时可以用来简化图的复杂性,减少问题的规模。在一些情况下,缩点操作可能会破坏某些图的属性,因此在使用时需要谨慎考虑。此外,缩点操作后的图可能不再是原始问题的精确表示,可能会导致问题的近似解。
二.缩短的作用
把一个环缩为一个超级点,可以由有环图-->DAG,从而更好的解决问题。
总之就是我们不想要环,直接缩为一个点,我们可以更好地解决问题,就就可以使用缩点法。
三.模板题
P3387 【模板】缩点 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)
四.思路
1.求点权之和最大,我们可以想到什么?最小生成树。
2.但这只需要解决一条路径的点权值最大,那可以怎么解决?拓扑+DP。
3.但是...拓扑只能解决DAG,这有环啊!!! 我们把环缩成一个超级点,然后再建一个新图不就行了吗?理论通过,实践开始!
五.实践
(1)tarjan缩点
主函数部分:
scanf("%d%d",&n,&m);for(int i=1;i<=n;i++){scanf("%d",&p[i]);}for(int i=1;i<=m;i++){int u,v;scanf("%d%d",&u,&v);add(u,v);}for(int i=1;i<=n;i++){if(!dfn[i]) tarjan(i);}
tarjan:
void tarjan(int u){dfn[u]=low[u]=++num;sta[++top]=u;ins[u]=1;for(int i=head[u];i;i=edge[i].next){int v=edge[i].v;if(!dfn[v]){tarjan(v);low[u]=min(low[u],low[v]);}else if(ins[v]){low[u]=min(low[u],dfn[v]);}}if(dfn[u]==low[u]){int j=0;while(1){j=sta[top--];ins[j]=0;h[j]=u; //j从此属于u if(j==u) break;p[u]+=p[j]; //点权值合并到第一个点(u点)上 }}
}
(2)重新建图
for(int i=1;i<=m;i++){int u=h[edge[i].u],v=h[edge[i].v];if(u!=v){ //不在一个环 add2(u,v);in[v]++; //入度++,拓扑用 }}
(3)拓扑排序+DP
int topu(){queue<int> q;for(int i=1;i<=n;i++){ if(!in[i] && i==h[i]){q.push(i); //这是这条路径的起点 dp[i]=p[i]; //记得赋值 } }//拓扑基础 while(!q.empty()){int k=q.front(); q.pop();for(int i=head2[k];i;i=ed[i].next){int v=ed[i].v;dp[v]=max(dp[v],dp[k]+p[v]);in[v]--;if(!in[v]) q.push(v);}}//找最大值,不一定n就最大,毕竟不止一条路 int ans=0;for(int i=1;i<=n;i++){ans=max(ans,dp[i]);}return ans;
}
六.参考代码(完整代码)
#include<bits/stdc++.h>
#define maxn 100005
using namespace std;
int n,m;
int p[maxn];
struct Edge{int u,v,next;
}edge[maxn],ed[maxn];
int head[maxn],head2[maxn],cnt,cnt2;
void add(int u,int v){edge[++cnt]=(Edge){u,v,head[u]}; head[u]=cnt;
}
void add2(int u,int v){ed[++cnt2]=(Edge){u,v,head2[u]}; head2[u]=cnt2;
}
int dfn[maxn],low[maxn],num;
int sta[maxn],ins[maxn],top;
int lg,h[maxn]; //环的个数,成员属于哪个环
void tarjan(int u){dfn[u]=low[u]=++num;sta[++top]=u;ins[u]=1;for(int i=head[u];i;i=edge[i].next){int v=edge[i].v;if(!dfn[v]){tarjan(v);low[u]=min(low[u],low[v]);}else if(ins[v]){low[u]=min(low[u],dfn[v]);}}if(dfn[u]==low[u]){int j=0;while(1){j=sta[top--];ins[j]=0;h[j]=u; //j从此属于u if(j==u) break;p[u]+=p[j]; //点权值合并到第一个点(u点)上 }}
}
int in[maxn],dp[maxn];
int topu(){queue<int> q;for(int i=1;i<=n;i++){ if(!in[i] && i==h[i]){q.push(i); //这是这条路径的起点 dp[i]=p[i]; //记得赋值 } }//拓扑基础 while(!q.empty()){int k=q.front(); q.pop();for(int i=head2[k];i;i=ed[i].next){int v=ed[i].v;dp[v]=max(dp[v],dp[k]+p[v]);in[v]--;if(!in[v]) q.push(v);}}//找最大值,不一定n就最大,毕竟不止一条路 int ans=0;for(int i=1;i<=n;i++){ans=max(ans,dp[i]);}return ans;
}
int main(){scanf("%d%d",&n,&m);for(int i=1;i<=n;i++){scanf("%d",&p[i]);}for(int i=1;i<=m;i++){int u,v;scanf("%d%d",&u,&v);add(u,v);}for(int i=1;i<=n;i++){if(!dfn[i]) tarjan(i);}for(int i=1;i<=m;i++){int u=h[edge[i].u],v=h[edge[i].v];if(u!=v){ //不在一个环 add2(u,v);in[v]++; //入度++,拓扑用 }}cout<<topu();return 0;
}