博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
树的最小支配集,最小点覆盖,最大独立集两种算法
阅读量:5251 次
发布时间:2019-06-14

本文共 4017 字,大约阅读时间需要 13 分钟。

1.基本概念

对图G=<V,E>,

最小支配集:从V中取尽量少的点组成一个集合,使得V中剩余的点都与取出来的点有边相连

最小点覆盖:从V中取尽量少的点组成一个集合,使得E中所有边都与取出来的点相连

最大独立集:从V中取尽量多的点组成一个集合,使得这些点之间没有边相连

 

2.贪心法求树的最小支配集,最小点覆盖,最大独立集模板

基本算法:

以最小支配集为例,首先选择一点为根,按照深度优先遍历得到遍历序列,按照所得序列的反向序列的顺序进行贪心,对于一个既不属于支配集也不与支配集中的点相连的点来说,如果他的父节点不属于支配集,将其父节点加入支配集。

struct Edge{    int v,next;}G[M];int fa[N];int vis[N];int pos[N],head[N];int now;int n,m;void DFS(int u){    pos[now++] = u;    for(int i=head[u];i!=-1;i=G[i].next)    {        int v = G[i].v;        if(!vis[v])        {            vis[v] = 1;            fa[v] = u;            DFS(v);        }     }}
int MDS(){    int s[N] = {
0}; int set[N] = {
0}; int ans = 0; for(int i=now-1;i>=0;i--) { int t = pos[i]; if(!s[t]) { if(!set[fa[t]]) { set[fa[t]] = 1; ans++; } s[t] = 1; s[fa[t]] = 1; s[fa[fa[t]]] = 1; } } return ans;}
最小支配集MDS
int MPC(){    int s[N] = {
0}; int set[N] = {
0}; int ans = 0; for(int i=now-1;i>=0;i--) { int t = pos[i]; if(!s[t] && !s[fa[t]]) { set[fa[t]] = 1; ans++; s[t] = 1; s[fa[t]] = 1; } } return ans;}
最小点覆盖MPC
int MIS(){    int s[N] = {
0}; int set[N] = {
0}; int ans = 0; for(int i=now-1;i>=0;i--) { int t = pos[i]; if(!s[t]) { set[t] = 1; ans++; s[t] = 1; s[fa[t]] = 1; } } return ans;}
最大独立集MIS
int main(){    //读入图信息    memset(vis,0,sizeof(vis));    now = 0;    vis[1] = 1;    fa[1] = 1;    DFS(1);    printf("%d\n",MIS()); //MDS | MPC    return 0;}

 

3.树形DP算法

(1)最小支配集:

定义:

dp[i][0]: 点i属于支配集,并且以点i为根的子树都被覆盖了的情况下支配集中包含的最少点数

dp[i][1]: 点i不属于支配集,且以i为根的子树都被覆盖,且i被其中不少于1个子节点覆盖的情况下支配集包含的最少点数

dp[i][2]: 点i不属于支配集,且以i为根的子树都被覆盖,且i没被子节点覆盖的情况下支配集包含的最少点数

则有:

dp[i][0] = SUM{min(dp[u][0],dp[u][1],dp[u][2])} (fa[u] = i)

dp[i][1]:

if(i没有子节点)  dp[i][1] = INF

else   dp[i][1] = SUM{min(dp[u][0],dp[u][1])} + inc  (fa[u] = i)

inc有:

if(上面式子中SUM{min(dp[u][0],dp[u][1])}包含某个dp[u][0])  inc = 0

else  inc = MIN(dp[u][0]-dp[u][1]) (fa[u] = i)

dp[i][2] = SUM(dp[u][1]) (fa[u] = i)

void MDS_DP(int u,int fa){    dp[u][2] = 0;    dp[u][0] = 1;    int s = 0;    int sum = 0;    int inc = Mod;    for(int i=head[u];i!=-1;i=G[i].next)    {        int v = G[i].v;        if(v == fa)            continue;        MDS_DP(v,u);        dp[u][0] += min(dp[v][0],min(dp[v][1],dp[v][2]));        if(dp[v][0] <= dp[v][1])        {            sum += dp[v][0];            s = 1;        }        else        {            sum += dp[v][1];            inc = min(inc,dp[v][0]-dp[v][1]);        }        if(dp[v][1] != Mod && dp[u][2] != Mod)            dp[u][2] += dp[v][1];        else            dp[u][2] = Mod;        if(inc == Mod && !s)  //i没有子节点            dp[u][1] = Mod;        else        {            dp[u][1] = sum;            if(!s)                dp[u][1] += inc;        }    }}
最小支配集MDS_DP

 

(2)最小点覆盖

定义:

dp[i][0]: 点i属于点覆盖,并且以点i为根的自述中所连接的边都被覆盖的情况下点覆盖集所包含的最少点数

dp[i][1]: 点i不属于点覆盖,并且以点i为根的自述中所连接的边都被覆盖的情况下点覆盖集所包含的最少点数

则有:

dp[i][0] = SUM{min(dp[u][0],dp[u][1])}+1 (fa[u] = i)

dp[i][1] = SUM(dp[u][0])  (fa[u] = i)

void MPC_DP(int u,int fa){    dp[u][0] = 1;    dp[u][1] = 0;    for(int i=head[u];i!=-1;i=G[i].next)    {        int v = G[i].next;        if(v = fa)            continue;        MPC_DP(v,u);        dp[u][0] += min(dp[v][0],dp[v][1]);        dp[u][1] += dp[v][0];    }}
最小点覆盖MPC_DP

 

(3)最大独立集

dp[i][0]: 点i属于独立集的情况下,最大独立集中点的个数

dp[i][1]: 点i属于独立集的情况下,最大独立集中点的个数

则有:

dp[i][0] = SUM(dp[u][1])+1  (fa[u] = i)

dp[i][1] = SUM{max(dp[u][0],dp[u][1])}  (fa[u] = i)

void MIS_DP(int u,int fa){    dp[u][0] = 1;    dp[u][1] = 0;    for(int i=head[u];i!=-1;i=G[i].next)    {        int v = G[i].v;        if(v == fa)            continue;        MIS_DP(v,u);        dp[u][0] += dp[v][1];        dp[u][1] += max(dp[v][0],dp[v][1]);    }}
最大独立集MIS_DP

 

例题:POJ 3398 

 

转载于:https://www.cnblogs.com/whatbeg/p/3776612.html

你可能感兴趣的文章
getopt_long
查看>>
TensorFlow MNIST CNN 代码
查看>>
javascript之Style物
查看>>
JSON跨域解决方案收集
查看>>
图的深度优先遍历
查看>>
C# 之 提高WebService性能大数据量网络传输处理
查看>>
[bzoj1004] [HNOI2008] Cards
查看>>
原生HttpClient详细使用示例
查看>>
几道面试题
查看>>
Factory Design Pattern
查看>>
python中贪婪与非贪婪
查看>>
guava API整理
查看>>
无锁编程笔记
查看>>
jquery mobile
查看>>
如何在vue单页应用中使用百度地图
查看>>
Springboot使用步骤
查看>>
Spring属性注入
查看>>
Springboot-配置文件
查看>>
Springboot-日志框架
查看>>
P1192-台阶问题
查看>>