Many textbooks will explain how the Euler product formula can be used as a basis for proving the infinite of primes:
$$\sum_{n}\frac{1}{n^{s}}=\prod_{p}\left(1-\frac{1}{p^{s}}\right)^{-1}$$
My concern is that this is a circular argument because the derivation of the product formula uses the fact there are infinite primes.
Maths history texts will suggest that Euler's derivation of the product formula was done by a process of "sieving the zeta function", outlined next:
We start with the zeta function, having proven (elsewhere) $s>1$ for convergence.
$$\zeta(s)=1+\frac{1}{2^{s}}+\frac{1}{3^{s}}+\frac{1}{4^{s}}+\frac{1}{5^{s}}+\frac{1}{6^{s}}+\ldots$$
We can divide this series by $2^{s}$.
$$\frac{1}{2^{s}}\zeta(s)=\frac{1}{2^{s}}+\frac{1}{4^{s}}+\frac{1}{6^{s}}+\frac{1}{8^{s}}+\frac{1}{10^{s}}+\frac{1}{12^{s}}\ldots$$
These denominators are multiples of $2^{s}$. By subtracting these terms from $\zeta(s)$, we sieve out terms with these multiples of $2^{s}$.
$$(1-\frac{1}{2^{s}})\cdot\zeta(s)=1+\frac{1}{3^{s}}+\frac{1}{5^{s}}+\frac{1}{7^{s}}+\frac{1}{9^{s}}+\frac{1}{11^{s}}+\ldots$$
Repeating for $3^s$ gives us:
$$(1-\frac{1}{3^{s}})\cdot(1-\frac{1}{2^{s}})\cdot\zeta(s)=1+\frac{1}{5^{s}}+\frac{1}{7^{s}}+\frac{1}{11^{s}}+\frac{1}{13^{s}}+\ldots$$
And again for $5^s$:
$$(1-\frac{1}{5^{s}})\cdot(1-\frac{1}{3^{s}})\cdot(1-\frac{1}{2^{s}})\cdot\zeta(s)=1+\frac{1}{7^{s}}+\frac{1}{11^{s}}+\frac{1}{13^{s}}+\ldots$$
Repeating the process, we can only remove multiples of successive primes. We can't remove multiple of a number C if it is non-prime, C=AB, because we would have already sieved them when we removed multiples of its factors A or B.
In the limit of this process (having removed multiples of an infinitude of primes), only 1 remains on the RHS, and on the left we have an infinite product over all primes:
$$\ldots\cdot(1-\frac{1}{11^{s}})\cdot(1-\frac{1}{7^{s}})(1-\frac{1}{5^{s}})\cdot(1-\frac{1}{3^{s}})\cdot(1-\frac{1}{2^{s}})\cdot\zeta(s)=1$$
Back to my question.
This derivation assumes there are an infinitude of primes in order to attain the single 1 on the RHS. Therefore this Euler product can't be used as basis to say that there is an infinitude of primes.
Where is the flaw in my logic?
Comment: I am aware of the other derivation which uses the FTA and doesn't assume an infinitude of primes, but this is the method used by Euler, and he used it to demonstrate the infinitude of primes, for the first time since Euclid, according to the maths history books.
I am not a university trained mathematician so would appreciate replies with minimal technical terminology.