dp专辑 G - Fibonacci [ 快速幂]

第一反映这一题也属于DP？？？是数学题呀，呵呵~不管了，开始做了......

题意：

求 斐波那契数列 的第n项，输出F_n mod 10000

思路：

快速幂

 模板：

long long qpow(int a, int b)
{
         long long c, d;
         c = 1;    //存a^b
         d = a;    //存a的倍幂
         while (b > 0)
         {
                 if (b & 1)     //或 if (b % 2 == 1)
                        c *= d;
                 b = b >> 1;   //或  b = b / 2
                 d = d * d;
          }
         return c;
}

稍微修改下就可以了~

//AC CODE：

#include<iostream>#include<cmath>#include<algorithm>#include<vector>#include<cstdio>#include<cstdlib>#include<cstring>#include<string>using namespace std;struct Mat{    int p[2][2];};Mat mul(Mat c, Mat b){    Mat result;    result.p[0][0] = (c.p[0][0]*b.p[0][0] + c.p[0][1]*b.p[1][0]) % 10000;    result.p[0][1] = (c.p[0][0]*b.p[0][1] + c.p[0][1]*b.p[1][1]) % 10000;    result.p[1][0] = (c.p[1][0]*b.p[0][0] + c.p[1][1]*b.p[1][0]) % 10000;    result.p[1][1] = (c.p[1][0]*b.p[0][1] + c.p[1][1]*b.p[1][1]) % 10000;    return result;}int qpow(int N){    Mat  c, d;    c.p[0][0] = 1;    c.p[0][1]=0;    c.p[1][0]=0;    c.p[1][1]=1;    //存a^b    d.p[0][0] = 1;    d.p[0][1]=1;    d.p[1][0]=1;    d.p[1][1]=0;    //存a的倍幂    while (N > 0)    {        if (N & 1)     //或 if (b % 2 == 1)            c = mul(c,d);        N = N >> 1;   //或  b = b / 2        d = mul(d,d);    }    return c.p[0][1];}int main(){    int N;    while(scanf("%d",&N)!=EOF&& N!=-1)    {        printf("%d\n",qpow(N));    }    return 0;}