1. 程式人生 > >輾轉相除的層數

輾轉相除的層數

nod long long iostream namespace pac fas printf 溢出 code

技術分享

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>
#include<cmath>
#define MOD 1000000007
using namespace std;
struct node{
    int  v[2][2]; //註意最好開 long log 
}ans,f,ans1;
long long k;
node ch(node a,node b)
{
    
    for(int i=0;i<=1
;i++) for(int j=0;j<=1;j++) { ans1.v[i][j]=0; for(int kk=0;kk<=1;kk++) ans1.v[i][j] = ( 1LL*ans1.v[i][j]+(1LL*a.v[i][kk]*b.v[kk][j])%MOD +MOD)%MOD; // 沒提前轉化成(long long )再 取模,溢出後再取模就錯了 } return ans1; } void fastlow() { while(k) { if
(k%2) ans=ch(ans,f); k/=2; f=ch(f,f); } } int main() { freopen("gcd.in","r",stdin); freopen("gcd.out","w",stdout); cin>>k; if(k==1) { printf("1 1"); return 0; } ans.v[0][0]=1;ans.v[0][1]=1;ans.v[1][0]=ans.v[1][1]=0; f.v[0][0]=1,f.v[0][1
]=1,f.v[1][0]=1;f.v[1][1]=0; fastlow(); printf("%d %d\n",(ans.v[0][1]+MOD)%MOD,(ans.v[0][0]+MOD)%MOD); return 0; }

輾轉相除的層數