hdu1757

添加链接描述

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll k,m;
struct mm{
	ll a[11][11];
}aa,bb;
mm init(){
	mm E;
	for(int i=0;i<10;i++){
		for(int j=0;j<10;j++){
			if(i==j) E.a[i][i]=1;
			else E.a[i][j]=0;
		}
	}
	return E;
}
mm operator *(mm x,mm y){
	mm c;
	memset(c.a,0,sizeof(mm));
	for(int i=0;i<10;i++){
		for(int j=0;j<10;j++){
			for(int k=0;k<10;k++)
				c.a[i][j]+=x.a[i][k]*y.a[k][j]%m;
				c.a[i][j]%=m;
		}
	}
	return c;
}
ll fast(mm x,mm y,ll k){
	mm c=init();
	while(k){
		if(k%2) c=c*x;
		x=x*x;
		k/=2;
	}
	c=c*y;
	return c.a[0][0];
}
int main()
{
	for(int i=0;i<10;i++){
		bb.a[i][0]=9-i;
		for(int j=1;j<10;j++)
			bb.a[i][j]=0;
	}
	while(scanf("%lld%lld",&k,&m)!=EOF){
		for(int i=0;i<10;i++){
			scanf("%lld",&aa.a[0][i]);
			if(i) aa.a[i][i]=0;
			if(i<9) aa.a[i+1][i]=1;
		}
		ll sum=0;
		if(k<10){
			sum=k;
		}
		else{
			sum=fast(aa,bb,k-9);
		}
		printf("%d\n",sum%m);
	}
	return 0;
}