$My$ $template$(持续更新)

树链剖分：（来源：树的统计）

#include<bits/stdc++.h>
#define rint register int
using namespace std;
inline void read(int &A)
{
    A=;int B=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-')B=-;ch=getchar();}
    while(ch>=''&&ch<=''){A=(A<<)+(A<<)+ch-'';ch=getchar();}
    A=A*B;return ;
}
int n,q,v[],w[],nxt[],first[],tot;
int f[],d[],siz[],son[];
int id[],rk[],cnt,tp[],ans;
char ch[];
struct node{int k,mx;}t[<<];
inline void build_line(int uu,int vv)
{
    v[++tot]=vv,nxt[tot]=first[uu];
    first[uu]=tot;
}
inline void dfs(int x,int fa,int dep)
{
    f[x]=fa,d[x]=dep,siz[x]=;
    for(rint i=first[x];i;i=nxt[i])
    {
        int y=v[i];if(y==fa)continue;
        dfs(y,x,dep+);siz[x]+=siz[y];
        if(siz[y]>siz[son[x]])son[x]=y;
    }
    return ;
}
inline void Dfs(int x,int t)
{
    tp[x]=t,id[x]=++cnt,rk[cnt]=x;
    if(!son[x])return ;Dfs(son[x],t);
    for(rint i=first[x];i;i=nxt[i])
        if(v[i]!=son[x]&&v[i]!=f[x])
            Dfs(v[i],v[i]);
    return ;
}
inline void update(int k)
{
    t[k].k=t[k<<].k+t[k<<|].k;
    t[k].mx=max(t[k<<].mx,t[k<<|].mx);
    return ;
}
inline void build(int k,int l,int r)
{
    if(l==r){t[k].k=t[k].mx=w[rk[l]];return ;}
    int mid=(l+r)>>;
    build(k<<,l,mid);build(k<<|,mid+,r);
    update(k);return ;
}
inline int Qsum(int k,int l,int r,int L,int R)
{
    if(L<=l&&r<=R)return t[k].k;
    int mid=(l+r)>>,res=;
    if(L<=mid)res+=Qsum(k<<,l,mid,L,R);
    if(R>mid)res+=Qsum(k<<|,mid+,r,L,R);
    return res;
}
inline int Qmax(int k,int l,int r,int L,int R)
{
    if(L<=l&&r<=R)return t[k].mx;
    int mid=(l+r)>>,res=-0x7fffffff;
    if(L<=mid)res=max(res,Qmax(k<<,l,mid,L,R));
    if(R>mid)res=max(res,Qmax(k<<|,mid+,r,L,R));
    return res;
}
inline void change(int k,int l,int r,int p,int dat)
{
    if(l==r){t[k].k=t[k].mx=dat;return ;}
    int mid=(l+r)>>;
    if(p<=mid)change(k<<,l,mid,p,dat);
    else change(k<<|,mid+,r,p,dat);
    update(k);return ;
}
inline int get_sum(int x,int y)
{
    int Sum=,fx=tp[x],fy=tp[y];
    while(fx!=fy)
    {
        if(d[fx]>=d[fy])
        {
            Sum+=Qsum(,,n,id[fx],id[x]);
            x=f[fx],fx=tp[x];
        }
        else
        {
            Sum+=Qsum(,,n,id[fy],id[y]);
            y=f[fy],fy=tp[y];
        }
    }
    if(id[x]<=id[y])Sum+=Qsum(,,n,id[x],id[y]);
    else Sum+=Qsum(,,n,id[y],id[x]);
    return Sum;
}
inline int get_max(int x,int y)
{
    int Max=-0x7fffffff,fx=tp[x],fy=tp[y];
    while(fx!=fy)
    {
        if(d[fx]>=d[fy])
        {
            Max=max(Max,Qmax(,,n,id[fx],id[x]));
            x=f[fx],fx=tp[x];
        }
        else
        {
            Max=max(Max,Qmax(,,n,id[fy],id[y]));
            y=f[fy],fy=tp[y];
        }
    }
    if(id[x]<=id[y])Max=max(Max,Qmax(,,n,id[x],id[y]));
    else Max=max(Max,Qmax(,,n,id[y],id[x]));
    return Max;
}
int main()
{
//    freopen("count1.in","r",stdin);
//    freopen("my.out","w",stdout);
    read(n);
    for(rint i=,ST,EN;i<n;++i)
    {
        read(ST),read(EN);
        build_line(ST,EN);
        build_line(EN,ST);
    }
    for(rint i=;i<=n;++i)read(w[i]);
    dfs(,,);Dfs(,);build(,,n);
    read(q);
    for(rint i=,uu,vv;i<=q;++i)
    {
        scanf("%s",ch);
        if(ch[]=='C')
        {
            read(uu),read(vv);
            change(,,n,id[uu],vv);
        }
        else
        {
            ans=;read(uu),read(vv);
            if(ch[]=='M')ans=get_max(uu,vv);
            else ans=get_sum(uu,vv);
            printf("%d\n",ans);
        }
    }
}

树链剖分模板

tarjan点双

inline void tarjan(int x)
{
    dfn[x]=low[x]=++cnt;s.push(x);
    if(x==root && firsty[x]==)
    {
        dcc[++totq].clear();
        dcc[totq].push_back(x);
        return ;
    }
    int flag=;
    for(rint i=firsty[x];i;i=edge[i].nxt)
    {
        int y=edge[i].v;
        if(!dfn[y])
        {
            tarjan(y);low[x]=min(low[x],low[y]);
            if(low[y]>=dfn[x])
            {
                flag++;
                if(x!=root||flag>)cut[x]=true;
                totq++;int z;
                dcc[totq].clear();
                do{
                    z=s.top();s.pop();
                    dcc[totq].push_back(z);
                }while(z!=y);
                dcc[totq].push_back(x);
            }
        }
        else low[x]=min(low[x],dfn[y]);
    }
}

tarjan点双模板

二分图匹配（带注释）

int point(int u)//这个函数的作用是寻找增广路和更新cx，xy数组，如果找到了增广路，函数返回1，找不到，函数返回0。
{
    for(int v=;v<=ny;v++)//依次遍历右边的所有顶点
    {
        if(e[u][v]&&!visited[v])//条件一：左边的u顶点和右边的v顶点有连通边，条件二：右边的v顶点在没有被访问过，这两个条件必须同时满足
        {
            visited[v]=;//将v顶点标记为访问过的
            if(cy[v]==-||point(cy[v]))//条件一：右边的v顶点没有左边对应的匹配的点，条件二：以v顶点在左边的匹配点为起点能够找到一条增广路（如果能够到达条件二，说明v顶点在左边一定有对应的匹配点）。
            {
                cx[u]=v;//更新cx，cy数组
                cy[v]=u;
                return ;
            }
        }
    }
    return ;//如果程序到达了这里，说明对右边所有的顶点都访问完了，没有满足条件的。
}

二分图匹配模板

中国剩余定理

LL exgcd(LL a,LL b,LL &x,LL &y)
{
    if(!b)
    {
        x=;y=;
        return a;
    }
    LL res=exgcd(b,a%b,x,y);
    LL t=x;
    x=y;
    y=t-a/b*x;
    return res;
}
LL China()
{
    LL m=;
    LL res=;
    LL x,y;
    for(rint i=;i<=n;i++)m*=a[i];
    for(rint i=;i<=n;i++)
    {
        LL tmp=m/a[i];
        exgcd(tmp,a[i],x,y);
        res=(res+tmp*x*b[i])%m;
    }
    res=(res+m)%m;
    return res;
}

中国剩余定理模板