HDU 1575 Tr A 矩陣快速冪

2021-07-04 08:33:19 字數 859 閱讀 8350

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^k)%9973。

2

2 21 0

0 13 99999999

1 2 3

4 5 6

7 8 9

2

2686

順被帶上快速冪演算法::

int quickpow(int m,int n,int k)

return b;

}

#include#include#includeusing namespace std;

const int maxn=15;

const int mod=9973;

int n;

struct node

;node a,b;

node mul(node a,node b)

} }return c;

}node pow_mod(node a,int m)

return b;

}int main()

}b=pow_mod(a,m);

ans=0;

for(i=1;i<=n;i++)

printf("%d\n",ans);

} }return 0;

}

hdu1575 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...

hdu 1575 TR A 矩陣快速冪

b tr a time limit 1000msmemory limit 32768kb64bit io format i64d i64u submit status description a為乙個方陣,則tr a表示a的跡 就是主對角線上各項的和 現要求tr a k 9973。input 資料的...