标签：des style blog http color java os io
4 2 5 5 4 4 5 4 0 0 4 2 5 5 1 3 1 5 6 3 1 2 3 0 3 0 2 3 4 4 3 2 2 5 5 0 5 0 3 4 5 1 1 0 5 3 2 3 3 2 3 1 5 4 5 2 0 0
14 56
题解及代码：
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
int l[1010][7];
int r[7][1010];
int nt[1010][1010];
int bt[1010][1010];
struct mat
{
int t[7][7];
void set()
{
memset(t,0,sizeof(t));
}
} a,b,c;
mat multiple(mat a,mat b,int n,int p)
{
int i,j,k;
mat temp;
temp.set();
for(i=0; i<n; i++)
for(j=0; j<n; j++)
{
if(a.t[i][j]!=0)
for(k=0; k<n; k++)
temp.t[i][k]=(temp.t[i][k]+a.t[i][j]*b.t[j][k]+p)%p;
}
return temp;
}
mat quick_mod(mat b,int n,int m,int p)
{
mat t;
t.set();
for(int i=0;i<n;i++) t.t[i][i]=1;
while(m)
{
if(m&1)
{
t=multiple(t,b,n,p);
}
m>>=1;
b=multiple(b,b,n,p);
}
return t;
}
int main()
{
int n,m;
while(scanf("%d%d",&n,&m)!=EOF)
{
if(n==0&&m==0) break;
for(int i=0;i<n;i++)
{
for(int j=0;j<m;j++)
scanf("%d",&l[i][j]);
}
for(int i=0;i<m;i++)
{
for(int j=0;j<n;j++)
scanf("%d",&r[i][j]);
}
a.set();
for(int i=0; i<m; i++)
for(int j=0; j<n; j++)
{
if(r[i][j]!=0)
for(int k=0; k<m; k++)
a.t[i][k]=(a.t[i][k]+r[i][j]*l[j][k])%6;
}
b=quick_mod(a,m,n*n-1,6);
/*for(int i=0;i<m;i++)
{
for(int j=0;j<m;j++)
printf("%d ",b.t[i][j]);
puts("");
}
puts("");*/
//memset(nt,0,sizeof(nt));
//memset(bt,0,sizeof(bt));
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
nt[i][j]=bt[i][j]=0;
for(int i=0; i<n; i++)
for(int j=0; j<m; j++)
{
if(l[i][j]!=0)
for(int k=0; k<m; k++)
nt[i][k]=(nt[i][k]+l[i][j]*b.t[j][k])%6;
}
for(int i=0; i<n; i++)
for(int j=0; j<m; j++)
{
if(nt[i][j]!=0)
for(int k=0; k<n; k++)
bt[i][k]=(bt[i][k]+nt[i][j]*r[j][k])%6;
}
int ans=0;
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
ans+=bt[i][j];
printf("%d\n",ans);
}
return 0;
}
hdu 4965 Fast Matrix Calculation,布布扣,bubuko.com
hdu 4965 Fast Matrix Calculation
标签：des style blog http color java os io
原文地址：http://blog.csdn.net/knight_kaka/article/details/38686737
