关闭

【NOIP2013模拟】归途与征程

标签: NOIP2013模拟kmp归途与征途暴力模拟
707人阅读 评论(0) 收藏 举报
分类:

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);
}
1
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:143905次
    • 积分:4817
    • 等级:
    • 排名:第6591名
    • 原创:335篇
    • 转载:4篇
    • 译文:0篇
    • 评论:77条
    博客公告
    我觉得我已经是一碗废蘑菇汤了……
    博客专栏
    文章分类
    最新评论