【矩阵】

9 篇文章 0 订阅

初学矩阵,建议看一下这个bloghttps://www.cnblogs.com/cmmdc/p/6936196.html讲得比较清楚简单。

codeforces450B

思路:很简单的水题,初学有几个注意点,矩阵幂一定要是方阵,填0凑成方阵;ans是单位矩阵,res要自己算,我根据上面博客最后求斐波那契数列的方式写的,一定要前后承接,多乘几个就能发现bug;负数的模运算,+mod凑成整数再模。

#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <stdio.h>
#include <iostream>
#include <cstdlib>
#include <cmath>
#include <cctype>
#include <string>
#include <cstring>
#include <algorithm>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <ctime>
#include <vector>
#include <fstream>
#include <list>
#include <iomanip>
#include <numeric>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
#define ms(s) memset(s, 0, sizeof(s))
const int inf = 0x3f3f3f3f;
const int mod =1000000007;
#define LOCAL
int x,y,n;
struct Maxtrix
{
	int a[2][2];
}ans,res;
Maxtrix Mul(Maxtrix A,Maxtrix B)
{
	Maxtrix temp;
	for(int i=0;i<2;++i)
		for(int j=0;j<2;++j)
			temp.a[i][j]=0;
	for(int i=0;i<2;++i)
		for(int j=0;j<2;++j)
			for(int k=0;k<2;++k)
				temp.a[i][j]+=(A.a[i][k]%mod*B.a[k][j]%mod)%mod;
	return temp;
}
void quick_power(int x,int y,int n)
{
	ans.a[0][0]=ans.a[1][1]=1;
	ans.a[0][1]=ans.a[1][0]=0;
	res.a[0][0]=res.a[1][0]=1;
	res.a[0][1]=-1;res.a[1][1]=0;
	while(n)
	{
		if(n&1)ans=Mul(res,ans);
		res=Mul(res,res);
/*		for(int i=0;i<2;++i)
			for(int j=0;j<2;++j)
				cout<<ans.a[i][j]<<" ";*/
		//cout<<endl;
		n=n>>1;
	}
}
int main() 
{
   // freopen("in.txt", "r", stdin);
    while(cin>>x>>y>>n)
    {
    	if(n==1){cout<<(x%mod+mod)%mod<<endl;}
    	else if(n==2){cout<<(y%mod+mod)%mod<<endl;}
    	else 
    	{
    		quick_power(x,y,n-2);
    		cout<<((y*ans.a[0][0]%mod+x*ans.a[0][1]%mod)%mod+mod)%mod<<endl;

    	}
    }
    return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值