for( j = i+1; j < M; ++j ) idx[ i ] += N - idx[ M - j + i ];