zoukankan      html  css  js  c++  java
  • [BJOI2018]求和

    link

    其实可以用$sum(i,j)$表示从$i$到$1$的$k$次方的值,然后就是$lca$的基本操作

    注意,能一起干的事情就一起搞,要不会超时

    #include<iostream>
    #include<cstring>
    #include<cstdio>
    #include<algorithm>
    #define int long long
    #define mod 998244353
    using namespace std;
    const int N=300001;
    inline int read(){
        int f=1,ans=0;char c=getchar();
        while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
        while(c>='0'&&c<='9'){ans=ans*10+c-'0';c=getchar();}
        return f*ans;
    }
    struct node{
        signed u,v,nex;
    }x[N<<1];
    signed n,m,cnt;
    int sum[N][51];
    signed head[N],deep[N],fa[N][20];
    void add(signed u,signed v){
        x[cnt].u=u,x[cnt].v=v,x[cnt].nex=head[u],head[u]=cnt++;
    }
    int ksm(int a,int b){
        if(a==1) return 1;
        int ans=1;
        a%=mod;
        while(b){
            if(b&1) ans*=a,ans%=mod;
            a*=a,a%=mod;
            b>>=1;
        }return ans;
    }
    void dfs(signed f,signed fath){
        deep[f]=deep[fath]+1,fa[f][0]=fath;
        for(signed i=1;i<=50;i++) sum[f][i]=sum[fath][i]+ksm(deep[f],i);
        for(signed i=1;(1<<i)<=deep[f];++i) fa[f][i]=fa[fa[f][i-1]][i-1];
        for(signed i=head[f];i!=-1;i=x[i].nex){
            if(x[i].v==fath) continue;
            dfs(x[i].v,f);
        }
    }
    int Log2[N];
    signed lca(signed u,signed v){
        if(deep[u]<deep[v]) swap(u,v);
        for(signed i=Log2[u];i>=0;--i)
            if(deep[u]-(1<<i)>=deep[v]) u=fa[u][i];
        if(u==v) return u;
        for(signed i=Log2[v];i>=0;--i){
            if(fa[u][i]==fa[v][i]) continue;
            u=fa[u][i],v=fa[v][i];
        }return fa[u][0];
    }
    signed q;
    bool ff;
    signed main(){
        memset(head,-1,sizeof(head));
        n=read();
        for(signed i=1;i<n;++i){
            int u=read(),v=read();
            add(u,v),add(v,u);
        }deep[0]=-1;Log2[0]=1;
        for(int i=1;i<=n;++i) Log2[i]=Log2[i>>1]+1;
        dfs(1,0);
        q=read();
        while(q--){
            int u=read(),v=read(),k=read();
            int ls=lca(u,v);
            printf("%d
    ",(((sum[u][k]+sum[v][k]-2*sum[ls][k]+ksm(deep[ls],k))%mod+mod)%mod));
        }
        return 0;
    }
    View Code
  • 相关阅读:
    JAVA微信公众号网页开发——获取公众号关注的所有用户
    删除mysl
    sql语言(mysql)
    mycat读写分离
    mysql双主双从技术
    实用的10个日志处理案例
    ansible基本操作
    MySQL改密
    mysql源码包安装
    ftp搭建mysql服务器
  • 原文地址:https://www.cnblogs.com/si-rui-yang/p/10088072.html
Copyright © 2011-2022 走看看