bzoj 2326: [HNOI2011]数学作业

 #include<cstdio>
 #include<iostream>
 #include<cstring>
 #include<cmath>
 #define ll long long
 ll b[],x[][],y[][],n,a[],ans,ans1,m;
 void cheng(ll q[][],ll q1[][])
 {
     ll q2[][];
     memset(q2,,sizeof(q2));
     for(int i=;i<;i++)
       for(int j=;j<;j++)
         for(int k=;k<;k++)
           q2[i][j]=(q2[i][j]+(q[i][k]*q1[k][j])%m)%m;
     for(int i=;i<;i++)
       for(int j=;j<;j++)
         q[i][j]=q2[i][j];
     return;
 }
 int main()
 {
     scanf("%lld%lld",&n,&m);
     b[]=;
     for(int i=;i<=;i++)
       {
         b[i]=(b[i-]*);
         if(n>=b[i])
           a[]=i;
       }
     a[]++;
     if(n<)
       a[]=n;
     else
       a[]=;
     for(int i=;i<a[];i++)
       a[i]=a[i-]*;
     if(n>=)
       a[a[]]=n-b[a[]-]+;
     for(int i=;i<=a[];i++)
       {
         memset(x,,sizeof(x));
         x[][]=x[][]=x[][]=x[][]=x[][]=;
         memset(y,,sizeof(y));
         x[][]=b[i]%m;
         y[][]=y[][]=y[][]=;
         ll k=a[i];
         for(;k;)
           {
             if(k%)
               cheng(y,x);
             cheng(x,x);
             k/=;
           }
         ans=(y[][]+(y[][]*ans1)%m+(y[][]*ans)%m)%m;
         ans1=(y[][]+(ans1*y[][])%m)%m;
       }
     printf("%lld\n",ans);
     return ;
 }

这是个矩阵乘法 对相同的位数构建一个矩阵。