日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 >

P7581-「RdOI R2」路径权值【长链剖分,dp】

發(fā)布時間:2023/12/3 28 豆豆
生活随笔 收集整理的這篇文章主要介紹了 P7581-「RdOI R2」路径权值【长链剖分,dp】 小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.

正題

題目鏈接:https://www.luogu.com.cn/problem/P7581


題目大意

給出nnn個點的有邊權有根樹,mmm次詢問一個節(jié)點xxx的所有kkk級兒子兩兩之間路徑長度。

1≤n,m≤1061\leq n,m\leq 10^61n,m106


解題思路

有根長剖,無根點分治。所以這題應該是長剖(?,先離線一下詢問

然后略微分析一下,兩兩的路徑長度所以需要合并兩棵子樹向上的路徑,合并的時候又需要記錄子樹的kkk級兒子到該節(jié)點的距離和,還有kkk級兒子個數(shù)。

所以要記錄三個東西,fi,jf_{i,j}fi,j?表示iii節(jié)點的jjj級兒子個數(shù),gi,jg_{i,j}gi,j?表示iii節(jié)點的jjj級兒子到根的距離和,hi,jh_{i,j}hi,j?表示iii節(jié)點的jjj級兒子兩兩之間的路徑。

然后這三個用長剖轉移就好了。

時間復雜度O(n)O(n)O(n)


code

#include<cstdio> #include<cstring> #include<algorithm> #include<vector> #define ll long long #define mp(x,y) make_pair(x,y) #define lowbit(x) (x&-x) using namespace std; const ll N=1e6+10,P=1e9+7; struct node{ll to,next,w; }a[N<<1]; ll n,m,tot,p[N],ls[N],ans[N],len[N],son[N],dis[N],F[N],G[N],H[N]; ll *now,*tmp,*buf,*f[N],*g[N],*h[N]; vector<pair<ll,ll> >v[N]; void addl(ll x,ll y,ll w){a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;a[tot].w=w;return; } void dfs(ll x,ll fa){for(ll i=ls[x];i;i=a[i].next){ll y=a[i].to;if(y==fa)continue;dis[y]=dis[x]+a[i].w;dfs(y,x);if(len[y]>len[son[x]])son[x]=y;}len[x]=len[son[x]]+1;return; } void solve(ll x,ll fa){f[x][0]=1;g[x][0]=dis[x];if(son[x]){f[son[x]]=f[x]+1;g[son[x]]=g[x]+1;h[son[x]]=h[x]+1;solve(son[x],x);}for(ll i=ls[x];i;i=a[i].next){ll y=a[i].to;if(y==fa||y==son[x])continue;f[y]=now;now+=len[y];g[y]=tmp;tmp+=len[y];h[y]=buf;buf+=len[y];solve(y,x);for(ll j=0;j<len[y];j++){ll t1=(g[x][j+1]-f[x][j+1]*dis[x])%P;ll t2=(g[y][j]-f[y][j]*dis[x])%P;(h[x][j+1]+=t1*f[y][j]%P+f[x][j+1]*t2%P)%=P;(h[x][j+1]+=h[y][j])%=P;(g[x][j+1]+=g[y][j])%=P;f[x][j+1]+=f[y][j];}}for(ll i=0;i<v[x].size();i++){ll k=v[x][i].first,id=v[x][i].second;if(k>=len[x])ans[id]=0;else ans[id]=h[x][k];}return; } signed main() {scanf("%lld%lld",&n,&m);for(ll i=1;i<n;i++){ll x,y,w;scanf("%lld%lld%lld",&x,&y,&w);addl(x,y,w);addl(y,x,w);}for(ll i=1;i<=m;i++){ll x,k;scanf("%lld%lld",&x,&k);v[x].push_back(mp(k,i));}dfs(1,1);now=f[1]=F;now+=len[1];tmp=g[1]=G;tmp+=len[1];buf=h[1]=H;buf+=len[1];solve(1,1);for(ll i=1;i<=m;i++)printf("%lld\n",(ans[i]+P)%P);return 0; }

總結

以上是生活随笔為你收集整理的P7581-「RdOI R2」路径权值【长链剖分,dp】的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內容還不錯,歡迎將生活随笔推薦給好友。