这道到是不用看题解,不过太经典了,早就被剧透一脸了

这道题很像ac自动机上的dp(其实就是)

然后注意到n很大,节点很小,于是就可以用矩阵快速幂优化了

时间复杂度为o(m^3 *log n);

蒟蒻kpm写得少,改了好久= =

CODE:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
int n,m,mod;
#define maxn 30
char c[maxn];
struct martix{
 int r,c;int f[maxn][maxn];
 int init(int x){
  r=c=x;
  for (int i=1;i<=x;i++) f[i][i]=1;
 }
}a,b;
martix operator * (const martix &x,const martix &y){
 martix ans;
 ans.r=x.r;
 ans.c=y.c;
 for (int i=1;i<=ans.r;i++)
  for (int j=1;j<=ans.c;j++){
   ans.f[i][j]=0;
   for (int k=1;k<=x.c;k++)
   ans.f[i][j]=((x.f[i][k]*y.f[k][j])%mod+ans.f[i][j])%mod;
  }
 return ans;
}
martix quick(int x,martix y){
 martix ans;
 ans.init(y.r);
 for (;x;x>>=1){
  if (x&1) ans=ans*y;
  y=y*y;
 }
 return ans;
}
int next[maxn];
int main(){
 scanf("%d%d%d",&n,&m,&mod);
 scanf("%s",c+1);
 int t=0;
 next[1]=0;
 for (int i=2;i<=m;i++) {
  while (t&&c[t+1]!=c[i]) t=next[t];
  next[i]=c[t+1]!=c[i]?0:++t;
 }
 a.r=a.c=m+1;
 for (int i=1;i<=m;i++) a.f[i][i+1]=1;
 a.f[1][1]=9;
 for (int i=1;i<m;i++) {
  int sum=1;
  {
   for (int j=2;j<=m+1;j++)
    if (a.f[next[i]+1][j]&&c[j-1]!=c[i+1]){
     a.f[i+1][j]=a.f[next[i]+1][j];sum+=a.f[i+1][j];
    }
  }
  a.f[i+1][1]=10-sum;
 }
 a=quick(n,a);
 b.r=1;b.c=1+m;
 b.f[1][1]=1;
 a=b*a;
 int ans=0;
 for (int i=1;i<=m;i++) (ans+=a.f[1][i])%=mod;
 printf("%d",ans);
 return 0;
}

