[BZOJ3924][Zjoi2015][点分树][暴力]幻想乡战略游戏

最新推荐文章于 2019-03-31 09:13:51 发布
最新推荐文章于 2019-03-31 09:13:51 发布
阅读量759 收藏
点赞数
分类专栏： 点分树 文章标签： 点分树
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Coldef/article/details/54882626
版权

年前的坑今天补……

题意

求一棵树的带权重心，支持修改权值。

动态树分治，也叫点分树。
就是把每层的重心连成一棵树，然后在这棵树上乱搞（具体网上教程多）。

不过第一次写这题暴力碾过去了…..好像还挺快的….

先讲暴力
假设上一次找到的重心在u,那么如果在某一点v增加了权值，那当前的重心一定是在u到v的相反方向上，只要沿着相反方向找就行了。

具体怎么找…可以这么想：
当前结点为x，y为与x相邻的结点，w[x]为x结点上的权值，cnt为总权值，那么如果cnt-w[y]< w[y]即cnt<2*w[y]时，往y点移动，直到不能移动为止。至于为什么……自己脑补一下

w[x]可以用dfs序加BIT维护。

#include <cstdio>
#include <iostream>
#include <string>
#include <cstring>
#define inf 1ll<<60
#define N 100010

using namespace std;

typedef long long ll;

int G[N],n,m,tt,rt,l[N],r[N],tc,lcA[N][25],dp[N];
ll B[N<<1],Ans,A[N],tot,ds[N];

struct edge{
    int w,t,nx;
}E[N<<1];

struct lef{
    int f,w,d;
}T[N];

inline char C(){
    static char buf[100000],*p1=buf,*p2=buf;
    if(p1==p2){
        p2=(p1=buf)+fread(buf,1,100000,stdin);
        if(p1==p2) return EOF;
    }
    return *p1++;
}

inline void reaD(int &x){
    char Ch=C();x=0;int f=1;
    for(;Ch>'9'||Ch<'0';Ch=C())if(Ch=='-')f=-1;
    for(;Ch>='0'&&Ch<='9';x=x*10+Ch-'0',Ch=C());x*=f;
}

inline void reaD(ll &x){
    char Ch=C();x=0;ll f=1;
    for(;Ch>'9'||Ch<'0';Ch=C())if(Ch=='-')f=-1;
    for(;Ch>='0'&&Ch<='9';x=x*10+Ch-'0',Ch=C());x*=f;
}

inline void InserT(int x,int y,int w){
    E[++tt].t=y;E[tt].nx=G[x];E[tt].w=w;G[x]=tt;
    E[++tt].t=x;E[tt].nx=G[y];E[tt].w=w;G[y]=tt;
}

inline void build(int g,int f,int d){
    T[g].f=f;l[g]=++tc;T[g].d=d;
    dp[g]=dp[f]+1;ds[g]=ds[f]+d;
    lcA[g][0]=f;
    for(int i=1;i<=20;i++) lcA[g][i]=lcA[lcA[g][i-1]][i-1];
    for(int i=G[g];i;i=E[i].nx)
        if(E[i].t!=f) build(E[i].t,g,E[i].w);
    r[g]=++tc;
}

inline void add(int x,int y){
    for(;x<=tc;x+=x&-x)
        B[x]+=y;
}

inline ll query(int x){
    ll res=0;
    for(;x;x-=x&-x) res+=B[x];
    return res;
}

inline ll Qlw(int x){
    return query(r[x])-query(l[x]);
}

int w[30],wt;

inline void Pt(ll x){
    if(!x){putchar(48);putchar('\n');return;}
    if(x<0){putchar('-');x=-x;};
    while(x)w[++wt]=x%10,x/=10;
    for(;wt;wt--)putchar(48+w[wt]);putchar('\n');
}

inline void swap(int &x,int &y){
    int z=x;x=y;y=z;
}

inline int clca(int x,int y){
    if(dp[x]<dp[y]) swap(x,y);
    int delt=dp[x]-dp[y],i;
    for(i=0;i<=20;i++)
        if(delt&(1<<i)) x=lcA[x][i];
    while(x!=y){
        for(i=-1;i;i++) if(lcA[x][i+1]==lcA[y][i+1]) break;
        if(i==-1) return lcA[x][0];
        x=lcA[x][i];y=lcA[y][i];
    }
    return x;
}

