程序员充电站（itcharge）

题目大意： 给定一个N*N的矩阵A和一个整数K，要求计算S = A + A^2 + A^3 + … + A^k。 思路： 分别用矩阵快速幂求出每一项的A^i，然后将每一项矩阵相加，考虑到k值很大，所有采用 二分求解。

const int MAXN = 110; struct Matrax { int m[MAXN][MAXN]; }a,per; int N,M; void Init() { for(int i = 0; i < N; ++i) for(int j = 0; j < N; ++j) { scanf("%d",&a.m[i]

题目大意： 定义Sn = [( a + √b )^n] % m，[x]表示x向上取整，比如[3.14] = 4。给你a b n m的值， 求Sn是多少。 思路： 这道题很经典，因为(a-1)^2< b < a^2，所以0 < |a-√(b)| < 1，所以 Sn = [( a + √b )^n] % m = ( [( a + √b )^n] + [( a - √b))^n] ) % m。 即右边其实是一个整数，如果将右边二项式展开的话，除了相互抵消的部分，剩下的部分全为 整数。这个式子设An = (a