斐波那契数列有一个公式: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;
}