Nonextensive statistics has attracted attention as a description of particle spectra in nuclear collisions at QCD energies. First, we construct the equation of state by incorporating Tsallis statistics based on the hadron resonance gas and parton gas models. Thermodynamic conditions are found to impose constraints on the q-parameter of Tsallis distribution. Next, we apply the equation of state to the relativistic hydrodynamic modeling of nuclear collisions. The Cooper-Frye prescription is consistently modified. Numerical demonstrations indicate that the model may describe charged particle spectra at Large Hadron Collider in the transverse momentum range up to 6-8 GeV. Elliptic flow, on the other hand, suggests a narrower range of applicability.