今天做的第二道矩陣快速冪題,因為是初次接觸,各種奇葩錯誤整整除錯了一下午。廢話不說,入正題。該題應該屬於矩陣快速冪的裸題了吧,知道快速冪原理(二進位制迭代法,非遞迴版)後,剩下的只是處理矩陣乘法的功夫了,我直接用個結構體來表示矩陣,確實能省去不少功夫(這裡一定要注意用單位矩陣來初次相乘,但不要把它放進建構函式中,我就是在這裡卡了好久
)。下面附上**:
1 #include2 #include3const
int mod= 9973;4
5struct
matrix
11void
identity()
15 matrix operator *(const
matrix m2)
23};
2425 matrix quick_mod(matrix &m, int
p)26
34return
ans;35}
3637
intmain()
3854
return0;
55 }
hdu1575 Tr A(矩陣快速冪)
今天做的第二道矩陣快速冪題,因為是初次接觸,各種奇葩錯誤整整除錯了一下午。廢話不說,入正題。該題應該屬於矩陣快速冪的裸題了吧,知道快速冪原理 二進位制迭代法,非遞迴版 後,剩下的只是處理矩陣乘法的功夫了,我直接用個結構體來表示矩陣,確實能省去不少功夫 這裡一定要注意用單位矩陣 來初次相乘,但不要把它...
HDU 1575 Tr A 矩陣快速冪
a為乙個方陣,則tr a表示a的跡 就是主對角線上各項的和 現要求tr a k 9973。資料的第一行是乙個t,表示有t組資料。每組資料的第一行有n 2 n 10 和k 2 k 10 9 兩個資料。接下來有n行,每行有n個資料,每個資料的範圍是 0,9 表示方陣a的內容。對應每組資料,輸出tr a ...
HDU 1575 Tr A(矩陣快速冪)
思路 簡單的矩陣快速冪 include include includeusing namespace std int n,mod struct mat mat mul mat a,mat b return t mat expo mat p,int k return e int add mat a i...