【JSOI2012】【BZOJ4327】玄武密码

Description

在美丽的玄武湖畔，鸡鸣寺边，鸡笼山前，有一块富饶而秀美的土地，人们唤作进香河。相传一日，一缕紫气从天而至，只一瞬间便消失在了进香河中。老人们说，这是玄武神灵将天书藏匿在此。
很多年后，人们终于在进香河地区发现了带有玄武密码的文字。更加神奇的是，这份带有玄武密码的文字，与玄武湖南岸台城的结构有微妙的关联。于是，漫长的破译工作开始了。
经过分析，我们可以用东南西北四个方向来描述台城城砖的摆放，不妨用一个长度为N的序列来描述，序列中的元素分别是‘E’，‘S’，‘W’，‘N’，代表了东南西北四向，我们称之为母串。而神秘的玄武密码是由四象的图案描述而成的M段文字。这里的四象，分别是东之青龙，西之白虎，南之朱雀，北之玄武，对东南西北四向相对应。
现在，考古工作者遇到了一个难题。对于每一段文字，其前缀在母串上的最大匹配长度是多少呢？
Input

第一行有两个整数，N和M，分别表示母串的长度和文字段的个数。
第二行是一个长度为N的字符串，所有字符都满足是E,S,W和N中的一个。
之后M行，每行有一个字符串，描述了一段带有玄武密码的文字。依然满足，所有字符都满足是E,S,W和N中的一个。
Output

输出有M行，对应M段文字。
每一行输出一个数，表示这一段文字的前缀与母串的最大匹配串长度。
Sample Input

7 3

SNNSSNS

NNSS

NNN

WSEE
Sample Output

4

2

0
HINT

对于100%的数据，N<=10^7，M<=10^5，每一段文字的长度<=100。

应上传者要求，此题不公开，如有异议，请提出.

Source

SB题
SAM+暴力匹配

#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<algorithm>#define MAXN 10000010using namespace std;int n,m,top;char ch[MAXN],t[110];struct sam{    int a[MAXN<<1][4],fa[MAXN<<1],len[MAXN<<1];    int p,q,np,nq,last,cnt;    sam()   {last=++cnt;}    void insert(int c)    {        p=last;np=last=++cnt;len[np]=len[p]+1;        while (!a[p][c]&&p) a[p][c]=np,p=fa[p];        if (!p) fa[np]=1;        else        {            q=a[p][c];            if (len[q]==len[p]+1)   fa[np]=q;            else            {                nq=++cnt;len[nq]=len[p]+1;memcpy(a[nq],a[q],sizeof(a[q]));                fa[nq]=fa[q];fa[np]=fa[q]=nq;                while (a[p][c]==q)  a[p][c]=nq,p=fa[p];            }        }    }}sam;int main(){    scanf("%d%d",&n,&m);    for (;top<n;)    {        char c=getchar();        if (c=='E') ch[++top]='A';        if (c=='N') ch[++top]='B';        if (c=='S') ch[++top]='C';        if (c=='W') ch[++top]='D';    }    for (int i=1;i<=n;i++)  sam.insert(ch[i]-'A');    for (int i=1;i<=m;i++)    {        scanf("%s",t);int l=strlen(t),x=1,maxn=0;        for (int j=0;j<l;j++)        {            if (t[j]=='E')  t[j]='A';            if (t[j]=='N')  t[j]='B';            if (t[j]=='S')  t[j]='C';            if (t[j]=='W')  t[j]='D';            if (sam.a[x][t[j]-'A']) maxn++,x=sam.a[x][t[j]-'A'];            else    break;        }        printf("%d\n",maxn);    }}