gcd

样例坑爹系列
本来输出的是对的
一看样例就倒过来了
(╯‵□′)╯︵┻━┻

打表找规律
(╯‵□′)╯︵┻━┻

没有提示
去研究斐波那契数列和这个题的关系
(╯‵□′)╯︵┻━┻

就是说
F_k+1和F_k+2
辗转相除的次数是k

但是对于k=1时
情况不一样
此时a=1,b=1

斐波那契数列的矩阵快速幂算法

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
#include<queue>
#include<vector>
#include<climits>
#include<string>
#include<cstdlib>
#include<ctime>
#define MOD 1000000007
#define LL long long
using namespace std;

struct node
{
    LL arr[3][3];
}a,b;

LL k;

node ch(node p1,node p2)
{
    node l;

    for(int i=1;i<=2;i++)
        for(int j=1;j<=2;j++)
        {
            l.arr[i][j]=0;

            for(int t=1;t<=2;t++)
                l.arr[i][j]=(l.arr[i][j]+(p1.arr[i][t]*p2.arr[t][j])%MOD)%MOD;          
        }

    return l;
}

void fastpow(LL p)
{
    while(p>1)
    {
        if(p%2) a=ch(a,b);

        p=p/2;
        b=ch(b,b);
    }

    a=ch(a,b);
}


int main()
{
    freopen("gcd.in","r",stdin);
    freopen("gcd.out","w",stdout);

    scanf("%lld",&k);

    if(k==1)
    {
        printf("1 1");
        return 0;
    }

    a.arr[1][1]=1;
    a.arr[1][2]=1;

    b.arr[1][1]=1;
    b.arr[1][2]=1;
    b.arr[2][1]=1;
    b.arr[2][2]=0;

    fastpow(k);

    printf("%lld %lld",a.arr[1][1]%MOD,a.arr[1][2]%MOD);

    return 0;
}