启发式合并(Peterson Polyglot ,cf 778C)

又抄了一波答案。。

先说下自己对树上dfs的时间复杂度的理解吧。

在很久很久以前,树形dp也不会,dfs序也不会,时间复杂度也不会分析,所以当时一遇到树上的问题就只会dfs里面套个dfs2,然后O(n^2)超时,从此以后便完全否定了这种方法,认为这样做一定妥妥的O(n^2),但是其实并不是这样,有些细节还挺有意思的。


dfs里面套个dfs2,大概有两种吧。。

第一种dfs2是可以往上dfs2的,这种就是铁定O(n^2)了,在最开始学习的时候就是这样瞎搞的,一般解决办法都是两遍dfs树形dp,第一遍往下dp,第二遍往上dp。

第二种dfs2是只可以往下dfs2的,这种在最优情况下是O(n),比如星,最坏情况下是O(n^2),比如链。这种方法一般也是不行的,出题人一般都会出一些长链的数据来卡你。解决办法同上。

值得注意的是出现了最优情况和最坏情况,你要知道哪怕是快速排序都有最优与最坏的时间复杂度。所以也不要太看低上面的解法了。

现在又发现了一种神奇的暴力,启发式合并,跟第二种dfs2有点像,这也是为什么不应该看低上述算法的原因,深入挖掘下还是有可能有所变化的。

启发式合并,听起来十分高大上,但是其实说起来挺简单的,就是要合并两个集合,那我们就不妨把小的往大的合并。

虽然说起来挺简单,但是分析起其中的时间复杂度,还是很有趣的。

我也不知道该怎么证明,只能把具体怎么做的,这种算法适用情况的特点,以及大概为什么对说一下把。


具体怎么做的:

就是说如果要合并k个集合(子树),那我们就把其他所有的集合(子树)往最大的那个集合(子树)里边合并就好了。

这样的话,我们就少遍历了一个最大的集合(子树)。


这种算法适用情况的特点:

大概就是不明确需要遍历所有子树的每个节点,但又需要获得某些子树或节点的信息吧。(比如合并?)


大概为什么对:

http://blog.csdn.net/zjznku/article/details/60605971

别人家的证明,我复制下来。


题意:给定一棵trie树,可以删除一层边,再将父边被删且父亲相同的结点对应的子树合并得到一棵新trie树,求新trie树的最小结点数。

分析:启发式合并。在合并结点u的子树时,选择将小子树合并到大子树里,这样总的合并的时间复杂度是O(nlgn)的。

        证明:合并的耗时来自于对小子树的遍历。设全体小子树的遍历总量为T,考虑每个结点u对T的贡献。设结点u可以作为小子树的第i1层、i2层、...、it层(i1<i2<...<it),则结点u对T的贡献为t。由合并的方式(小子树合并到大子树)知对应的子树tree1,tree2,...,treet有size(tree1)<size(tree2)/2<...<size(treet)/2^t。所以t不超过lgn。所以T为nlgn级别的量。

        程序实现上可以先正着将结点u的小子树合并到u的最大子树里,再倒着删除,然后考虑u的子结点。


这个证明挺好的,我看明白了,就讲一些自己卡住的细节。

为什么size(tree1)<size(tree2)/2呢?(i1<i2)

首先由于(i1<i2)所以 tree1是tree2的子树。

如果tree1是tree2最大的那个子树,那么此时u就不做贡献。(最大子树不参与合并)

否则一定满足size(tree1)<size(tree2)/2,而且此时u做出了贡献。


抄别人的代码真是爽,可惜没什么乱用,就当学习启发式合并了吧。

#include<bits/stdc++.h>
using namespace std;
const int maxn = 600010;
const int maxst = 26;

int n,N;
int ch[maxn][maxst];
int ans[maxn];
int cnt;

int idx(char c)
{
    return c-'a';
}

void Read()
{
    int u,v;
    char c;
    scanf("%d",&n);
    for(int i=1;i<n;i++)
    {
        scanf("%d %d %c",&u,&v,&c);
        ch[u][idx(c)]=v;
    }
}

int merge(int x,int y)
{
    if(!x) return y;
    if(!y) return x;
    int rt=++N;cnt++;
    for(int i=0;i<maxst;i++) ch[rt][i]=merge(ch[x][i],ch[y][i]);
    return rt;
}

void dfs(int u,int d)
{
    N=n;
    int rt=++N;
    cnt=0;
    for(int i=0;i<maxst;i++) if(ch[u][i]) rt=merge(rt,ch[u][i]);
    ans[d]+=cnt;
    for(int i=0;i<maxst;i++) if(ch[u][i]) dfs(ch[u][i],d+1);
}

int main()
{
    Read();
    N=n;
    dfs(1,1);
    int p=-1;
    for(int i=1;i<=n;i++) if(p==-1||ans[p]<ans[i]) p=i;
    printf("%d\n%d\n",n-ans[p],p);
}



  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值