【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 10000010
using 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);
    }
}