The divmod() function of Numpy returns the element-wise quotient and remainder simultaneously. You can then use the zip() function to get the result in your desired format.
Here is an example:
>>> import numpy as np>>> x=[12,13,14,15]>>> y=[2,3,4,5]>>> q,r=np.divmod(x,y)>>> qarray([6, 4, 3, 3])>>> rarray([0, 1, 2, 0])>>> list(zip(q,r))[(6, 0), (4, 1), (3, 2), (3, 0)]