If we focus on the integrand, then the substitution u=3x+1 seems suitable, since we then get u which we can integrate. There is also no risk involved in using a linear substitution such as u=3x+1, because the relation between dx and du will be a constant factor,