description


analysis

  • 考虑矩阵乘法

  • 设初始\(m×m\)矩阵上\(i\)行\(j\)列的数字表示该矩阵第\(j\)位上\(f[i]\)的指数

  • 那么一开始表示\(f[1..k]\)的矩阵就长这个样子,举样例\(k=4\)的例子

\[\left(
\begin{matrix}
1,0,0,0\\
0,1,0,0\\
0,0,1,0\\
0,0,0,1\\
\end{matrix}
\right)\]

  • 也就是\(f[1]=f[1]^1,f[2]=f[2]^1\)等等

  • 可知\(f[5]=f[4]^{b[1]}f[3]^{b[2]}f[2]^{b[3]}f[1]^{b[4]}\),那表示\(f[2..k+1]\)的矩阵就是

\[\left(
\begin{matrix}
0,0,0,b[4]\\
1,0,0,b[3]\\
0,1,0,b[2]\\
0,0,1,b[1]\\
\end{matrix}
\right)
\]

  • 不懂可以把这个矩阵的各项拆出来,同一列从上往下做\(f\)的次幂再相乘就可以分别得到\(f[2..k+1]\)

  • 由于第一个矩阵相当于矩阵意义的\(1\),所以转移矩阵就是第二个矩阵

  • 好像这就没了,但是要知道指数的矩乘不能直接取模,比如\(3^{15}\mod 7≠3^{15\mod 7}\)

  • 费马小定理告诉你\(a^{p-1}≡1(\mod p)\),也就是说每\(p-1\)个\(a\)相乘的积模\(p\)等于\(1\)

  • 于是矩乘里的模数取原来的模数\(-1\)就可以了

  • 我在考场上最后十分钟推出第二个矩阵对然后我™就不知道那个就是转移矩阵然后傻逼地对着转移矩阵发呆


code

#pragma GCC optimize("O3")
#pragma G++ optimize("O3")
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define MAXK 205
#define mod 998244353
#define MOD 998244352
#define ll long long
#define reg register ll
#define fo(i,a,b) for (reg i=a;i<=b;++i)
#define fd(i,a,b) for (reg i=a;i>=b;--i) using namespace std; ll b[MAXK],f[MAXK];
ll n,m,ans; struct matrix
{
ll f[MAXK][MAXK],n,m;
matrix(){memset(f,0,sizeof(f));}
}tmp;
inline ll read()
{
ll x=0,f=1;char ch=getchar();
while (ch<'0' || '9'<ch){if (ch=='-')f=-1;ch=getchar();}
while ('0'<=ch && ch<='9')x=x*10+ch-'0',ch=getchar();
return x*f;
}
inline matrix operator*(matrix a,matrix b)
{
matrix c;
fo(i,1,m)fo(j,1,m)fo(k,1,m)
c.f[i][j]=(c.f[i][j]+a.f[i][k]*b.f[k][j])%MOD;
return c;
}
inline matrix pow(matrix x,ll y)
{
matrix z;
fo(i,1,m)z.f[i][i]=1;
if (y==0)return z;
while (y)
{
if (y&1)z=z*x;
x=x*x,y>>=1;
}
return z;
}
inline ll ksm(ll x,ll y)
{
ll z=1;
while (y)
{
if (y&1)z=z*x%mod;
x=x*x%mod,y>>=1;
}
return z;
}
int main()
{
freopen("T1.in","r",stdin);
//freopen("seq.in","r",stdin);
//freopen("seq.out","w",stdout);
n=read(),m=read();
fo(i,1,m)b[i]=read();
fo(i,1,m)f[i]=read();
if (n<=m){printf("%lld\n",f[n]);return 0;}
fo(i,2,m)tmp.f[i][i-1]=1;
fo(i,1,m)tmp.f[i][m]=b[m-i+1];
tmp=pow(tmp,n-m),ans=1ll;
fo(i,1,m)ans=(ans*(ksm(f[i],tmp.f[i][m])))%mod;
printf("%lld\n",ans);
return 0;
}

