We present p-adaptive hybridized flux reconstruction schemes to reduce the computational cost of implicit discretizations using a nondimensional vorticity indicator. At each adaptation level, we apply a projection operation to the new space based on the element-wise projected solution and the transmission conditions. We validate our implementation and analyze performance via numerical examples. Specifically, we show that p-adaptation of hybridizable flux reconstruction (HFR) methods results in comparable numerical error to standard p-adaptive and p-uniform standard flux reconstruction (FR) discretizations with a fraction of degrees of freedom. Performance results for a cylinder at Re = 150 showcase speedup factors in excess of 6 for hybridized methods in comparison with p-adaptive standard FR schemes, and up to 40 compared to p-uniform FR discretizations. Similarly, results for a NACA 0012 airfoil at Re = 10, 000 demonstrate speedup factors close to 6 compared to p-adaptive FR discretizations, and up to 33 compared to p-uniform conventional FR. Hence, combinining hybridization with adaptation yields a significant reduction in computatonal cost compared to conventional implicit discretizations.