In dilute solutions of rodlike polyelectrolytes some counterions are distributed far from polyions while others are located in their vicinity in the regions of cylindrical symmetry of the electrostatic potential. For these cylindrical regions around rodlike polyelectrolytes we find an exact solution of the nonlinear Poisson-Boltzmann equation for the case of nonzero net charge in these regions. This exact solution implies three qualitatively different phases of counterion distribution around the polyions with second order phase transitions between these phases.