Well, let's try to use a bit of Fourier Analysis.
As $f$ is continuous, 1-periodic, we can then calculate its Fourier Transform, which is defined by
$$ \hat{f}(n) = \int_0^1 f(t)e^{-2 \pi i t n} dt $$
Using this in the given formula, we have that
$$ \hat{f}(n) = (\hat{f}(n))^2 $$
Which implies that $\hat{f}(n) \in \{0,1\}$. But $f$ is continuous $\Rightarrow f \in L^2(0,1)$, which, by the Plancherel's Theorem, implies that
$$ \sum_n |\hat{f}(n)|^2 < \infty $$
But this is a sum of zeroes and ones. So, there must be an $n_0$ such that $|N| \ge n_0 \Rightarrow \hat{f}(N) = 0 $.
This shows that the solutions to the given problem are exactly the functions of the form
$$ f(x) = \sum_{k \in \mathbb{Z}} \varepsilon_k e^{2 \pi i k x} $$
Where $\varepsilon_k \in\{0,1\}$, and, for $|k|$ sufficiently large, $\varepsilon_k = 0$.