您的位置:首页 > 其它

[hdu2243][ac自动机]考研路茫茫——单词情结

2018-02-07 10:57 387 查看
考研路茫茫——单词情结

Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)

Total Submission(s): 6436 Accepted Submission(s): 2209

Problem Description

背单词,始终是复习英语的重要环节。在荒废了3年大学生涯后,Lele也终于要开始背单词了。

一天,Lele在某本单词书上看到了一个根据词根来背单词的方法。比如”ab”,放在单词前一般表示”相反,变坏,离去”等。

于是Lele想,如果背了N个词根,那这些词根到底会不会在单词里出现呢。更确切的描述是:长度不超过L,只由小写字母组成的,至少包含一个词根的单词,一共可能有多少个呢?这里就不考虑单词是否有实际意义。

比如一共有2个词根 aa 和 ab ,则可能存在104个长度不超过3的单词,分别为

(2个) aa,ab,

(26个)aaa,aab,aac…aaz,

(26个)aba,abb,abc…abz,

(25个)baa,caa,daa…zaa,

(25个)bab,cab,dab…zab。

这个只是很小的情况。而对于其他复杂点的情况,Lele实在是数不出来了,现在就请你帮帮他。

Input

本题目包含多组数据,请处理到文件结束。

每组数据占两行。

第一行有两个正整数N和L。(0

sol::

在我前面做了那些ac自动机的基础上,这题考的就是矩乘了,这题坑点是模数,而且m<2^31,如果你+1的话就ll了

#include<cstdio>
#include<algorithm>
#include<string>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
using namespace std;
typedef unsigned long long ll;
int n,m;
inline int read()
{
char c;
int res,flag=0;
while((c=getchar())>'9'||c<'0') if(c=='-')flag=1;
res=c-'0';
while((c=getchar())>='0'&&c<='9') res=(res<<3)+(res<<1)+c-'0';
return flag?-res:res;
}
const int N=40;

int tot,l,top
[26],fail
;
int go
[26];
char sr[11];
bool mark
;
inline void insert(int &x,int y)
{
if(!x)
{
x=++tot;
fail[x]=0;
mark[x]=0;
for(int i=0;i<=25;++i) top[x][i]=go[x][i]=0;
}
if(y>l)
{
mark[x]=1;
return;
}
insert(go[x][sr[y]-'a'],y+1);
}
int q
;
inline void build()
{
int t=0,w,u,v;
q[w=1]=1;
for(int i=0;i<=25;++i) top[0][i]=go[0][i]=1;
while(t<w)
{
u=q[++t];
mark[u]=mark[u]|mark[fail[u]];

for(int ch=0;ch<=25;++ch)
if(go[u][ch])
{
top[u][ch]=go[u][ch];
q[++w]=go[u][ch];
v=fail[u];
while(!go[v][ch]) v=fail[v];
v=go[v][ch];
fail[go[u][ch]]=v;
}
else top[u][ch]=top[fail[u]][ch];
}
}
int rt;
ll tmp

,g

,f

;
inline void mul(ll a

,ll b

)
{
memset(tmp,0,sizeof(tmp));
for(int i=1;i<=tot;++i)
for(int j=1;j<=tot;++j)
for(int k=1;k<=tot;++k)
tmp[i][j]+=a[i][k]*b[k][j];
for(int i=1;i<=tot;++i)
for(int j=1;j<=tot;++j)
a[i][j]=tmp[i][j];
}
inline void ksm(ll t)
{
while(t)
{
if(t&1) mul(f,g);
mul(g,g);
t>>=1;
}
}
int main()
{
//    freopen("2243.in","r",stdin);
while(scanf("%d %d",&n,&m)!=EOF)
{
tot=rt=0;
for(int i=1;i<=n;++i)
{
scanf("%s",sr+1);
l=strlen(sr+1);
insert(rt,1);
}
build();
for(int i=1;i<=tot+1;++i)
for(int j=1;j<=tot+1;++j)
f[i][j]=g[i][j]=0;
f[1][1]=1;
for(int j=1;j<=tot;++j)
if(!mark[j])
for(int ch=0;ch<=25;++ch)
g[j][top[j][ch]]++;
else g[j][j]+=26,g[j][tot+1]++;
++tot;
g[tot][tot]=1;
ksm((ll)m+1);
printf("%I64u\n",f[1][tot]);
}
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: 
相关文章推荐