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

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

CF932F-Escape Through Leaf【树上启发式合并,CDQ分治,斜率优化dp】

發布時間:2023/12/3 编程问答 24 豆豆
生活随笔 收集整理的這篇文章主要介紹了 CF932F-Escape Through Leaf【树上启发式合并,CDQ分治,斜率优化dp】 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

正題

題面鏈接:https://www.luogu.com.cn/problem/CF932F


題目大意

nnn個點的一棵樹,從xxx跳到yyy(要求yyyxxx的子樹中)會產生Ax?ByA_x*B_yAx??By?的代價,求每個節點出發跳到某個葉節點的最小代價。


解題思路

考慮dpdpdp的話,那么有fx=fy+Ax?Byf_x=f_y+A_x*B_yfx?=fy?+Ax??By?,這個式子可以考慮斜率優化,若y1y_1y1?y2y_2y2?優,那么有fy1?fy2By1?By2≥Ax\frac{f_{y_1}-f_{y_2}}{B_{y_1}-B_{y_2}}\geq A_xBy1???By2??fy1???fy2???Ax?
也就是我們對于每個節點要維護一個子樹里所有點構成的一個下凸殼。

考慮樹上啟發式合并+CDQ+CDQ+CDQ,我們要求一個序列使得被貢獻的點排在貢獻點的后面。維護一個序列,每次我們保留重子樹的序列,然后再加入其它輕子樹的序列,當到一個節點的頭頂上是一條輕邊時我們就對這個序列跑一次CDQCDQCDQ來維護凸殼,然后清空序列。需要注意的是對于二次掃描輕子樹的節點需要標記不能在CDQCDQCDQ分治中被修改答案。

時間復雜度O(nlog?2n)O(n\log^2 n)O(nlog2n)


codecodecode

#include<cstdio> #include<cstring> #include<algorithm> #include<stack> #define ll long long using namespace std; const ll N=1e5+10; struct node{ll to,next; }e[N*2]; ll n,tot,ls[N],a[N],b[N],f[N],siz[N],son[N]; ll cnt,q[N],p[N],v[N],st[N],rfn[N]; void addl(ll x,ll y){e[++tot].to=y;e[tot].next=ls[x];ls[x]=tot;return; } void dfs(ll x,ll fa){siz[x]=1;for(ll i=ls[x];i;i=e[i].next){ll y=e[i].to;if(y==fa)continue;dfs(y,x);siz[x]+=siz[y];if(siz[y]>siz[son[x]])son[x]=y;}return; } bool cmp(ll x,ll y) {return a[x]<a[y];} bool cMp(ll x,ll y) {return (b[x]==b[y])?(f[x]<f[y]):(b[x]<b[y]);} double slope(ll x,ll y) {return (double)(f[x]-f[y])/(b[x]-b[y]);} void cdq(ll l,ll r){if(l==r)return;ll mid=(l+r)>>1,cnt1=l-1,cnt2=mid;for(ll i=l;i<=r;i++)if(rfn[p[i]]<=mid)q[++cnt1]=p[i];else q[++cnt2]=p[i];for(ll i=l;i<=r;i++)p[i]=q[i];cdq(l,mid);ll tot=0;for(ll i=l;i<=mid;i++){if(b[p[i]]==b[p[i-1]]&&i!=l)continue;while(tot>1&&slope(st[tot-1],st[tot])>slope(st[tot-1],p[i]))tot--;st[++tot]=p[i];}for(ll i=mid+1;i<=r;i++){if(v[p[i]])continue;while(tot>1&&slope(st[tot-1],st[tot])>-a[p[i]])tot--;ll x=p[i],y=st[tot];f[x]=min(f[x],f[y]+a[x]*b[y]);}cdq(mid+1,r);sort(p+l,p+1+r,cMp);return; } void calc(ll x,ll fa){p[++cnt]=x;rfn[x]=cnt;v[x]=1;for(ll i=ls[x];i;i=e[i].next){ll y=e[i].to;if(y==fa)continue;calc(y,x);} } void solve(ll x,ll fa,ll top){for(ll i=ls[x];i;i=e[i].next){ll y=e[i].to;if(y==fa||y==son[x])continue;solve(y,x,y);}if(son[x])solve(son[x],x,top);else f[x]=0;for(ll i=ls[x];i;i=e[i].next){ll y=e[i].to;if(y==fa||y==son[x])continue;calc(y,x);}p[++cnt]=x;rfn[x]=cnt;v[x]=(!son[x]);if(x==top){sort(p+1,p+1+cnt,cmp);cdq(1,cnt);cnt=0;}return; } int main() {scanf("%lld",&n);for(ll i=1;i<=n;i++)scanf("%lld",&a[i]);for(ll i=1;i<=n;i++)scanf("%lld",&b[i]);for(ll i=1;i<n;i++){ll x,y;scanf("%lld%lld",&x,&y);addl(x,y);addl(y,x);}memset(f,0x3f,sizeof(f));dfs(1,1);solve(1,1,1);for(ll i=1;i<=n;i++)printf("%lld\n",f[i]);return 0; }

總結

以上是生活随笔為你收集整理的CF932F-Escape Through Leaf【树上启发式合并,CDQ分治,斜率优化dp】的全部內容,希望文章能夠幫你解決所遇到的問題。

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