bzoj 2754: [SCOI2012]喵星球上的点名

Description

a180285幸运地被选做了地球到喵星球的留学生。他发现喵星人在上课前的点名现象非常有趣。   假设课堂上有N个喵星人，每个喵星人的名字由姓和名构成。喵星球上的老师会选择M个串来点名，每次读出一个串的时候，如果这个串是一个喵星人的姓或名的子串，那么这个喵星人就必须答到。 然而，由于喵星人的字码过于古怪，以至于不能用ASCII码来表示。为了方便描述，a180285决定用数串来表示喵星人的名字。

现在你能帮助a180285统计每次点名的时候有多少喵星人答到，以及M次点名结束后每个喵星人答到多少次吗？  

Input

 

现在定义喵星球上的字符串给定方法：

先给出一个正整数L，表示字符串的长度，接下来L个整数表示字符串的每个字符。

输入的第一行是两个整数N和M。

接下来有N行，每行包含第i 个喵星人的姓和名两个串。姓和名都是标准的喵星球上的

字符串。

接下来有M行，每行包含一个喵星球上的字符串，表示老师点名的串。

Output

 

对于每个老师点名的串输出有多少个喵星人应该答到。

然后在最后一行输出每个喵星人被点到多少次。

Sample Input

2 3
6 8 25 0 24 14 8 6 18 0 10 20 24 0
7 14 17 8 7 0 17 0 5 8 25 0 24 0
4 8 25 0 24
4 7 0 17 0
4 17 0 8 25

Sample Output

2
1
0
1 2
【提示】
事实上样例给出的数据如果翻译成地球上的语言可以这样来看
2 3
izayoi sakuya
orihara izaya
izay
hara
raiz

HINT

【数据范围】 

 对于30%的数据，保证： 

1<=N,M<=1000，喵星人的名字总长不超过4000，点名串的总长不超过2000。

对于100%的数据，保证：

1<=N<=20000，1<=M<=50000，喵星人的名字总长和点名串的总长分别不超过100000，保证喵星人的字符串中作为字符存在的数不超过10000。

Source

 

这个题是叉姐考试的一道AC自动机题，当场这题有n种做法都可以AC

这题我是看了洛谷标签上有莫队才做的。。。

首先还是要拼成一个串(把询问串也拼起来)，然后记录每个后缀的开头属于哪个串，对于询问串还要记下串的开头在串中的位置st[i];

对于询问的话，如果询问串是别的串的子串的话，那么以询问串的开头为开头的后缀与以别的串的字母开头的后缀的lcp>=该串长度；

因为靠的越近lcp越大，所以我们只需要从rank[st[i]]，向前后扩展直到height<le[i]为止，那么中间的l--r的后缀就是满足条件的后缀，具体实现和TJOI 单词一样

第一问相当于询问l--r中有多少个不同的数，是莫队的基础套路题，开个桶即可；

第二问相当于l--r的区间加法，因为不能重复加，所以一开始是想线段树直接区间加法，但一个串每个询问只加一次不好保证。

其实这个东西可以直接在莫队的过程中实现：

每新出现一个数，就默认它在后面的询问区间都出现，即加上(m-i+1);

每有一个数消失，就默认它在后面的询问区间都不出现，然后把原来多加的减去，即减去(m-i+1);

 

几个傻逼错误：1.block=sqrt(len)写成=sqrt(n)了；

2.从rank[st[i]]扩展的时候，左边可能扩展到0,右边如果是while(height[r]>=le[i]) r++;最后的r是不满足条件的，要最后r--；

代码能力真的越来越不行了。。。

(upd:我真的是石乐志。。。都while暴跳左右两边了，就是裸暴力的复杂度了，我就说我的莫队怎么跑得这么慢,跟个暴力似的QAQ。。找l和r的话二分长度即可logn找完。。。)

