题目链接

题目大意

$$F[i]=F[i-1]+F[i-2]\ (\ F[1]=1\ ,\ F[2]=1\ )$$

$$T[i]=F[1]+2F[2]+3F[3]+...+nF[n]$$

求$T[n]\ mod\ m$

$n,m<=2^{31}-1$

这题的递推式推导有点神仙,完全想不到多用两个数组来形成递推式。研究了一下一本通上面的两个辅助数组的用途然后才会推出来这个转移矩阵

还是太菜了

题解

由题目可得

 $$F[i]=F[i-1]+F[i-2]\ (\ F[1]=1\ ,\ F[2]=1\ )$$

$$T[i]=F[1]+2F[2]+3F[3]+...+nF[n]$$

然而$T[i]$并不可以形成递推式,所以可以设

$$nS[i]=nF[1]+nF[2]+nF[3]+...+nF[n]$$

\begin{align*}
p[i]&=nS[i]-T[i]\\
&=(n-1)F[1]+(n-2)F[2]+...+F[n]\\
&=p[i-1]+s[i-1]
\end{align*}

这个推导可以把$p[i-1]$和$S[i-1]$以及$p[i]$展开来,就可以很好理解为什么可以这样子化成递推式了

于是我们就得到了一个递推式

 $$p[i]=p[i-1]+s[i-1]$$

而答案即为$T[n]=nS[n]-p[N]$

有了递推式就可以使用矩阵乘法优化:

转移矩阵随便推推就出来了,也不难推

$$\left[
\begin {matrix}
p[i]\\
S[i]\\
F[i]\\
F[i-1]
\end {matrix}
\right]
*
\left[
\begin{matrix}
1&1&0&0\\
0&1&1&0\\
0&0&1&1\\
0&0&1&0
\end{matrix}
\right]
=
\left[
\begin{matrix}
p[i+1]\\
S[i+1]\\
F[i+1]\\
F[i]
\end{matrix}
\right]$$

然后就是套板子了

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstdlib>
#include <iostream>
#include <vector>
#include <cmath>
#include <set>
#include <map> using namespace std ; #define I_int int
inline I_int read() {
I_int x = , f = ; char c = getchar() ;
while( c < '' || c > '' ) {
if( c == '-' ) f = - ;
c = getchar() ;
}
while( c >= '' && c <= '' ) {
x = x * + c -'' ;
c = getchar() ;
}
return x * f ;
}
#undef I_int #define ll long long
#define inf 0x3f3f3f3f
#define in(a) a = read()
#define out(a) printf( "%d " , a )
#define outn(a) printf( "%d\n" , a )
#define N 100010 ll mod ;
struct matrix {
ll m[ ][ ] ;
matrix() { memset( m , , sizeof( m ) ) ; }
ll *operator[] ( ll a ) { return m[ a ] ; }
matrix operator * ( matrix &x ) {
matrix ans ;
memset( ans.m , , sizeof( ans.m ) ) ;
for( int i = ; i < ; i ++ ) {
for( int j = ; j < ; j ++ ) {
for( int k = ; k < ; k ++ ) {
ans[ i ][ j ] = ( ans[ i ][ j ] + m[ i ][ k ] * x[ k ][ j ] % mod ) % mod ;
}
}
}
return ans ;
}
} A , B ; int n = read() ; void init() {
A[][]=A[][]=A[][]=A[][]=A[][]=A[][]=A[][]=;
B[][]=B[][]=B[][]=;
} matrix power( ll p ) {
matrix ans , base = A ;
for( int i = ; i < ; i ++ ) ans[ i ][ i ] = ;
while( p ) {
if( p & ) ans = ans * base ;
base = base * base ;
p >>= ;
}
return ans ;
} int main() {
scanf( "%lld" , &mod ) ;
init() ;
A = power( n - ) ;
matrix ans = A * B ;
printf( "%lld\n" , (n*ans[][]-ans[][]+mod)%mod ) ;
}

