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

定义:

最小支配集:对于图G = (V, E) 来说,最小支配集指的是从 V 中取尽量少的点组成一个集合, 使得 V 中剩余的点都与取出来的点有边相连.也就是说,设 V' 是图的一个支配集,则对于图中的任意一个顶点 u ,要么属于集合 V', 要么与 V' 中的顶点相邻. 在 V' 中除去任何元素后 V' 不再是支配集, 则支配集 V' 是极小支配集.称G 的所有支配集中顶点个数最少的支配集为最小支配集,最小支配集中的顶点个数称为支配数.

最小点覆盖:对于图G = (V, E) 来说,最小点覆盖指的是从 V 中取尽量少的点组成一个集合, 使得 E 中所有边都与取出来的点相连.也就是说设 V' 是图 G 的一个顶点覆盖,则对于图中任意一条边(u, v), 要么 u 属于集合 V', 要么 v 属于集合 V'. 在集合 V' 中除去任何元素后 V' 不再是顶点覆盖, 则 V' 是极小点覆盖. 称 G 的所有顶点覆盖中顶点个数最小的覆盖为最小点覆盖.

最大独立集: 对于图G = (V, E) 来说,最大独立集指的是从 V 中取尽量多的点组成一个集合,使得这些点之间没有边相连.也就是说设 V' 是图 G 的一个独立集,则对于图中任意一条边(u, v), u 和 v 不能同时属于集合 V', 甚至可以 u 和 v 都不属于集合 V'. 在 V' 中添加任何不属于 V' 元素后 V' 不再是独立集,则 V' 是极大独立集.称 G 的所有顶点独立集中顶点个数最多的独立即为最大独立集.

求解:

　　对于任意图 G 来说,最小支配集, 最小点覆盖, 最大独立集问题是不存在多项式时间的解法,即属于NP问题.但是这里所要求的图 G 是一棵树.可以在多项式的时间内给予解决.在这里利用贪心的思想来解决树上的这三个问题.

(1)最小支配集

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

伪代码:

　　第一步:以根节点深度优先遍历整棵树,求出每个点在深度优先遍历序列中的编号和每个点的父节点编号.

　　第二步:按照深度优先遍历的反向顺序检查每个点,如果当前点不属于支配集也不与支配集的点相连,且它的父节点不属于支配集,将其父节点加入到支配集,支配集中点的个数加 1, 标记当前节点, 当前节点的父节点, 当前节点的父节点的父节点,因为这些节点要么属于支配集(当前点的父节点),要么与支配集中的点相连(当前节点 和 当前节点的父节点的父节点).

具体实现:

　　采用链式前向星存储整棵树.整形数组newpos[i] 表示深度优先遍历序列的第 i 个点是哪个点, now 表示当前深度优先遍历序列已经有多少个点了. bool形数组visit[]用于深度优先遍历的判重,整形pre[i]表示点 i 的父节点编号,  bool型数组s[i]如果为 true, 表示第 i 个点被覆盖, bool型数组set[i]如果为 true,表示点 i 属于要求的点的集合.

代码:

#include <bits/stdc++.h>

using namespace std;
const int maxn = 1000;
int pre[maxn];//存储父节点
bool visit[maxn];//DFS标记数组
int newpos[maxn];//遍历序列
int now;
int n, m;

int head[maxn];//链式前向星
struct Node
{
    int to;
    int next;
};
Node edge[maxn];

void DFS(int x)
{
    newpos[now ++] = x;//记录遍历序列
    for(int k = head[x]; k != -1; k = edge[k].next)
    {
        if(!visit[ edge[k].to ])
        {
            visit[ edge[k].to ] = true;
            pre[edge[k].to] = x;//记录父节点
            DFS(edge[k].to);
        }
    }
}

int MDS()
{
    bool s[maxn] = {0};
    bool set[maxn] = {0};
    int ans = 0;
    for(int i = n - 1; i >= 0; i--)  //逆序进行贪心
    {
        int t = newpos[i];
        if(!s[t])   //如果当前点没被覆盖
        {
            if(! set[ pre[t] ])  //当前点的父节点不属于支配集
            {
                set[ pre[t] ] = true;//当前点的父节点加入支配集
                ans ++;  //支配集节点个数加 1
            }
            s[t] = true; //标记当前点已被覆盖
            s[ pre[t] ] = true;// 标记当前点的父节点被覆盖
            s[ pre[ pre[t] ] ] = true;//标记当前点的父节点的父节点被覆盖
        }
    }
    return ans;
}

