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

    考虑每一次要求的是一条路径上深度的(k)次方和,但是(k leq 50),所以可以先把(k)的所有情况预处理出来,然后直接树上差分就行了。

    /*
      mail: mleautomaton@foxmail.com
      author: MLEAutoMaton
      This Code is made by MLEAutoMaton
    */
    #include<stdio.h>
    #include<stdlib.h>
    #include<string.h>
    #include<math.h>
    #include<algorithm>
    #include<queue>
    #include<set>
    #include<map>
    #include<iostream>
    using namespace std;
    #define ll long long
    #define re register
    #define file(a) freopen(a".in","r",stdin);freopen(a".out","w",stdout)
    inline int gi(){
    	int f=1,sum=0;char ch=getchar();
    	while(ch>'9' || ch<'0'){if(ch=='-')f=-1;ch=getchar();}
    	while(ch>='0' && ch<='9'){sum=(sum<<3)+(sum<<1)+ch-'0';ch=getchar();}
    	return f*sum;
    }
    const int N=300010,Mod=998244353;
    int n,front[N],cnt,f[N][23],dep[N],dis[N][60];
    struct node{int to,nxt;}e[N<<1];
    int qpow(int a,int b){int ret=1;while(b){if(b&1)ret=1ll*ret*a%Mod;b>>=1;a=1ll*a*a%Mod;}return ret;}
    void Add(int u,int v){e[++cnt]=(node){v,front[u]};front[u]=cnt;}
    void dfs(int u,int ff){
    	dep[u]=dep[ff]+1;f[u][0]=ff;
    	for(int k=1;k<=50;k++)dis[u][k]=(dis[ff][k]+qpow(dep[u],k))%Mod;
    	for(int i=front[u];i;i=e[i].nxt){
    		int v=e[i].to;if(v==ff)continue;
    		dfs(v,u);
    	}
    }
    int getlca(int x,int y){
    	if(dep[x]<dep[y])swap(x,y);
    	for(int i=20;~i;i--)
    		if(dep[x]-(1<<i)>=dep[y])x=f[x][i];
    	if(x==y)return x;
    	for(int i=20;~i;i--)
    		if(f[x][i]!=f[y][i])x=f[x][i],y=f[y][i];
    	return f[x][0];
    }
    int solve(int x,int y,int k){
    	int lca=getlca(x,y);
    	return (((dis[x][k]+dis[y][k])%Mod-dis[lca][k]+Mod)%Mod-dis[f[lca][0]][k]+Mod)%Mod;
    }
    int main(){
    	n=gi();
    	for(int i=1;i<n;i++){
    		int u=gi(),v=gi();
    		Add(u,v);Add(v,u);
    	}
    	dep[1]=-1;dfs(1,1);
    	f[1][0]=0;
    	for(int j=1;j<=20;j++)
    		for(int i=1;i<=n;i++)
    			f[i][j]=f[f[i][j-1]][j-1];
    	int m=gi();
    	while(m--){
    		int u=gi(),v=gi(),k=gi();
    		printf("%d
    ",solve(u,v,k));
    	}
    	return 0;
    }
    
  • 相关阅读:
    sql server 跨库操作
    FPGA使用技巧
    MATLAB学习(3)
    FPGA统计摄像头输出-基于MD9T112
    zedboard VmodCAM 图像采集 HDMI输出显示
    VmodCAM图像采集 VGA显示
    EDK中如何使用ISE中生成的IP
    如何将HDL文件实例化到XPS中
    ubuntu下安装 Source insight
    VmodCAM 初始化
  • 原文地址:https://www.cnblogs.com/mleautomaton/p/11146010.html
Copyright © 2011-2022 走看看