2021-05-06

4 篇文章 0 订阅
2 篇文章 0 订阅
#include<iostream>
#include<algorithm>
#include<vector>
using namespace std; 
int t,n,p;
//矩阵乘法
vector<vector<int> > f1(vector<vector<int> >a, vector<vector<int> > b){
	int n=a.size(), k=a[0].size(), m=b[0].size();
	vector<vector<int> > c(n,vector<int>(m,0));
	for(int i=0;i<n;++i){
		for(int j=0;j<m;++j){
			for(int x=0;x<k;++x){
				c[i][j]=(c[i][j]+a[i][x]*b[x][j])%p;
			}
		}
	}
	return c;
}
//矩阵快速幂
vector<vector<int> > f2(vector<vector<int> >a,int x){
	vector<vector<int> >b(2);
	b[0].push_back(1);
	b[0].push_back(0);
	b[1].push_back(0);
	b[1].push_back(1); 
	if(x<0) return b;
	while(x){
		if(x&1) b=f1(b,a);
		a=f1(a,a);
		x>>=1;
	}
	return b;
}
//斐波那契数列
int main(){
	cin>>t;
	vector<vector<int> > b(2);
	b[0].push_back(1);
	b[0].push_back(1);
	b[1].push_back(1);
	b[1].push_back(0);
	vector<vector<int> > e(2);
 	e[0].push_back(1);
	e[1].push_back(1);
	
	while(t--){
		cin>>n>>p;
		vector<vector<int> > c=f2(b,n-1);
		vector<vector<int> > d=f1(c,e);
		cout<<d[0][0]%p<<endl;	
	}
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值