标签:
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
#include<bits/stdc++.h> using namespace std; #define ll __int64 #define esp 0.00000000001 const int N=3e5+10,M=1e6+10,inf=1e9,mod=1e9+7; struct is { ll a[10][10]; }; is juzhenmul(is a,is b,ll hang ,ll lie) { int i,t,j; is ans; memset(ans.a,0,sizeof(ans.a)); for(i=1;i<=hang;i++) for(t=1;t<=lie;t++) for(j=1;j<=lie;j++) { ans.a[i][t]+=(a.a[i][j]*b.a[j][t]); ans.a[i][t]%=mod; } return ans; } is quickpow(is ans,is a,ll x) { while(x) { if(x&1) ans=juzhenmul(ans,a,2,2); a=juzhenmul(a,a,2,2); x>>=1; } return ans; } ll getans(ll x) { if(x==1) return 1; if(x==2) return 2; is ans,base; memset(ans.a,0,sizeof(ans.a)); ans.a[1][1]=1; ans.a[2][2]=1; base.a[1][1]=2; base.a[1][2]=0; base.a[2][1]=1; base.a[2][2]=2; ans=quickpow(ans,base,x-2); return (ans.a[1][1]*2+ans.a[2][1])%mod; } int main() { ll x,y,z,i,t; int T; scanf("%d",&T); while(T--) { scanf("%I64d%I64d",&x,&y); if(x>=y) printf("%I64d\n",getans(x-y+1)); else printf("0\n"); } return 0; }
标签:
原文地址:http://www.cnblogs.com/jhz033/p/5655877.html