用矩陣快速冪求解遞推數(shù)列的第n項(xiàng)

斐波那契數(shù)列

  • 每一項(xiàng)的數(shù)字等于前兩項(xiàng)的和, f_0=0, f_1 =1, f_2=1
  1. 遞推公式
    f_i = f_{i-1} + f_{i-2}
  2. 矩陣形式
    [f_i,f_{i-1}] = [f_{i-1}, f_{i-2}] * \begin{bmatrix} 1&1 \\ 1& 0 \end{bmatrix} = [1, 0] * \begin{bmatrix} 1&1 \\ 1& 0 \end{bmatrix} ^ {i-1}

擊鼓傳花

  • 題目描述:
    一共有k個(gè)人, 每一次可以把自己手中的花傳給其它任何一個(gè)人, 初始時(shí)在自己手中, 求M次后有多少種方式回到自己手中
  1. 遞推公式
    a_i表示第i次在自己手中, b_i表示第i次不在自己手中, a_0=1, b_0=0
  2. 矩陣形式
    [a_i,b_i] = [a_{i-1}, b_{i-1}] * \begin{bmatrix} 0&k-1 \\ 1& k-2 \end{bmatrix} = [1, 0] * \begin{bmatrix} 0&k-1 \\ 1& k-2 \end{bmatrix} ^ i

求數(shù)列的第n項(xiàng)

  1. 遞推公式
    f_i = a f_{i-1} + bf_{i-2} + c f_{i-3} + 2 i^2 - i + 32767

  2. 矩陣形式
    [f_i, f_{i-1}, f_{i-2}, (i+1)^2, i+1, 1] = [f_{i-1}, f_{i-2}, f_{i-3}, i^2, i, 1] * \begin {bmatrix} a & 1 & 0 & 0 & 0 & 0\\ b & 0 & 1 & 0 & 0 & 0\\ c & 0 & 0 & 0 & 0 & 0\\ 2 & 0 & 0 & 1 & 0 & 0\\ -1 & 0 & 0 & 2 & 1 & 0\\ 32767 & 0 & 0 & 1 & 1 & 1 \end {bmatrix}
    首先將向量改寫成concat([f_{i-1}, f_{i-2}, \cdots, f_0],[i^t, i ^{t-1}, \cdots, i, 1])的形式, 即向量中只有系數(shù)為1數(shù)列的前幾項(xiàng)或者i的多項(xiàng)式, 然后待定系數(shù)求轉(zhuǎn)移矩陣

矩陣快速冪的參考代碼


public class MatQuickPower {

    public static final int MOD = 1_000_000_007;

    public static int[][] matQuickPower(int [][]a, int k) {
        if (k == 1) return a;
        int [][]tmp = matQuickPower(a, k >> 1);
        if ((k & 1) == 0) return matmul(tmp, tmp);
        else return matmul(matmul(tmp, tmp), a);
    }

    public static int[][] matmul(int [][]a, int [][]b) {
        int [][]c = new int[a.length][b[0].length];
        for (int i = 0; i < a.length; i++) {
            for (int j = 0; j < b[0].length; j++) {
                for (int k = 0; k < a[i].length; k ++){
                    c[i][j] = (int) ((c[i][j] + (long)(a[i][k] % MOD) * (long)(b[k][j] % MOD)) % MOD);
                }
            }
        }
        return c;
    }

}

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個(gè)濱河市补疑,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌莲组,老刑警劉巖,帶你破解...
    沈念sama閱讀 206,968評(píng)論 6 482
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件锹杈,死亡現(xiàn)場(chǎng)離奇詭異,居然都是意外死亡竭望,警方通過查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,601評(píng)論 2 382
  • 文/潘曉璐 我一進(jìn)店門闭专,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人影钉,你說我怎么就攤上這事∑轿” “怎么了?”我有些...
    開封第一講書人閱讀 153,220評(píng)論 0 344
  • 文/不壞的土叔 我叫張陵廉赔,是天一觀的道長。 經(jīng)常有香客問我昂勉,道長,這世上最難降的妖魔是什么岗照? 我笑而不...
    開封第一講書人閱讀 55,416評(píng)論 1 279
  • 正文 為了忘掉前任,我火速辦了婚禮攒至,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘迫吐。我一直安慰自己,他們只是感情好志膀,可當(dāng)我...
    茶點(diǎn)故事閱讀 64,425評(píng)論 5 374
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著烫止,像睡著了一般。 火紅的嫁衣襯著肌膚如雪戳稽。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 49,144評(píng)論 1 285
  • 那天互躬,我揣著相機(jī)與錄音,去河邊找鬼吼渡。 笑死,一個(gè)胖子當(dāng)著我的面吹牛诞吱,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播竭缝,決...
    沈念sama閱讀 38,432評(píng)論 3 401
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場(chǎng)噩夢(mèng)啊……” “哼咙俩!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起阿趁,我...
    開封第一講書人閱讀 37,088評(píng)論 0 261
  • 序言:老撾萬榮一對(duì)情侶失蹤,失蹤者是張志新(化名)和其女友劉穎坛猪,沒想到半個(gè)月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體墅茉,經(jīng)...
    沈念sama閱讀 43,586評(píng)論 1 300
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡呜呐,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 36,028評(píng)論 2 325
  • 正文 我和宋清朗相戀三年悍募,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片坠宴。...
    茶點(diǎn)故事閱讀 38,137評(píng)論 1 334
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡,死狀恐怖喜鼓,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情颠通,我是刑警寧澤膀懈,帶...
    沈念sama閱讀 33,783評(píng)論 4 324
  • 正文 年R本政府宣布,位于F島的核電站启搂,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏胳赌。R本人自食惡果不足惜牢撼,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,343評(píng)論 3 307
  • 文/蒙蒙 一熏版、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧撼短,春花似錦、人聲如沸挺勿。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,333評(píng)論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至蚊丐,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間麦备,已是汗流浹背讲竿。 一陣腳步聲響...
    開封第一講書人閱讀 31,559評(píng)論 1 262
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留题禀,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 45,595評(píng)論 2 355
  • 正文 我出身青樓迈嘹,卻偏偏與公主長得像,于是被迫代替她去往敵國和親秀仲。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 42,901評(píng)論 2 345

推薦閱讀更多精彩內(nèi)容