A general automatic method for optimal construction of matrix product operators using bipartite graph theory