#include<bits/stdc++.h>
using namespace std;
const int mn=5005,mod=1e9+7;
int f[mn][mn];
int sum[mn];
int n,a,b,k;
int main()
{
    scanf("%d%d%d%d",&n,&a,&b,&k);
    for(int i=1;i<=n;i++)
    {
        if(i!=b)f[i][0]=1;
        sum[i]=sum[i-1]+f[i][0];
    }
    for(int i=1;i<=k;i++)
    {
        for(int j=1;j<=n;j++)
        {
            int dis=abs(j-b)-1;
            f[j][i]=((((sum[min(j+dis,n)]-sum[max(j-dis,1)-1])%mod+mod)%mod-f[j][i-1])%mod+mod)%mod;
        }
        for(int j=1;j<=n;j++)
        {
            sum[j]=(sum[j-1]+f[j][i])%mod;
        }
    }
    printf("%d",f[a][k]);
}