又写次矩阵模乘,表示看了解题报告后,才知道是矩阵题。练得题少的缘故。应该看懂题就能分析出是什么题。勤加练习,加油!吖飒~~~~
hdu2276,数据比较水,很快就AC了,nyoj一直TLE,后来用结构体写AC啦。大概是初始化和函数的调用浪费的时间。
题意:
一个环形排列的灯,灯的左侧灯状态为1,改变该灯的状态。否则不改变。n轮过后,所有灯的状态。
递推矩阵为:
1 1 0 0 0 0 0
0 1 1 0 0 0 0
0 0 1 1 0 0 0
0 0 0 1 1 0 0
0 0 0 0 1 1 0
1 0 0 0 0 1 1
乘的结果对2取余
hdu AC的代码:
#include<stdio.h>
#include<string.h>
char num[101];
int initmatrix[1][101];//初始矩阵
int matrix[101][101]={0};//递推矩阵
int temp[101][101];
int len;//矩阵行/列数
void calculate(int a[][101],int b[][101],int m,int p)//m*len...len*p
{
int i,j,k;
int res[101][101];
for(i=0;i<m;i++)
{
for(j=0;j<p;j++)
{
res[i][j]=0;
for(k=0;k<len;k++)
res[i][j]+=a[i][k]*b[k][j];
res[i][j]%=2;
}
}
for(i=0;i<m;i++)
for(j=0;j<p;j++)
a[i][j]=res[i][j];
}
void multiply(int temp[][101],int x)
{
if(x==1) return ;
multiply(temp,x>>1);
calculate(temp,temp,len,len);
if(x&1)
calculate(temp,matrix,len,len);
}
int main()
{
//freopen("Input.txt","r",stdin);
//freopen("1.txt","w",stdout);
int n;int i,j;
//初始化递推矩阵
for(i=0;i<101;i++)
{
for(j=i;j<i+2;j++)
{
matrix[i][j]=1;
}
}
while(~scanf("%d",&n))
{
scanf("%s",num);
len=strlen(num);
//初始化初始矩阵
for(i=0;i<len;i++)
{
initmatrix[0][i]=num[i]-'0';
}
//临时递推矩阵
for(i=0;i<len;i++)
{
for(j=0;j<len;j++)
{
temp[i][j]=matrix[i][j];
}
}
temp[len-1][0]=1;matrix[len-1][0]=1;
multiply(temp,n);//结果修改到matrix
calculate(initmatrix,temp,1,len);//结果修改到matrix
for(i=0;i<len;i++)
printf("%d",initmatrix[0][i]);
printf("\n");
matrix[len-1][0]=0;
}
return 0;
}
nyoj的AC代码:、【可以当模板】
#include<iostream>
#include<string.h>
#include<stdio.h>
using namespace std;
#define N 101
#define CLR(data,val) memset(data,val,sizeof(data))
char str[101];
struct Matrix
{
Matrix(){}
void init(int n)
{
row=n;col=n;CLR(data,0);
for(int i=0;i<n;i++)
{
data[i][i]=1;
}
}
void init(int m, int n)
{
row=m;col=n;CLR(data,0);
}
int row,col;
bool data[N][N];
const Matrix &pow(int n);
};
const Matrix &operator *(const Matrix & a,const Matrix &b)
{
Matrix re;
re.init(a.row,b.col);
int i,j,k;
for(i=0;i<a.row;i++)
{
for(j=0;j<a.col;j++)
{
if(a.data[i][j])
{
for(k=0;k<b.col;k++)
{
if(b.data[j][k])
re.data[i][k]=re.data[i][k]^(a.data[i][j]&b.data[j][k]);
}
}
}
}
return re;
}
const Matrix &Matrix::pow(int n)
{
Matrix temp,pre;
temp=*this;pre.init(this->row);
while(n)
{
if(n&1) pre=temp*pre;
temp=temp*temp;
n>>=1;
}
return pre;
}
int main()
{
//freopen("Input.txt","r",stdin);
//freopen("1.txt","w",stdout);
int t,i;
Matrix a,b,res;
while(~scanf("%d",&t))
{
scanf("%s",str);
int len=strlen(str);
a.init(1,len);
for(i=0;i<len;i++)
a.data[0][i]=str[i]-'0';
b.init(len);
for(i=0;i<len;i++)
{
b.data[i][i]=1;
b.data[i][(i+1)%len]=1;
}
res=a*b.pow(t);
for(i=0;i<len;i++)
printf("%d",res.data[0][i]);
printf("\n");
}
return 0;
}