The answer can be expressed with the formula:

\[\large{1 + \sum_{b=2}^{\lfloor \log_2 N \rfloor} -\mu(b) \cdot (\lfloor \sqrt[b]N \rfloor - 1)}\]

where \(\mu(b)\) represents the Möbius function. The explanation is similar to the first official editorial.