2019獨(dú)角獸企業(yè)重金招聘Python工程師標(biāo)準(zhǔn)>>>

? 寫出來也被它的速度嚇了一跳，程序最坑的地方應(yīng)該就是防止溢出了。 #include <stdio.h> #include <stdlib.h> #define MOD (19999997) struct matrix {unsigned long long a;unsigned long long b; //{a,b}unsigned long long c; //{c,d}unsigned long long d; }; struct matrix mul_matrix(struct matrix a,struct matrix b) //return a*b {struct matrix temp;temp.a=(a.a*b.a+a.b*b.c)%MOD;temp.b=(a.c*b.a+a.d*b.b)%MOD;temp.c=(a.a*b.b+a.b*b.d)%MOD;temp.d=(a.c*b.b+a.d*b.d)%MOD;return temp; } struct matrix pow_matrix(struct matrix ori) {return mul_matrix(ori,ori); } int main(int argc, char*argv[]) {unsigned long max=atoi(argv[1]);struct matrix ms[32]; //pow(M,n)int i=1;ms[0].a=0;ms[0].b=1;ms[0].c=1;ms[0].d=1;for(i=1;i<32;i++) ms[i]=pow_matrix(ms[i-1]);struct matrix answer;answer.a=1;answer.b=0;answer.c=0;answer.d=1;// E matrixi=0;while(max){if(max&0x01) answer=mul_matrix(answer,ms[i]);i++;max>>=1;}printf("%lld\n",(answer.d)%MOD);return 0;

}

轉(zhuǎn)載于:https://my.oschina.net/000quanwei/blog/402844

總結(jié)

以上是生活随笔為你收集整理的编程之美-快速费波拉契数列的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。