大致题意: 一棵树,每个节点有一个人,他打水需要(T_i)的时间,每次询问两点之间所有人去打水的最小等待时间。
伪·强制在线
这题看似强制在线,但实际上,(pre mod 2)只能为(0)或(1),因此只要将两种情况下的答案都求出来,最后视情况输出即可。
这样就可以用离线算法乱搞了。
树上莫队+树状数组
其实,这道题是可以用树上莫队来做的。
考虑当前已有若干人要去打水,现在新加入一个人,他的打水时间为(x),求改变的贡献值。
显然,根据贪心的思想,我们应让打水时间越久的人越早打水。
则设有(rk)个人打水时间(>x),(val)为打水时间(le x)的人打水时间总和,则这个人所需等待时间为((rk+1)*x+val)。
不难发现,(rk)和(val)两个值是可以用树状数组来进行维护的。
这样这题就做完了。
代码
#include<bits/stdc++.h>
#define N 50000
#define LL long long
#define add(x,y) (e[++ee].nxt=lnk[x],e[lnk[x]=ee].to=y)
#define swap(x,y) (x^=y^=x^=y)
using namespace std;
int n,query_tot,key,ee,a[N+5],lnk[N+5];
struct edge
{
int to,nxt;
}e[(N<<1)+5];
class Class_FIO
{
private:
#define Fsize 100000
#define tc() (A==B&&(B=(A=Fin)+fread(Fin,1,Fsize,stdin),A==B)?EOF:*A++)
#define pc(ch) (void)(FoutSize<Fsize?Fout[FoutSize++]=ch:(fwrite(Fout,1,Fsize,stdout),Fout[(FoutSize=0)++]=ch))
int Top,FoutSize;char ch,*A,*B,Fin[Fsize],Fout[Fsize],Stack[Fsize];
public:
Class_FIO() {A=B=Fin;}
inline void read(int &x) {x=0;while(!isdigit(ch=tc()));while(x=(x<<3)+(x<<1)+(ch&15),isdigit(ch=tc()));}
inline void readc(char &x) {while(isspace(x=tc()));}
inline void writeln(LL x) {if(!x) return pc('0'),pc('
');while(x) Stack[++Top]=x%10+48,x/=10;while(Top) pc(Stack[Top--]);pc('
');}
inline void clear() {fwrite(Fout,1,FoutSize,stdout),FoutSize=0;}
}F;
class Class_CaptainMotao_on_Tree//树上莫队
{
private:
#define bp(x) (((x)-1)/S+1)
#define F5(x) ((v=Di.get_val(a[t=x]),op[t]^=1)?(res+=1LL*(T.QSize(Di.cnt)-T.QSize(v)+1)*a[t]+T.QSum(v),T.Add(v,1,a[t]),0):(T.Add(v,-1,-a[t]),res-=1LL*(T.QSize(Di.cnt)-T.QSize(v)+1)*a[t]+T.QSum(v)))
int Q,S,op[N+5];LL ans[N+5][2];char opt[N+5];
class Class_Dfser//DFS预处理
{
private:
#define LogN 16
int cnt,Depth[N+5],fa[N+5][LogN+5];
public:
int s[(N<<1)+5],I[N+5],O[N+5];
inline void Init(int x=1,int lst=0)
{
register int i;
for(s[I[x]=++cnt]=x,i=1;i<=LogN;++i) fa[x][i]=fa[fa[x][i-1]][i-1];
for(i=lnk[x];i;i=e[i].nxt) e[i].to^lst&&(Depth[e[i].to]=Depth[fa[e[i].to][0]=x]+1,Init(e[i].to,x),0);
s[O[x]=++cnt]=x;
}
inline int LCA(int x,int y)
{
register int i;
for(Depth[x]<Depth[y]&&swap(x,y),i=0;Depth[x]^Depth[y];++i) ((Depth[x]^Depth[y])&(1<<i))&&(x=fa[x][i]);
if(!(x^y)) return x;
for(i=LogN;~i;--i) fa[x][i]^fa[y][i]&&(x=fa[x][i],y=fa[y][i]);
return fa[x][0];
}
}D;
class Class_Discretization//离散化
{
private:
int data[N+5];
public:
int cnt;
inline void Init(int len,int *num)
{
for(register int i=1;i<=len;++i) data[i]=num[i];
sort(data+1,data+len+1),cnt=unique(data+1,data+len+1)-data-1;
}
inline int get_val(int x)
{
register int l=1,r=cnt,mid;
for(mid=l+r>>1;l<=r;mid=l+r>>1) data[mid]<x?l=mid+1:r=mid-1;
return l;
}
}Di;
class Class_TreeArray//树状数组
{
private:
#define lowbit(x) (x&-x)
int Size[N+5];LL Sum[N+5];
public:
int len;
inline void Add(int pos,int val1,int val2) {while(pos<=len) Size[pos]+=val1,Sum[pos]+=val2,pos+=lowbit(pos);}
inline int QSize(int pos,int res=0) {while(pos) res+=Size[pos],pos-=lowbit(pos);return res;}
inline LL QSum(int pos,LL res=0) {while(pos) res+=Sum[pos],pos-=lowbit(pos);return res;}
}T;
struct Query
{
int l,r,pos,op,bl,flag;
Query(int x=0,int y=0,int p=0,int o=0,int b=0,int f=0):l(x),r(y),pos(p),op(o),bl(b),flag(f){}
inline friend bool operator < (Query x,Query y) {return x.bl^y.bl?x.bl<y.bl:(x.bl&1?x.r<y.r:x.r>y.r);}
}q[(N<<1)+5];
public:
inline void Solve()
{
int i,x,y,z,k,s=1,t,v,L=1,R=0;LL res=0,lst_ans=0;
for(D.Init(),Di.Init(n,a),T.len=Di.cnt,S=sqrt(n),i=1;i<=query_tot;++i)//每个询问分pre%2的值存储两个
{
if(F.readc(opt[i]),F.read(k),opt[i]^'Q') {s=k;continue;}
D.I[x=k%n+1]>D.I[y=s]&&swap(x,y),q[++Q]=(z=D.LCA(x,y))^x?Query(D.O[x],D.I[y],i,0,bp(D.O[x]),z):Query(D.I[x],D.I[y],i,0,bp(D.I[x]),0);
D.I[x=(k+key)%n+1]>D.I[y=s]&&swap(x,y),q[++Q]=(z=D.LCA(x,y))^x?Query(D.O[x],D.I[y],i,1,bp(D.O[x]),z):Query(D.I[x],D.I[y],i,1,bp(D.I[x]),0);
}
for(sort(q+1,q+Q+1),i=1;i<=Q;++i)//处理询问
{
while(R<q[i].r) F5(D.s[++R]);while(L>q[i].l) F5(D.s[--L]);while(R>q[i].r) F5(D.s[R--]);while(L<q[i].l) F5(D.s[L++]);
q[i].flag&&F5(q[i].flag),ans[q[i].pos][q[i].op]=res,q[i].flag&&F5(q[i].flag);
}
for(i=1;i<=query_tot;++i) opt[i]^'C'&&(F.writeln(lst_ans=ans[i][lst_ans&1]),0);//输出答案
}
}C;
int main()
{
register int i,x,y;register char op;
for(F.read(n),F.read(query_tot),F.read(key),i=1;i<=n;++i) F.read(a[i]);
for(i=1;i<=n;++i) F.read(x),add(x,i);
return C.Solve(),F.clear(),0;
}