The aim of this study is to find a generic method for generating a path of the solution of a given stochastic differential equation which is more efficient than the standard Euler-Maruyama scheme with Gaussian increments. First we characterize the asymptotic distribution of pathwise error in the Euler-Maruyama scheme with a general partition of time interval and then, show that the error is reduced by a factor (d+2)/d when using a partition associated with the hitting times of sphere for the driving d-dimensional Brownian motion. This reduction ratio is the best possible in a symmetric class of partitions. Next we show that a reduction which is close to the best possible is achieved by using the hitting time of a moving sphere which is easier to implement.