[BJWC2011]禁忌
IX.[BJWC2011]禁忌
老套路,建出AC自动机,并求出每个节点的ok。
之后呢?
考虑dp。
令fi,j表示在一个长度为i的串,当前在点j的可能性。
则有:
对于任意k∈[0,alphabet),
如果t[t[j].ch[k]].ok为true,那么fi+1,k+=fi,j/alphabet;
否则,在统计答案的变量中加上fi,j/alphabet,同时fi+1,1+=fi,j/alphabet(匹配成功,默认跳到根节点)。
但是!!!这里的len巨大,我们不能这么跑,否则复杂度是O(size∗len∗alphabet)的。
考虑一下,这每次的转移的步骤,都与i和具体的fi,j的值无关。因此可以考虑矩阵快速幂。
构建转移矩阵base。
for(int i=1;i<=cnt;i++){
for(int j=0;j<m;j++){
if(t[t[i].ch[j]].ok)base.num[i][t[i].ch[j]]+=1;
else base.num[i][1]+=1,base.num[i][cnt+1]+=1;
}
for(int i=1;i<=base.N;i++)for(int j=1;j<=base.N;j++)base.num[i][j]/=m;
base.num[base.N][base.N]=1.0;
}
(m即为alphabet。)
该矩阵的最后一列为统计答案列。但是,该矩阵的右下角处必需强制赋成1。为什么呢?
看到矩阵乘法x∗y的过程中什么时候要用到右下角:
or(int i=1;i<=x.N;i++)for(int j=1;j<=y.N;j++)for(int k=1;k<=y.N;k++)z.num[i][j]+=x.num[i][k]*y.num[k][j];
则在做最后一行和最后一列的乘法时,它调用了右下角。
最后一行全为0(因为原本应该是一个size∗size的方阵,但是我们强制补了最后一列),可以忽略。
但是,计算最后一列时,这一点是必须的。
然后我们就可以暴力矩阵快速幂了。
最后答案为矩阵的右上角。
因为最右一列是答案列,而最上一行为根节点行。
代码:
#include<bits/stdc++.h>
using namespace std;
int n,len,m,S,cnt=1;
struct AC_Automaton{
int ch[26],fail;
bool ok;
}t[80];
char s[20];
void ins(){
int x=1;
for(int i=0;i<S;i++){
if(!t[x].ch[s[i]-'a'])t[x].ch[s[i]-'a']=++cnt,t[cnt].ok=true;
x=t[x].ch[s[i]-'a'];
}
t[x].ok=false;
}
queue<int>q;
void build(){
for(int i=0;i<m;i++){
if(t[1].ch[i])t[t[1].ch[i]].fail=1,q.push(t[1].ch[i]);
else t[1].ch[i]=1;
}
while(!q.empty()){
int x=q.front();q.pop();
for(int i=0;i<m;i++){
if(t[x].ch[i])t[t[x].ch[i]].fail=t[t[x].fail].ch[i],q.push(t[x].ch[i]);
else t[x].ch[i]=t[t[x].fail].ch[i];
}
t[x].ok&=t[t[x].fail].ok;
}
}
struct Matrix{
long double num[80][80];
int N;
Matrix(){
N=0;
memset(num,0,sizeof(num));
}
friend Matrix operator *(const Matrix &x,const Matrix &y){
Matrix z;
z.N=x.N;
for(int i=1;i<=x.N;i++)for(int j=1;j<=y.N;j++)for(int k=1;k<=y.N;k++)z.num[i][j]+=x.num[i][k]*y.num[k][j];
return z;
}
}base,I;
Matrix pov(Matrix x,int y){
if(!y)return I;
Matrix tmp=pov(x,y>>1);
return tmp*tmp*(y&1?x:I);
}
int main(){
scanf("%d%d%d",&n,&len,&m),t[1].ok=true;
for(int i=1;i<=n;i++)scanf("%s",s),S=strlen(s),ins();
build();
for(int i=1;i<=cnt;i++){
for(int j=0;j<m;j++){
if(t[t[i].ch[j]].ok)base.num[i][t[i].ch[j]]+=1;
else base.num[i][1]+=1,base.num[i][cnt+1]+=1;
}
}
base.N=I.N=cnt+1;
for(int i=1;i<=base.N;i++)for(int j=1;j<=base.N;j++)base.num[i][j]/=m;
base.num[base.N][base.N]=1.0;
for(int i=1;i<=I.N;i++)I.num[i][i]=1;
base=pov(base,len);
printf("%.8Lf\n",base.num[1][base.N]);
return 0;
}
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· .NET 9 new features-C#13新的锁类型和语义
· Linux系统下SQL Server数据库镜像配置全流程详解
· 现代计算机视觉入门之:什么是视频
· 你所不知道的 C/C++ 宏知识
· 聊一聊 操作系统蓝屏 c0000102 的故障分析
· 不到万不得已,千万不要去外包
· C# WebAPI 插件热插拔(持续更新中)
· 会议真的有必要吗?我们产品开发9年了,但从来没开过会
· 如何打造一个高并发系统?
· 《SpringBoot》EasyExcel实现百万数据的导入导出