原题链接
题目的意思还是非常明了的,但是如果从l开始一个个地加到r势必会超时,所以就需要找到sn=sum( T(k) )0<=k<=n的方法,也就是找出前n项和的方法。发现这里前n项和有这样的规律那就是Sn=Sn-1 + Sn-2 +Sn-3 S0=1;S1=2;S2=3;
那么当矩阵有n-2次那么多的时候,得到的c矩阵的第一行的第一列的元素就是我们要求的那一项了
//http://acm.xidian.edu.cn/problem.php?id=1145
#include <algorithm>
#include <iostream>
#include <utility>
#include <sstream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <queue>
#include <stack>
#include <cmath>
#include <map>
#include <set>
using namespace std;
typedef long long ll;
const int MOD = int(1e9) + 7;
//int MOD = 99990001;
const int INF = 0x3f3f3f3f;
const ll INFF = 0x3f3f3f3f3f3f3f3fLL;
const double EPS = 1e-9;
const double OO = 1e20;
const double PI = acos(-1.0); //M_PI;
const int fx[] = {-1, 1, 0, 0};
const int fy[] = {0, 0, -1, 1};
const int maxn=3;
struct matrix{
ll mat[maxn][maxn];
};
ll fastmultiply(ll x,ll y){
ll res=0;
while(y){
if(y&1) res=(res+x)%MOD;
x=(x*2)%MOD;
y>>=1;
}
return res;
}
matrix matrix_multiply(matrix a,matrix b,int n){
matrix c;
memset(c.mat,0,sizeof(c.mat));
for(int i=0;i<n;i++)
for(int k=0;k<n;k++){
if(a.mat[i][k]!=0){
for(int j=0;j<n;j++){
if(b.mat[k][j]==0) continue;
c.mat[i][j]=(c.mat[i][j] + fastmultiply(a.mat[i][k],b.mat[k][j]))%MOD;
}
}
}
return c;
}
matrix matrix_index(matrix a,ll mi,int n){
matrix t;
memset(t.mat,0,sizeof(t.mat));
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
t.mat[i][j]=(i==j);
while(mi){
if(mi&1) t=matrix_multiply(t,a,n);
a=matrix_multiply(a,a,n);
mi>>=1;
}
return t;
}
//F(x)是T(x)的前n项和
ll F(ll x){
if(x==0 || x==1 || x==2) return x+1;
matrix a,t,tt,c;
memset(a.mat,0,sizeof(a.mat));memset(t.mat,0,sizeof(t.mat));
memset(tt.mat,0,sizeof(tt.mat));memset(c.mat,0,sizeof(c.mat));
a.mat[0][0]=3;a.mat[0][1]=2;a.mat[0][2]=1;
t.mat[0][0]=1;t.mat[1][0]=1;t.mat[2][0]=1;t.mat[0][1]=1;t.mat[1][2]=1;
tt=matrix_index(t,x-2,3);
c=matrix_multiply(a,tt,3);
return c.mat[0][0]%MOD;
}
int main(){
//这个题一定要找出求和公式中前n项和的公式,否则就会超时,其通项也相当于是一个泰波拉契数列
ll l,r;
while(scanf("%lld%lld",&l,&r)==2){
if(l>=1)
printf("%lld\n",(F(r)-F(l-1)+MOD)%MOD);//这里不能用减法,因为模后前者可能会比后者小,所以需要加一个mod以实现将模后的负数转化成正数的效果
else
printf("%lld\n",F(r)%MOD);
}
return 0;
}