#pragma GCC optimize("Ofast")
#pragma GCC target("avx2")

char*mmap();
#define RD(v) int v=0;{int _c;while(_c=*rp++-48,_c>=0)v=v*10+_c;}

int wbuf[1<<26];

#define MD 1000000007

long modpow(long a, long b, long m){
  long r=1;
  while(b){
    if(b&1){
      r=r*a%m;
    }
    b>>=1;
    a=a*a%m;
  }
  return r;
}

int inverse(int a){
  int b=MD;
  int u=1;
  int v=0;
  while(b){
    int q=a/b,t;
	t=b, b=a-q*b, a=t;
	t=v, v=u-q*v, u=t;
  }
  if(u<0){
    u+=MD;
  }
  return u;
}

int* f2(wp,v0,v1,i,e)
	int*wp;
{
	int b[6];
	sprintf(b,"%11d\n%11d\n",v0,v1);
	int b0=b[0];
	int b1=b[1];
	int b2=b[2];
	int b3=b[3];
	int b4=b[4];
	int b5=b[5];
	for(i+=i&1;i<e;i+=2){
		*wp++=b0;
		*wp++=b1;
		*wp++=b2;
		*wp++=b3;
		*wp++=b4;
		*wp++=b5;
	}
	return wp;
}

int* f1(wp,v,i,e)
	int*wp;
{
	int b[3];
	sprintf(b,"%11d\n",v);
	int b0=b[0];
	int b1=b[1];
	int b2=b[2];
	for(i+=i&1;i<e;i+=2){
		*wp++=0x0a302020;
		*wp++=b0;
		*wp++=b1;
		*wp++=b2;
	}
	return wp;
}

main(){
	char*rp=mmap(0l,1l<<28,1,2,0,0ll);
	RD(n);
	RD(k);
	int r=modpow(2,k-1,n);
	int q=(modpow(2,k-1,MD)-r+MD)%MD*inverse(n)*(n&1?2:4)%MD;
	int d=modpow((MD+1)/2,k,MD);
	q=1l*q*d%MD;
	r*=2;
	int*wp=wbuf;
	if(n&1){
		if(n<r){
			r-=n;
			q=(q+d)%MD;
			int q0=q;
			int q1=(q+d)%MD;
			if(r<n-r){
				wp=f2(wp,q1,q0,0,r);
				wp=f2(wp,q0,q0,r,n-r);
				wp=f2(wp,q0,q1,n-r,n);
			}else{
				wp=f2(wp,q1,q0,0,n-r);
				wp=f2(wp,q1,q1,n-r,r);
				wp=f2(wp,q0,q1,r,n);
			}
		}else{
			int q0=q;
			int q1=(q+d)%MD;
			if(r<n-r){
				wp=f2(wp,q0,q1,0,r);
				wp=f2(wp,q0,q0,r,n-r);
				wp=f2(wp,q1,q0,n-r,n);
			}else{
				wp=f2(wp,q0,q1,0,n-r);
				wp=f2(wp,q1,q1,n-r,r);
				wp=f2(wp,q1,q0,r,n);
			}
		}
		wp-=3;
	}else{
		if(n<r){
			r-=n;
			q=(q+d)%MD;
			q=(q+d)%MD;
		}
		int q0=q;
		int q1=(q0+d)%MD;
		int q2=(q1+d)%MD;
		if(r<n-r){
			wp=f1(wp,q1,0,r);
			wp=f1(wp,q0,r,n-r);
			wp=f1(wp,q1,n-r,n);
		}else{
			wp=f1(wp,q1,0,n-r);
			wp=f1(wp,q2,n-r,r);
			wp=f1(wp,q1,r,n);
		}
	}
	write(1,wbuf,(wp-wbuf)*sizeof*wbuf);
	_exit(0);
}