目录

题目链接

loj#2128. 「HAOI2015」数字串拆分

题解

\(f(s)\)对于\(f(i) = \sum_{j = i - m}^{i - 1}f(j)\)

这个可以用转移矩阵通过矩阵乘法处理出来

预处理出\(A[i][j]\)表示数S为\(j * 10 ^ i\)的转移矩阵

对于g的转移

\(g(i) = \sum_{j = 0}^{i - 1}g(j) * D(j + 1,i)\)

D[i][j]表示第i位到底j位构成的数的f,(转移矩阵

对于g的转移也是需要矩阵的

g(0) = 1也就是{0,0,0,.....,0,1}

代码


#include<map>
#include<vector>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define gc getchar()
#define pc putchar
#define LL long long
inline int read() {
int x = 0,f = 1;
char c = gc;
while(c < '0' || c > '9')c = gc;
while(c <= '9' && c >= '0') x = x * 10 + c - '0',c = gc;
return x * f;
}
void print(int x) {
if(x < 0) {
pc('-');
x = -x;
}
if(x >= 10) print(x / 10);
pc(x % 10 + '0');
}
const int maxn = 507;
int n,m;
const int mod = 998244353;
struct Ma {
int a[6][6];
Ma() { memset(a,0,sizeof a); }
Ma operator * (const Ma & t)const {
Ma ret;
memset(ret.a,0,sizeof ret.a);
for(int i = 1;i <= m;++ i)
for(int j = 1;j <= m;++ j)
for(int k = 1;k <= m;++ k)
(ret.a[i][j] += 1ll * a[i][k] * t.a[k][j] % mod) %= mod;
return ret;
}
Ma operator + (const Ma &k)const {
Ma ret;
memset(ret.a,0,sizeof ret.a);
for(int i = 1;i <= m;++ i)
for(int j = 1;j <= m;++ j)
ret.a[i][j] = (a[i][j] + k.a[i][j]) % mod;
return ret;
}
} A[maxn][10],f[maxn];
char s[maxn];
int main() {
scanf("%s",s + 1);
n = strlen(s + 1);
m = read();
for(int i = 1;i <= m;++ i)
A[0][0].a[i][i] = 1 ,
A[0][1].a[i][m] = 1;
for(int i = 1;i < m;++ i)
A[0][1].a[i + 1][i] = 1;
for(int i = 2;i <= 9;++ i)
A[0][i] = A[0][i - 1] * A[0][1];
for(int i = 1;i <= n;++ i) { //十进制快速幂
A[i][0] = A[0][0];
A[i][1] = A[i - 1][9] * A[i - 1][1];
for(int j = 2;j <= 9;++ j)
A[i][j] = A[i][j - 1] * A[i][1];
}
f[0].a[1][m] = 1;
Ma now;
for(int i = 1;i <= n;++ i) {
now = A[0] [s[i] - '0'];
for(int j = i - 1;j >= 0;-- j) {
f[i] = f[i] + (f[j] * now);
if(j) now = A[i - j][s[j] - '0'] * now;
}
}
print(f[n].a[1][m]);
pc('\n');
return 0;
}

loj#2128. 「HAOI2015」数字串拆分 矩阵乘法的更多相关文章

  1. 【LOJ】#2128. 「HAOI2015」数字串拆分

    题解 题中给的函数可以用矩阵快速幂递推 我们记一个数组dp[i](这个数组每个元素是一个矩阵)表示从1到i所有的数字经过拆分矩阵递推的加和 转移方法是 \(dp[i] = \sum_{j = 0}^{ ...

  2. loj#2002. 「SDOI2017」序列计数(dp 矩阵乘法)

    题意 题目链接 Sol 质数的限制并没有什么卵用,直接容斥一下:答案 = 忽略质数总的方案 - 没有质数的方案 那么直接dp,设\(f[i][j]\)表示到第i个位置,当前和为j的方案数 \(f[i ...

  3. loj #2006. 「SCOI2015」小凸玩矩阵

    #2006. 「SCOI2015」小凸玩矩阵   题目描述 小凸和小方是好朋友,小方给小凸一个 N×M N \times MN×M(N≤M N \leq MN≤M)的矩阵 A AA,要求小凸从其中选出 ...

  4. LOJ#2127「HAOI2015」按位或

    用$ Min-Max$容斥之后要推的东西少了好多 无耻的用实数快读抢了BZOJ.Luogu.LOJ三个$ OJ$的Rank 1 即将update:被STO TXC OTZ超了QAQ 题意:集合$ [0 ...

  5. LOJ#2083. 「NOI2016」优秀的拆分

    $n \leq 30000$的字符串,问其所有子串的所有AABB形式的拆分有多少种.$t \leq 10$组询问. $n^3$过80,$n^2$过95,鬼去写正解.. $n^2$:先枚举一次算每个位置 ...

  6. 「BZOJ2510」弱题(矩阵乘法,降维)

    有M个球,一开始每个球均有一个初始标号,标号范围为1-N且为整数,标号为i的球有ai个,并保证Σai = M. 每次操作等概率取出一个球(即取出每个球的概率均为1/M),若这个球标号为k(k < ...

  7. loj#2013. 「SCOI2016」幸运数字 点分治/线性基

    题目链接 loj#2013. 「SCOI2016」幸运数字 题解 和树上路径有管...点分治吧 把询问挂到点上 求出重心后,求出重心到每个点路径上的数的线性基 对于重心为lca的合并寻味,否则标记下传 ...

  8. Loj #3089. 「BJOI2019」奥术神杖

    Loj #3089. 「BJOI2019」奥术神杖 题目描述 Bezorath 大陆抵抗地灾军团入侵的战争进入了僵持的阶段,世世代代生活在 Bezorath 这片大陆的精灵们开始寻找远古时代诸神遗留的 ...

  9. Loj #3055. 「HNOI2019」JOJO

    Loj #3055. 「HNOI2019」JOJO JOJO 的奇幻冒险是一部非常火的漫画.漫画中的男主角经常喜欢连续喊很多的「欧拉」或者「木大」. 为了防止字太多挡住漫画内容,现在打算在新的漫画中用 ...

随机推荐

  1. Linux内核驱动之GPIO子系统(一)GPIO的使用【转】

    转自:http://blog.csdn.net/tommy_wxie/article/details/9427047 一 概述 Linux内核中gpio是最简单,最常用的资源(和 interrupt  ...

  2. oracle删除表字段和oracle表增加字段

    这篇文章主要介绍了oracle表增加字段.删除表字段修改表字段的使用方法,大家参考使用吧   添加字段的语法:alter table tablename add (column datatype [d ...

  3. MPI 在Windows10 上安装,使用VS2013编译生成可执行程序

    原文地址:http://www.cnblogs.com/leijin0211/p/6851789.html 参考博客: http://www.cnblogs.com/shixiangwan/p/662 ...

  4. transfer pdf to png

    #! /bin/bash # # transfer pdf to png if [ $# != 1 ] ; then echo "USAGE: $0 PDF FILE ABSOLUTELY ...

  5. vue系列之项目结构

    参考地址:链接 build webpack配置相关 config webpack配置相关 node_modules  npm install 安装的依赖代码库 src  存放项目源码 static 存 ...

  6. ssh-keygen -t rsa -b 4096 -C "邮箱"

    ssh-keygen -t rsa -b 4096 -C "邮箱":这条命令的目的是为了让本地机器ssh登录远程机器上的GitHub账户无需输入密码.将这条命令分解: 1.ssh- ...

  7. 5个php实例,细致说明传值与传引用的区别

    传值:是把实参的值赋值给行参 ,那么对行参的修改,不会影响实参的值 传引用 :真正的以地址的方式传递参数传递以后,行参和实参都是同一个对象,只是他们名字不同而已对行参的修改将影响实参的值 说明: 传值 ...

  8. 【mysql】autocommit=0后,commit, rollback无效

    之前在[mysql]MySQLdb中的事务处理中用autocommit和commit()以及rollback()实现了事务处理. 但后来,用同样的代码在另一个数据库中运行却失败了.找了一个下午的原因. ...

  9. hdu3308

    区间合并比较模板的题,就是求一个区间的LCIS 线段树维护左最大LCIS,右最大LCIS,区间LCIS 看代码就行 #include<iostream> #include<cstri ...

  10. 步步为营-36-ADO.Net简介

    与数据库进行连接交互 方法一 #region 01连接对象 //01 连接字符串 string connstr = "server=.;uid=sa;pwd=sa;database=Demo ...