【Ybt OJ】[数学基础 第1章] 矩阵快速幂

$「$数学基础$」$第$1$章 矩阵快速幂
目录：

A.序列的第k个数
B.斐波那契数列
C.行为方案
D.矩阵求和
E.最短路径

$A.$ 例题$1$ 序列的第$k$个数

分析：

等差 等比数列通项就行了

CODE：

#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<iostream>
using namespace std;
typedef long long ll;
const int Mod=200907;
int T;
ll ans;
ll ksm(ll x,ll y)
{
	ll res=1;
	while(y)
	{
		if(y&1) res=(res*x)%Mod;
		y>>=1;
		x=x*x%Mod;
	}
	return res;
}
void query(ll a,ll b,ll c,ll k)
{
	ans=0;
	if(c-b==b-a)
	{
		ll sub=(b-a)%Mod;
		ans=(a+(k-1)*sub)%Mod;
	}
	else
	{
		ll div=(b/a)%Mod;
		ans=a*ksm(div,k-1)%Mod;
	}
}
int main()
{
	scanf("%d",&T);
	while(T--)
	{
		ll a,b,c,k;
		scanf("%lld%lld%lld%lld",&a,&b,&c,&k);
		query(a,b,c,k);
		printf("%lld\n",ans);
	}
	
	return 0;
}

$B.$ 例题$2$ 斐波那契数列

洛谷$link$

分析：

矩阵快速幂模板题了 $blog$_$link$

$C.$ 例题$3$ 行为方案

洛谷$link$

分析：

邻接矩阵$A^k$的$i$行$j$列表示 $i$到$j$经过$k$步的方案总数

将 邻接矩阵建出 矩阵快速幂出$t$次方 最后统计${\sum }_{i=1}^n{A}_{1,i}$
原地停留 可以认为存在自环 就不用管它了

自爆 可以将自爆看成一个城市 编号为$0$ 将所有城市与它连边 但它不往外连 就满足无法恢复状态了
最后答案$={\sum }_{i=0}^n{A}_{1,i}$

CODE：

#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<iostream>
using namespace std;
typedef long long ll;
const int Mod=2017,N=105;
int n,m,t,ans;
struct Matrix{
	int a[N][N];
}G;
Matrix operator *(Matrix A,Matrix B)
{
	Matrix res;
	memset(&res,0,sizeof(res));
	for(int i=0;i<=n;i++)
		for(int j=0;j<=n;j++)
			for(int k=0;k<=n;k++)
				res.a[i][j]=(res.a[i][j]+A.a[i][k]*B.a[k][j])%Mod;
	return res;
}
Matrix ksm(Matrix A,int k)
{
	Matrix res;
	memset(&res,0,sizeof(res));
	for(int i=0;i<=n;i++)
		res.a[i][i]=1;
	while(k)
	{
		if(k&1) res=res*A;
		k>>=1;
		A=A*A;
	}
	return res;
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++)
	{
		int u,v;
		scanf("%d%d",&u,&v);
		G.a[u][v]=G.a[v][u]=1;
	}
	for(int i=0;i<=n;i++)
		G.a[i][i]=G.a[i][0]=1;
	scanf("%d",&t);
	Matrix Ans=ksm(G,t);
	for(int i=0;i<=n;i++)
		ans+=Ans.a[1][i];
	printf("%d",ans%Mod);	
	
	return 0;
}

$D.$ 例题$4$ 矩阵求和

$Poj$_$link$

分析：

矩阵快速幂经典题了
不过$ybt$上面数据不知道有啥问题 只有$90pts$ 建议去$poj$做
$blog$_$link$

$E.$ 例题$5$ 最短路径

洛谷$link$

分析：

将矩阵乘法变成矩阵加法 然后做最短路
$blog$_$link$