题目描述

瞎做

显然是AC自动机dp。
小的直接dp大的矩阵乘法。

#include<cstdio>
#include<algorithm>
#include<cstring>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
typedef long long ll;
const int maxn=200+10,mo=1000000007;
char a[maxn][maxn],b[maxn][maxn];
int l1[maxn],l2[maxn];
int g[maxn][26],fail[maxn],next[maxn][26],trans[maxn][maxn],f[maxn][maxn],dl[maxn],dy[maxn],sta[80];
int o[maxn][maxn],an[maxn][maxn],dis[maxn][maxn],wd[maxn][maxn];
bool bz[maxn];
int i,j,k,l,t,n,m,len,tot,top,head,tail,root,last,now,ans;
int insert(int x,int c){
    if (!g[x][c]) g[x][c]=++tot;
    return g[x][c];
}
void work(){
    fo(i,1,tot)
        fo(j,1,n)
            if (trans[i][j]!=-1){                t=trans[i][j];
                if (l1[j]==1){                    //dis[i][t]++;
                    dis[i+tot][t+tot]++;
                }
                else dis[i][t+tot]++;
            }
    fo(i,1,tot) dis[i+tot][i]++;
    wd[1][tot+1]=1;
    while (len){        sta[++top]=len%2;
        len/=2;
    }
    fo(i,1,2*tot) an[i][i]=1;
    while (top){        fo(i,1,2*tot)
            fo(j,1,2*tot)
                o[i][j]=0;
        fo(k,1,2*tot)
            fo(i,1,2*tot)
                fo(j,1,2*tot)
                    o[i][j]=(o[i][j]+(ll)an[i][k]*an[k][j]%mo)%mo;
        fo(i,1,2*tot)
            fo(j,1,2*tot)
                an[i][j]=o[i][j];
        if (sta[top]){            fo(i,1,2*tot)
                fo(j,1,2*tot)
                    o[i][j]=0;
            fo(k,1,2*tot)
                fo(i,1,2*tot)
                    fo(j,1,2*tot)
                        o[i][j]=(o[i][j]+(ll)an[i][k]*dis[k][j]%mo)%mo;
            fo(i,1,2*tot)
                fo(j,1,2*tot)
                    an[i][j]=o[i][j];
        }
        top--;
    }
    fo(i,1,2*tot)
        fo(j,1,2*tot)
            o[i][j]=0;
    fo(k,1,2*tot)
        fo(i,1,2*tot)
            fo(j,1,2*tot)
                o[i][j]=(o[i][j]+(ll)wd[i][k]*an[k][j]%mo)%mo;
    fo(i,1,2*tot)
        fo(j,1,2*tot)
            an[i][j]=o[i][j];
    fo(j,tot+1,2*tot) ans=(ans+an[1][j])%mo;
    printf("%d\n",ans);
}
int main(){
    freopen("sorcery.in","r",stdin);freopen("sorcery.out","w",stdout);
    scanf("%d%d%d",&n,&m,&len);
    fo(i,1,n)
        scanf("%s",a[i]+1),l1[i]=strlen(a[i]+1);
    fo(i,1,m)
        scanf("%s",b[i]+1),l2[i]=strlen(b[i]+1);
    tot=root=1;
    fo(i,1,m){        last=1;
        fo(j,1,l2[i]) last=insert(last,b[i][j]-'a');
        dy[i]=last;
        bz[last]=1;
    }
    fo(i,0,25)
        if (g[1][i]) next[1][i]=g[1][i];else next[1][i]=1;
    fo(i,0,25)
        if (g[1][i]){            dl[++tail]=g[1][i];
            fail[g[1][i]]=1;
            j=g[1][i];
            fo(k,0,25)
                if (g[j][k]) next[j][k]=g[j][k];else next[j][k]=next[1][k];
        }
    while (head<tail){        now=dl[++head];
        bz[now]|=bz[fail[now]];
        fo(i,0,25)
            if (g[now][i]){                j=fail[now];
                while (j>1&&!g[j][i]) j=fail[j];
                if (g[j][i]) j=g[j][i];
                fail[g[now][i]]=j;
                j=g[now][i];
                fo(k,0,25)
                    if (g[j][k]) next[j][k]=g[j][k];else next[j][k]=next[fail[j]][k];
                dl[++tail]=g[now][i];
            }
    }
    fo(i,1,tot){        if (bz[i]) continue;
        fo(j,1,n){            last=i;
            fo(k,1,l1[j]+1){                if (k>l1[j]) break;
                last=next[last][a[j][k]-'a'];
                if (bz[last]) break;
            }
            if (k<=l1[j]) trans[i][j]=-1;else trans[i][j]=last;
        }
    }
    if (len>100){        work();
        return 0;
    }
    f[0][1]=1;
    fo(i,0,len-1)
        fo(j,1,tot)
            if (f[i][j])
                fo(k,1,n)
                    if (i+l1[k]<=len&&trans[j][k]!=-1) (f[i+l1[k]][trans[j][k]]+=f[i][j])%=mo;
    fo(j,1,tot) (ans+=f[len][j])%=mo;
    printf("%d\n",ans);
}

