矩阵快速幂---斐波那契数列（最大公约数）

题目描述

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

Update：加入了一组数据。

输入输出格式

输入格式：

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

注意：数据很大

输出格式：

Fn和Fm的最大公约数。

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

输入输出样例

输入样例#1： 复制

4 7

输出样例#1： 复制

1

说明

用递归&递推会超时

用通项公式也会超时

注意：：公式：gcd（f(n),f(m)）=f(gcd(n,m))------------很重要

#include<cstdio>
#include<iostream>
#include<string.h>
#include<algorithm>
using namespace std;
//定义矩阵类型的结构体
struct Maxtrix{
	long long m[2][2];
}res,ans;
//定义矩阵乘法
Maxtrix mul(Maxtrix a,Maxtrix b,int n){
	Maxtrix temp;
	memset(temp.m,0,sizeof(temp.m));
	for(long long i=0;i<n;i++){
		for(long long j=0;j<n;j++){
			for(long long k=0;k<n;k++){
				temp.m[i][j]=(temp.m[i][j]+(a.m[i][k]*b.m[k][j]))%100000000;
			}
		}
	}
	return temp;
}

int main(){
	long long x,y;
	scanf("%lld %lld",&x,&y);
	long long g=__gcd(x,y);
	if(g==1||g==2){
		printf("1\n");
		return 0;
	}
//初始化单位矩阵和关系矩阵
	ans.m[0][0]=ans.m[0][1]=ans.m[1][0]=1;
	ans.m[1][1]=0;
	res.m[0][0]=res.m[1][1]=1;
	res.m[0][1]=res.m[1][0]=1;
	long long ss=g-3;
	while(ss){
		if(ss&1){
			res=mul(res,ans,2);
		}
		ans=mul(ans,ans,2);
		ss=ss>>1;
	}
	long long aa=res.m[0][0]+res.m[0][1];
	printf("%lld\n",aa%100000000);
	return 0;
}