觀察上式可以得出下面等價方程組
很明顯可以構造乙個矩陣
構造矩陣為
將左邊的矩陣乘n-1次
#include using namespace std;
typedef long long llt;
const int cube_size = 3;///矩陣大小
struct cube;
//單位矩陣
cube _unitcube = ;
///矩陣乘機
cube _multiply(cube a,cube b,llt mod)
}return _tmpcube;
}///矩陣快速冪
cube power_cube(cube a,llt n,llt mod) //矩陣快速冪
return _tmpcube;
}int main();
//單位矩陣
cube _unitcube = ;
///矩陣乘機
cube _multiply(cube a,cube b,llt mod)
}return _tmpcube;
}///矩陣快速冪
cube power_cube(cube a,llt n,llt mod) //矩陣快速冪
return _tmpcube;
}int main();
//單位矩陣
cube _unitcube = ;
cube tmp = ;
///矩陣乘機
cube _multiply(cube a,cube b,llt mod)
}return _tmpcube;
}///矩陣快速冪
cube power_cube(cube a,llt n,llt mod) //矩陣快速冪
return _tmpcube;
}int main()
arr = power_cube(tmp,k - 9ll,mod);
llt ans = 0;
for (int i = 0;i < 10;++i)
cout <
矩陣快速冪(模板 例題)
矩陣快速冪推導過程 矩陣快速冪求解數列第n項的關鍵在於計算係數矩陣a。之後就是套模板了。模板 求解斐波那契數列第n項 include include define ll long long define mod 2147493647 using namespace std struct mat ma...
矩陣快速冪模板 例題
矩陣快速冪這個演算法,理解起來很容易,但是我之前自己寫的 有bug,也是因為上課不聽課,對形參和實參沒理解,平常用的都是全域性變數,是不是實參影響不大,這次定義乙個結構體的矩陣,矩陣需要初始化為0,然後,因為形參和實參沒怎麼理解,導致輸出的答案差異很大,前提是矩陣快速冪,矩陣需要初始化陣列,並不會預...
矩陣快速冪模板 例題(快速入門)
本來早在乙個月前就應該把這給掌握的,硬是得等到卡題了,才想到來補知識 傳送門 矩陣快速冪基礎講解 單純的求某一矩陣的快速冪是沒有什麼意義的,主要是由於所求的遞推式的數比較大時,如果使用暴力遞推,必定會超時,這時矩陣快速冪的作用就完美體現了 根據個人碼風和編碼習慣,敲了乙份模板,如果覺得有用,大佬們不...