#include<bits/stdc++.h>

using namespace std;

const int N = 5010, mod = 1e9+7;

typedef long long ll;
int n,a,b,k;
ll d[N<<2],P[N<<2];
void update(int x,int s,int t,int p,int c){
	if(s==t){P[p]=c;return ;}
	int mid=(s+t)>>1;
	if(x<=mid) update(x,s,mid,p*2,c);
	else update(x,mid+1,t,p*2+1,c);
	P[p]=(P[p*2]+P[p*2+1])%mod;
}
ll get(int l,int r,int s,int t,int p){
	if(l>r) return 0;
	if(l<=s&&t<=r) return d[p];
	int mid=(s+t)>>1;ll re=0;
	if(l<=mid) re=get(l,r,s,mid,p*2);
	if(r>mid) re+=get(l,r,mid+1,t,p*2+1);
	return re%mod;
}

int main(){
	ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
	
	cin>>n>>a>>b>>k;
	update(a,1,n,1,1);
	for(int i=1;i<=n*4;i++)
		d[i]=P[i];
	while(k--){
		for(int i=1;i<=n;i++){
			if(i==b) continue;
			update(i,1,n,1,-get(i,i,1,n,1));
			int dis=abs(i-b);
//			if(dis/2==0) continue;
//			cout<<i<<' '<<1<<' '<<b-dis/2-1<<' '<<b+dis/2+1<<' '<<n<<'\n';
			if(i<b)
				update(i,1,n,1,(get(1,b-dis/2-1,1,n,1)-get(i,i,1,n,1))%mod);
			else
				update(i,1,n,1,(get(b+dis/2+1,n,1,n,1)-get(i,i,1,n,1))%mod);
		}
		for(int i=1;i<=n*4;i++)
			d[i]=P[i];
//		for(int i=1;i<=n;i++)
//			cout<<get(i,i,1,n,1)<<' ';
//		cout<<'\n';
	}
	
	cout<<(get(1,n,1,n,1)+mod)%mod;
	
	return 0;
}