逛了一圈。。。觉得这篇讲的比较清楚:传送门~
简要概括:
1、线性代数的知识,单位矩阵的利用;(如果不知道单位矩阵的,先去补习一下线代,做几题行列式就会了);
2、然后构造好矩阵以后,直接做M次乘积运算,然后利用一个[1, 0 , ... , 0 ]矩阵提取答案就好了,对,那个[1,0,..., 0 ] 就是获取答案的作用。
PS
以前的矩阵快速幂总是这样写的:先构造一个矩阵,然后跑一发矩阵快速幂,然后这个矩阵a乘以另一个矩阵b就会得出答案,而答案矩阵往往会包括前一个答案的值。所以自己就有点不知所措。【这里写自己的智障困惑。。】
#include <iostream>
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
typedef long long LL;
struct asd{
LL a[110][110];
};
int n;
asd mul(asd a, asd b)
{
asd ans;
memset(ans.a,0,sizeof(ans.a));
for(int k=0;k<=n;k++)
for(int i=0;i<=n;i++)
if(a.a[i][k])
for(int j=0;j<=n;j++)
ans.a[i][j]+=a.a[i][k]*b.a[k][j];
return ans;
}
asd quickmul(int g,asd x)
{
asd ans;
memset(ans.a,0,sizeof(ans.a));
for(int i=0;i<=n;i++)
ans.a[i][i]=1;
while(g)
{
if(g&1) ans=mul(ans,x);
x=mul(x,x);
g>>=1;
}
return ans;
}
asd init(int k)
{
asd tmp;
memset(tmp.a,0,sizeof(tmp.a));
for(int i=0;i<=n;++i)
tmp.a[i][i]=1;
char x[5];
int a,b;
while(k--)
{
scanf("%s",x);
if(x[0]=='g')
{
scanf("%d",&a);
tmp.a[0][a]++;
}
else if(x[0]=='s')
{
scanf("%d%d",&a,&b);
for(int i=0;i<=n;i++)
swap(tmp.a[i][a],tmp.a[i][b]);
}
else
{
scanf("%d",&a);
for(int i=0;i<=n;i++)
tmp.a[i][a]=0;
}
}
return tmp;
}
int main()
{
int m,k;
asd ans,tmp;
while(~scanf("%d%d%d",&n,&m,&k))
{
if((!n)&&(!m)&&(!k))
break;
tmp=init(k);
ans=quickmul(m,tmp);
memset(tmp.a,0,sizeof(tmp.a));
tmp.a[0][0]=1;
ans=mul(tmp,ans);
for(int i=1;i<=n;i++)
{
if(i!=1) printf(" ");
printf("%lld",ans.a[0][i]);
}
puts("");
}
return 0;
}