倍增和LCA

倍增和LCA

现在终于可以写这篇文章了，拖了很久没有学的倍增在去NOIP2015打酱油之前终于初步学习完。再次特别声明十分感谢fye的帮助，生命中第一篇倍增和LCA就是以她的代码为模板的，并且应该会一直延续下去。。。= =
【倍增是什么】
顾名思义，倍增就是根据已经得到的信息，将考虑的范围扩大一倍，从而加速操作的一种思想。
【有关倍增的算法】
使用了倍增思想的算法有很多，包括归并排序、快速幂、基于ST表的RMQ算法和树上倍增找LCA等，还有FFT、后缀数组等高级算法（wok…），而我们只是初级入门，学习了树上倍增找LCA。
【树上倍增找LCA】
算法的基本思路是：求两个点的LCA，先把深度较深的点移到和另一个点相同的深度，然后一起向上移动，直到重合，就找到了LCA，倍增其实是加速区间操作的一种手段，因为相当于把每个数进行二进制分解，倍增的复杂度在logN级别。
在这里不做详细讲解，只给出普通查找和使用倍增优化的代码模板，详情可以百度“倍增算法 高天宇”，SLYZGTY学长为你详细讲解。= =
【普通查找LCA】

#include
#include
#include
#include
using namespace std;
int n,a,b,x,y,ans;
int h[105],f[105][105],fa[105];
bool bb[105];
void build(int t,int dep){
    h[t]=dep;
    int r;
    for (r=1;r<=n;++r)
      if (f[t][r]&&!bb[r]){
        bb[r]=true;
        fa[r]=t;
        build(r,dep+1);
      }
}
int lca(int x,int y){
    if (h[x]int k=h[x]-h[y];
    for (int i=1;i<=k;++i)
      x=fa[x];
    while (fa[x]!=fa[y]){
        x=fa[x];y=fa[y];
    }
    return fa[x];
}
int main(){
    scanf("%d%d%d",&n,&a,&b);
    for (int i=1;iscanf("%d%d",&x,&y);
        f[x][y]=f[y][x]=1;
    }
    bb[1]=true;
    build(1,1);
    ans=lca(a,b);
    printf("%d",ans);
    return 0;
}

【倍增优化】

#include
#include
#include
#include
#define sz 8
using namespace std;
int n,a,b,i,t,ans,x,y;
int v[105],next[105],point[105],h[105],mi[sz],f[105][sz],root[105];
inline void add(int x,int y,int i){
    v[i]=y; next[i]=point[x]; point[x]=i;
}
inline void build(int x,int fa,int dep){
    h[x]=dep;
    for (i=1;i
        if (h[x]-mi[i]<1) break;
        f[x][i]=f[f[x][i-1]][i-1];
    }
    for (int i=point[x];i;i=next[i])
      if (v[i]!=fa){
        f[v[i]][0]=x;
        build(v[i],x,dep+1);
      }     
}
inline int lca(int x,int y){
    if (h[x]
    for (int i=0;i
      if ((k>>i)&1)
        x=f[x][i];
    if (x==y) return x;
    for (i=sz-1;i>=0;--i)
      if (f[x][i]!=f[y][i]){x=f[x][i];y=f[y][i];}
    return f[x][0];
}
int main()
{
    mi[0]=1; for (i=1;i
    scanf("%d%d%d",&n,&a,&b);
    for (i=1;i
        scanf("%d%d",&x,&y);
        add(x,y,i<<1);
        add(y,x,i<<1|1);
    }
    build(1,0,1);
    ans=lca(a,b);
    printf("%d",ans);
    return 0;
}

其实上面两个程序的数据是很弱的，暴力查找LCA都可以承受，只是为了练练手。。。
使用倍增优化用到了next数组，因为如果使用倍增时数据很大，其他存储方式无法存储，next数组是优先考虑的存储形式。
f[x][i]表示结点x向上蹦2^i个点到达的结点是哪个，从小到大枚举i，转移方程为：f[x][i]=f[f[x-1][i-1]][i-1]。
如果在预处理时和查找LCA时记录一些别的信息，我们还可以处理很多维护和查询问题。
算法精妙，用心体会。。。