int main(){
    freopen("tree.in","r",stdin);
    freopen("tree.out","w",stdout);
    reaD(n);reaD(m);
    for(int i=1,x,y,w;i<n;i++)reaD(x),reaD(y),reaD(w),InserT(x,y,w);
    /*int ok=0;
    for(int i=1;i<=n;i++) if(du[i]>2){ok=1;break;}
    if(!ok){linktime();return 0;}*/
    build(1,0,0);reaD(rt);
    reaD(A[rt]);tot+=A[rt];
    add(r[rt],A[rt]);
    Pt(Ans=0);
    for(int i=1,x,y,j,lca;i<m;i++){
        reaD(x);reaD(y);
        add(r[x],y);A[x]+=y;tot+=y;
        lca=clca(x,rt);Ans+=1ll*(ds[x]+ds[rt]-2*ds[lca])*y;
        //S(rt,0,Ans);
        while(1){
            if(T[rt].f&&tot>2ll*Qlw(rt)){Ans-=1ll*(tot-2*Qlw(rt))*T[rt].d;rt=T[rt].f;continue;}
            for(j=G[rt];j;j=E[j].nx)
                if(T[rt].f!=E[j].t&&tot<2ll*Qlw(E[j].t)){
                    Ans-=1ll*(2*Qlw(E[j].t)-tot)*E[j].w;
                    rt=E[j].t;
                    break;
                }
            if(!j) break;
        }
        Pt(Ans);
    }
    return 0;
}

复杂度是 O(λnlogn) , λ 为一个常数（由数据决定），由于随机数据，所以 λ 较小，所以复杂度竟然比标算小……

正确做法是点分树。
建出点分树，每次找只要从根节点开始分治地找就行了。

复杂度 O(nlog2n) 

#include <cstdio>
#define N 100010

typedef long long ll;

int n,m,maxs,root,sizz,trot;
ll Anst,Ans,nAns;
int nG[N],G[N],cnt,V[N],p[N],dfslt[N<<2],lca[N<<2][20],dept[N],pst[N],tw[N<<2],ben[N];
ll w[N],dist[N],subd[N],dis2[N];

struct edge{
    int t,nx,t1;
    ll w;
}nE[N<<2],E[N<<2];

inline char C(){
    static char buf[100000],*p1=buf,*p2=buf;
    if(p1==p2){
        p2=(p1=buf)+fread(buf,1,100000,stdin);
        if(p1==p2) return EOF;
    }
    return *p1++;
}

inline void reaD(int &x){
    char Ch=C();x=0;int f=1;
    for(;Ch>'9'||Ch<'0';Ch=C())if(Ch=='-')f=-1;
    for(;Ch>='0'&&Ch<='9';x=x*10+Ch-'0',Ch=C());x*=f;
}

inline void InsEdge(int u,int v,int w){
    nE[++cnt].t=v;nE[cnt].nx=nG[u];nE[cnt].w=w;nG[u]=cnt;
    nE[++cnt].t=u;nE[cnt].nx=nG[v];nE[cnt].w=w;nG[v]=cnt;
}

inline int max(const int &a,const int &b){
    return a<b?b:a;
}

void Getsz(int x,int f){
    sizz++;
    for(int i=nG[x];i;i=nE[i].nx)
        if(!V[nE[i].t]&&nE[i].t!=f) Getsz(nE[i].t,x);
}

int GetRoot(int x,int f){
    int mx=0,sz=1,nsz;
    for(int i=nG[x];i;i=nE[i].nx){
        int t=nE[i].t;
        if(V[t]||t==f) continue;
        nsz=GetRoot(t,x);
        mx=max(mx,nsz);
        sz+=nsz;
    }
    mx=max(mx,sizz-sz);
    if(mx<maxs) maxs=mx,root=x;
    return sz;
}

inline void AddEdge(int x,int y,int z){
    E[++cnt].t=y;E[cnt].nx=G[x];E[cnt].t1=z;G[x]=cnt;
}

