#include <bits/stdc++.h>
const int P = 1000000007;
using namespace std;
int f[2][200005], sum[200005];
int main() {
  int a, b, k, t;
  scanf("%d%d%d%d", &a, &b, &k, &t);
  int i, j;
  f[0][k * t + 1] = 1;
  for (i = 1; i <= t; i++) {
    int cur = i & 1, pre = cur ^ 1;
    sum[0] = 0;
    for (j = 1; j <= t * k * 2 + 1; j++) sum[j] = (sum[j - 1] + f[pre][j]) % P;
    for (j = -k * i; j <= k * i; j++)
      f[cur][j + k * t + 1] = (sum[min(k * t * 2 + 1, j + k * t + 1 + k)] -
                               sum[max(0, j + k * t + 1 - k - 1)] + P) %
                              P;
  }
  sum[0] = 0;
  for (j = 1; j <= t * k * 2 + 1; j++) sum[j] = (sum[j - 1] + f[t & 1][j]) % P;
  int ans = 0;
  for (i = -k * t; i <= k * t; i++) {
    ans += 1LL * f[t & 1][i + k * t + 1] *
           ((sum[max(0, min(k * t * 2 + 1, i + a - b + k * t))] - sum[0] + P) %
            P) %
           P;
    if (ans > P) ans -= P;
  }
  printf("%d\n", ans);
  return 0;
}