// MADE BY QT666
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<iostream>
#include<cstring>
#define RG register
using namespace std;
typedef long long ll;
const int N=1000050;
const int Inf=19260817;
int gi(){
    int x=0,flag=1;
    char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-') flag=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*flag;
}
int sa[N],y[N],rk,rnk[N],height[N],bl[N],len,a[N],n,m,st[N],le[N],cnt[N],ans[N],ans2[N],pos[N],block;
struct data{
    int fir,sec,id;
}x[N];
struct ac{
    int l,r,id;
}q[N];
bool cmp(const data &a,const data &b){
    if(a.fir==b.fir) return a.sec<b.sec;
    else return a.fir<b.fir;
}
void work2(){
    rk=1;y[x[1].id]=rk;
    for(RG int i=2;i<=len;i++){
        if(x[i-1].fir!=x[i].fir||x[i-1].sec!=x[i].sec) rk++;
        y[x[i].id]=rk;
    }
}
void work(){
    sort(x+1,x+1+len,cmp);work2();
    for(RG int i=1;i<=len;i<<=1){
        for(RG int j=1;j+i<=len;j++) x[j].fir=y[j],x[j].sec=y[j+i],x[j].id=j;
        for(RG int j=len-i+1;j<=len;j++) x[j].fir=y[j],x[j].sec=0,x[j].id=j;
        sort(x+1,x+1+len,cmp);work2();
        if(rk==len) break;
    }
}
void get_height(){
    int kk=0;for(RG int i=1;i<=len;i++) rnk[sa[i]]=i;
    for(RG int i=1;i<=len;i++){
        if(kk) kk--;
        int j=sa[rnk[i]-1];
        while(a[i+kk]==a[j+kk]) kk++;
        height[rnk[i]]=kk;
    }
}
bool cmp2(ac a,ac b){
    if(pos[a.l]==pos[b.l]) return a.r<b.r;
    return pos[a.l]<pos[b.l];
}
void Modui(){
    int num=0;
    for(RG int i=1,l=1,r=0;i<=m;i++){
        while(l>q[i].l) {l--,cnt[bl[sa[l]]]++;if(cnt[bl[sa[l]]]==1&&bl[sa[l]]<=n&&bl[sa[l]]) num++,ans2[bl[sa[l]]]+=(m-i+1);}
        while(r<q[i].r) {r++,cnt[bl[sa[r]]]++;if(cnt[bl[sa[r]]]==1&&bl[sa[r]]<=n&&bl[sa[r]]) num++,ans2[bl[sa[r]]]+=(m-i+1);}
        while(l<q[i].l) {cnt[bl[sa[l]]]--;if(!cnt[bl[sa[l]]]&&bl[sa[l]]<=n&&bl[sa[l]]) num--,ans2[bl[sa[l]]]-=(m-i+1);l++;}
        while(r>q[i].r) {cnt[bl[sa[r]]]--;if(!cnt[bl[sa[r]]]&&bl[sa[r]]<=n&&bl[sa[r]]) num--,ans2[bl[sa[r]]]-=(m-i+1);r--;}
        ans[q[i].id]=num;
    }
}
int main(){
    n=gi(),m=gi();
    for(RG int i=1;i<=n;i++){
        int l=gi();
        for(int j=1;j<=l;j++) a[++len]=gi(),bl[len]=i;
        int l2=gi();
        a[++len]=Inf+i;bl[len]=N-1;
        for(int j=1;j<=l2;j++) a[++len]=gi(),bl[len]=i;
        a[++len]=Inf;bl[len]=N-1;
    }
    for(RG int i=1;i<=m;i++){
        int l=gi();st[i]=len+1;le[i]=l;
        for(int j=1;j<=l;j++) a[++len]=gi(),bl[len]=n+i;
        a[++len]=Inf;bl[len]=N-1;
    }
    for(RG int i=1;i<=len;i++){
        x[i].id=i;x[i].fir=x[i].sec=a[i];
    }
    work();for(int i=1;i<=len;i++) sa[y[i]]=i;
    get_height();
    for(RG int i=1;i<=m;i++){
        int g=rnk[st[i]],l=g,r=g;
        while(height[l]>=le[i]) l--;
        while(height[r+1]>=le[i]) r++;
        q[i].l=l,q[i].r=r;q[i].id=i;
    }
    block=(int)sqrt(len);
    for(RG int i=1;i<=len;i++) pos[i]=(i-1)/block+1;
    sort(q+1,q+1+m,cmp2);Modui();
    for(RG int i=1;i<=m;i++) printf("%d\n",ans[i]);
    for(RG int i=1;i<n;i++) printf("%d ",ans2[i]);
    printf("%d",ans2[n]);
    return 0;
}