bzoj[3881]Divljak(dfs序+树状数组+fail树)
这道题利用了fail树的神奇性质————父节点为其子节点的前缀
先对Alice的集合建一个fail树，
Bob每插入一个串，都将串在自动机上经过的点在树上打上标记(+1)
每次查询的答案就是询问串的结束节点的子树的贡献
所以还需要用到树状数组来维护dfs序
因为Bob的一个串至多只能对Alice的某一个串做出1的贡献，所以不能单纯加和
要对树链取并
可以将一个Bob串在fail树上经过的所有点按入栈顺序排序，相邻每两个点的lca贡献-1
可以证明，这样会使每次某个点的子树内没有重复贡献
1 #include<queue> 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 using namespace std; 6 int n,m,cnt,tot,tp,num; 7 char s[2000005]; 8 bool ed[2000005]; 9 int vis[2000005]; 10 int ed_pos[100005]; 11 int head[2000005]; 12 int fa[2000005]; 13 int dep[2000005]; 14 int grand[2000005]; 15 int siz[2000005]; 16 int son[2000005]; 17 int ad[4000005]; 18 int stk[4000005]; 19 int tmp[2000005]; 20 struct Trie{ 21 int son[26]; 22 int fail; 23 }tr[2000005]; 24 struct Edge{ 25 int fr; 26 int to; 27 int nxt; 28 }edge[2000005]; 29 struct node{ 30 int in; 31 int out; 32 }pos[2000005]; 33 int lowbit(int x){ 34 return x&(-x); 35 } 36 void init(){ 37 memset(head,-1,sizeof(head)); 38 } 39 void addx(int x,int v){ 40 while(x<=2*tot){ 41 ad[x]+=v; 42 x+=lowbit(x); 43 } 44 } 45 int getsum(int x){ 46 int ret=0; 47 while(x){ 48 ret+=ad[x]; 49 x-=lowbit(x); 50 } 51 return ret; 52 } 53 void addedge(int u,int v){ 54 edge[cnt].fr=u; 55 edge[cnt].to=v; 56 edge[cnt].nxt=head[u]; 57 head[u]=cnt++; 58 } 59 void insert(char h[],int mrk){ 60 int now=0; 61 int len=strlen(h+1); 62 for(int i=1;i<=len;i++){ 63 int k=h[i]-'a'; 64 if(!tr[now].son[k])tr[now].son[k]=++tot; 65 now=tr[now].son[k]; 66 } 67 ed[now]=true; 68 ed_pos[mrk]=now; 69 } 70 void getfail(){ 71 queue<int>que; 72 for(int i=0;i<26;i++){ 73 if(tr[0].son[i]){ 74 addedge(0,tr[0].son[i]); 75 que.push(tr[0].son[i]); 76 } 77 } 78 while(!que.empty()){ 79 int u=que.front(); 80 que.pop(); 81 for(int i=0;i<26;i++){ 82 if(tr[u].son[i]){ 83 tr[tr[u].son[i]].fail=tr[tr[u].fail].son[i]; 84 addedge(tr[tr[u].fail].son[i],tr[u].son[i]); 85 que.push(tr[u].son[i]); 86 } 87 else tr[u].son[i]=tr[tr[u].fail].son[i]; 88 } 89 } 90 } 91 void dfs1(int u){ 92 siz[u]=1;stk[++tp]=u;pos[u].in=tp; 93 for(int i=head[u];i!=-1;i=edge[i].nxt){ 94 int v=edge[i].to; 95 if(v==fa[u])continue; 96 fa[v]=u;dep[v]=dep[u]+1; 97 dfs1(v);siz[u]+=siz[v]; 98 if(siz[v]>siz[son[u]]||(!son[u]))son[u]=v; 99 } 100 stk[++tp]=u;pos[u].out=tp; 101 } 102 void dfs2(int u){ 103 if(u!=son[fa[u]])grand[u]=u; 104 else grand[u]=grand[fa[u]]; 105 for(int i=head[u];i!=-1;i=edge[i].nxt){ 106 int v=edge[i].to; 107 if(v==fa[u])continue; 108 dfs2(v); 109 } 110 } 111 int lca(int x,int y){ 112 while(grand[x]!=grand[y]){ 113 if(dep[grand[x]]<dep[grand[y]])swap(x,y); 114 x=fa[grand[x]]; 115 } 116 if(dep[x]>dep[y])swap(x,y); 117 return x; 118 } 119 int cmp(int a,int b){ 120 return pos[a].in<pos[b].in; 121 } 122 void match(char b[],int tim){ 123 int now=0; 124 int len=strlen(s+1); 125 for(int i=1;i<=len;i++){ 126 int k=b[i]-'a'; 127 now=tr[now].son[k]; 128 if(vis[now]==tim)continue; 129 vis[now]=tim;tmp[++num]=now; 130 } 131 sort(tmp+1,tmp+num+1,cmp); 132 for(int i=1;i<=num;i++)addx(pos[tmp[i]].in,1); 133 for(int i=1;i<num;i++){ 134 int f=lca(tmp[i],tmp[i+1]); 135 addx(pos[f].in,-1); 136 } 137 } 138 int main(){ 139 init(); 140 scanf("%d",&n); 141 for(int i=1;i<=n;i++){ 142 scanf("%s",s+1); 143 insert(s,i); 144 } 145 getfail();dfs1(0);dfs2(0); 146 scanf("%d",&m); 147 for(int i=1;i<=m;i++){ 148 int ac; 149 scanf("%d",&ac); 150 if(ac==1){ 151 scanf("%s",s+1); 152 num=0;match(s,i); 153 }else{ 154 int p; 155 scanf("%d",&p); 156 int ans=getsum(pos[ed_pos[p]].out)-getsum(pos[ed_pos[p]].in-1); 157 printf("%d\n",ans); 158 } 159 } 160 return 0; 161 }