LOJ #10222. 「一本通 6.5 例 4」佳佳的 Fibonacci的更多相关文章

  1. LOJ #10222. 「一本通 6.5 例 4」佳佳的 Fibonacci 题解

    题目传送门 如果之前推过斐波那契数列前缀和就更好做(所以题目中给出了). 斐波那契数列前缀和题目链接 先来推一下斐波那契数列前缀和: \[\sum\limits_{i=1}^nf(i) \] 其中 \ ...

  2. LOJ#10064. 「一本通 3.1 例 1」黑暗城堡

    LOJ#10064. 「一本通 3.1 例 1」黑暗城堡 题目描述 你知道黑暗城堡有$N$个房间,$M$条可以制造的双向通道,以及每条通道的长度. 城堡是树形的并且满足下面的条件: 设$D_i$为如果 ...

  3. LOJ #10131 「一本通 4.4 例 2」暗的连锁

    LOJ #10131 「一本通 4.4 例 2」暗的连锁 给一棵 \(n\) 个点的树加上 \(m\) 条非树边 , 现在需要断开一条树边和一条非树边使得图不连通 , 求方案数 . $n \le 10 ...

  4. Loj 10115 「一本通 4.1 例 3」校门外的树 (树状数组)

    题目链接:https://loj.ac/problem/10115 题目描述 原题来自:Vijos P1448 校门外有很多树,学校决定在某个时刻在某一段种上一种树,保证任一时刻不会出现两段相同种类的 ...

  5. LOJ#10065. 「一本通 3.1 例 2」北极通讯网络

    题目链接:https://loj.ac/problem/10065 题目描述 原题来自:Waterloo University 2002 北极的某区域共有 nnn 座村庄,每座村庄的坐标用一对整数 ( ...

  6. LOJ#10106. 「一本通 3.7 例 2」单词游戏

    题目链接:https://loj.ac/problem/10106 题目描述 来自 ICPC CERC 1999/2000,有改动. 有 NNN 个盘子,每个盘子上写着一个仅由小写字母组成的英文单词. ...

  7. LOJ #10132. 「一本通 4.4 例 3」异象石

    题目地址 LOJ 题解 神仙思路.思路参考自<算法竞赛进阶指南>. 考虑维护dfs序中相邻两个石头的距离,那么每次?的答案就是sum/2(首尾算相邻) 然后维护一下拿个平衡树/set维护一 ...

  8. LOJ 10138 -「一本通 4.5 例 1」树的统计

    树链剖分模板题,详见这篇博客.

  9. LOJ 10155 - 「一本通 5.2 例 3」数字转换

    前言 从现在开始,这个博客要写一些题解了.起初,开这个博客只是好玩一样,没事就写写CSS.JS,然后把博客前端搞成了现在这个样子.以前博客只是偶尔记录一些东西,刷题也从来不记录,最近受一些学长的影响, ...

随机推荐

  1. 【笔记】mysql两条数据的某个属性值互换

    update groupuser as g1 join groupuser as g2 on (g1.user_id=1 and g2.user_id = 2) or(g1.user_id = 2 a ...

  2. Cosh.3

    查壳.没有 拖 OD 查找字符串 找到有用的东西   顺线往上看       找到  大概算入口处  下断 就从这里开始  单步下去吧 name的变幻 Serial的变幻 直接翻译出来了   继续看看 ...

  3. 批量修改java文件的包名

    需求:我复制进批量的java文件,但是包名需要用现在创建的包名 导入进来时,由于包名不一致会报错 解决办法,点击包名鼠标右键>Refactor>Rename 结果如下ok:

  4. Xcode的代码片段快捷方式-Code Snippet Library(代码片段库)

    最近换了新电脑,装上Xcode敲代码发现很多以前攒的Code Snippet忘记备份了,总结了一下Code Snippet的设置方法,且行且添加,慢慢积累吧. 如下图:   Title - Code ...

  5. POJ1032 Parliament(数论)

    New convocation of The Fool Land's Parliament consists of N delegates. According to the present regu ...

  6. 【cogs 597】【dp】交错匹配

    597. 交错匹配 ★☆ 输入文件:crossa.in 输出文件:crossa.out 简单对照 时间限制:1 s 内存限制:128 MB [问题描写叙述] 有两行自然数. UP[1..N] . DO ...

  7. angularjs 控制器、作用域、广播详解

    一.控制器 首先列出几种我们平常使用控制器时的几种误区: 我们知道angualrJs中一个控制器时可以对应不同的视图模板的,但这种实现方式存在的问题是: 如果视图1和视图2根本没有任何逻辑关系,这样& ...

  8. php导出超大csv导出方法,读取超大文件或者接受超大数组,防止内存溢出

    基本思路就是,知道总数之后分割成2万一个数组进行查询,最后独立写入csv,避免数据过大导致溢出 速度还不错,在php7下,机器I5 8G内存,128G,SSD,52W多条,大概也就30秒,出来整个文件 ...

  9. Ubuntu18.04 安装tomcat9

    1.官网下载 2.移动到/usr/local/tomcat 3.解压 4.修改权限,否则在idea中不能正常使用

  10. memory prefix pre,pro,penta,pseudo out _p 1

    1● pre 前的,预先   2● pro 在前,很多,赞同       3● penta 5,五     4● pseaudo   伪,假 pseudo