Our previously developed analytical infinite order sudden (IOS) quantum theory of triatomic photodissociation is extended to describe indirect photodissociation processes through a real or virtual intermediate state. The theory uses the IOS approximation for the dynamics in the final dissociative channels and an Airy function approximation for the continuum states. These approximations enable us to evaluate the multi-dimensional non-separable transition amplitudes analytically (as one-dimensional quadratures), despite the different natural coordinates for the initial bound, the intermediate resonant, and the final dissociative states. The fragment internal energy distributions are described as a function of the initial and final quantum states and the photon excitation energy. The theory readily permits the evaluation of rotational distributions for high values of the total angular momentum J in the initial bound molecular state, a feature that would be very difficult with close-coupled methods. In paper II we apply the theory to describe the photofragment yield spectrum of NOCl in the region of the T-1(1(3)A ")<--S-0(1(1)A’) transition.