BZOJ 3924: [Zjoi2015]幻想乡战略游戏(动态点分治)
2018-08-15 09:43
344 查看
题目描述
傲娇少女幽香正在玩一个非常有趣的战略类游戏,本来这个游戏的地图其实还不算太大,幽香还能管得过来,但是不知道为什么现在的网游厂商把游戏的地图越做越大,以至于幽香一眼根本看不过来,更别说和别人打仗了。在打仗之前,幽香现在面临一个非常基本的管理问题需要解决。 整个地图是一个树结构,一共有n块空地,这些空地被n-1条带权边连接起来,使得每两个点之间有一条唯一的路径将它们连接起来。
在游戏中,幽香可能在空地上增加或者减少一些军队。同时,幽香可以在一个空地上放置一个补给站。 如果补给站在点u上,并且空地v上有dv个单位的军队,那么幽香每天就要花费dv*dist(u,v)的金钱来补给这些军队。
由于幽香需要补给所有的军队,因此幽香总共就要花费为Sigma(Dv*dist(u,v),其中1<=V<=N)的代价。其中dist(u,v)表示u个v在树上的距离(唯一路径的权和)。
因为游戏的规定,幽香只能选择一个空地作为补给站。在游戏的过程中,幽香可能会在某些空地上制造一些军队,也可能会减少某些空地上的军队,进行了这样的操作以后,出于经济上的考虑,幽香往往可以移动他的补给站从而省一些钱。
但是由于这个游戏的地图是在太大了,幽香无法轻易的进行最优的安排,你能帮帮她吗? 你可以假定一开始所有空地上都没有军队。
输入输出格式
输入格式:第一行两个数n和Q分别表示树的点数和幽香操作的个数,其中点从1到n标号。 接下来n-1行,每行三个正整数a,b,c,表示a和b之间有一条边权为c的边。 接下来Q行,每行两个数u,e,表示幽香在点u上放了e单位个军队(如果e<0,就相当于是幽香在u上减少了|e|单位个军队,说白了就是du←du+e)。数据保证任何时刻每个点上的军队数量都是非负的。
输出格式:
对于幽香的每个操作,输出操作完成以后,每天的最小花费,也即如果幽香选择最优的补给点进行补给时的花费。
输入输出样例
输入样例#1: 复制10 5 1 2 1 2 3 1 2 4 1 1 5 1 2 6 1 2 7 1 5 8 1 7 9 1 1 10 1 3 1 2 1 8 1 3 1 4 1
输出样例#1: 复制
0 1 4 5 6
说明
对于所有数据,1<=c<=1000, 0<=|e|<=1000, n<=10^5, Q<=10^5 非常神奇的是,对于所有数据,这棵树上的点的度数都不超过20,且N,Q>=1[b]题解[/b]
噩梦……真的噩梦……我快被幽香给玩死了……
用这道题来理解动态淀粉质点分治的(因为当初捉迷藏那题直接学岛娘的括号序列……就没去再打一遍……),然而题解大佬们都直接默认我们已经会动态点分然后直接上……结果我昨天看了一个晚上才弄懂……然后今天又花了一个早上码完,交上去竟然1A了不可思议
先说一下什么是动态点分治吧。对于一般的点分治,因为我们每一次都找出重心,所以每一次递归往左右子树找的时候深度不会超过$O(log n)$层。但是如果有了修改操作怎么办呢?每修改一次点分一次么?那怕是得T飞
我们来看一下这道题目,它的每一次修改,更改的只有点权,树的结构是没有变化的(如果有的话怕是只能上LCT了……),也就是说,每一次点分的时候找到的重心是不会改变的。那么,我们可不可以把点分治每一层的重心给连成一棵树呢?因为点分的递归层数只有$O(log n)$层,所以这棵点分树的深度也是$O(log n)$的
那么考虑一下在每一个点分树的节点维护什么东西,对于每一个节点,我们维护他的子树中的所有信息,也就是在原树中它被选为重心时的那个子树的所有信息。修改的时候,只要从一个点开始在点分树里往上跳,并不断更新信息即可。查询一个点的时候,点分树里跳,不断考虑与父亲之间的贡献就好了。
上面那一段看不懂也没关系,因为我只是在口胡,假装自己已经很懂动态点分的样子
那么我们具体来分析一下这道题目我们要维护什么
题目要求使$\sum d_v*dis(u,v)$最小,其中$d_v$为$v$点的点权,$dis(u,v)$为原树中$u,v$两点的距离。题目要求就是求带权重心。假设我们当前已经选定了点$v$为答案,那么考虑它的一个子节点$w$,如果把补给站从点$v$转移到$w$,那么$w$的所有子树内的点到补给站的距离少了$dis(v,w)$,而其他所有点到补给站的距离多了$dis(v,w)$。我们假设$sum_v[v]$为$v$的子树内的点权和,那么答案总共的变化量是$$dis(u,v)*(sum_v[v]-sum_v[w]-sum_v[w])$$
不难发现,当$sum_v[w]*2>sum_v[v]$的时候,答案的变化量是小于零的,也就是说代价减小,可以变得更优。而且,满足这样条件的点$w$最多只有一个
那么我们可以每一次选定一个点,然后看看往他的哪个子树走更优,如果没有说明他自己就已经是最优的了。这样不断下去肯定能找到答案。
但是由于原图可能是一条链,要怎么做才能保证复杂度呢?我们选择在点分树上走,每一次都跳到它下一层的重心,这样可以保证层数最多只有$O(log n)$
然后考虑如何维护答案。不难发现,对于点分树上一个点$v$,它的子树中的点就是在点分治时它被选为重心时的那棵树上的点。考虑点分树上的一对父子$u,v$,我们设$sum_a[v]$表示$v$的子树内的所有点到他的代价之和,$sum_b[v]$为$v$的子树内的所有点到$v$点父亲(也就是点$u$)的距离之和,$sum_v[v]$还是表示子树的点权之和。那么我们设答案已经选定为点$v$,加上$u$的不包括$v$的子树后答案就是$sum_a[v]-sum_b[v]+sum_a[u]+sum_v[u]*dis(u,v)$。于是只要在点分树上不断找父亲并合并,就可以知道答案了
然后我们只要能在修改时维护好这三个数组就可以了!!!
至于修改时如何维护呢?我们修改一个点之后,然后不断在点分树上往父节点跳,并不断更新即可。查询的时候也是,不断跳并合并答案
然后又新学会了一招,用$RMQ O(1)$查询$LCA$(只会倍增和树剖的我瑟瑟发抖),总时间复杂度$O(nlog^2n)$
//minamoto #include<cstdio> #include<iostream> #include<cstring> #define ll long long #define N 100005 #define inf 0x3f3f3f3f #define rint register int using namespace std; #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++) char buf[1<<21],*p1=buf,*p2=buf; template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;} inline int read(){ #define num ch-'0' char ch;bool flag=0;int res; while(!isdigit(ch=getc())) (ch=='-')&&(flag=true); for(res=num;isdigit(ch=getc());res=res*10+num); (flag)&&(res=-res); #undef num return res; } char sr[1<<21],z[20];int C=-1,Z; inline void Ot(){fwrite(sr,1,C+1,stdout),C=-1;} inline void print(ll x){ if(C>1<<20)Ot();if(x<0)sr[++C]=45,x=-x; while(z[++Z]=x%10+48,x/=10); while(sr[++C]=z[Z],--Z);sr[++C]='\n'; } struct G{ int head ,Next[N<<1],edge[N<<1],ver[N<<1],tot; G(){tot=0;memset(head,0,sizeof(head));} inline void add(int u,int v,int e){ ver[++tot]=v,Next[tot]=head[u],head[u]=tot,edge[tot]=e; } }T1,T2; int n,q,st[N<<1][18],logn[N<<1],bin[25],tp; ll sum,ans,d ,dis1 ,dis2 ,sumv ; int dfn ,num; void dfs1(int u,int fa){ st[dfn[u]=++num][0]=d[u]; for(int i=T1.head[u];i;i=T1.Next[i]){ int v=T1.ver[i]; if(v==fa) continue; d[v]=d[u]+T1.edge[i],dfs1(v,u),st[++num][0]=d[u]; } } inline ll LCA(int a,int b){ if(dfn[a]>dfn[b]) a^=b^=a^=b; int k=logn[dfn[b]-dfn[a]+1]; return min(st[dfn[a]][k],st[dfn[b]-bin[k]+1][k])<<1; } inline ll dis(int a,int b){return d[a]+d[b]-LCA(a,b);} int sz ,son ,size,rt,fa ;bool vis ; void dfs2(int u,int fa){ sz[u]=1,son[u]=0; for(int i=T1.head[u];i;i=T1.Next[i]){ int v=T1.ver[i]; if(vis[v]||v==fa) continue; dfs2(v,u),sz[u]+=sz[v],cmax(son[u],sz[v]); } cmax(son[u],size-sz[u]); if(son[u]<son[rt]) rt=u; } void dfs3(int u){ vis[u]=true; for(int i=T1.head[u];i;i=T1.Next[i]){ int v=T1.ver[i]; if(vis[v]) continue; rt=0,size=sz[v],son[0]=n+1; dfs2(v,0),T2.add(u,rt,v),fa[rt]=u,dfs3(rt); } } inline void update(int u,int val){ sumv[u]+=val; for(int p=u;fa[p];p=fa[p]){ ll dist=dis(fa[p],u)*val; dis1[fa[p]]+=dist; dis2[p]+=dist; sumv[fa[p]]+=val; } } inline ll calc(int u){ ll ans=dis1[u]; for(int p=u;fa[p];p=fa[p]){ ll dist=dis(fa[p],u); ans+=dis1[fa[p]]-dis2[p]; ans+=dist*(sumv[fa[p]]-sumv[p]); } return ans; } ll query(int u){ ll ans=calc(u); for(int i=T2.head[u];i;i=T2.Next[i]){ ll tmp=calc(T2.edge[i]); if(tmp<ans) return query(T2.ver[i]); } return ans; } void init(){ n=read(),q=read(); bin[0]=1,logn[0]=-1; for(rint i=1;i<=20;++i) bin[i]=bin[i-1]<<1; while(bin[tp+1]<=(n<<1)) ++tp; for(rint i=1;i<=(n<<1);++i) logn[i]=logn[i>>1]+1; for(rint i=1;i<n;++i){ rint u=read(),v=read(),e=read(); T1.add(u,v,e),T1.add(v,u,e); } dfs1(1,0),rt=0,son[0]=n+1,size=n,dfs2(1,0); for(rint j=1;j<=tp;++j) for(rint i=1;i+bin[j]-1<=(n<<1);++i) st[i][j]=min(st[i][j-1],st[i+bin[j-1]][j-1]); } int main(){ init(); int LastOrder=rt;dfs3(rt); while(q--){ int x=read(),y=read();update(x,y); print(query(LastOrder)); } Ot(); return 0; }
相关文章推荐
- [BZOJ3924][Zjoi2015]幻想乡战略游戏(动态点分治)
- [BZOJ3924][ZJOI2015]幻想乡战略游戏(动态点分治)
- 【bzoj3924】[Zjoi2015]幻想乡战略游戏 动态点分治
- BZOJ 3924: [Zjoi2015]幻想乡战略游戏(动态点分治)
- bzoj3924 [Zjoi2015]幻想乡战略游戏(动态点分治)
- [ZJOI2015][bzoj3924] 幻想乡战略游戏 [动态点分治]
- 【BZOJ3924】【Zjoi2015】幻想乡战略游戏(树链剖分+点分治)
- bzoj 3924: [Zjoi2015]幻想乡战略游戏
- bzoj 3924: [Zjoi2015]幻想乡战略游戏 (树链剖分)
- BZOJ 3924 Zjoi2015 幻想乡战略游戏 动态树分治
- [点分树] BZOJ 3924 [Zjoi2015]幻想乡战略游戏
- 洛谷P3345:[ZJOI2015]幻想乡战略游戏 (动态点分治)
- BZOJ 3924: [Zjoi2015]幻想乡战略游戏 动态树分治
- BZOJ.3924.[ZJOI2015]幻想乡战略游戏(动态点分治)
- 【BZOJ 3924】[Zjoi2015]幻想乡战略游戏
- bzoj3924 [Zjoi2015]幻想乡战略游戏
- BZOJ3924 : [Zjoi2015]幻想乡战略游戏
- BZOJ3924 : [Zjoi2015]幻想乡战略游戏
- bzoj3924: [Zjoi2015]幻想乡战略游戏 //动态点分治
- bzoj 3924: [Zjoi2015]幻想乡战略游戏 动态树分治