#include <bits/stdc++.h>
const long long MODULO = 1000 * 1000 * 1000 + 7;
const long long MIN_SOMME = -100 * 1000;
const long long MAX_SOMME = +100 * 1000;
const long long NB_SOMMES = MAX_SOMME - MIN_SOMME + 1;
long long scoreA, scoreB, taille, nbTours;
long long nbPossibilites[2][NB_SOMMES];
inline long long mod(long long a) { return (a % MODULO + MODULO) % MODULO; }
void calculerNbPossibilites() {
  for (long long somme = 0; somme < NB_SOMMES; somme++)
    nbPossibilites[0][somme] = 0;
  nbPossibilites[0][0 - MIN_SOMME] = 1;
  for (long long tour = 0; tour < nbTours; tour++) {
    long long sommePossibilites = 0;
    for (long long delta = 0; delta <= taille; delta++)
      if (delta < NB_SOMMES)
        sommePossibilites =
            mod(sommePossibilites + nbPossibilites[tour % 2][delta]);
    for (long long somme = 0; somme < NB_SOMMES; somme++) {
      nbPossibilites[(tour + 1) % 2][somme] = sommePossibilites;
      if (0 <= somme - taille)
        sommePossibilites =
            mod(sommePossibilites - nbPossibilites[tour % 2][somme - taille]);
      if (somme + taille + 1 < NB_SOMMES)
        sommePossibilites = mod(sommePossibilites +
                                nbPossibilites[tour % 2][somme + taille + 1]);
    }
  }
}
int main() {
  scanf("%lld%lld%lld%lld", &scoreA, &scoreB, &taille, &nbTours);
  calculerNbPossibilites();
  long long totalPossibilites = 0;
  long long sommeB = MIN_SOMME;
  long long sommePossibilites = 0;
  for (long long sommeA = MIN_SOMME; sommeA <= MAX_SOMME; sommeA++) {
    while (sommeB <= MAX_SOMME && scoreA + sommeA > scoreB + sommeB) {
      sommePossibilites = mod(sommePossibilites +
                              nbPossibilites[nbTours % 2][sommeB - MIN_SOMME]);
      sommeB++;
    }
    totalPossibilites = mod(totalPossibilites +
                            nbPossibilites[nbTours % 2][sommeA - MIN_SOMME] *
                                sommePossibilites);
  }
  printf("%lld", totalPossibilites);
  return 0;
}