int main()
{
    /* read Graph message*/ //建图
    memset(visit, false, sizeof(visit));//初始化
    now = 0;
    visit[1] = true;
    pre[1] = 1;
    DFS(1);//从根节点开始寻摘遍历序列
    MDS();
    return 0;
}

 (2)最小点覆盖

贪心策略:同样需要按照反方向的深度优先遍历序列来进行贪心.每检查一个结点,如果当前点和当前点的父节点都不属于顶点覆盖集合,则将父节点加入到顶点覆盖集合,并标记当前节点和其父节点都被覆盖.注意此贪心策略不适用于根节点,所以要把根节点排除在外.

具体实现:实现上基本和求最小支配集差不多,仅仅需要改动贪心部分即可.

代码:

#include <bits/stdc++.h>

using namespace std;
const int maxn = 1000;
int pre[maxn];//存储父节点
bool visit[maxn];//DFS标记数组
int newpos[maxn];//遍历序列
int now;
int n, m;

int head[maxn];//链式前向星
struct Node
{
    int to;
    int next;
};
Node edge[maxn];

void DFS(int x)
{
    newpos[now ++] = x;//记录遍历序列
    for(int k = head[x]; k != -1; k = edge[k].next)
    {
        if(!visit[ edge[k].to ])
        {
            visit[ edge[k].to ] = true;
            pre[edge[k].to] = x;//记录父节点
            DFS(edge[k].to);
        }
    }
}

int MVC()
{
    bool s[maxn] = {0};
    bool set[maxn] = {0};
    int ans = 0;
    for(int i = n - 1; i >= 1; i--)  //逆序进行贪心,排除掉其根节点
    {
        int t = newpos[i];
        if(!s[t] && !s[ pre[t] ])  //如果当前节点和其父节点都不属于顶点覆盖集合
        {
            set[ pre[t] ] = true;//把其父节点加入到顶点覆盖集合
            ans ++; //集合内顶点个数加 1
            s[t] = true;//标记当前节点被覆盖
            s[ pre[t] ] = true;//标记其父节点被覆盖
        }
    }
    return ans;
}

int main()
{
    /* read Graph message*/ //建图
    memset(visit, false, sizeof(visit));//初始化
    now = 0;
    visit[1] = true;
    pre[1] = 1;
    DFS(1);//从第一个根节点开始寻找遍历序列
    MDS();
    return 0;
}

(3)最大独立集

贪心策略:同样和以上两个贪心问题的贪心方法差不多,需要反向遍历DFS的遍历序列,检查每一个点,如果当前节点没有被覆盖,则将当前节点加入独立集,并标记当前点和其父节点都被覆盖.

代码:

#include <bits/stdc++.h>

using namespace std;
const int maxn = 1000;
int pre[maxn];//存储父节点
bool visit[maxn];//DFS标记数组
int newpos[maxn];//遍历序列
int now;
int n, m;

int head[maxn];//链式前向星
struct Node
{
    int to;
    int next;
};
Node edge[maxn];

void DFS(int x)
{
    newpos[now ++] = x;//记录遍历序列
    for(int k = head[x]; k != -1; k = edge[k].next)
    {
        if(!visit[ edge[k].to ])
        {
            visit[ edge[k].to ] = true;
            pre[edge[k].to] = x;//记录父节点
            DFS(edge[k].to);
        }
    }
}

int MIS()
{
    bool s[maxn] = {0};
    bool set[maxn] = {0};
    int ans = 0;
    for(int i = n - 1; i >= 0; i--)  //按照DFS遍历序列的逆序进行贪心
    {
        int t = newpos[i];
        if(!s[t])  //如果当前节点没有被覆盖
        {
            set[t] = true;//把当前节点加入到独立集
            ans ++;//独立集中点的个数加 1
            s[t] = true;//标记当前点已经被覆盖
            s[ pre[t] ] = true;//标记当前点的父节点已经被覆盖
        }
    }
    return ans;
}

int main()
{
    /* read Graph message*/ //建图
    memset(visit, false, sizeof(visit));//初始化
    now = 0;
    visit[1] = true;
    pre[1] = 1;
    DFS(1);//从第一个根节点开始寻找遍历序列
    MDS();
    return 0;
}

 

参考书籍<<图论及应用>>以及网上的其他资料.