比赛 20120217 评测结果 AAAAAAAAAA
题目名称 线性递推式 最终得分 100
用户昵称 Czb。 运行时间 0.000 s
代码语言 C++ 内存使用 0.00 MiB
提交时间 2012-02-17 22:16:26
显示代码纯文本
#include<stdio.h>

long long n,m,ans,a[12][12],tmp[12][12],h[100],t[100],f[12],g[12],top;

int main()
{
	freopen("recursion.in","r",stdin);
	freopen("recursion.out","w",stdout);
	long long i,j,k,x;
	scanf("%lld%lld",&n,&m);
	m=m-n+1;
	for(i=1;i<=n;i++)
	{
		a[i+1][i]=1;
	}
	for(i=1;i<=n;i++)
	{
		scanf("%lld",&a[i][n]);
	}
	a[n+1][n+1]=1;
	scanf("%lld",&ans);
	f[n+1]=ans;
	for(i=1;i<=n;i++)
	{
		scanf("%lld",&f[i]);
	}
	h[0]=1;
	for(i=1;i<=60;i++)
	{
		h[i]=h[i-1]*2;
	}
	for(i=60;i>=0;i--)
	{
		if(m>=h[i])
		{
			top++;
			t[top]=i;
			m-=h[i];
		}
	}
	for(x=0;x<=60;x++)
	{
		if(t[top]==x)
		{
			for(i=1;i<=n+1;i++)
			{
				g[i]=f[i];
				f[i]=0;
			}
			for(i=1;i<=n+1;i++)
			{
				for(j=1;j<=n+1;j++)
				{
					f[i]=(f[i]+g[j]*a[j][i])%9973;
				}
			}
			top--;
		}
		for(i=1;i<=n+1;i++)
		{
			for(j=1;j<=n+1;j++)
			{
				tmp[i][j]=a[i][j];
				a[i][j]=0;
			}
		}
		for(i=1;i<=n+1;i++)
		{
			for(j=1;j<=n+1;j++)
			{
				for(k=1;k<=n+1;k++)
				{
					a[i][j]=(a[i][j]+tmp[i][k]*tmp[k][j])%9973;
				}
			}
		}
	}
	printf("%lld\n",f[n]%9973);
	return 0;
}