Formulas for evaluating analytic energy gradient are derived for combined time-dependent density functional theory (TDDFT) and polarizable force field methods that incorporate dipole polarizability tensors and linearly induced point dipoles. The Z-vector method for determining relaxed one-particle difference density matrix in regular TDDFT methods is extended to include induced dipoles. The analytic gradient of the mutual polarization energy of the force field and the TDDFT excited state can be formulated by using the TDDFT difference density-induced dipoles and the transition state density-induced dipoles. All the forces and torques involving induced dipoles can be efficiently evaluated using standard electrostatic formulas as if the induced dipoles were permanent dipoles. The formulas are given in the most general form and are applicable to various flavors of polarizable force fields. Implementation and tests with a polarizable five-point water model show that the formulas are rigorous. The carbonyl vibration modes and infrared spectrum intensities of a cluster formed by acetone and two water molecules are studied.