Tr A hdu 1575

/* 其实跟自己以前写的是貌似没有区别的 */ #include <iostream> #include <cstdio> #include <cstring> using namespace std; int N; const int mod = 9973; struct Mat { int num[11][11]; Mat() { for(int i = 0; i < 11; i++) for(int j = 0; j < 11; j++) num[i][j] = 0; } }; Mat init, unit; Mat mul(Mat a, Mat b) {// 听大牛说这样可以优化，但是这样一来感觉理解起来有点要绕弯，还是习惯自己传统的写法 Mat ans; for(int i = 0; i < N; i++) { for(int k = 0; k < N; k++) { //ans.num[i][j] = 0; if(a.num[i][k] == 0) continue; for(int j = 0; j < N; j++) { if(b.num[k][j] == 0) continue; ans.num[i][j] += a.num[i][k] * b.num[k][j]; if(ans.num[i][j] >= mod) ans.num[i][j] %= mod; } //ans.num[i][j] %= mod; } } return ans; } Mat mal(int n) { Mat p, q; q = unit; p = init; while(n) { //if(n % 2 == 1) if(n & 1) { p =mul(p, q); n--; } else { q = mul(q, q); //n /= 2; n >>= 1; } } return p; //return mul(p, q); } int main() { int t; scanf("%d", &t); while(t--) { int n, k; scanf("%d %d", &n, &k); for(int i = 0; i < n; i++) for(int j = 0; j < n; j++) { scanf("%d", &unit.num[i][j]); init.num[i][j] = (i == j); } N = n; Mat ans; ans = mal(k); int sum = 0; for(int i = 0; i < n; i++) sum += ans.num[i][i]; printf("%d/n", sum % mod); } }

http://acm.hdu.edu.cn/showproblem.php?pid=1575

#include<iostream>//典型的矩阵二分，只要记住模板，稍微修改就OK了
#include<cstdio>
#include<string>
using namespace std;
int MAX;
const int mod = 9973;

struct Mat
{
  int df[11][11];
}init, unit;

Mat Mul(Mat a,Mat b)
{
 int i, j, k;
 Mat c;
 for(i = 0; i < MAX; i++)
  for(j = 0; j < MAX; j++)
    {
     c.df[i][j] = 0;
     for(k = 0; k < MAX; k++)
     c.df[i][j] += a.df[i][k] * b.df[k][j];
     if(c.df[i][j] >= mod)
       c.df[i][j] %= mod;
                 
    }
 return c;
}

Mat cal(int t)
{
 Mat p, q;
 p = init;
 q = unit;
 while(t != 1)
 {
  if(t&1)
  {
   p = Mul(p, q);
   t--;
  }
  else
  {
   q = Mul(q, q);
   t >>= 1;
  }
 }
 p = Mul(p, q);
 return p;
}

int main()
{
 int i, j, k, sum, t;
 Mat tt;
 scanf("%d", &t);
 while(t--)
 {
  scanf("%d%d",&MAX, &k);
  //memset(unit,0,sizeof(unit));
  memset(unit.df,0,sizeof(unit.df));
  memset(init.df,0,sizeof(init.df));
  for(i = 0; i < MAX; i++)
   for(j = 0; j < MAX; j++)
  {
   scanf("%d",&unit.df[i][j]);
   init.df[i][i] = 1;
   
  }
  tt = cal(k);
  sum = 0;
  for(i = 0; i < MAX ; i++)
   for(j = 0; j < MAX; j++)
   if(i == j)
   sum += tt.df[i][j] % mod;
  printf("%d/n", sum % mod);
 }
 return 0;
}