时光真疯狂, 我一路执迷于匆忙.

那么,你是想放手一搏,还是等到年华老去,心中充满遗憾,孤独地迈向黄泉路?...

【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);
}
阅读更多
版权声明:既然是蒟蒻写的文,那么各位大爷就将就着看吧~跑的比西方记者慢多了233 https://blog.csdn.net/alan_cty/article/details/51540529
个人分类: 暴力 kmp
想对作者说点什么? 我来说一句

jzoj 3427 归途征程

jzoj 3427

ssl_lyy ssl_lyy

2017-08-18 21:56:44

阅读数:126

没有更多推荐了,返回首页

不良信息举报

【NOIP2013模拟】归途与征程

最多只允许输入30个字

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