HDU 2276 小小矩阵乘法

这些天刷了蛮多的矩阵题，感觉对思维的训练还是很大的，一向不喜欢做题的我这几天被惊奇郭虐菜了...

果断跟进做题！

题目大意：

有一排数列用0,1构成，若当前数的左边为1，则当前数改变状态这样转一圈。为一次。

那我们就可以这么构造矩阵了。

先看这个矩阵：

|010|

|001|=B

|100|

用A=[101]去乘，可见第一列的1与A最后的1有关系，也就是左边为1状态改变的意思。

再多多看看就是了。

然后由于变化与不变化，与自身有关系只要B+E就可以了。

然后矩阵N次幂..... soso.....

Time:62ms 还不错啦...

#include<iostream>
#include<string.h>
#include<cstdio>
using namespace std;

struct node
{	int s[111]; };
struct nodem
{	int m[111][111]; }res,temp;

int n,len;

nodem matriXmult( nodem a,nodem b )
{	nodem c;
	memset( c.m,0,sizeof(c.m) );
	for( int i=0;i<len;i++ )
	for( int k=0;k<len;k++ )
	if( a.m[i][k] )
	for( int j=0;j<len;j++ )
	 	 c.m[i][j]+=a.m[i][k]&b.m[k][j];
 	for( int i=0;i<len;i++ )
	for( int j=0;j<len;j++ )
		 c.m[i][j]=c.m[i][j]&1;
	return c;
}

void matrix_Power()
{	memset( res.m,0,sizeof(res.m) );
 	memset( temp.m,0,sizeof(temp.m) );
 	
 	for( int i=0;i<len;i++ )
 	{	temp.m[i][i]=res.m[i][i]=1,
  		temp.m[i][(i+1)%len]=1;
	}//out(res),out(temp);
 	for( int i=0;n>=(1<<i);i++ )
	{	if( n&(1<<i) )
			res=matriXmult(res,temp);
  	  	temp=matriXmult(temp,temp);
	}
}
  	  	
int main()
{	while( scanf("%d",&n)!=EOF )
 	{	char str[111];
 		scanf( "%s",&str );
	 	len=strlen(str);
	 	node src,dst;
		 for( int i=0;i<len;i++ )
   		 	  src.s[i]=(str[i]=='1')?1:0;
   		 matrix_Power();
   		 //out(res);
   		 memset( dst.s,0,sizeof(dst.s) );
   		 for( int j=0;j<len;j++ )
   		 for( int k=0;k<len;k++ )
   		 	  dst.s[j]+=src.s[k]*res.m[k][j];
   		 for( int i=0;i<len;i++ )
			  printf( "%d",dst.s[i]&1 );
		 printf( "\n" );
  	}
 	return 0;
}