矩阵快速幂

在这里插入图片描述
在这里插入图片描述
解题思路:原来是这么构造的,之前想了一个晚上都没有想明白,,哎。。。数组不开long long 也会wa哦

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<map>
#include<vector>
#include<queue>
#include<set>
#define IL inline
#define x first
#define y second
typedef long long ll;
using namespace std;
const	int N=10;
ll f[N][N];
ll a[N];
ll g[N][N]; 
ll k,m;
void mul(ll a[][N],ll b[][N],ll c[][N])
{
	ll temp[N][N];
	memset(temp,0,sizeof temp);
	for(int i=0;i<N;i++)
	for(int j=0;j<N;j++)
	for(int k=0;k<N;k++)
	temp[i][j]=(temp[i][j] + b[i][k]*c[k][j])%m;
	memcpy(a,temp,sizeof temp);
}
int main()
{
	while(cin >> k >> m)
	{
		for(int i=0; i <10;i++)
		cin >>a[i];
		memset(g,0,sizeof g);
		memset(f,0,sizeof f);
		for(int i=0;i<10;i++)
		g[0][i] = a[i];
		for(int i=1;i<10;i++)
		g[i][i-1] = 1 ;
		for(int i = 0;i<10;i++)
		f[i][0] =  9 - i;
		if(k<=9)
		cout<<k%m<<endl;
		else
		{
			k -= 9;
			while(k)
			{
				if(k&1)	mul(f , g ,f);
				k>>=1;
				mul(g,g,g);
			 }
			 cout << f[0][0] << endl; 
		}
	}
    return 0;
}




  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值