Placing Bishops on a Chessboard¶
Find the number of ways to place $K$ bishops on an $N \times N$ chessboard so that no two bishops attack each other.
Algorithm¶
This problem can be solved using dynamic programming.
Let's enumerate the diagonals of the chessboard as follows: black diagonals have odd indices, white diagonals have even indices, and the diagonals are numbered in non-decreasing order of the number of squares in them. Here is an example for a $5 \times 5$ chessboard.
Let D[i][j] denote the number of ways to place j bishops on diagonals with indices up to i which have the same color as diagonal i.
Then i = 1...2N-1 and j = 0...K.
We can calculate D[i][j] using only values of D[i-2] (we subtract 2 because we only consider diagonals of the same color as $i$).
There are two ways to get D[i][j].
Either we place all j bishops on previous diagonals: then there are D[i-2][j] ways to achieve this.
Or we place one bishop on diagonal i and j-1 bishops on previous diagonals.
The number of ways to do this equals the number of squares in diagonal i minus j-1, because each of j-1 bishops placed on previous diagonals will block one square on the current diagonal.
The number of squares in diagonal i can be calculated as follows:
int squares (int i) {
if (i & 1)
return i / 4 * 2 + 1;
else
return (i - 1) / 4 * 2 + 2;
}
The base case is simple: D[i][0] = 1, D[1][1] = 1.
Once we have calculated all values of D[i][j], the answer can be obtained as follows:
consider all possible numbers of bishops placed on black diagonals i=0...K, with corresponding numbers of bishops on white diagonals K-i.
The bishops placed on black and white diagonals never attack each other, so the placements can be done independently.
The index of the last black diagonal is 2N-1, the last white one is 2N-2.
For each i we add D[2N-1][i] * D[2N-2][K-i] to the answer.
Implementation¶
int bishop_placements(int N, int K)
{
if (K > 2 * N - 1)
return 0;
vector<vector<int>> D(N * 2, vector<int>(K + 1));
for (int i = 0; i < N * 2; ++i)
D[i][0] = 1;
D[1][1] = 1;
for (int i = 2; i < N * 2; ++i)
for (int j = 1; j <= K; ++j)
D[i][j] = D[i-2][j] + D[i-2][j-1] * (squares(i) - j + 1);
int ans = 0;
for (int i = 0; i <= K; ++i)
ans += D[N*2-1][i] * D[N*2-2][K-i];
return ans;
}