P1306 斐波那契公约数（矩阵快速幂+公式推导）

斐波那契数列有一个公式：gcd(F(n) , F(m)) = F(gcd(n,m))，正好解决这道例题，再用上我们的矩阵快速幂递推斐波那契数列就可以得出答案。

 

 

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int mod=100000000 ;
int n,m,k;
int a[4][4];
int b[4][4];
int b1[4][4];
void jzc(int c[4][4],int d[4][4])
{
	//c*d
	memset(b1,0,sizeof(b1));
	for(int i=0;i<2;i++)
	{
		for(int j=0;j<2;j++)
		{
			for(int k=0;k<2;k++)
			{
				b1[i][j]=(b1[i][j]+c[i][k]*d[k][j])%mod;	
			}
		}
	}
}
int gcd(int a,int b)
{
	while(b>0)
	{
		int t;
		t=a%b;
		a=b;
		b=t;
	}
	return a;
}
signed main()
{
	cin>>n>>m;
	a[0][0]=1;
	a[0][1]=1;
	a[1][0]=1;
	a[1][1]=0;
	b[0][0]=1;
	b[1][1]=1;
	k=gcd(n,m);
	//cout<<k<<endl;
	if(k<=2)
	{
		cout<<1<<endl;
		return 0;
	}
	k-=2;
	//快速幂 
	while(k)
	{
		if(k&1)
		{
			jzc(b,a);
			for(int i=0;i<2;i++)
			{
				for(int j=0;j<2;j++)
				{
					b[i][j]=b1[i][j];
				}
			}
		}
		jzc(a,a);
		for(int i=0;i<2;i++)
		{
			for(int j=0;j<2;j++)
			{
				a[i][j]=b1[i][j];
			}
		}
		k>>=1;
	}
	cout<<(b[0][0]+b[1][0])%mod<<endl;
}