Just to show another approach to the answer already given by C.Blatter, let's start and consider
$$
\left\{ \matrix{
0 \le b_{\,j} \in \mathbb Z \hfill \cr
b_{\,1} + b_2 + b_{\,3} + \cdots + b_{\,k} \le n \hfill \cr} \right.
$$
Then the number of solutions is the integral volume $V(k,n)$ of k-D simplex, with edge segments $(0,n)$, thus of length $n+1$.
It is not difficult to demonstrate that
$$V(k,n)= \binom{n+k}{n}= \binom{n+k}{k}=V(n,k)$$
then the hypothesis is already true in this case.
We can then partition the above solutions into those that contains $j$ variables
null and the remaining $k-j$ having values greater than $0$. We can assign the null variables in $\binom{k}{j}$ ways.
That corresponds to
$$
V(k,n)= \binom{n+k}{n}= \binom{n+k}{k}
= \sum\limits_{0\, \le \,j\, \le \,k} {
\binom{n}{k-j} \binom{k}{j}
}
$$
Coming to our present case, for each non null variable $b_j$ we can assign two values to
the corresponding $a_j$, and just one to the null variables. So in this case the volume will be
$$
V_{\,a} \,(k,n) = \sum\limits_{\left( {0\, \le } \right)\,j\,\left( { \le \,k} \right)} {\binom{n}{k-j} \binom{k}{j} 2^{\,k - j} }
= \sum\limits_{\left( {0\, \le } \right)\,l\,\left( { \le \,k} \right)} {\binom{n}{l} \binom{k}{l} 2^l }
= V_{\,a} \,(n,k)
$$
Therefore the hypothesis is demonstrated.