码迷,mamicode.com
首页 > 其他好文 > 详细

(模板)hdoj5977 Garden of Eden(点分治)

时间:2019-08-24 13:20:09      阅读:79      评论:0      收藏:0      [点我收藏+]

标签:problem   return   ems   min   oid   roo   void   printf   names   

题目链接:https://vjudge.net/problem/HDU-5977

题意:给一颗树,每个结点上有一个权值a[i],a[i]<=10,求有多少条路径满足这条路径上所有权值的结点都出现了。

思路:

  首先利用二进制的思想,将a[i]转化为1<<(a[i]-1)。我们在子树中,计算出结点到重心的路径,用二进制表示,比如011表示该路径中权值3没有出现、权值1和2出现。因为k最大为10,那么我们在计算结果时把所有可能枚举一遍,也就1024,如果枚举的i和当前路径取或后=(1<<k)-1,那么该路径满足要求,加上即可。具体实现时用桶记录信息,mine[i]表示权值为i的路径的个数。

  另外,题目规定不同的路径仅当起点终点均不同,所以(1,2)和(2,1)是两个合法解,我的处理是先getdis一遍得到桶的信息,处理子结点,现将该子结点的子树的信息清除掉,即change函数,dfs之后再恢复回来。这种处理很重要,所以把这题当作模板记录一下。

  不得不说,写点分治时要非常细心,我总是半小时代码,1小时改bug,老是一些简单错误。

AC代码:

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;

typedef long long LL;
const int maxn=50005;
const int inf=0x3f3f3f3f;
struct node{
    int v,nex;
}edge[maxn<<1];

int n,k,cnt,head[maxn],a[maxn],sz[maxn],mson[maxn],Min,size,root;
int vis[maxn],flag;
LL ans,mine[1050];

void adde(int u,int v){
    edge[++cnt].v=v;
    edge[cnt].nex=head[u];
    head[u]=cnt;
}

void getroot(int u,int fa){
    sz[u]=1,mson[u]=0;
    for(int i=head[u];i;i=edge[i].nex){
        int v=edge[i].v;
        if(vis[v]||v==fa) continue;
        getroot(v,u);
        sz[u]+=sz[v];
        mson[u]=max(mson[u],sz[v]);
    }
    mson[u]=max(mson[u],size-sz[u]);
    if(mson[u]<Min) Min=mson[u],root=u;
}

void getdis(int u,int fa,int len){
    ++mine[len];
    for(int i=head[u];i;i=edge[i].nex){
        int v=edge[i].v;
        if(vis[v]||v==fa) continue;
        getdis(v,u,len|a[v]);
    }
}

void change(int u,int fa,int len,int f){
    mine[len]+=f;
    for(int i=head[u];i;i=edge[i].nex){
        int v=edge[i].v;
        if(vis[v]||v==fa) continue;
        change(v,u,len|a[v],f);
    }
}

void dfs(int u,int fa,int len){
    for(int i=0;i<(1<<k);++i){
        if((i|len)!=flag) continue;
        if(!mine[i]) continue;
        ans+=mine[i];
    }
    for(int i=head[u];i;i=edge[i].nex){
        int v=edge[i].v;
        if(vis[v]||v==fa) continue;
        dfs(v,u,len|a[v]);
    }
}

void solve(int u){
    getdis(u,0,a[u]);
    for(int i=0;i<(1<<k);++i){
        if((i|a[u])!=flag) continue;
        if(!mine[i]) continue;
        ans+=mine[i];
    }
    for(int i=head[u];i;i=edge[i].nex){
        int v=edge[i].v;
        if(vis[v]) continue;
        change(v,u,a[u]|a[v],-1);
        dfs(v,u,a[u]|a[v]);    
        change(v,u,a[u]|a[v],1);
    }
    memset(mine,0,sizeof(mine));
}

void fenzhi(int u,int ssize){
    vis[u]=1;
    solve(u);
    for(int i=head[u];i;i=edge[i].nex){
        int v=edge[i].v;
        if(vis[v]) continue;
        Min=inf,root=0;
        size=sz[v]<sz[u]?sz[v]:ssize-sz[u];
        getroot(v,0);
        fenzhi(root,size);
    }
}

int main(){
    while(~scanf("%d%d",&n,&k)){
        cnt=0;
        ans=0;
        flag=(1<<k)-1;
        for(int i=0;i<=n;++i)
            head[i]=vis[i]=0;
        memset(mine,0,sizeof(mine));
        for(int i=1;i<=n;++i){
            scanf("%d",&a[i]);
            a[i]=1<<(a[i]-1);
        }
        for(int i=1;i<n;++i){
            int u,v;
            scanf("%d%d",&u,&v);
            adde(u,v);
            adde(v,u);
        }
        Min=inf,root=0,size=n;
        getroot(1,0);
        fenzhi(root,n);
        printf("%lld\n",ans);
    }
    return 0;
}

 

(模板)hdoj5977 Garden of Eden(点分治)

标签:problem   return   ems   min   oid   roo   void   printf   names   

原文地址:https://www.cnblogs.com/FrankChen831X/p/11404151.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!