矩阵

Q1：何为矩阵？
一个mn的由整数构成的图案：
如：32的矩阵：
[ 1 2 ]
[ 3 4 ]
[ 5 6 ]
矩阵的加减法：

1 2   11 12   12 14
3 4 + 13 14=  16 18
5 6   15 16   20 22

结论：矩阵加or减法：
对于两个分别为A[n][m]和B[n][m]的矩阵相加：
答案C[i][j]=A[i][j]+B[i][j]。
减法则为:C[i][j]=A[i][j]-B[i][j]
矩阵矩阵：
要求:
A:mn B:n*k

1 2  
3 4 * 8 9 10 11
5 6   12 13 14 15
=
3行2列*2行4列=3行4列
32 35 38 41
72 79 86 93
112 123 134 145
如：32=1*8+2*12
	35=1*9+2*13
A*B=C:
for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
		for(int k=1;k<=n;k++)
			C[i][j]=a[i][k]*b[k][j]

矩阵的乘法：
1.不满足交换律:A * B != B * A
2.满足结合律:A * B * C= A * ( B * C )
3.满足分配律: ( A + B ) * C = A * C + B * C
例题：求a的n次方：原理：分治
a的n次
a的1/2次 a的1/2次

快速幂：
3的11次:
(11)10=(1011)2
11=2的立方+2的一次方+2的零次方

斐波那契数列：

f(n-1)=f(n-2)* 0 1
			   1 1
f(n)=f(1)*0 1*0 1*-----(n-1)个
		  1 1 1 1

练习题：
I.矩阵快速幂
两种思路:
1.分治：如下：

		a的n次
a的1/2次	a的1/2次

递归求解
代码:

#include <bits/stdc++.h>
using namespace std;
long long a,n,k;
long long helper(long long a,long long n){
	if(n==1)
		return a%k;
	long long tmp=helper(a,n/2)%k;
	if(n%2==0)
		return tmp%k*tmp%k;
	else
		return tmp%k*tmp%k*a%k;
}
int main(){
	cin>>a>>n>>k;
	cout<<a<<"^"<<n<<" "<<"mod"<<" "<<k<<"="<<helper(a,n)%k;
	return 0;	
}

2.利用位运算优化上述过程：
代码:

#include <bits/stdc++.h>
using namespace std;
long long a,n,k;
long long helper(long long a,long long n){
	long long base=a;
    long long ans=1;
    while(n){
    	if(n&1)
    		ans=ans*base%k;
    	base=base%k*base%k;
    	n>>=1;
	}
	return ans%k;
}
int main(){
	cin>>a>>n>>k;
	cout<<a<<"^"<<n<<" "<<"mod"<<" "<<k<<"="<<helper(a,n)%k;
	return 0;	
}

斐波那契数列：
分析：对于这道题，第一次考虑到的方法自然是运用三个变量a，b，c滚动着来模拟斐波那契数列的过程（40分）：
代码：

#include<bits/stdc++.h>
using namespace std;
const int mod=1000000000+7;
long long a,b,c,n;
long long read(){
	long long x=0,f=1;
	char ch=getchar();
	while(ch<'0'||ch>'9'){
		if(ch=='-')
			f=-1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9'){
		x=(x<<1)+(x<<3)+(ch^48);
		ch=getchar();
	}
	return x*f;
}
int main(){
	n=read();
	a=1;
	b=1;
	for(long long i=3;i<=n;i++){
		c=a+b;
		c%=mod;
		a=b%mod;
		b=c%mod;
	}
	printf("%lld",c%mod);
	return 0;
}

思路2：运用矩阵乘法，模拟斐波那契数列的递推过程：
[ f(n-1) ，f(n)]* [?] = [f(n)，f(n+1)]
已知f(n-1)+f(n)=f(n+1)
且可以得出需要的是一个2*2的矩阵，所以可推得该矩阵为
0 1
1 1
代码:

#include<iostream>
#include<cstdio>
#include<cstring>
const int maxn=100+10;
#define int long long
const int b=1000000000+7;
using namespace std;
int k,a[maxn][maxn],s[maxn][maxn],tmp[maxn][maxn],n=2;
signed main(){
    cin>>k;
    if(k<=2){
    	cout<<1;
    	return 0;
	}
	k-=2;
    s[1][1]=a[1][1]=0;
    s[1][2]=a[1][2]=1;
    s[2][1]=a[2][1]=1;
    s[2][2]=a[2][2]=1;
    while(k){
        if(k%2){
            memset(tmp,0,sizeof(tmp));
            for(int i=1;i<=n;i++)
                for(int j=1;j<=n;j++)					
                    for(int k=1;k<=n;k++){
                        tmp[i][j]+=s[i][k]*a[k][j];					
                        tmp[i][j]%=b;
                    }
            for(int i=1;i<=n;i++)
                for(int j=1;j<=n;j++){
                    s[i][j]=tmp[i][j];
                }
        }
        memset(tmp,0,sizeof(tmp));
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)					
                for(int k=1;k<=n;k++){
                    tmp[i][j]+=a[i][k]*a[k][j];					
                    tmp[i][j]%=b;
                }
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++){
                a[i][j]=tmp[i][j]%b;
            }
        k/=2;	
    }
    cout<<s[2][2];
    return 0;
}