bzoj3052: [wc2013]糖果公园
又是一代神题。
uoj测速rank10,bzoj测速rank26(截止当前2016.5.30 12:58)
带修改的树上莫队。
修改很少,块的大小随便定都能A
然而我一开始把开3次根写成了pow(blabla,1/3)
我一副见了鬼的样子.jpg
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <algorithm>
#define ll long long
#define N 200005
#define M 200005 using namespace std;
inline int read(){
int ret=0;char ch=getchar();
while (ch<'0'||ch>'9') ch=getchar();
while ('0'<=ch&&ch<='9'){
ret=ret*10-48+ch;
ch=getchar();
}
return ret;
} struct edge{
int adj,next;
edge(){}
edge(int _adj,int _next):adj(_adj),next(_next){}
} e[M];
int n,g[N],m;
void AddEdge(int u,int v){
e[++m]=edge(v,g[u]);g[u]=m;
e[++m]=edge(u,g[v]);g[v]=m;
}
int fa[N],deep[N],dfn[N],stamp;
void dfs(int u){
dfn[u]=++stamp;
deep[u]=deep[fa[u]]+1;
for (int i=g[u];i;i=e[i].next){
int v=e[i].adj;
if (v==fa[u]) continue;
fa[v]=u;
dfs(v);
}
++stamp;
}
int anc[N][19];
void precompute(){
fa[1]=fa[0]=deep[0]=0;stamp=0;
dfs(1);
memset(anc[0],0,sizeof(anc[0]));
for (int i=1;i<=n;++i) anc[i][0]=fa[i];
for (int k=1;k<17;++k)
for (int i=1;i<=n;++i)
anc[i][k]=anc[anc[i][k-1]][k-1];
}
int qlca(int u,int v){
if (deep[u]<deep[v]) swap(u,v);
for (int k=16;k>=0;--k) if (deep[anc[u][k]]>=deep[v]) u=anc[u][k];
for (int k=16;k>=0;--k) if (anc[u][k]!=anc[v][k]) u=anc[u][k],v=anc[v][k];
return u==v?u:fa[u];
} struct Query{
int u,v,tm,id;
} q[N];
int bl[N];
inline bool operator <(const Query &x,const Query &y){
if (bl[dfn[x.u]]!=bl[dfn[y.u]]) return bl[dfn[x.u]]<bl[dfn[y.u]];
else if (bl[dfn[x.v]]!=bl[dfn[y.v]]) return bl[dfn[x.v]]<bl[dfn[y.v]];
else return ((x.tm<y.tm)^(bl[dfn[x.v]]&1));
}
int Q; int cdy,c[N],val[N],kk[N];
int cnt[N];
bool exist[N];
ll nowans;
void change(int u){
if (exist[u]) nowans-=(ll)val[c[u]]*kk[cnt[c[u]]--];
else nowans+=(ll)val[c[u]]*kk[++cnt[c[u]]];
exist[u]^=1;
}
void moveTo(int u,int v){
int w=qlca(u,v);
for (;u!=w;u=fa[u]) change(u);
for (;v!=w;v=fa[v]) change(v);
}
void modify(int u,int to){
if (exist[u]) nowans-=(ll)val[c[u]]*kk[cnt[c[u]]--];
c[u]=to;
if (exist[u]) nowans+=(ll)val[c[u]]*kk[++cnt[c[u]]];
} int mdf,pt[N],fr[N],to[N];
int tmpc[N];
ll ans[N];
int main(){
// freopen("candyland.in","r",stdin);
// freopen("candyland.out","w",stdout);
n=read();cdy=read();Q=read();
for (int i=1;i<=cdy;++i) val[i]=read();
for (int i=1;i<=n;++i) kk[i]=read();
memset(g,0,sizeof(g));m=1;
for (int i=1;i<n;++i) AddEdge(read(),read());
for (int i=1;i<=n;++i) tmpc[i]=c[i]=read();
precompute();
mdf=0;
for (int i=1;i<=Q;++i)if (read()){
q[i].u=read();q[i].v=read();
if (dfn[q[i].u]>dfn[q[i].v]) swap(q[i].u,q[i].v);
q[i].tm=mdf;
q[i].id=i;
}
else{
--Q;--i;++mdf;
pt[mdf]=read();
fr[mdf]=c[pt[mdf]];
to[mdf]=c[pt[mdf]]=read();
}
int K=max(sqrt(n),pow((double)4*n*n*mdf/Q,1.0/3));
for (int i=1;i<=2*n;++i) bl[i]=i/K;
sort(q+1,q+Q+1);
for (int i=1;i<=n;++i) c[i]=tmpc[i];
memset(cnt,0,sizeof(cnt));
memset(exist,0,sizeof(exist));
q[0].u=q[0].v=1;q[0].tm=0;
nowans=0;
int nowtm=0,w;
for (int i=1;i<=Q;++i){
moveTo(q[i-1].u,q[i].u);
moveTo(q[i-1].v,q[i].v);
change(w=qlca(q[i].u,q[i].v));
for (;nowtm<q[i].tm;++nowtm) modify(pt[nowtm+1],to[nowtm+1]);
for (;nowtm>q[i].tm;--nowtm) modify(pt[nowtm],fr[nowtm]);
ans[q[i].id]=nowans;
change(w);
} for (int i=1;i<=Q;++i) printf("%lld\n",ans[i]);
return 0;
}
bzoj3052: [wc2013]糖果公园的更多相关文章
- BZOJ3052:[WC2013]糖果公园(树上莫队)
Description Input Output Sample Input 4 3 51 9 27 6 5 12 33 13 41 2 3 21 1 21 4 20 2 11 1 21 4 2 Sam ...
- BZOJ3052: [wc2013]糖果公园【树上带修莫队】
Description Input Output Sample Input Sample Input Sample Output 84 131 27 84 HINT 思路 非常模板的树上带修莫队 真的 ...
- BZOJ3052 [wc2013] 糖果公园 【树上莫队】
树上莫队和普通的序列莫队很像,我们把树进行dfs,然后存一个长度为2n的括号序列,就是一个点进去当作左括号,出来当作右括号,然后如果访问从u到v路径,我们可以转化成括号序列的区间,记录x进去的时候编号 ...
- 【树上莫队】【带修莫队】bzoj3052 [wc2013]糖果公园
#include<cstdio> #include<cstring> #include<algorithm> #include<cmath> using ...
- 【BZOJ3052】[wc2013]糖果公园 带修改的树上莫队
[BZOJ3052][wc2013]糖果公园 Description Input Output Sample Input Sample Input Sample Output 84 131 27 84 ...
- [BZOJ3052][UOJ#58][WC2013]糖果公园
[BZOJ3052][UOJ#58][WC2013]糖果公园 试题描述 Candyland 有一座糖果公园,公园里不仅有美丽的风景.好玩的游乐项目,还有许多免费糖果的发放点,这引来了许多贪吃的小朋友来 ...
- 【BZOJ-3052】糖果公园 树上带修莫队算法
3052: [wc2013]糖果公园 Time Limit: 200 Sec Memory Limit: 512 MBSubmit: 883 Solved: 419[Submit][Status] ...
- bzoj 3052: [wc2013]糖果公园 带修改莫队
3052: [wc2013]糖果公园 Time Limit: 250 Sec Memory Limit: 512 MBSubmit: 506 Solved: 189[Submit][Status] ...
- 洛谷 P4074 [WC2013]糖果公园 解题报告
P4074 [WC2013]糖果公园 糖果公园 树上待修莫队 注意一个思想,dfn序处理链的方法,必须可以根据类似异或的东西,然后根据lca分两种情况讨论 注意细节 Code: #include &l ...
随机推荐
- SSH远程会话管理工具 - screen使用教程
一.screen命令是什么? Screen是一个可以在多个进程之间多路复用一个物理终端的全屏窗口管理器.Screen中有会话的概念,用户可以在一个screen会话中创建多个screen窗口,在每一个s ...
- arcgis api for js共享干货系列之一自写算法实现地图量算工具
众所周知,使用arcgis api for js实现地图的量算工具功能,无非是调用arcgisserver的Geometry服务(http://localhost:6080/arcgis/rest/s ...
- 为什么我会认为SAP是世界上最好用最牛逼的ERP系统,没有之一?
为什么我认为SAP是世界上最好用最牛逼的ERP系统,没有之一?玩过QAD.Tiptop.用友等产品,深深觉得SAP是贵的有道理! 一套好的ERP系统,不仅能够最大程度承接适配企业的管理和业务流程,在技 ...
- 提交本地项目到github服务器
已经完成的本地项目 提交到github 并不是按照先在github上创建一个仓库 然后clone下来的顺序 1.在github上创建仓库 2.在本地项目初始化git仓库 $ git init 3.添加 ...
- 阶段一:AsyncTask的三个属性值和四个步骤
“阶段一”是指我第一次系统地学习Android开发.这主要是对我的学习过程作个记录. 最近学到用AsyncTask来处理有关网络的操作.虽然代码看上去不是很复杂,但仍有很多地方有疑惑.所以研读了一下A ...
- 表单中Readonly和Disabled的区别
1.readonly是要锁定这个控件,通过在界面上无法修改他(但是通过javascript可以修改他). 2.disabled和readonly有相同的地方也是可以锁定这个控件用户不能改变他的值,但是 ...
- css权重计算方法浅谈
在这之前只知道css权重的皮毛,比如说:行内权重比头部权重高,头部比外部样式权重高----工作中才知道真正理解css权重重要性.理解权重了才能写出来最优css选择器来.对后面学习less,scss有很 ...
- java遍历给定目录,树形结构输出所有文件,包括子目录中的文件
(转自:http://blog.csdn.net/gangwazi0525/article/details/7569701) import java.io.File; public class Rea ...
- windows中,端口查看&关闭进程及Kill使用
测试过程中遇到的问题,杂记一: 1.netstat -ano | findstr "8001" 查看端口8001被哪个进程占用:由下图可以看出,被进程为3736的占用 2.查 ...
- S5PV210_时钟系统
1.S5PV210的时钟获得:外部晶振+内部时钟发生器+内部PLL产生高频时钟+内部分频器分频 S5PV210外部有4个W晶振接口,可以根据需要来决定在哪里接晶振.接了晶振之后上电相应的模块就能产生振 ...