【JZOJ6275】小L的数列的更多相关文章

  1. NC15979 小q的数列

    NC15979 小q的数列 题目 题目描述 小q最近迷上了各种好玩的数列,这天,他发现了一个有趣的数列,其递推公式如下: \[f[i] = \left\{ \begin{array}{l} 0 ,&a ...

  2. 洛谷U4727小L的二叉树[树转序列 LIS]

    题目背景 勤奋又善于思考的小L接触了信息学竞赛,开始的学习十分顺利.但是,小L对数据结构的掌握实在十分渣渣. 所以,小L当时卡在了二叉树. 题目描述 在计算机科学中,二叉树是每个结点最多有两个子结点的 ...

  3. ACM D的小L

    D的小L 时间限制:4000 ms  |  内存限制:65535 KB 难度:2   描述       一天TC的匡匡找ACM的小L玩三国杀,但是这会小L忙着哩,不想和匡匡玩但又怕匡匡生气,这时小L给 ...

  4. nyoj-366-D的小L(求全排列)

    D的小L 时间限制:4000 ms  |  内存限制:65535 KB 难度:2 描述       一天TC的匡匡找ACM的小L玩三国杀,但是这会小L忙着哩,不想和匡匡玩但又怕匡匡生气,这时小L给匡匡 ...

  5. 洛谷U4727 小L 的二叉树

    U4727 小L 的二叉树 题目背景 勤奋又善于思考的小L接触了信息学竞赛,开始的学习十分顺利.但是,小L对数据结构的掌握实在十分渣渣. 所以,小L当时卡在了二叉树. 题目描述 在计算机科学中,二叉树 ...

  6. 【BZOJ4030】[HEOI2015]小L的白日梦

    [BZOJ4030][HEOI2015]小L的白日梦 题面 BZOJ 洛谷 题解 要求的是最小的不开心连续段的期望. 然后发现自己就不会做了. 然后就可以来抄题解啦. 首先来猜性质: 第一个,一定是按 ...

  7. BZOJ 4030: [HEOI2015]小L的白日梦

    4030: [HEOI2015]小L的白日梦 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 172  Solved: 39[Submit][Statu ...

  8. 小L记单词

    题目描述 小L最近在努力学习英语,但是对一些词组总是记不住,小L小把这些词组中每一个单词的首字母都记一下,这样形成词组的缩写,通过这种方式小L的学习效率明显提高. 输入 输入有多行,每组测试数据占一行 ...

  9. 小L的试卷

    题目描述 小L期末考试结束,高高兴兴放假回家了,可是那么多试卷,老师还要加班批改,有n份试卷由k个老师批改,n份试卷进行了密封编号,由于试卷上的做题情况和书写的规范程序不一样,批改不同的试卷用时也可能 ...

随机推荐

  1. 10-MySQL-Ubuntu-数据表中数据的查询(三)

    数据的查询(select) (1)查询整个表的数据: select  * from 表名; (2)查询给定条件的数据: select  * from 表名 where 条件; (3)查询表中某些字段: ...

  2. 1、如何在列表,字典,集合种根据条件筛选数据?2、如何为元组中的每个元素命名,提高程序的可读性3、如何统计出序列中元素出现的频度4、如何根据字典中value的大小,对字典的key进行排序

    一.数据筛选: 处理方式: 1.filter函数在py3,返回的是个生成式. from random import randint data = [randint(-100,100) for i in ...

  3. C语言结构体嵌套

    #include <stdio.h> int main() { /*************************************************** *结构体嵌套:结构 ...

  4. Vue番外篇-路由进阶(一)

    Vue的router默认是 export default new Router({ mode: 'history', routes: [ { path: '/', name: 'HelloWorld' ...

  5. 驾驶心率和呼吸,疲劳检测系统,通过安全带,坐垫等内置sensor

    HARKEN - Heart and respiration in car embedded nonintrusive sensors 2 https://www.youtube.com/watch? ...

  6. 两台linux服务器传文件方法

    scp -r /当前服务器的文件的路径 第二台服务器的用户名@第二台服务器的ip:/第二台服务器的存放文件的路径

  7. thinkphp 动态配置

    之前的方式都是通过预先定义配置文件的方式,而在具体的操作方法里面,我们仍然可以对某些参数进行动态配置(或者增加新的配置),主要是指那些还没有被使用的参数. 设置新的值: C('参数名称','新的参数值 ...

  8. 17 win7 sp1 x64/VS2015下配置creo4.0二次开发环境——调用了众多开源库(ceres-solver,PCL1.8.0,office 2016COM接口,MySql数据库等)

    0 引言 本次开发环境的配置是在综合考虑了开源库的版本.VS版本以及CREO4.0的版本,同时针对甲方需求选择了win7 sp1 x64系统. 配置的过程中遇到了形形色色的问题,但是一一解决了.通过这 ...

  9. pytorch实现kaggle猫狗识别

    参考:https://blog.csdn.net/weixin_37813036/article/details/90718310 kaggle是一个为开发商和数据科学家提供举办机器学习竞赛.托管数据 ...

  10. NX二次开发-UFUN移动工程图视图到另一个图纸页UF_DRAW_move_view_to_drawing

    #include <uf.h> #include <uf_draw.h> #include <uf_obj.h> UF_initialize(); //找名字获取视 ...