Problem:
Score : 1600 points

Problem Statement

There are N zeros and M ones written on a blackboard.
Starting from this state, we will repeat the following operation: select K of the rational numbers written on the blackboard and erase them, then write a new number on the blackboard that is equal to the arithmetic mean of those K numbers.
Here, assume that N + M - 1 is divisible by K - 1.

Then, if we repeat this operation until it is no longer applicable, there will be eventually one rational number left on the blackboard.

Find the number of the different possible values taken by this rational number, modulo 10^9 + 7.

Constraints

1 ≦ N, M ≦ 2000

2 ≦ K ≦ 2000

N + M - 1 is divisible by K - 1.

Input

The input is given from Standard Input in the following format:

N M K

Output

Print the number of the different possible values taken by the rational number that will be eventually left on the blackboard, modulo 10^9 + 7.

Sample Input 1

2 2 2

Sample Output 1

5

There are five possible values for the number that will be eventually left on the blackboard: \frac{1}{4}, \frac{3}{8}, \frac{1}{2}, \frac{5}{8} and \frac{3}{4}.

For example, \frac{3}{8} can be eventually left if we:

Erase 0 and 1, then write \frac{1}{2}.

Erase \frac{1}{2} and 1, then write \frac{3}{4}.

Erase 0 and \frac{3}{4}, then write \frac{3}{8}.

Sample Input 2

3 4 3

Sample Output 2

9

Sample Input 3

150 150 14

Sample Output 3

937426930

Solution in C++14 (GCC 5.4.1) (judge verdict: Accepted)

#include <iostream>
#include <fstream>
#include <set>
#include <map>
#include <string>
#include <vector>
#include <bitset>
#include <algorithm>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <cassert>
#include <queue>

#define mp make_pair
#define pb push_back


typedef long long ll;
typedef long double ld;

using namespace std;
const ll MOD = 1e9 + 7;
const ll MX = 4010;

ll dp[MX];
ll dp2[MX];

int main() {
	int n, m, k;
	cin >> n >> m >> k;
	int k2 = k - 1;
	ll ans = 0;
	dp[0] = 1;
	int it = (n + m - 1) / k2;
	for (int i = 1; i <= it; ++i) {
		for (int j = 0; j <= m; ++j)
			dp2[j] = 0;
		for (int j = 0; j <= m; ++j) {
			for (int l = 0; l < k && j + l <= m; ++l)
				dp2[j + l] = (dp2[j + l] + dp[j]) % MOD;
		}
		for (int j = 0; j <= m; ++j)
			dp[j] = dp2[j];
		for (int j = 1; j <= m; ++j) {
			int x = j;
			int y = i * k2 - x;
			++y;
			if (m >= x && (m - x) % k2 == 0 && n == y) {
				ans = (ans + dp[j]) % MOD;
			}
		}
	}
	cout << ans << "\n";
	return 0;
}