傳送門

洛谷1306傳送門

?

-----------------------------------------------------------------------------------------------------------------------------------------------------

內網題題面

?

題目描述

????煩神看上了一個妹子，煩神想找她約會，可是妹子出了一道數學題來考驗煩神，煩神只有做對了，妹子才會跟他去約會，題目是這樣的：

????給出兩個正整數A和B，要求求出斐波那契數列的第A項和第B項的最大公約數，結果對19990208取模（煩神約的妹子的生日）。

????煩神當然可以秒殺啦，但是他沒有時間做，于是他把這個問題留給了你們，他還要節省下時間去約下一個妹子，你能幫幫他嗎?

?

輸入

第一行一個正整數T，表示有T組數據,

接下來T行，每行兩個正整數A和B。

輸出

輸出有T行，第 i 行對應第 i 組數據的答案.

樣例輸入

2 1 1 36 48

樣例輸出

1 144

提示

?

0 < A, B <= 10^9.

1 <= T <= 100000.

?

-----------------------------------------------------------------------------------------------------------------------------------------------------

?

有一個很有意思的點

　　　　gcd(F[n],F[m])=F[gcd(n,m)]

?

于是再+矩陣快速冪就好了

?

#include<bits/stdc++.h> #define il inline #define ll long long #define mem(p) memset(&p,0,sizeof(p)) using namespace std; const ll mod=1e8; ll n,m; struct mat{ll a[3][3],r,c;}; il mat mul(mat x,mat y) {mat p;mem(p);for(int i=0;i<x.r;i++)for(int j=0;j<y.c;j++)for(int k=0;k<x.c;k++)p.a[i][j]=(p.a[i][j]+x.a[i][k]*y.a[k][j])%mod;p.r=x.r,p.c=y.c;return p; } il void fast(ll k) {mat p,ans;mem(p),mem(ans);p.r=p.c=2;p.a[0][0]=p.a[0][1]=p.a[1][0]=1;ans.r=1,ans.c=2;ans.a[0][0]=ans.a[0][1]=1;while(k){if(k&1)ans=mul(ans,p);p=mul(p,p);k>>=1;}cout<<ans.a[0][0]; } il ll gcd(ll a,ll b){return b?gcd(b,a%b):a;} int main() {ios::sync_with_stdio(0);cin>>n>>m;n=gcd(n,m);if(n<=2)cout<<1;else fast(n-2);return 0; }

?

#include<cstdio> #include<algorithm> #include<cstring> #define ll long long using namespace std; const ll mod = 19990208; ll u,f[100005]; struct mat {ll m[4][4],l,r; }; inline mat mul(mat x,mat y) {mat p;memset(&p,0,sizeof(p));for(int i = 1;i <= x.l;i++)for(int j = 1;j <= y.r;j++)for(int k = 1;k <= x.r;k++)p.m[i][j] = (p.m[i][j] + x.m[i][k] *y.m[k][j]%mod)% mod;p.l = x.l,p.r = y.r;return p; } void fast(ll k) {mat p,ans;memset(&p,0,sizeof(p));memset(&ans,0,sizeof(ans));p.m[1][1] = p.m[1][2] = p.m[2][1] = 1;p.l = 2;p.r = 2;ans.m[1][1] = ans.m[1][2] = 1;ans.l = 1,ans.r = 2;while(k){if(k & 1)ans = mul(ans,p);p = mul(p,p);k >>= 1;}printf("%lld\n",ans.m[1][1]); } int main() {int t;scanf("%d",&t);ll a,b,maxn = -10;f[1] = f[2] = 1;for(int i = 1;i <= t;i++){scanf("%lld%lld",&a,&b);u = __gcd(a,b);if(u <= 2)printf("1\n");elsefast(u - 2);}return 0; }

emm...

兩次都因為struct結構體后面忘加“；”而gg

好在

還是一次a了

?

轉載于:https://www.cnblogs.com/darlingroot/p/10650094.html

總結

以上是生活随笔為你收集整理的烦神的斐波那契洛谷-1306-斐波那契公约数的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。