BZOJ 1009 :[HNOI2008]GT考试(KPM算法+dp+矩阵快速幂)的更多相关文章

  1. bzoj 1009: [HNOI2008]GT考试【kmp+dp+矩阵快速幂】

    看n和k的范围长得就很像矩阵乘法了 设f[i][j]表示到第i个位置的后缀最长匹配目标串的j位.转移的话显然是枚举0~9,然后选择f[i+1]中能被他转移的加起来,需要用到next数组.然后构造矩阵的 ...

  2. BZOJ 1009: [HNOI2008]GT考试(kmp+dp+矩阵优化)

    http://www.lydsy.com/JudgeOnline/problem.php?id=1009 题意: 思路:真的是好题啊! 对于这种题目,很有可能就是dp,$f[i][j]$表示分析到第 ...

  3. 【BZOJ】4861: [Beijing2017]魔法咒语 AC自动机+DP+矩阵快速幂

    [题意]给定n个原串和m个禁忌串,要求用原串集合能拼出的不含禁忌串且长度为L的串的数量.(60%)n,m<=50,L<=100.(40%)原串长度为1或2,L<=10^18. [算法 ...

  4. [BZOJ 1009] [HNOI2008] GT考试 【AC自动机 + 矩阵乘法优化DP】

    题目链接:BZOJ - 1009 题目分析 题目要求求出不包含给定字符串的长度为 n 的字符串的数量. 既然这样,应该就是 KMP + DP ,用 f[i][j] 表示长度为 i ,匹配到模式串第 j ...

  5. BZOJ 1009 HNOI2008 GT考试 KMP算法+矩阵乘法

    标题效果:给定的长度m数字字符串s.求不包括子s长度n数字串的数目 n<=10^9 看这个O(n)它与 我们不认为这 令f[i][j]长度i号码的最后的字符串j位和s前者j数字匹配方案 例如,当 ...

  6. [BZOJ1009] [HNOI2008] GT考试(KMP+dp+矩阵快速幂)

    [BZOJ1009] [HNOI2008] GT考试(KMP+dp+矩阵快速幂) 题面 阿申准备报名参加GT考试,准考证号为N位数X1X2-.Xn,他不希望准考证号上出现不吉利的数字.他的不吉利数学A ...

  7. 【BZOJ1009】GT考试(KMP算法,矩阵快速幂,动态规划)

    [BZOJ1009]GT考试(KMP算法,矩阵快速幂,动态规划) 题面 BZOJ 题解 看到这个题目 化简一下题意 长度为\(n\)的,由\(0-9\)组成的字符串中 不含串\(s\)的串的数量有几个 ...

  8. 【BZOJ】2004: [Hnoi2010]Bus 公交线路 状压DP+矩阵快速幂

    [题意]n个点等距排列在长度为n-1的直线上,初始点1~k都有一辆公车,每辆公车都需要一些停靠点,每个点至多只能被一辆公车停靠,且每辆公车相邻两个停靠点的距离至多为p,所有公车最后会停在n-k+1~n ...

  9. BZOJ5298 CQOI2018 交错序列 【DP+矩阵快速幂优化】*

    BZOJ5298 CQOI2018 交错序列 [DP+矩阵快速幂优化] Description 我们称一个仅由0.1构成的序列为"交错序列",当且仅当序列中没有相邻的1(可以有相邻 ...

随机推荐

  1. 安卓 handler解析

    参考: 1.http://www.2cto.com/kf/201302/190591.html(安卓更新UI的两种方式,handler,Activity.runOnUIThread()) 2.http ...

  2. ceentos5.5 配置samba服务&用户&组

    准备 Change Root Password passwd root 在提示下建立新密码 静态IP vi /etc/sysconfig/network-scripts/ifcfg-eth0  #网络 ...

  3. leetcode[149]Max Points on a Line

    Given n points on a 2D plane, find the maximum number of points that lie on the same straight line. ...

  4. TortoiseSVN使用简介

    TortoiseSVN使用简介 2009-04-24 来源: dev.idv.tw 1.安装及下载client 端 2.什么是SVN(Subversion)? 3.为甚么要用SVN? 4.怎么样在Wi ...

  5. C#综合揭秘——细说进程、应用程序域与上下文之间的关系

    引言 本文主要是介绍进程(Process).应用程序域(AppDomain)..NET上下文(Context)的概念与操作.虽然在一般的开发当中这三者并不常用,但熟悉三者的关系,深入了解其作用,对提高 ...

  6. Cocos2d-x 多分辨率支持

    最近遇到多分辨率支持问题,所以查了一些资料.将一些收获共享一下,以便自己和其他需要的朋友日后参考. 如果我要建立一个cocos2d-x项目,我的目标是支持iphone3G( 480, 320 ),ip ...

  7. Python内置数据类型总结

    python的核心数据类型:(很多语言之提供了数字,字符串,文件数据类型,其他形式的数据类型都以标准库的形式表示 也就是用之前需要import ) ,但是python有很多都是内置的,不需要impor ...

  8. java基础概略总结

    /*************************************************/ String b=""; String c=""; // ...

  9. 个人对现在大众对perl的偏见的一些见解

    最近我都在一些论坛和交流社区学习,发现很多人对perl有很大的偏见. 以我学习那么久的Perl来说吧,也算是有一些小经验了,所以我总结了大家对perl的偏见的原因,无非就是是下面两个的两个原因: 1. ...

  10. 条形码--JsBarcode

    介绍一下在GitHub生成条形码的js插件→JsBarcode 支持的条形码: 条码支持的有: CODE128  CODE128 (自动模式切换)  CODE128 A/B/C (强制模式)EAN  ...