[bzoj4861]魔法咒语相关推荐

【AC自动机】【矩阵加速】BZOJ4861魔法咒语
分析: 巨恶心的一题双代码题. 对于前半部分,把禁止出现的字符串建一颗AC自动机. 然后枚举每个位置用了某个模板串后转移到哪里. 然后直接DP即可. 但是对于后半部分数据,则必须写一个矩阵加速... ...
BZOJ4861 魔法咒语【AC自动机】【DP】【矩阵快速幂】
题意: 用n个基本词汇拼出长度为L的字符串,每个基本词汇可以用0次,1次或多次,要求拼出的字符串不能出现m个禁忌词汇中的任何一个,求方案数 mod 1e9+7 60% :1<=N,M<=5 ...
【BZOJ4861】[Beijing2017]魔法咒语矩阵乘法+AC自动机+DP
[BZOJ4861][Beijing2017]魔法咒语题意:别看BZ的题面了,去看LOJ的题面吧~ 题解:显然,数据范围明显的分成了两部分:一个是L很小,每个基本词汇长度未知:一个是L很大,每个基本 ...
bzoj4861 / P3715 [BJOI2017]魔法咒语
P3715 [BJOI2017]魔法咒语 AC自动机+dp+矩阵乘法常规思路是按基本串建立AC自动机然鹅这题是按禁忌串建立AC自动机对后缀是禁忌的点以及它的失配点做上标记$(a[i].ed)$, ...
【BZOJ4861】【BJOI2017】—魔法咒语（AC自动机+矩阵快速幂优化dp）
传送门当 l ≤ 100 l\le 100 l≤100时显然的自动机上 d p dp dp就完了当 l ≤ 1 e 8 l\le1e8 l≤1e8时直接 d p dp dp显然是不行的但是发现 ...
[BZOJ4861][Beijing2017]魔法咒语 AC自动机+动态规划+矩阵快速幂
对忌讳词语构建AC自动机 Fi,j F_{i,j}表示长度为 i i,匹配到AC自动机第jj位的合法串方案数当状态数少的时候用矩阵转移 #include <bits/stdc++.h> ...
编程语言居然是魔法咒语！
作为一个本职是程序员的魔幻爱好人士,THE DAY,我终于意识到了! 编程语言,本质上就是一种现代魔法咒语! 如果你不懂魔法的麻瓜,没关系,让我一点点来揭开这个世纪大奥妙. 首先我们从最简单的HTML ...
【记录】ChatGPT｜图片预览魔法咒语魔改，使用 ChatGPT 返回大量可以跳转的链接
很早的时候,我已经留意到 ChatGPT 会以返回图片的 markdown 格式来显示图片,很可能拥有一定的图片上传功能,但是它往往会显示得有些问题.一些代码图片之类的或者风景图什么的都不是很会. 但 ...
P3715 [BJOI2017]魔法咒语
P3715 [BJOI2017]魔法咒语用基本词汇组成$L$长度的单词,其中不能包含禁忌词汇用禁忌词汇建强大的$tire$图解决: 分类讨论,$L<=100$用普通dp暴力在\ ...

[bzoj4861]魔法咒语

题目描述

瞎做

[bzoj4861]魔法咒语相关推荐

最新文章

热门文章