ENH: implement `digitize` with `PyArray_SearchSorted`

Closes #2656 and replaces #4184.
Maybe add "In numpy versions < 1.10 the array must be 1-D."