#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
int n,a,b,k,dp[5005][5005],tt[5005][5005];
int main(){
	scanf("%d%d%d%d",&n,&a,&b,&k);
	dp[0][a]=1;
	for(int i=1;i<=n;i++)tt[0][i]=tt[0][i-1]+dp[0][i];
	for(int i=1;i<=k;i++){
		for(int j=1;j<=n;j++){
			int t=j;
			if(j==b);
			else if(j<b){
				if((t&1)==(b&1))t--;
				int p=(t+b-1)/2;
				dp[i][j]=(tt[i-1][p]+mod-dp[i-1][j])%mod;
			}else{
				if((t&1)==(b&1))t++;
				int p=(t+b+1)/2;
				dp[i][j]=((tt[i-1][n]+mod-tt[i-1][p-1])%mod+mod-dp[i-1][j])%mod;
			}
			tt[i][j]=(tt[i][j-1]+dp[i][j])%mod;
		}
	}
	printf("%d\n",tt[k][n]); 
}