【51nod 1538】 一道难题 - n阶常系数线性递推

原创 2017年01月04日 08:46:14

  题意:给定a,求

(ni=1aibi)=m(ni=1bi)!ni=1(bi!)

  m1018,ai[1,23333],bi0

【奇怪的背景】
  不知道多少天前。。。
  Q神:hgr你有 nlognlogk 求 n 阶线性递推第 k 项的板子吗
  我:好像没有啊。。。
  Q神:这个辣鸡51nod出了个裸题,写不动了
  于是看了波题。。。然后一脸懵逼???这怎么就线性递推了???
  然后这题就鸽了

  最近似乎比较有空。。。冷静一下之后推了一波这个题发现。。。怎么出个裸题啊
  这个式子看起来比较吓人。。其实是在搞笑。。。
  设Fm=(ni=1aibi)=m(ni=1bi)!ni=1(bi!)
  我们稍微加上个东西然后一波xjb变换
  

Fm=(ni=1aibi)=m(ni=1bi)!ni=1(bi!)nj=1bj[bj>0]ni=1bi  =j=1n(ni=1aibi)=m[bj>0](ni=1bi1)!ni=1(bi!)bj  =j=1n(ni=1aibi)=m[bj>0](ni=1bi1)!(bj1)!ij(bi!)  =j=1n(ni=1aibi)=m[bj10](ni=1bi1)!(bj1)!ij(bi!)  =j=1n(ni=1aibi)=maj[bj0](ni=1bi)!ni=1(bi!)

  因为bi0,所以直接有
  
Fm=i=1nFmai

  于是直接套n阶常系数齐次线性递推就行了。
  具体如何实现的话。。。叉姐论文和picks老师的博客有讲。。。这东西要完整地写出来不太好写。。。式子太多了= =
  而且这玩意写起代码来也是蛋疼的一b。。。
  
  有个东西我还不太清楚。。。比如现在递推式是fm=ni=1cifmi。。。然后要先求初值f0...fn对吧。。。然而似乎除了直接用生成函数搞个F(z)=11C(z)以外我没有想到别的复杂度更好的办法。。。(虽然这里多项式求逆和分治FFT速度好像没啥区别了hhhhh)
  然而这里一点都不瓶颈的样子→_→
  
  代码太丑了不贴了其实是TLE到现在弃疗了
  UPD:终于过了。。。辣鸡数据错了卡我这么久。。。有两个点的a_1是会大于23333的,要把他忽略掉。。。GTMD。。。

版权声明:本文为博主原创文章,未经博主允许不得转载。

矩阵快速幂在常系数线性递推关系中的应用

先引入一下,知乎上有一个问题 关于斐波拉契数列的一个低级问题 。题主询问了关于求解斐波拉契数列第n项对10007取模的结果。而这个n,可以达到10610^6甚至10910^9 。 解法已经在排名第一...
  • ArrowLLL
  • ArrowLLL
  • 2016年09月09日 00:34
  • 690

常系数线性递推式的快速求单项值方法

常系数线性递推式的快速求单项值方法前导知识 常系数线性递推式 形如f(n)=∑ki=1aif(n−i)f(n)=\sum_{i=1}^{k} a_if(n-i),其中kk、aia_i为常数的递推式,...
  • Le_ballon_rouge
  • Le_ballon_rouge
  • 2016年05月14日 20:31
  • 1161

[51nod 1538] 一道难题 -- 特征多项式,多项式取模

齐次线性递推式 (FFT优化)
  • whzzt
  • whzzt
  • 2017年05月05日 19:59
  • 470

【51nod 1538】 一道难题 - n阶常系数线性递推

辣鸡数据毁我青春 【奇怪的背景】   不知道多少天前。。。   Q神:hgr你有 nlognlogk 求 n 阶线性递推第 k 项的板子吗   我:好像没有啊。。。   ~~Q神:这个辣鸡51nod出...
  • GEOTCBRL
  • GEOTCBRL
  • 2017年01月04日 08:46
  • 673

矩阵乘幂优化k阶常系数线性递推关系

  • 2012年02月23日 20:58
  • 202KB
  • 下载

常系数齐次线性递推优化矩阵快速幂

一般矩阵快速幂的形式 : f(n)=∑i=1k1aif(n−i)+∑i=1k2big(n−i)+bf(n)=\sum_{i=1}^{k_1} a_if(n-i)+\sum_{i=1}^{k_2}b_...
  • qq_35649707
  • qq_35649707
  • 2017年12月01日 16:41
  • 128

bzoj 4161: Shlw loves matrixI 常系数线性齐次递推

题意给定数列 {hn}前k项,其后每一项满足 hn = a1*h(n-1) + a2*h(n-2) + … + ak*h(n-k) 其中 a1,a2…ak 为给定数列。请计算 h(n),并将结果对...
  • qq_33229466
  • qq_33229466
  • 2017年12月29日 17:00
  • 168

【BZOJ4944】【NOI2017】泳池 概率DP 常系数线性递推 特征多项式 多项式取模

题目大意   有一个1001×n1001\times n的的网格,每个格子有qq的概率是安全的,1−q1-q的概率是危险的。   定义一个矩形是合法的当且仅当: 这个矩形中每个格子都...
  • ez_yww
  • ez_yww
  • 2017年11月30日 18:51
  • 362

一种写程序快速计算常系数线性齐次递推关系的指定项的方法

上次学习到了利用Fibonacci恒等式在log(n)时间内计算出F(n)的方法。这种方法也适用于任何常系数线性齐次递推数列的计算。在这里把能想到的东西小结一下。例一:扩展Euclid算法。两个整数a...
  • cclsoft
  • cclsoft
  • 2010年05月13日 12:38
  • 2668

[齐次常数线性递推式][多项式取模] LOJ #6017. Shlw loves matrix I

Solution" role="presentation" style="position: relative;">SolutionSolutionSolution 记得以前数学竞赛的同桌大佬告诉我...
  • Vectorxj
  • Vectorxj
  • 2018年01月24日 14:35
  • 77
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:【51nod 1538】 一道难题 - n阶常系数线性递推
举报原因:
原因补充:

(最多只允许输入30个字)