ZOJ Alice's Sequence II //矩阵乘法

矩阵乘法果然犀利！  

#include<iostream>
#include<stdio.h>
#include<string>
using namespace std;
const int man=35;
const int mod=10000007;
char kk[10]="rdasit";
#define ll long long 
class node
{
public:
	ll a[man][man];
};
node date,matrix[12];
char  step[12];
int b[man];
int n,p,m;
int judge(char a[])
{
	for(int i=0;i<6;i++)
		if(a[0]==kk[i]) return i;
	return 0;
}
void init(node &a)
{
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
			a.a  [i][j]=(i==j);
}
void do1(int k)
{
	int temp ;
	scanf("%d",&temp);
	matrix[k].a [temp-1][temp-1]=0;
}
void do2(int k)
{
	int temp;
	scanf("%d",&temp);
	matrix[k].a [temp-1][temp-1]=2;
}
void do3(int k)
{
	int t1,t2;
	scanf("%d%d",&t1,&t2);
	if(t1==t2) matrix[k].a [t1-1][t1-1]=2;
	else matrix[k].a [t2-1][t1-1]=1;
}
void do4(int k)
{
	int t1,t2;
	scanf("%d%d",&t1,&t2);
	matrix[k].a [t1-1][t1-1]=0;
	matrix[k].a [t2-1][t2-1]=0;
	matrix[k].a [t2-1][t1-1]=1;
	matrix[k].a [t1-1][t2-1]=1;
}
void do5(int k)
{
	int t1,t2,i,j;
	scanf("%d%d",&t1,&t2);
	for(j=t2,i=t1;i<=t2;i++,j--)
	{
		matrix[k].a [i-1][i-1]=0;
		matrix[k].a [j-1][i-1]=1;
	}

}
void do6(int k)
{
	int i,j,c,d;
	scanf("%d%d%d%d",&c,&i,&d,&j);
	matrix[k].a [i-1][i-1]=c;
	matrix[k].a [j-1][i-1]=d;
	matrix[k].a [j-1][j-1]=0;
}
void input()
{
	int i;
	//memset(date.a,0,sizeof(date.a));
	for(int i=0;i<n;i++)
		scanf("%lld",&date.a [0][i]);
	scanf("%d",&m);
	for(int i=0;i<m;i++)
	{
		scanf("%s",step);
		int k=judge(step );
		init(matrix[i] );
		switch(k)
		{
		case 0:
			do1(i);
			break;
		case 1:
			do2(i);
			break;
		case 2:
			do3(i);
			break;
		case 3:
			do4(i);
			break;
		case 4:
			do5(i);
			break;
		case 5:
			do6(i);
			break;
		}
	}
	scanf("%d",&p);
}
node mutli(node a,node b)
{
	node re;
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++) re.a [i][j]=0;
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
		{
			int t1;
			for(t1=0;t1<n;t1++)
				re.a [i][j]=(re.a [i][j]+((a.a [i][t1]*b.a[t1][j])%mod))%mod;
		}
	return re;
}
void solve()
{
	node re;
	init(re);
	for(int i=0;i<m;i++)
		re=mutli(re,matrix[i]);
	node ans ;
	init(ans);
	int t=p/m;
	while(t)
	{
		if(t&1) 
			ans= mutli(ans,re);
		t=t>>1;
		re= mutli(re,re);
	}
	t= p%m;
	for(int i=0;i<t;i++)
	{
		ans= mutli(ans,matrix[i] );
	}
	ans =mutli(date,ans);
	for(int i=0;i<n-1;i++)
		cout<<ans.a [0][i]<<" ";
	cout<<ans.a [0][n-1]<<endl;
}
int main()
{
	while(scanf("%d",&n)!=EOF)
	{
		input();
		solve();
	}
	return 0;
}