BZOJ4460 [Jsoi2013]广告计划

题目

蛤蛤题，数据范围这么小，于是怎么暴力怎么来

既然保证有解，那么答案不会超过\(L\)，于是我们枚举一个\(k\)，看其能否成为答案

有了\(k\)还不够，再枚举一个\(j,j<k\)，表示这个单词不在一整列的前缀的长度为多少，有了\(j,k\)我们就能确定每一行的子串是什么了，于是往SAM上跑一跑，看看\(\rm endpos\)有没有交就好了

复杂度大概是\(O(L^3+\frac{L^3n}{w})\)，非常不满

代码

#include<bits/stdc++.h>
#define re register
int m,n,lim,num;
int lst,cnt;
char S[205];
std::bitset<100> ans,vis[20005];
struct E{int v,nxt;}e[20005];
int len[20005],fa[20005],head[20005],son[20005][26];
inline void ins(int c,int o) {
    int p=++cnt,f=lst;lst=p;
    len[p]=len[f]+1,vis[p][o]=1;
    while(f&&!son[f][c]) son[f][c]=p,f=fa[f];
    if(!f) {fa[p]=1;return;}
    int x=son[f][c];
    if(len[f]+1==len[x]) {fa[p]=x;return;}
    int y=++cnt;len[y]=len[f]+1,fa[y]=fa[x],fa[x]=fa[p]=y;
    for(re int i=0;i<26;i++) son[y][i]=son[x][i];
    while(f&&son[f][c]==x) son[f][c]=y,f=fa[f];
}
inline void add(int x,int y){e[++num].v=y;e[num].nxt=head[x];head[x]=num;}
void dfs(int x) {
    for(re int i=head[x];i;i=e[i].nxt) dfs(e[i].v),vis[x]|=vis[e[i].v];
}
int main() {
    scanf("%d%d",&m,&n);cnt=1;
    for(re int i=1;i<=m;i++) {
        scanf("%s",S+1);lst=1;
        for(re int j=1;j<=n;j++) ins(S[j]-'a',j-1);
    }
    for(re int i=2;i<=cnt;i++) add(fa[i],i);dfs(1);
    scanf("%s",S+1);lim=strlen(S+1);
    for(re int k=1;k<=lim;k++) {
        int flag=0;
        for(re int j=0;j<k;++j) {
            for(re int i=0;i<n;i++) ans[i]=1;
            int ml=0;
            for(re int t=j+1;t<=k;t++) {
                int nw=1,cnt=0;
                for(re int p=t;p<=lim;p+=k)nw=son[nw][S[p]-'a'],++cnt;
                if(!ml) ml=cnt;if(ml==cnt) ans&=vis[nw];
                if(cnt<ml) ans&=(vis[nw]<<1);if(!(ans.any())) break;
            }
            for(re int t=1;t<=j;++t) {
                int cnt=0,nw=1;
                for(re int i=t;i<=lim;i+=k)++cnt,nw=son[nw][S[i]-'a'];
                if(cnt>ml) ans&=vis[nw];if(cnt==ml) ans&=(vis[nw]<<1);
                if(!(ans.any())) break;
            }
            flag=ans.any();
            if(flag) break;
        }
        if(flag) {printf("%d\n",k);break;}
    }
}