void swap(int &x,int &y){
    int z=x;x=y;y=z;
}

int LCA(int x,int y){
    int a=pst[x],b=pst[y];
    if(a>b) swap(a,b);
    int t=tw[b-a+1];
    return dept[lca[a][t]]<dept[lca[b-(1<<t)+1][t]]?lca[a][t]:lca[b-(1<<t)+1][t];
}

int divont(int x,int f){
    sizz=0,maxs=1<<30,Getsz(x,0),GetRoot(x,0);
    int nRoot=root,nxRoot;V[nRoot]=1;p[nRoot]=f;
    for(int i=nG[nRoot];i;i=nE[i].nx)
        if(nE[i].t!=f&&!V[nE[i].t]){
            nxRoot=divont(nE[i].t,nRoot);
            AddEdge(nRoot,nxRoot,nE[i].t);
        }
    return nRoot;
}

void dfs(int x,int f){
    dept[x]=dept[f]+1;dfslt[pst[x]=++cnt]=x;
    for(int i=nG[x];i;i=nE[i].nx)
        if(nE[i].t!=f)dist[nE[i].t]=dist[x]+nE[i].w,dfs(nE[i].t,x),dfslt[++cnt]=x;
}

void Prelca(){
    for(int i=1;i<=cnt;i++) tw[i]=tw[i-1]+((1<<tw[i-1]+1)==i);
    for(int i=1;i<=cnt;i++) lca[i][0]=dfslt[i];
    for(int k=1;k<=tw[cnt];k++)
        for(int i=1;i+(1<<k)-1<=cnt;i++)
        lca[i][k]=dept[lca[i][k-1]]<dept[lca[i+(1<<k-1)][k-1]]?lca[i][k-1]:lca[i+(1<<k-1)][k-1];
}

ll dis(int x,int y){
    return dist[x]+dist[y]-2*dist[LCA(x,y)];
}

inline void Addtr(int x,int y){
    w[x]+=y;
    for(int i=x,j;p[i];i=p[i]){
        w[p[i]]+=y;
        subd[p[i]]+=dis(x,p[i])*y;
        dis2[i]+=dis(x,p[i])*y;
    }
}

inline ll min(const ll &a,const ll &b){
    return a<b?a:b;
}

inline ll disf(int x){
    ll re=subd[x];
    for(int i=x;p[i];i=p[i]){
        ll disr=dis(x,p[i]);
        re+=subd[p[i]]-dis2[i];
        re+=(w[p[i]]-w[i])*disr;
    }
    return re;
}

inline ll query(int x){
    bool flg=1;
    ll tot=disf(x);
    for(int i=G[x];i;i=E[i].nx){
        ll cost=disf(E[i].t1);
        if(cost<tot) return query(E[i].t);
    }
    return tot;
}

int main(){
    freopen("tree.in","r",stdin);
    freopen("tree.out","w",stdout);
    reaD(n);reaD(m);
    for(int i=1,u,v,w;i<n;i++)reaD(u),reaD(v),reaD(w),InsEdge(u,v,w);
    cnt=0,dfs(1,0),Prelca(),
    cnt=0,trot=divont(1,0),cnt=0;
    for(int i=1;i<=m;i++){
        int x,y;
        reaD(x);reaD(y);
        Ans=(Anst+=y*dis(x,trot));
        Addtr(x,y),cnt+=y;
        printf("%lld\n",query(trot));
    }
}

事实证明如果考场上想不出标算或没时间写标算，这种信仰暴力还是可以接受的……

