CF1106F Lunar New Year and a Recursive Sequence
题目链接:CF1106F Lunar New Year and a Recursive Sequence
大意:已知\(f_1,f_2,\cdots,f_{k-1}\)和\(b_1,b_2,\cdots,b_k\),且有递推关系
\]
对于所有\(i>k\)均成立,给出\(f_n=m\),求\(f_k\),\(p=998244353\)
分析:
数论板子大集合
首先很显然有\(f_n=f_k^{q}\),于是考虑先求出这个\(q\)
看到\(k\leq100\)一般就会想到矩阵乘法之类的
但是这个式子的\(b\)处在乘方的位置,无法直接使用矩阵
注意到模数为998244353,其原根为3
那么我们就可以使用原根来改写这个转移
\]
这个是可以利用矩阵乘法的,初值\(g_i=0(i<k),g_k=1\)
我们于是可以利用矩乘找到这样的关系式:\(f_k^{q}\equiv f_n(mod\ p),q=g_n\)
这个玩意似乎可以类比的被定义为k次剩余,不过我们显然不会这个玩意
于是我们继续使用原根
我们知道\(f_n\)可以被写作\(3^t\),那么我们设\(f_k\)可以被写作为\(3^x\)
那么我们就是求满足这个式子的x值:\(3^{xq}\equiv3^t(mod\ p)\)
那么对于指数我们可以直接用exgcd,因为此时有\(xq\equiv t(mod\ p-1)\)
将m转为\(3^t\)时可以使用BSGS
问题就得到了解决QAQ
注意模数会在\(p-1\)和\(p\)之间切换,指数的模数为\(p-1\)(包括上面的矩乘),而在BSGS中的模数则为\(p\)
#include<iostream>
#include<string.h>
#include<string>
#include<stdio.h>
#include<algorithm>
#include<math.h>
#include<vector>
#include<queue>
#include<map>
using namespace std;
#define int long long
const int maxd=998244353,N=100000;
const double pi=acos(-1.0);
typedef long long ll;
int n,K,b[120],m;
struct matrix{
ll x[120][120];
}ans,sum;
matrix operator *(matrix a,matrix b)
{
matrix c;
memset(c.x,0,sizeof(c.x));
int i,j,k;
for (i=0;i<K;i++)
{
for (j=0;j<K;j++)
{
for (k=0;k<K;k++)
{
c.x[i][j]=(c.x[i][j]+a.x[i][k]*b.x[k][j])%(maxd-1);
}
}
}
return c;
}
void qpow(int tim)
{
while (tim)
{
int tmp=tim&1;tim/=2;
if (tmp) ans=ans*sum;
sum=sum*sum;
}
}
ll qpow(ll x,ll y)
{
ll ans=1,sum=x;
while (y)
{
int tmp=y&1;y/=2;
if (tmp) ans=(ans*sum)%maxd;
sum=(sum*sum)%maxd;
}
return ans;
}
map<ll,ll> mp;
ll bsgs(int a,int b)
{
int siz=sqrt(maxd-1)+1,i;ll sum=b;
for (i=0;i<=siz;i++)
{
mp[sum]=i;
sum=(sum*a)%maxd;
}
ll tmp=qpow(a,siz);sum=1;
for (i=1;i<=siz;i++)
{
sum=(sum*tmp)%maxd;
if (mp[sum]) return i*siz-mp[sum];
}
}
ll gcd(ll x,ll y)
{
if (!y) return x; else return gcd(y,x%y);
}
void exgcd(ll a,ll b,ll &g,ll &x,ll &y)
{
if (!b) {x=1;y=0;g=a;return;}
else
{
ll tmpx,tmpy;
exgcd(b,a%b,g,tmpx,tmpy);
x=tmpy;y=tmpx-a/b*tmpy;
}
}
int read()
{
int x=0,f=1;char ch=getchar();
while ((ch<'0') || (ch>'9')) {if (ch=='-') f=-1;ch=getchar();}
while ((ch>='0') && (ch<='9')) {x=x*10+(ch-'0');ch=getchar();}
return x*f;
}
signed main()
{
K=read();
int i,j;
for (i=0;i<K;i++) {ans.x[i][i]=1;sum.x[0][i]=read();}
for (i=0;i<K-1;i++) sum.x[i+1][i]=1;
ans.x[0][0]=1;
n=read();m=read();
qpow(n-K);
ll q=ans.x[0][0],p=bsgs(3,m);
//cout << q << " " << p << endl;
ll tmp=gcd(q,maxd-1);
if (p%tmp) {printf("-1");return 0;}
ll x,y,g;
exgcd(q,maxd-1,g,x,y);
x=(p/g*x)%(maxd-1);
if (x<0) x+=(maxd-1);
printf("%lld",qpow(3,x));
return 0;
}
/*
3
2 3 5
4 16
*/
CF1106F Lunar New Year and a Recursive Sequence的更多相关文章
- CF1106F Lunar New Year and a Recursive Sequence 原根、矩阵快速幂、BSGS
传送门 好久没写数论题了写一次调了1h 首先发现递推式是一个乘方的形式,线性递推和矩阵快速幂似乎都做不了,那么是否能够把乘方运算变成加法运算和乘法运算呢? 使用原根!学过\(NTT\)的都知道\(99 ...
- CF1106F Lunar New Year and a Recursive Sequence 线性递推 + k次剩余
已知\(f_i = \prod \limits_{j = 1}^k f_{i - j}^{b_j}\;mod\;998244353\),并且\(f_1, f_2, ..., f_{k - 1} = 1 ...
- CF1106F Lunar New Year and a Recursive Sequence(矩阵快速幂+bsgs+exgcd)
题面 传送门 前置芝士 \(BSGS\) 什么?你不会\(BSGS\)?百度啊 原根 对于素数\(p\)和自然数\(a\),如果满足\(a^x\equiv 1\pmod{p}\)的最小的\(x\)为\ ...
- CF1106F Lunar New Year and a Recursive Sequence——矩阵快速幂&&bsgs
题意 设 $$f_i = \left\{\begin{matrix}1 , \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ i < k\\ ...
- Codeforces 1106F Lunar New Year and a Recursive Sequence | BSGS/exgcd/矩阵乘法
我诈尸啦! 高三退役选手好不容易抛弃天利和金考卷打场CF,结果打得和shi一样--还因为queue太长而unrated了!一个学期不敲代码实在是忘干净了-- 没分该没分,考题还是要订正的 =v= 欢迎 ...
- Codeforces 1106F Lunar New Year and a Recursive Sequence (数学、线性代数、线性递推、数论、BSGS、扩展欧几里得算法)
哎呀大水题..我写了一个多小时..好没救啊.. 数论板子X合一? 注意: 本文中变量名称区分大小写. 题意: 给一个\(n\)阶递推序列\(f_k=\prod^{n}_{i=1} f_{k-i}b_i ...
- @codeforces - 1106F@ Lunar New Year and a Recursive Sequence
目录 @description@ @solution@ @accepted code@ @details@ @description@ 定义递推数列 f: (1)f[1] = f[2] = ... f ...
- hdu 5950 Recursive sequence 矩阵快速幂
Recursive sequence Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Other ...
- HDU 5950 Recursive sequence 【递推+矩阵快速幂】 (2016ACM/ICPC亚洲区沈阳站)
Recursive sequence Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Other ...
随机推荐
- docker容器与宿主交互数据
1.查看容器 [root@localhost ~]# docker ps CONTAINER ID IMAGE COMMAND CREATED STATUS PORTS NAMES cd6957191 ...
- 小P的金字塔
题目描述 小P感到自己前几天太作了,于是非常有远见的决定为自己建立一座金字塔. 现在他有n种标准长方体的石头,每种石头只有两个,第i种石头是长宽高分别为Xi,Yi,Zi的长方体.由于整个工程只有小P一 ...
- 微信开发 提示 Redirect_uri(错误10003)
情景: 搭建完成一个网站,使用微信打开链接地址,结果报错1003 完整的错误信息: 出现这种情况一般有两种原因: 1.没有配置网页授权 我们可以根据微信的开发者文档http://mp.weixin. ...
- Azure系列2.1.14 —— CopyState
(小弟自学Azure,文中有不正确之处,请路过各位大神指正.) 网上azure的资料较少,尤其是API,全是英文的,中文资料更是少之又少.这次由于公司项目需要使用Azure,所以对Azure的一些学习 ...
- [转帖]浅谈程序中的text段、data段和bss段
作者:百问科技链接:https://zhuanlan.zhihu.com/p/28659560来源:知乎著作权归作者所有.商业转载请联系作者获得授权,非商业转载请注明出处. 一般情况,一个程序本质上都 ...
- 【学亮IT手记】PL/SQL游标编程
游标提供了一种从表中检索数据并进行操作的灵活手段,主要用在服务器上,处理由客户端发送给服务器端的sql语句,或者是批处理.存储过程.触发器中的数据处理请求. 显式游标 是由用户声明和操作的一种游标,通 ...
- 关于wordpress更新提示无法创建目录问题
说说自己的看法和解决办法 看法: 网上很多人说:是权限问题,那么将文件目录权限设置为777就可以解决.恩,没错,是可以解决更新问题,可是却带来了更大的问题——安全.像他们这个设置后,网站被攻破,数据被 ...
- DNS_PROBE_FINISHED_NXDOMAIN & MacOS
DNS_PROBE_FINISHED_NXDOMAIN 内网 DNS bug 8.8.8.8 8.8.4.4 # new inner Wi-Fi 10.1.3.10 10.1.3.13 Windows ...
- 常见IT工具软件总结
1. 阿里云在线迁移服务 2.智能媒体管理 格式转换 业务域名管理 1. 每个业务有一个英文单词, 1. 每个 git 的命名应该是 chgg-业务英文-种类 2. 例如 chgg-plant-api ...
- tornado web框架简介
https://www.cnblogs.com/aylin/p/5702994.html