矩阵模板之一
一個大數的k次方? 斐波那契數
#include <stdio.h> #include <string.h> #include <algorithm> using namespace std; struct app {int p[30][30]; }; app mult(app a, app b) {app c;int i , j, k;for( i =0; i<2; i++)for(j =0; j<2; j++){c.p[i][j]=0; //注意清零的操作for(k =0; k<2; k++){c.p[i][j]=(c.p[i][j]+a.p[i][k]*b.p[k][j])%10000; //取每個元素的后四位}}return c; }app pow(app a, int k) {app b; //用來返回值。memset(b.p, 0, sizeof(b.p)); //注意清空的類型是數組時要確定到數組上。for(int i=0;i<2; i++){b.p[i][i]=1;}while(k){if(k&1) //所在為的二進制是一的時候b=mult(b,a);a=mult(a,a); //二進制是零的情況也要繼續相乘k>>=1; //向右一位。}return b; } int main() {int n,t;scanf("%d", &t);while(t--){ scanf("%d", &n);app a;a.p={{1,1},{1,0}}; //f1 的值app d=pow(a, n);printf("%d\n",d.p[0][1]); //d.p[0][0]是f[n+1]的值}return 0; }
?
總結
- 上一篇: NOJI 148 fibonacci数
- 下一篇: 循环节模板 NOJ427Number