LowestJN
关注
专栏目录
bzoj4876: [Zjoi2017]线段树（树上差分）
苟为蒟蒻又何妨
06-27 182
传送门 对于这种广义线段树的题，想要定位区间[l,r][l,r][l,r]可以转化成定位节点[l−1,l−1],[r+1,r+1][l-1,l-1],[r+1,r+1][l−1,l−1],[r+1,r+1]路径上的一些儿子节点。 设lca=tlca=tlca=t 那么路径[l−1,l−1]→t[l-1,l-1]\rightarrow t[l−1,l−1]→t的右儿子和路径[r+1,r+1]→t[r...
bzoj5515: [Zjoi2019]线段树（线段树）
苟为蒟蒻又何妨
06-27 214
传送门 每次等价于一部分操作一部分不操作。 显然一个点被标记只可能是被直接标记或者祖先被标记，于是我们对于每一个节点 维护它被标记的概率与祖先被标记的概率。 把每个区间分类讨论更新即可。 有如下几种情况： 1.覆盖祖先但未覆盖自己 2. 直接覆盖当前点 3. 经过当前点继续在子树中修改 4. 修改到父亲处向兄弟子树走 5. 修改到非父亲祖先向兄弟子树走 然后用线段树维护即可。 详见代码 代码： #...
zjoi2015.zip
02-23
浙江省选 2015 zjoi 数据 题目zjoi2015.zip 2015浙江省省赛数据。
【BZOJ3924】【Zjoi2015】幻想乡战略游戏（树链剖分+点分治）
(■-■)
03-29 518
Description傲娇少女幽香正在玩一个非常有趣的战略类游戏，本来这个游戏的地图其实还不算太大，幽香还能管得过来，但是不知道为什么现在的网游厂商把游戏的地图越做越大，以至于幽香一眼根本看不过来，更别说和别人打仗了。 在打仗之前，幽香现在面临一个非常基本的管理问题需要解决。 整个地图是一个树结构，一共有n块空地，这些空地被n-1条带权边连接起来，使得每两个点之间有一条唯一的路径将它们连接起来。在游
BZOJ3924: [Zjoi2015]幻想乡战略游戏
Don't fake it till you made it
03-15 1048
SB题目卡我常数 最后寄刀片成功！时限提高到100s A了！！#include<cstdio> #include<cstring> #include<iostream> #include<cstdlib> using namespace std; #define ll long long #define youhua __attribute__((optimize("O3"))) char c;
[zjoi2015]幻想乡战略游戏
zhouyuheng2003的博客
03-31 547
前言 略略略 题目相关 链接 题目大意 给出一棵树，每次修改一个点的权值，维护一个带权重心 啥是带权重心？ 设点iii的值为ViV_iVi​我们要选一个点uuu，每个点对应一个值： ∑v=1nVv∗dis(u,v)\sum_{v=1}^n V_v*dis(u,v)v=1∑n​Vv​∗dis(u,v) 我们要选一个uuu使得这个值最小 输出最小的值 数据范围 n,q≤105n,q\le10^5n,q...
[ZJOI2015]幻想乡战略游戏
待成熟的葡萄
05-13 2093
题目大意一个nn个节点的树，点带有权值dd，边也带有权值（小于等于10001000的正整数）。初始时点权为00。 定义dist(u,v)dist(u,v)为uu与vv两点之间简单路径的边权和。 会有QQ次操作，形如： ∙u e\bullet u\ e，表示dud_u加上e(0≤|e|≤1000)e(0\le|e|\le1000) 每次操作之后，你都需要找到一个点xx，最小化 w=∑i=1n
BZOJ 1057 ZJOI 2007 棋盘制作 DP+悬线法
(╯°口°)╯(┴—┴
01-14 915
题目大意：给出一个由01形成的矩阵，问这个矩阵中最大面积的正方形和矩形，其中任意一个方块相邻的都是不同的格子。 思路：其实吧所有(i + j)&1的位置上的数字异或一下，就变成都是0或者都是1的最大正方形和矩形了。第一问就是水DP，第二问可以单调栈或者悬线。都很好写。 CODE： #include #include #include #include #defin
BZOJ 4027 HEOI2015 兔子与樱花 树形贪心
世界
05-05 1749
题目大意：给定一棵有根树，每个点上有一些樱花，现在要求删除一些节点，删除节点的樱花和子节点都会连到父节点上，要求每个节点的樱花数+子节点数不超过mm，求最多删多少个节点这数据范围也只能贪心了吧= = 令fif_i为以节点ii为根的子树中能删除的最多节点(ii节点不删)，gig_i为删除最多节点的情况下ii号节点的最小负重 那么首先对于每个节点我们对于所有的子节点为根的子树尽量删，然后考虑如何删除
BZOJ 3924: [Zjoi2015]幻想乡战略游戏
Flere825的博客
06-03 602
题意：动态维护支持修改点权的树上带权重心 Sol: 看见东方就点进来了233 膜拜下陈老师出的题 首先考虑如何求一个点到树上其他点的带权距离和，动态点分治可以做到单词O(logn) 现在考虑暴力，发现如果求出了每个点到其他点的带权距离和记为权值，那么以带权重心为根，权值从上到下递增，换句话说对于任何一条链，权值都是单峰的。 暴力做法可以随便选一个点枚举出边，计算出边指向的点的权值的最小
【ZJOI 2015 幻想乡战略游戏】【动态点分治】
Kumii
12-18 666
【ZJOI 2015 幻想乡战略游戏】【动态点分治】
3924: [Zjoi2015]幻想乡战略游戏(点分树)
少主大人殿下的博客
04-03 335
题目 链接 分析 点分树裸题！！！ 其实是不想写题解 大概就是要新建立点分树，基本上是一个新树了。要维护以每个结点为根的子树信息来辅助查询。 查询的时候就是找子树中比当前结点优的点然后。这里跳到那个结点的根结点去。 还有就是会被菊花卡成O(n)的查询，但是这题有限制子树不超过20个。 代码 #include&lt;bits/stdc++.h&gt;//注意距离都要开LL ...
bzoj3924 [Zjoi2015]幻想乡战略游戏
Elijahqi
02-24 277
Description 傲娇少女幽香正在玩一个非常有趣的战略类游戏，本来这个游戏的地图其实还不算太大，幽香还能管得过来，但是不知道为什么现在的网游厂商把游戏的地图越做越大，以至于幽香一眼根本看不过来，更别说和别人打仗了。 在打仗之前，幽香现在面临一个非常基本的管理问题需要解决。 整个地图是一个树结构，一共有n块空地，这些空地被n-1条带权边连接起来，使得每两个点之间有一条唯一的路径将它们连接起来...
[点分树] BZOJ 3924 [Zjoi2015]幻想乡战略游戏
雯舞
12-27 489
%%%CLJ 首先动态维护每个点的答案？动态树分治就好了 怎么找重心 对于分治结构的每一个点，我们枚举它的出边 如果某条出边连向的点的距离之和小于当前点，那么答案一定在那条出边指向的子树中，分治做下去就行了 这样是两个log再乘20的  实际上最优可以一个log 跑的竟然还没有当年不懂事打的暴力快 #include #include #include us
[ZJOI2015] 幻想乡战略游戏
INCINCIBLE的博客
04-25 678
P3254【ZJOI2015 Day1】幻想乡战略游戏 时间限制 : 120000 MS   空间限制 : 512000 KB 问题描述         傲娇少女幽香正在玩一个非常有趣的战略类游戏，本来这个游戏的地图其实还不算太大，幽香还能管得过来，但是不知道为什么现在的网游厂商把游戏的地图越做越大，以至于幽香一眼根本看不过来，
bzoj 3924: [Zjoi2015]幻想乡战略游戏
lethalboy
12-21 2844
3924: [Zjoi2015]幻想乡战略游戏 Time Limit: 100 Sec  Memory Limit: 256 MB Submit: 630  Solved: 291 [Submit][Status][Discuss] Description  傲娇少女幽香正在玩一个非常有趣的战略类游戏，本来这个游戏的地图其实还不算太大，幽香还能管得过来，但是不知道为什么现在的网游厂商
洛谷P3345：[ZJOI2015]幻想乡战略游戏 （动态点分治）
KsCla
03-20 798
题目传送门：https://www.luogu.org/problemnew/show/P3345 题目分析：这题简直是噩梦，调了我一个下午。 首先考虑暴力怎么做：在树上随机找一个点，然后计算以该点为补给点的答案，并算出它周围的点的答案。如果它周围存在一个点答案更优，就往该点走，直到当前点为最优点就停止。这样一次询问就是n2dn2dn^2d的，其中d=20d=20d=20，因为每个点的度数...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值