題目描述

對於Fibonacci數列：1,1,2,3,5,8,13......大家應該很熟悉吧~~~但是現在有一個很“簡單”問題：第n項和第m項的最大公約數是多少？

Update：加入了一組數據。

輸入輸出格式

輸入格式：

兩個正整數n和m。（n,m<=10^9）

註意：數據很大

輸出格式：

Fn和Fm的最大公約數。

由於看了大數字就頭暈，所以只要輸出最後的8位數字就可以了。

輸入輸出樣例

4 7

1

說明

用遞歸&遞推會超時

用通項公式也會超時

code

#include<stdio.h>
#include<string
 
.h>
#include<algorithm>
#define ll long long 
using namespace std;
const int P=1e8;
struct node {
    ll mp[3][3];
    int h,l;
}p,ans;

int gcd(int x,int y)
{
    if(x<y) swap(x,y);
    if(x%y==0) return y;
    else return gcd(y,x%y);
}

node mul(node x,node y) {
    node tep;
    memset(
 
&tep,0,sizeof(tep));
    for(int i=0;i<x.h;++i) 
        for(int j=0;j<y.l;++j) 
            for(int k=0;k<x.l;++k) 
                tep.mp[i][j]=(tep.mp[i][j]+x.mp[i][k]*y.mp[k][j])%P;
    tep.h=x.h,tep.l=y.l;
    return tep;
}

int juc(ll k)
{
    ans.mp[0][0]=ans.mp[0][1]=1;
    ans.h=1,ans.l=2
 
;
    p.mp[0][0]=p.mp[0][1]=p.mp[1][0]=1;
    p.h=p.l=2;
    while(k) {
        if(k&1) ans=mul(ans,p);
        p=mul(p,p);
        k>>=1;
    }
    return ans.mp[0][0];
}

int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    int d=gcd(n,m);
    if(d<=2) printf("1");
    else printf("%d",juc(d-2));
    return 0;
}

斐波那契公約數（luogu 1306）