P1306 斐波那契公约数

题目描述

对于Fibonacci数列：1,1,2,3,5,8,13......大家应该很熟悉吧~~~但是现在有一个很“简单”问题：第n项和第m项的最大公约数是多少？

输入输出格式

输入格式：

两个正整数n和m。（n,m<=10^9）

注意：数据很大

输出格式：

Fn和Fm的最大公约数。

由于看了大数字就头晕，所以只要输出最后的8位数字就可以了。

输入输出样例

输入样例#1： 

4 7

输出样例#1： 

1

说明

用递归&递推会超时

用通项公式也会超时

斐波那契数列有一个很高级的推论：gcd(f(n),f(m))==f(gcd(n,m))

于是，矩阵快速幂+暴力gcd

附代码：

#include<iostream>
#include<algorithm>
#include<cstdio>
#define MOD 100000000
using namespace std;
int n,m,k;
struct node{
	long long a[3][3];
}ans,base,s;
inline int read(){
	int date=0,w=1;char c=0;
	while(c<'0'||c>'9'){if(c=='-')w=-1;c=getchar();}
	while(c>='0'&&c<='9'){date=date*10+c-'0';c=getchar();}
	return date*w;
}
int gcd(int x,int y){
	if(y==0)return x;
	return gcd(y,x%y);
}
node operator *(const node &x,const node &y){
	node ret;
	for(int i=1;i<=2;i++)
		for(int j=1;j<=2;j++){
			ret.a[i][j]=0;
			for(int k=1;k<=2;k++){
				ret.a[i][j]+=x.a[i][k]*y.a[k][j];
				ret.a[i][j]%=MOD;
			}
		}
	return ret;
}
void mexp(int k){
	for(int i=1;i<=2;i++)s.a[i][i]=1;
	base.a[1][2]=base.a[2][1]=base.a[2][2]=1;
	while(k){
		if(k&1)s=s*base;
		base=base*base;
		k>>=1;
	}
}
int main(){
	n=read();m=read();
	k=gcd(n,m);
	if(k==1||k==2){
		printf("1\n");
		return 0;
	}
	mexp(k-2);
	ans.a[1][1]=ans.a[1][2]=1;
	ans=ans*s;
	printf("%lld\n",ans.a[1][2]);
	return 0;
}