很不错的一道题,除了矩阵快乘以外,还有线性代数和代数结构的相关知识。
/************************************************
**文件名:
**Copyright (c) 2010-2020 OrdinaryCrazy
**创建人:OrdinaryCrazy
**日期:20170804
**描述:百炼-2786参考答案
**版本:1.0
*************************************************/
#include <stdio.h>
#include <stdlib.h>
#define divisor 32767
/************************************************
和斐波那契数列一样,我们只需要将递推矩阵计算k-2次方即可
好吧,我们遇到了一个很尴尬的问题,就是在这个矩阵中a11的增长比
2的k-2次方还快一点,如果用置的数据类型,对于k>64的情况就没办法计算了
所以我们还要考虑这样一个问题:
a(n)=x*2+y
x=m*32767+a
y=n*32767+b
a(n)%32767=(2*a+b)%32767 所以为了防止矩阵越界,每次计算矩阵都只保留模32767的值
*************************************************/
int n,k;
struct matrix
{