做了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))=
刚开始就直接用了矩阵连乘的 算法 但对线性代数的一些知识不了解 错了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);
}
}