loj#6517. 「雅礼集训 2018 Day11」字符串(回滚莫队)
模拟赛的时候纯暴力竟然骗了\(70\)分……
首先对于一堆\(g\)怎么计算概率应该很好想,用总的区间数减去不合法的区间数就行了,简而言之对\(g\)排个序,每一段长为\(d\)的连续序列的区间有\(\frac{d(d+1)}{2}\),那么对于每一个\([g_{i-1}+1,g_{i}-1]\)的区间,把它能贡献的子区间减去就行了
其次它是询问,而且不强制在线,那么很容易想到莫队
然而单纯的莫队不可做,因为很有可能加一个字符串会导致好几个长度满足条件,这样一次转移就不是\(O(1)\)的了
这个问题也不难解决,我们把所有的字符串首尾相接看做一个字符串,那么每一次转移只会加上一个前缀,转移就可以做到\(O(1)\)了。对于\(g\)开一个\(set\),加入数字的时候查一下前驱\(pre\)和后继\(next\),把\([next+1,pre-1]\)的贡献删掉,同时把\([next+1,g-1]\)和\([g+1,pre-1]\)的贡献加上去。删除数字的时候同理
\(ok\)到了现在交上去发现还是\(T\),因为\(set\)有一个\(\log\),还是不够
我们发现\(set\)可以用链表优化,链表的删除是\(O(1)\)的,撤销也是\(O(1)\)的
但是莫队除了删除还有插入,这个时候就需要用到一个叫做回滚莫队的黑科技了
简单来说,我们在对询问排序的时候先按左端点所在块从小到大排序,再按右端点(不是所在块!)从大到小排序。对于左端点在同一块内的询问,右端点是递减的,那么每一次删除的复杂度是\(O(1)\),暴力删除右端点的贡献是\(O(n)\),这一部分的复杂度是\(O(n\sqrt{n})\)。对于左端点,每一次都从该块的最左端一直删除到左端点,每一次复杂度为\(O(\sqrt{n})\),这一部分的府再度也是\(O(n\sqrt{n})\)
所以最后的复杂度是\(O(n\sqrt{n})\)
//minamoto
#include<bits/stdc++.h>
#define R register
#define inf 0x3f3f3f3f
#define ll long long
#define fp(i,a,b) for(R int i=a,I=b+1;i<I;++i)
#define fd(i,a,b) for(R int i=a,I=b-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
template<class T>inline bool cmin(T&a,const T&b){return a>b?a=b,1:0;}
using namespace std;
//char buf[1<<21],*p1=buf,*p2=buf;
//inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;}
#define getc() getchar()
int read(){
R int res,f=1;R char ch;
while((ch=getc())>'9'||ch<'0'){(ch=='-')&&(f=-1);if(ch==EOF)return EOF;}
for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');
return res*f;
}
int read(char *s){
R int len=0;R char ch;while((ch=getc())>'z'||ch<'a');
while(s[++len]=ch,((ch=getc())>='a'&&ch<='z'));return s[len+1]=0,len;
}
char sr[1<<21],z[20];int K=-1,Z=0;
inline void Ot(){fwrite(sr,1,K+1,stdout),K=-1;}
void print(R ll x){
if(K>1<<20)Ot();if(x<0)sr[++K]='-',x=-x;
while(z[++Z]=x%10+48,x/=10);
while(sr[++K]=z[Z],--Z);sr[++K]='\n';
}
const int N=3e5+5;
int a[N],len[N],ch[N][26],pos[N],v[N],dep[N];
int bl[N],g[N],ls[N],rs[N],st[N],cnt[N];
struct node{
int l,r,id;
inline bool operator <(const node &b)const{return bl[l]==bl[b.l]?r>b.r:bl[l]<bl[b.l];}
}q[N];
char s[N];ll F[N],ans[N],res,d,qwq;
int n,m,tot,top,num,A,B,C,mx,l,r,S;
inline bool ck(R ll f,R int len){return f?f*B+1ll*A*len>=C:0;}
inline ll calc(R int x){return 1ll*x*(x+1)>>1;}
void ins(char *s,int vva,int len){
int p=0;
fp(i,1,len){
int c=s[i]-'a';
if(!ch[p][c])ch[p][c]=++num,dep[num]=dep[p]+1;
p=ch[p][c],pos[++tot]=p,v[tot]=vva;
}
}
void del(R int i){
int x=pos[i];
if(ck(F[x],dep[x])&&!ck(F[x]-v[i],dep[x])&&!--cnt[g[dep[x]]]){
int k=g[dep[x]],l=ls[k],r=rs[k];
ls[r]=l,rs[l]=r,res-=calc(r-k-1),res-=calc(k-l-1),res+=calc(r-l-1);
}F[x]-=v[i];
}
void add(R int i){
int x=pos[i];
if(!ck(F[x],dep[x])&&ck(F[x]+v[i],dep[x])&&!cnt[g[dep[x]]]++){
int k=g[dep[x]],l=ls[k],r=rs[k];
ls[r]=rs[l]=k,res-=calc(r-l-1),res+=calc(r-k-1),res+=calc(k-l-1);
}F[x]+=v[i];
}
int main(){
// freopen("testdata.in","r",stdin);
// freopen("string.in","r",stdin);
// freopen("string.out","w",stdout);
n=read(),A=read(),B=read(),C=read();fp(i,1,n)a[i]=read();
fp(i,1,n)cmax(mx,len[i]=read(s)),ins(s,a[i],len[i]),len[i]+=len[i-1];
S=sqrt(tot);fp(i,1,tot)bl[i]=(i-1)/S+1;
fp(i,1,mx)g[i]=read();
m=read(),qwq=calc(mx);
fp(i,1,m){
l=read(),r=read();
q[i].l=len[l-1]+1,q[i].r=len[r],q[i].id=i;
}
sort(q+1,q+1+m);
for(R int i=1,x=pos[i];i<=tot;F[x]+=v[i],x=pos[++i])
if(!ck(F[x],dep[x])&&ck(F[x]+v[i],dep[x])&&!cnt[g[dep[x]]]++)
st[++top]=g[dep[x]];
st[++top]=0,st[++top]=mx+1,sort(st+1,st+1+top);
fp(i,1,top)ls[st[i]]=st[i-1],rs[st[i]]=st[i+1];
fp(i,2,top)res+=calc(st[i]-st[i-1]-1);
for(R int i=1,j=S,t=1,l=i,r=tot;t<=m&&i<=tot;i+=S,j+=S){
while(t<=m&&q[t].l>=i&&q[t].l<=j){
while(r>q[t].r)del(r--);while(l<q[t].l)del(l++);
ans[q[t].id]=qwq-res;while(l>i)add(--l);++t;
}while(r<tot)add(++r);while(l<=j)del(l++);
}fp(i,1,m){
d=__gcd(qwq,ans[i]);
print(ans[i]/d),sr[K]='/',print(qwq/d);
}
return Ot(),0;
}