需要用dp的思想去找状态,我们先维护出从一端开始的情况,还是可以分成两种情况,一种是回到原来位置,另一种是不是,
很明显要想最后回到原来位置,每次都要往前步进,计数为2的长度幂次,下面就是不回到原点了,一种是直接下一步把当前位置再走掉然后状态这里就开始切分,得到了一个F(n-1),其中F函数其实就是我们本来在维护的含义,即f(n,0)+f(n,1).
那么可以得出:f(n,1)=2*(f(n-1,1)+2*f(n-2,1)).
下面就是从中间部分开始的情况了,首先我们不可能下一步还在同一个列位置的,不然会走不回来,那么就两部分计数,往左走和往右走,而这次一定要回来所以我们用f(i,0)函数,而右边的走法又是F状态的了,另一种情况一样.详见代码.
#include<bits/stdc++.h>
using namespace std;
#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long
#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))
#define pii pair<int,int>
#define fi first
#define se second
#define mk(x,y) make_pair(x,y)
const int mod=1e9+7;
const int maxn=1e4+100;
const int ub=1e6;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){
if(y==0) return x;
return gcd(y,x%y);
}
ll f[maxn][2];
ll solve(int x){
if(x==0) return 1LL;
if(x==1) return 2LL;
return (f[x][0]+f[x][1])%mod;
}
int n;
int main(){
int t;scanf("%d",&t);
f[0][1]=1,f[1][0]=2,f[1][1]=0;
rep(i,2,maxn){
f[i][0]=f[i-1][0]*2LL%mod;
f[i][1]=(solve(i-1)%mod+2LL*solve(i-2)%mod)%mod;
f[i][1]=f[i][1]*2LL%mod;
}
while(t--){
scanf("%d",&n);
if(n==1) {puts("2");continue;}
ll ans=solve(n)*2%mod;
rep(i,2,n){
ans=(ans+powmod(2LL,i)*solve(n-i)%mod)%mod;
ans=(ans+powmod(2LL,n-i+1)*solve(i-1)%mod)%mod;
}
printf("%lld\n",ans);
}
return 0;
}