Tr A

Time Limit: 1000/1000 MS (Java/Others)????Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 6348????Accepted Submission(s): 4707


Problem Description A為一個(gè)方陣，則Tr A表示A的跡（就是主對(duì)角線上各項(xiàng)的和），現(xiàn)要求Tr(A^k)%9973。

Input 數(shù)據(jù)的第一行是一個(gè)T，表示有T組數(shù)據(jù)。
每組數(shù)據(jù)的第一行有n(2 <= n <= 10)和k(2 <= k < 10^9)兩個(gè)數(shù)據(jù)。接下來(lái)有n行，每行有n個(gè)數(shù)據(jù)，每個(gè)數(shù)據(jù)的范圍是[0,9]，表示方陣A的內(nèi)容。

Output 對(duì)應(yīng)每組數(shù)據(jù)，輸出Tr(A^k)%9973。
Sample Input 2 2 2 1 0 0 1 3 99999999 1 2 3 4 5 6 7 8 9
Sample Output 2 2686
Author xhd
Source HDU 2007-1 Programming Contest
Recommend

linle???|???We have carefully selected several similar problems for you:??1757?1588?2256?2604?2254?

#include <iostream> #include <cstdio> #include <cstring> //#include <bits/stdc++.h> using namespace std; int n; long long k,mod=9973; struct matrix {long long x[15][15]; }a,e; void init() {memset(e.x,0,sizeof(e.x));for(int i=0;i<n;++i)e.x[i][i]=1; } matrix POw(matrix xx,matrix yy) {matrix sum;memset(sum.x,0,sizeof(sum.x));for(int i=0;i<n;++i){for(int j=0;j<n;++j){for(int k=0;k<n;++k){if(xx.x[i][k]&&yy.x[k][j])sum.x[i][j]+=(xx.x[i][k]%mod*yy.x[k][j]%mod)%mod;//sum.x[i][j]+=((xx.x[i][k]%mod)*yy.x[k][j]%mod)%mod;}}}return sum; } void martix_pow(long long k) {init();while(k){if(k&1){e=POw(e,a);}a=POw(a,a);k>>=1;}/* for(int i=0;i<n;++i){for(int j=0;j<n;++j)cout<<e.x[i][j]<<' ';cout<<endl;}*/long long ans=0;for(int i=0;i<n;++i)ans+=e.x[i][i]%mod;cout<<(ans%mod)<<endl; } int main() {int t;//freopen("in.txt","r",stdin);// freopen("out.txt","w",stdout);cin>>t;while(t--){cin>>n>>k;for(int i=0;i<n;++i)for(int j=0;j<n;++j){scanf("%d",&a.x[i][j]);}martix_pow(k);//<<endl;}return 0; }

總結(jié)

以上是生活随笔為你收集整理的hdu 1575Tr A的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。