【NOIP2013模拟】归途与征程

Description

给出两个字符串，A和B，|A|=n，|B|=m，并且A中有许多*——通配符。这个东西可以匹配任意的字符串，包括空串。现在问你，B的循环同构的字符串中，有多少个可以和A匹配？
n<=100,m<=10^5

Solution

什么东东？
还以为是扩k呢。。。
实际上只是一个暴力。
我们把A串分成一段一段的，然后枚举起点，暴力往后跑，用KMP判断第一个能匹配某段的位置，然后就能拿到80分。
再然后，100分算法只需要用一个类似于前缀和的东西，处理处每个位值后面第一个能个某个串匹配的位置在哪里，然后每次跑的复杂度就降到O(n)了。

Code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define N 100005
#define M 105
using namespace std;
char s[N],str[N*2],p[M][M],pre[M],suf[M];
int n,m,st,en,S,next[M][M],suc[N*2][M],len[M],ans,lenp,lens,cnt;
bool bz[M][N*2];
int main() {
    scanf("%s",s+1);scanf("%s",str+1);
    n=strlen(s+1);m=strlen(str+1);
    fo(i,m+1,2*m) str[i]=str[i-m];
    fo(i,1,n) {
        if (s[i]=='*') {st=i+1;break;}
        pre[++lenp]=s[i];
    }
    fd(i,n,1) {
        if (s[i]=='*') {en=i-1;break;}
        suf[++lens]=s[i];
    }
    fo(i,1,lens/2) swap(suf[i],suf[lens-i+1]);cnt=1;
    fo(i,st,en) {
        if (s[i]=='*') {
            if (len[cnt]) cnt++;
            continue;
        }
        p[cnt][++len[cnt]]=s[i]; 
    }
    fo(i,1,cnt) {
        int k=0;
        fo(j,2,len[i]) {
            while (k&&p[i][k+1]!=p[i][j]) k=next[i][k];
            next[i][j]=k+=(p[i][k+1]==p[i][j]);
        }
    }
    fo(i,1,cnt) {
        int k=0;
        fo(j,1,m*2) {
            while (k&&p[i][k+1]!=str[j]) k=next[i][k];
            k+=(p[i][k+1]==str[j]);
            if (k==len[i]) bz[i][j-len[i]+1]=1,k=next[i][k];
        }
    }
    fd(i,m*2-1,1) {
        fo(j,1,cnt) suc[i][j]=suc[i+1][j];
        fo(j,1,cnt) if (bz[j][i]) suc[i][j]=i;
    }
    fo(i,1,m) {
        bool pd=1;int k=i;
        fo(j,1,lenp) {
            if (pre[j]!=str[k]) {pd=0;break;}
            k++;
        }
        if (!pd) continue;
        if (k==i+m) {ans++;continue;}
        S=k;k=i+m-1;
        fd(j,lens,1) {
            if (suf[j]!=str[k]) {pd=0;break;}
            k--;
        }
        if (!pd) continue;int tot=1;
        while (tot<=cnt) {
            if (!suc[S][tot]||suc[S][tot]>k) {pd=0;break;}
            S=suc[S][tot]+len[tot++];
        }
        if (!pd) continue;
        if (tot>cnt&&S-1<=k) ans++;
    }
    printf("%d",ans);
}