兩個矩陣間的乘法
struct node
;node jz
(int n,node a,node b)
return ans1;
}
矩陣的快速冪
如果快速冪不懂得可以點 快速冪取模詳解
node ksqm
(int n,node a,
int cf)
return ans;
}
a為乙個方陣,則tr a表示a的跡(就是主對角線上各項的和),現要求tr(a^k)%9973。
input
資料的第一行是乙個t,表示有t組資料。
每組資料的第一行有n(2 <= n <= 10)和k(2 <= k < 10^9)兩個資料。接下來有n行,每行有n個資料,每個資料的範圍是[0,9],表示方陣a的內容。
output
對應每組資料,輸出tr(a^k)%9973。
sample input
22 2
1 00 1
3 99999999
1 2 3
4 5 6
7 8 9
sample output
22686
#include
using
namespace std;
#define long long ll;
const
int mod=
9973
;struct node
;node jz
(int n,node a,node b)
return ans1;
}node ksqm
(int n,node a,
int cf)
return ans;
}int
main()
an=ksqm
(n,an,k)
;int jl=0;
for(
int i=
1;i<=n;i++
) cout<}return0;
}
矩陣快速冪 例題 模板
觀察上式可以得出下面等價方程組 很明顯可以構造乙個矩陣 構造矩陣為 將左邊的矩陣乘n 1次 include using namespace std typedef long long llt const int cube size 3 矩陣大小 struct cube 單位矩陣 cube unitc...
矩陣快速冪(模板 例題)
矩陣快速冪推導過程 矩陣快速冪求解數列第n項的關鍵在於計算係數矩陣a。之後就是套模板了。模板 求解斐波那契數列第n項 include include define ll long long define mod 2147493647 using namespace std struct mat ma...
矩陣快速冪模板 例題
矩陣快速冪這個演算法,理解起來很容易,但是我之前自己寫的 有bug,也是因為上課不聽課,對形參和實參沒理解,平常用的都是全域性變數,是不是實參影響不大,這次定義乙個結構體的矩陣,矩陣需要初始化為0,然後,因為形參和實參沒怎麼理解,導致輸出的答案差異很大,前提是矩陣快速冪,矩陣需要初始化陣列,並不會預...