hdu A Simple Math Problem

做了2天终于AC了，真是可喜可贺
刚开始就直接用了矩阵连乘的 算法 但对线性代数的一些知识不了解 错了N次
方法：(f(n),f(n-1),f(n-2),f(n-3),f(n-4),f(n-5),f(n-6),f(n-7),f(n-8),f(n-9))=

(f(9),f(8),f(7),f(6),f(5),f(4),f(3),f(2),f(1),f(0))*

^(n-9);

f(n)=f(9)*a0+f(8)*a1+f(7)*a2+f(6)*a3+f(5)*a4+f(4)*a5+f(3)*a6+f(2)*a7+f(1)*a8+f(0)*a9;(a0到a9都是变化后的值)
代码:

#include<stdio.h>
#include<math.h>
#include<string.h>
#include<stdlib.h>
int main()
{
long long ju[11][11],zh[11][11],chu[11][11],f[10]={9,8,7,6,5,4,3,2,1,0},zhan[1000],len;
   long long n,a,i,j,x,y,z,m,k,ji,result,sum,sign;
   while(scanf("%I64d%I64d",&n,&m)==2)
   {
       sum=0;
    memset(ju,0,sizeof(ju));
    memset(chu,0,sizeof(chu));
    for(i=0;i<=9;i++)
    {
    scanf("%I64d",&ju[i][0]);
    chu[i][0]=ju[i][0];
    }
    if(n==10)
    {
        for(i=0;i<=9;i++)
        sum+=ju[i][0]*f[i];
        printf("%I64d\n",sum%m);
        continue;
    }
    if(n<9)
    {
        printf("%I64d\n",n%m);
    }
    n-=9;
    for(i=0,j=1;j<=9;i++,j++)
    ju[i][j]=chu[i][j]=1;
    memset(zhan,0,sizeof(zhan));
     for(len=0;n>1;)
    {
        if(n%2==1)
        zhan[len++]=1;
        else
        zhan[len++]=2;
        n/=2;
    }
    for(len--;len>=0;len--)
    {
        memset(zh,0,sizeof(zh));
    for(j=0;j<=9;j++)
        {
            for(x=0;x<=9;x++)
            {
                sign=0;
            for(y=0;y<=9;y++)
                {
                    sign+=ju[j][y]*ju[y][x];
                }
                zh[j][x]=sign%m;
            }
        }
        for(k=0;k<=9;k++)
        for(z=0;z<=9;z++)
        ju[k][z]=zh[k][z];
        memset(zh,0,sizeof(zh));
    if(zhan[len]==1)
    {
        for(x=0;x<=9;x++)
        {
        for(i=0;i<=9;i++)
            {
                sign=0;
                for(j=0;j<=9;j++)
                {
                    sign+=ju[x][j]*chu[j][i];
                }
                zh[x][i]=sign%m;
            }
        }
    }
    if(zhan[len]==1)
    for(k=0;k<=9;k++)
        for(z=0;z<=9;z++)
        ju[k][z]=zh[k][z];
    }
        sign=0;
        for(j=0;j<=9;j++)
        {
            sign+=f[j]*ju[j][0];
        }
        sum+=sign;
    printf("%I64d\n",sum%m);
   }
}