poj3070_斐波那契数列(Fibonacci)

用矩阵求斐波那契数列，快速幂log(n),只用求最后4位(加和乘的运算中前面的位数无用)

 #include <stdio.h>

 #include <stdlib.h>

 int main()

 {

     /*

     (x y)(x y)=  (x*x+y*s x*y+y*t) =(x*x+y*s  y*(x+t))

     (s t)(s t)   (s*x+t*s s*y+t*t)  (s*(x+t)  t*t+y*s)

     (x y)(a b)= (x*a+y*c x*b+y*d)

     (s t)(c d)  (s*a+t*c s*b+t*d)

     (a b)(f1=1)=(a+b)

     (c d)(f2=1) (c+d)

     */

     //(a,b)(c,d)一开始为E2

     long n,x,y,s,t,xx,yy,ss,tt,a,b,c,d,aa,bb,cc,dd;

     while (scanf("%ld",&n))

     {

         if (n==)

         {

             printf("0\n");

             continue;

         }

         else if (n==-)

             break;

         n-=;

         x=;

         y=;

         s=;

         t=;

         a=;

         b=;

         c=;

         d=;

         while (n)

         {

             if (n%==)

             {

                 aa=a;

                 bb=b;

                 cc=c;

                 dd=d;

                 a=(x*aa+y*cc)%;

                 b=(x*bb+y*dd)%;

                 c=(s*aa+t*cc)%;

                 d=(s*bb+t*dd)%;

             }

             xx=x;

             yy=y;

             ss=s;

             tt=t;

             x=(xx*xx+yy*ss)%;

             y=(yy*(xx+tt))%;

             s=(ss*(xx+tt))%;

             t=(tt*tt+yy*ss)%;

             n=n/;

         }

         printf("%ld\n",(c+d)%);

     }

     return ;

 }

 #include <stdio.h>

 #include <stdlib.h>

 int main()

 {

     /*

     (x y)(x y)=  (x*x+y*s x*y+y*t) =(x*x+y*s  y*(x+t))

     (s t)(s t)   (s*x+t*s s*y+t*t)  (s*(x+t)  t*t+y*s)

     (x y)(a b)= (x*a+y*c x*b+y*d)

     (s t)(c d)  (s*a+t*c s*b+t*d)

     (a b)(f1=1)=(a+b)

     (c d)(f2=1) (c+d)

     */

     //(a,b)(c,d)一开始为E2

     long n,x[],y[],s[],t[],a[],b[],c[],d[],w,r,i;

     x[]=;

     y[]=;

     s[]=;

     t[]=;

     a[]=;

     b[]=;

     c[]=;

     d[]=;

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

     {

         x[i]=(x[i-]*x[i-]+y[i-]*s[i-])%;

         y[i]=(y[i-]*(x[i-]+t[i-]))%;

         s[i]=(s[i-]*(x[i-]+t[i-]))%;

         t[i]=(t[i-]*t[i-]+y[i-]*s[i-])%;

     }

     while (scanf("%ld",&n))

     {

         if (n==)

         {

             printf("0\n");

             continue;

         }

         else if (n==-)

             break;

         n-=;

         w=;

         r=;

         while (n)

         {

             if (n%==)

             {

                 a[r+]=(x[w]*a[r]+y[w]*c[r])%;

                 b[r+]=(x[w]*b[r]+y[w]*d[r])%;

                 c[r+]=(s[w]*a[r]+t[w]*c[r])%;

                 d[r+]=(s[w]*b[r]+t[w]*d[r])%;

                 r++;

             }

             w++;

             n=n/;

         }

         printf("%ld\n",(c[r]+d[r])%);

     }

     return ;

 }