P1306 斐波那契公约数详细题解

题意

题目链接
求斐波拉契第n项与第m项的gcd的后8位。

题解

首先，对于这题有一个结论，这还是一次华东师范oj上遇到的，就是gcd(F(n) , F(m)) = F(gcd(n,m))，当时只知道这个结论，不知如何来的>_<。
证明
假设 F(n) = a , F(n+1) = b.
F(n+2) = a + b , F(n+3) = a + 2b , F(n+4) = 2a + 3*b…
F(m) = F(m-n-1)*a + F(m-n)*b ==> F(m) = F(m-n-1)*F(n) + F(m-n)*F(n+1)
所以求gcd( F(n) , F(m) ) = gcd( F(n) , F(m-n-1)*F(n) + F(m-n)*F(n+1) )
由于F(n) | F(m-n-1)*F(n) 和 更相减损法
① gcd( F(n) , F(m-n-1)*F(n) + F(m-n)*F(n+1) ) = gcd( F(n) , F(m-n)*F(n+1) )
对于斐波拉契相邻2项一定互质，即：gcd( F(n) , F(n + 1) ) = 1
:: ( F(n) ,F(n-1) + F(n) ) = ( F(n) ,F(n-1) ) = (F(n-2) ,F(n-1) )… = (F(1) ,F(2) ) = 1

所以对于① = gcd( F(n) , F(m-n) )
gcd( F(n) ,F(m%n)) = gcd( F(n%(m%n)) , F(m%n) ) … = gcd(F(n1) , F(0) ) = F(n1)
最终就是再求n与m的gcd 的 F();
所以结论得证！

然后接下来就只需求F(gcd),用矩阵快速幂解。

AC代码(14ms)

#pragma GCC optimize(2)
#include<bits/stdc++.h>
using namespace std;
const int man = 2e5+10;
template<typename T>T gcd(T a, T b) {return b==0?a:gcd(b, a%b);}
template<typename T>T exgcd(T a,T b,T &g,T &x,T &y){if(!b){g = a,x = 1,y = 0;}else {exgcd(b,a%b,g,y,x);y -= x*(a/b);}}
#ifndef ONLINE_JUDGE
#define debug(fmt, ...) {printf("debug ");printf(fmt,##__VA_ARGS__);puts("");}
#else
#define debug(fmt, ...)
#endif
typedef long long ll;
const int mod = 1e8;
struct matrx{
    ll a[2][2];
    friend matrx operator *(matrx A,matrx B){
        matrx ans;
        memset(ans.a,0,sizeof(ans.a));
        for(int i = 0;i < 2;i++){
            for(int j = 0;j < 2;j++){
                //ans.a[i][j] = 0;
                for(int k = 0;k < 2;k++){
                    ans.a[i][j] = (ans.a[i][j] + A.a[i][k] * B.a[k][j])%mod;
                }
            }
        }
        return ans; 
    }
};

int quick(matrx A,int b){
    matrx ans;
    ans.a[0][1] = ans.a[1][0] = 0;
    ans.a[0][0] = ans.a[1][1] = 1;
    while(b){
        if(b&1)ans = ans * A;
        A = A*A;
        b /= 2;
    }
    return (ans.a[0][0] + ans.a[0][1])%mod;
}
  
int main() {
    #ifndef ONLINE_JUDGE
        //freopen("in.txt", "r", stdin);
        //freopen("out.txt","w",stdout);
    #endif
    int n,m;
    scanf("%d%d",&n,&m);
    int d = gcd(n,m);
    //cout << d << endl;
    if(d<=2)cout << 1 << endl;
    else{
        matrx ans;
        ans.a[1][1] = 0;
        ans.a[0][0] = ans.a[1][0] = ans.a[0][1] = 1;
        int res = quick(ans,d-2);
        cout << res << endl;
    }
    return 0;
}