UVA10870遞推關係（矩陣乘法）

給以個遞推f(n) = a1 f(n - 1) + a2 f(n - 2) + a3 f(n - 3) + ... + ad f(n - d), for n > d.，給你n,d,a1,a2..ad ,f[1],f[2]..f[d],讓你求f[n]%m.

思路：比較基礎的矩陣題目，每次都構造乙個d*d的矩陣，然後用快速冪求出來它的n-1次冪，然後在求出乘積就行了，簡單構造，沒有什麼坑點。

#include

#include

typedef struct

mat;

long long n ,mod ,d;

mat mm(mat a ,mat b)

mat quick(mat a ,long long b)

return c;

}int main ()

for(i = 1 ;i <= d ;i ++) 

if(n <= d)

memset(a.mat ,0 ,sizeof(a.mat));

int x = 2 ,y = 1;

for(i = 2 ;i <= d ;i ++)

for(i = 1 ;i <= d ;i ++)

a.mat[i][d] = d[d-i+1];

a = quick(a ,n - 1);

long long ans = 0;

for(i = 1 ;i <= d ;i ++)

printf("%lld\n" ,ans);

}return 0;