While solving a problem I came across the following identity, which holds by numerical evidence: $$ \sum_{k=1}^i\frac1k\binom{i}{k-1}\binom kj{B_{k-j}}=\delta_{ij}. $$ where $B$ are the Bernoulli numbers.
I have no experience with Bernoulli numbers, so any hint for proving the equality will be appreciated.
Allow me to change the notation as to reserve (as far as possible) to $i$ it's common meaning, and use $k,j,l$ as indices.
First let's simplify the sum working on the properties of the binomial coefficients $$ \eqalign{ & S(n,m) = \sum\limits_{1 \le \,k\, \le \,n} {{1 \over k}\left( \matrix{ n \cr k - 1 \cr} \right)\left( \matrix{ k \cr m \cr} \right)B_{\,k - m} } = \cr & = {1 \over {n + 1}}\sum\limits_{1 \le \,k\, \le \,n} {\left( \matrix{ n + 1 \cr k \cr} \right)\left( \matrix{ k \cr m \cr} \right)B_{\,k - m} } = \cr & = {1 \over {n + 1}}\left( \matrix{ n + 1 \cr m \cr} \right)\sum\limits_{1 \le \,k\, \le \,n} {\left( \matrix{ n + 1 - m \cr k - m \cr} \right)B_{\,k - m} } = \cr & = {1 \over {n + 1}}\left( \matrix{ n + 1 \cr m \cr} \right)\sum\limits_{\max \left( {1 - m,0} \right) \le \,j\, \le \,n - m} {\left( \matrix{ n + 1 - m \cr j \cr} \right)B_{\,j} } \cr} $$
where:
- in the first step, we use the "absorption" indentity;
- in the second step, we use the "trinomial revision" indentity;
- in the third step, we changed the summation index.
Now, assuming $1 \le m$, we can use the fundamental reursive identity of Bernoulli numbers (the "standard" definition $B_{\,j} ^ -$) $$ \sum\limits_{0 \le \,j\, \le \,n - m} {\left( \matrix{ n + 1 - m \cr j \cr} \right)B_{\,j} } = \delta _{\,n - m,\,0} $$
to get $$ S(n,m) = \sum\limits_{1 \le \,k\, \le \,n} {{1 \over k}\binom{n}{k-1} \binom{k}{m} B_{\,k - m}^{\, - } } = \delta _{\,n,\;m} \quad \left| {\;1 \le n,m} \right. $$