Simulating overtaking interactions (OIs) of dust-acoustic N-shock waves (DANSWs) in dust plasmas with opposite polarities and q-distributed electrons and Boltzmann distributed ions are examined in this work. Reductive perturbation theory is used to obtain a Burgers’ equation (BE), and dust-acoustic N-shock wave solutions are derived using the Cole-Hopf transformation and exponential function. Numerical simulations show that the evolution and OIs of DANSWs and their electric fields are modified by the effects of negative/positive dust fluid kinematic viscosity, the ratio of negative dust to positive dust mass, and the nonextensive q-distribution of electrons. This work may help explain the properties of dust-acoustic shock waves in cometary tails and Jupiter's magnetosphere that support the propagation of nonlinear waves.