BZOJ 4316 小C的独立集 仙人掌DP

25 篇文章 1 订阅

题目大意:给定一棵仙人掌,求最大点独立集

老办法,每个环新建一个点,连接环上所有点,然后把环上边都拆掉,就变成了一棵树
然后跑树形DP,对于一个点如果是普通节点就正常DP,如果是环上节点参考BZOJ1040

#include <vector>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 120100
using namespace std;

int n,m;

namespace Tree{
    int n;
    struct abcd{
        int to,next;
    }table[M];
    int head[M],tot=1;
    int f[M][2];
    void Add(int x,int y)
    {
        table[++tot].to=y;
        table[tot].next=head[x];
        head[x]=tot;
    }
    void Tree_DP(int x,int from)
    {
        int i;
        if(x<=::n)
        {
            f[x][0]=0;f[x][1]=1;
            for(i=head[x];i;i=table[i].next)
                if(table[i].to!=from)
                {
                    Tree_DP(table[i].to,x);
                    if(table[i].to>::n)
                        continue;
                    f[x][0]+=max(f[table[i].to][0],f[table[i].to][1]);
                    f[x][1]+=f[table[i].to][0];
                }
        }
        else
        {
            static int g[M][2];int top=0;
            for(i=head[x];i;i=table[i].next)
                if(table[i].to!=from)
                    Tree_DP(table[i].to,x);
            for(i=head[x];i;i=table[i].next)
            {
                g[++top][0]=f[table[i].to][0];
                g[  top][1]=f[table[i].to][1];
            }
            for(i=top-1;i;i--)
            {
                g[i][0]+=max(g[i+1][0],g[i+1][1]);
                g[i][1]+=g[i+1][0];
            }
            f[from][0]=g[1][0];
            for(i=1;i<top;i++)
            {
                g[i][0]-=max(g[i+1][0],g[i+1][1]);
                g[i][1]-=g[i+1][0];
            }
            g[top][1]=0xefefefef;
            for(i=top-1;i;i--)
            {
                g[i][0]+=max(g[i+1][0],g[i+1][1]);
                g[i][1]+=g[i+1][0];
            }
            f[from][1]=g[1][1];
        }
    }
}

namespace Cactus{
    struct abcd{
        int to,next;
    }table[M];
    int head[M],tot=1;
    int fa[M],dfn[M],T;
    bool on_ring[M];
    void Add(int x,int y)
    {
        table[++tot].to=y;
        table[tot].next=head[x];
        head[x]=tot;
    }
    void DFS(int x,int from)
    {
        int i;
        dfn[x]=++T;
        for(i=head[x];i;i=table[i].next)
            if(i^from^1)
            {
                if(!dfn[table[i].to])
                {
                    fa[table[i].to]=x;
                    on_ring[x]=false;
                    DFS(table[i].to,i);
                    if(!on_ring[x])
                    {
                        Tree::Add(x,table[i].to);
                        Tree::Add(table[i].to,x);
                    }
                }
                else
                {
                    if(dfn[table[i].to]>dfn[x])
                        continue;
                    int temp=x;++Tree::n;
                    while(1)
                    {
                        Tree::Add(Tree::n,temp);
                        Tree::Add(temp,Tree::n);
                        on_ring[temp]=true;
                        if(temp==table[i].to)
                            break;
                        temp=fa[temp];
                    }
                }
            }

    }
}

int main()
{
    int i,x,y;
    cin>>n>>m;Tree::n=n;
    for(i=1;i<=m;i++)
    {
        scanf("%d%d",&x,&y);
        Cactus::Add(x,y);
        Cactus::Add(y,x);
    }
    Cactus::DFS(1,0);
    Tree::Tree_DP(1,0);
    cout<<max(Tree::f[1][0],Tree::f[1][1])<<endl;
}
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值