\section{Introduction}
Let me begin with an example which explains the problem I consider in this paper.
Consider the elliptic curve $E:y^2+y=x^3-x/{\mathbb Q}$. This is an elliptic curve with conductor $37$. Let $p$ be a prime of good reduction, here and at all relevant places I will always assume $p$ is a prime of good reduction and call such a prime a good prime (for $E$). Let $N_p(E)$ be the number of points on this elliptic curve modulo $p$.
It is a well-known theorem of Hasse and Weil \citep{silverman-arithmetic} that
\begin{equation}\label{eq:hw-bound}
p+1-2\sqrt{p}\leq \np{p} \leq p+1+2\sqrt{p}.
\end{equation}

A simple computation with \citep{sage} and \citep{pari} reveals that following equalities hold:
\begin{eqnarray}\label{eq:elliptic-examples}
\np{2}\cdot\np{13}\cdot\np{43}=& 3360 &=
\np{3}\cdot\np{5}\cdot\np{67}\\
\np{5}\cdot\np{43}\cdot\np{73}=&25200&=\np{17}\cdot\np{19}\cdot\np{61}\\
\np{101}\cdot\np{107}\cdot\np{251}=&3107520&=\np{113}\cdot\np{127}\cdot\np{167}\ {\rm or}\\
99\cdot 120\cdot 254 =&3107520 &= 132 \cdot  127\cdot  180.\\
\np{1009}\cdot\np{1181}\cdot\np{1601}=&1988217000&=\np{1063}\cdot\np{1283}\cdot\np{1399}\ {\rm or}\\
1057\cdot 1125\cdot 1648 = & 1988217000 & = 1057 \cdot 1320\cdot 1425.
\end{eqnarray}
All  the triples of primes are all pairwise distinct; moreover the six numbers which enter these equalities are typically distinct but occasionally not.

More generally, for a fixed $k\geq 3$, any elliptic curve $E/{\mathbb Q}$, and a given  natural number $n\geq 1$, I consider the equation (in pair wise distinct primes $p_1,\ldots,p_k$)
\begin{equation}\label{eq:erdos}
\np{p_1}\cdots \np{p_k} = n.
\end{equation}
Let $G_k(E,n)$ (for $k\geq 3$) be the number of solutions to \eqref{eq:erdos} in pairwise distinct primes $p_1,\ldots,p_k$. The number $G_k(E,n)$ is always finite. Indeed for any prime $p$ one has
\begin{equation}
\frac{p+1-2\sqrt{p}}{p}\geq\begin{cases}
	1-\frac{2}{\sqrt{p}}\geq \frac{1}{100} & {\rm if}\ p\geq 5\\
	1+\frac{1}{3}-\frac{2}{\sqrt{3}}\geq \frac{1}{100} & {\rm if}\ p=3\\
	1+\frac{1}{2}-\frac{2}{\sqrt{2}}\geq \frac{1}{100} & {\rm if}\ p=2\\
\end{cases}
\end{equation}
Hence
$N_p(E)\geq p/100$ for all primes $p$ and so  $$n\geq  \frac{p_1\cdots p_k}{100^k}.$$ Thus the primes which contribute to $G_k(E,n)$ are finite in number and hence $G_k(E,n)$ is finite for any $n$. The purpose of this paper is to prove  the following theorem which shows that  these equalities \eqref{eq:elliptic-examples} are not small numerical accidents.
\begin{theorem}\label{th:main}
Let $E/\Q$ be an elliptic curve. If $E$ does not have complex multiplication assume GRH holds for $E$.\
Then for any $k\geq 3$
$$
	\varlimsup_{n\to\infty} G_k(E,n)=\infty.
$$
More precisely, for every $k\geq 3$ and for all sufficiently large $x$ there exist  integers $n\leq x$ with $$G_k(E,n)\geq (\log x)^{\delta}$$
for some $\delta=\delta(k)>0$.
\end{theorem}
It is unlikely that the growth of $G_k(E,n)$ given by the proof of Theorem~\ref{th:main} is optimal. Let me point of that  in \citep{erdos36}, Paul Erdos considered the question of solutions (in triples of pairwise distinct primes $p_1,p_2,p_3$) to the equation
\begin{equation}
(p_1-1)(p_2-1)(p_3-1)=n,\end{equation}
and proved the corresponding result in this case
and that \textit{loc. cit.} served as an inspiration to this note. Let me also point out that one can also replace the assumption that GRH holds for $E$ by other types of hypothesis such as  $\theta$-quasi-GRH etc., however since I make no claim of providing sharpest error terms nor best possible constants, the choice of stronger or weaker hypothesis is insignificant (to me).
While \citep{erdos36} served as inspiration for this paper, my approach  differs from that of \textit{loc. cit.} in several important points. Notably I use the fourth moment to prove Lemma~\ref{le:convergence} instead of the Hardy-Ramanujan Theorem. I also prove Theorem~\ref{th:erdos-product}, which was proved in \citep{erdos36} for $k=3$, a bit differently. Note that Theorem~\ref{th:main} assumes $k\geq 3$.  But in fact I expect that the result also holds for $1\leq k\leq 2$.  I hope to pursue $k=2$  in a separate paper.

The case $k=1$ seems to be interesting in its own right: I conjecture that the result holds for $k=1$. In other words I conjecture that there exists a sequence of integers $n\to\infty$ and a sequence of primes $p_1<p_2<\cdots<p_m$ with $m\to\infty$ (with $n$) such that $$\np{p_1}=\np{p_2}=\cdots=\np{p_m}=n.$$ I call this an elliptic progression of primes and thus the assertion for $k=1$ is equivalent to existence of arbitrarily long sequences of primes in elliptic progressions (for any given elliptic curve $E/{\mathbb Q}$).  This seems difficult at the moment. Even the weakest assertion: that for any $E/{\mathbb Q}$ there exists infinitely many pairs of primes $\np{p}=\np{q}$ seems difficult at the moment. See Section~\ref{re:conj} for a detailed discussion of what I expect to be true  for $k=1$ and numerical examples. Let me also point out that in \citep{silverman11} it has been conjectured that that there exists infinitely many pairs of primes $p,q$ such that $\np{p}=q$ and $\np{q}=p$ (I am grateful to Joseph Silverman for providing this reference).

I thank M.~Ram Murty for his comments which have improved the readability of this paper and for pointing out that in \citep{murty84a} and \citep{murty84b}   estimation of higher moments (used here and calculated by the sieving method of \citep{granville07}) is carried using a different method.

\section{Basic Estimates}
Throughout this paper $\vare$ will be  a positive number which will be sufficiently small.  Symbols $c_1,c_2,\ldots$ will be positive constants as will $c_1(\vare),c_2(\vare),\ldots$ which will depend  on $\vare$. In either of the cases the  numerical values of these constants will be immaterial and I caution the reader that in some situations the same constant may be denoted by different symbols. The letter $x\geq 1$ will denote a real number and symbols $x_0,x_0(\vare),x_1,x_1(\vare),\ldots$ will also be real numbers whose values will not be important to us except for the fact that such numbers will exist (in the contexts where they appear) and again such symbols may occur in multiple contexts (and may differ from the ones which appearing in other contexts). Hopefully there will be no confusion caused by this notational conflation which I will indulge in through out this paper.

For an integer $n$ let $\omega(n)$ be the number of distinct prime factors of $n$ and let $\Omega(n)$ be the number of prime factors of $n$, counted with multiplicity. Let $E$ be an elliptic curve over ${\mathbb Q}$. Let $E[m]$ be the $m$-torsion of the elliptic curve $E$ for an integer $m\geq 1$. I will say that \textit{GRH holds for $E/{\mathbb Q}$} if  the Generalized Riemann Hypothesis holds for the Dedekind zeta function of the finite Galois extension ${\mathbb Q}(E[m])/{\mathbb Q}$ for every integer $m\geq 1$. A prime $p$ will be called a \textit{good prime for} $E$ if $E$ has good reduction at $p$.

Let $d$ be a square-free integer. Let $\pi_E(x,d)$ be the number of good primes $p\leq x$ such that $d|\np{p}$. Write
\begin{equation}
\pi_E(x,d)=\delta(d){\rm Li}(x)+r_d,
\end{equation}
where ${\rm Li}(x)=\int_2^x\frac{dt}{\log(t)}$, $\delta(d)$ is a positive fraction  (whose explicit form is not needed for the moment) and $r_d$ is the ''error term.''
If we assume that GRH holds for $E$ then the precise form of error term can be made explicit (this the only point where one needs to assume GRH). This is done by means of the explicit Chebotarev density theorem of \citep{murty88,serre82} and hence one has the following estimate for $r_d$ due to \citep{steuding05,steuding05b}:
\begin{equation}
|r_d|=O\left(d^{3/2}x^{1/2}\log(dx) \right),
\end{equation}
and the implied constant is dependent on $E$ (specifically the conductor of $E$).

The following result of \citep{liu05,liu06} will be crucial in the proof of Theorem~\ref{th:main}.
\begin{theorem}\label{th:liu}
Let $E/\Q$ be an elliptic curve. If $E$ does not have complex multiplication assume GRH holds for $E$.\ Then there exists a positive constant $c_i$ such that for all $x$ sufficiently large one has
$$\sump x (\omega(\np{p})-\log\log x)^2=c_{\arabic{constantCounter}\stepcounter{constantCounter}} x\log\log x(1+o(1)).$$
\end{theorem}

Theorem~\ref{th:liu} says that the numbers $\omega(\np{p})$ follow a normal distribution with mean $\log\log x$ and variance $(\log\log x)^{1/2}$. The next lemma computes the fourth moment of this distribution.
\begin{lemma}\label{le:fourth-moment}Let $E/\Q$ be an elliptic curve. If $E$ does not have complex multiplication assume GRH holds for $E$.\ There is a constant $c_i$ such that
for all $x$ sufficiently large, one has
$$
\sump x \left(\omega(\np{p})-\log\log x\right)^4=c_{\arabic{constantCounter}\stepcounter{constantCounter}} \pi(x)(\log\log x)^2+o(\pi(x)(\log\log x)^2).
$$
\end{lemma}
\begin{proof}
This proved by the sieving argument of \citep{granville06}. An alternative method to this using \citep{murty84a} and \citep{murty84b} was suggested to me by M.~Ram Murty.  Let me write $z=x^\beta$ where $\beta<1$ is a constant. At the end of this proof one is led to the choice  $\beta=\frac{1}{29}$. Let $P=\prod_{q<z} q$.   Define $$\omega_z(n) = \sum_{q|n, q<z}1,$$ so $\omega_z(n)$  counts the number of distinct prime divisors  $q|n$ such that $q<z$. As $z=x^\beta$, any $n\leq x$ has at most a bounded number of prime divisors $q\geq z$. So one has
\begin{eqnarray}
\omega(\np{p})-\log\log x & =& (\omega_z(\np{p})-\log\log x)
					  + (\omega(\np{p})-\omega_z(\np{p})) \\
				  & = &  (\omega_z(\np{p})-\log\log x) + O(1).
\end{eqnarray}
Thus
\begin{eqnarray}
(\omega(\np{p})-\log\log x)^4 & = & \left((\omega_z(\np{p})-\log\log x) + (\omega(\np{p})-\omega_z(\np{p}))\right)^4\\
& = & ((\omega_z(\np{p})-\log\log x) + O(1))^4.
\end{eqnarray}
So one has
\newcommand{\operatorname{O}\nolimits}{\operatorname{O}\nolimits}
\begin{eqnarray}
\sump x (\omega(\np{p})-\log\log x)^4 & = & \sump x (\omega_z(\np{p})-\log\log x)^4+\\
			& & \qquad+ \operatorname{O}\nolimits\left( \sump x (\omega_z(\np{p})-\log\log x)^3 \right)+ \\
			&  & \qquad+ \operatorname{O}\nolimits\left( \sump x (\omega_z(\np{p})-\log\log x)^2 \right)+ \\
			&  & \qquad+ \operatorname{O}\nolimits\left( \sump x (\omega_z(\np{p})-\log\log x)^x \right)+ \\
			&  & \qquad+ \operatorname{O}\nolimits\left(\pi(x)\right). \\
\end{eqnarray}
By \citep[Proposition 3]{granville06} one has
\begin{eqnarray}
\sump x (\omega_z(\np{p})-\log\log x)^4 &=& c_{\arabic{constantCounter}\stepcounter{constantCounter}} \pi(x)(\log\log x)^2\left(1+\operatorname{O}\nolimits\left(\frac{1}{\log\log x}\right)\right)\\
                               && +\operatorname{O}\nolimits\left((\log\log x)^4\sum_{d\in D_4} |r_d|\right),
\end{eqnarray}
where in the last sum  ``$d\in D_4$'' is short for sum over  sum over all square-free products  $d=q_1\cdots q_4$ of at most four primes $q_i$ each $q_i<z$.
Similarly
\begin{eqnarray}
\sump x (\omega_z(\np{p})-\log\log x)^3 &\leq& c_{\arabic{constantCounter}\stepcounter{constantCounter}} \pi(x)\log\log x\left(1+\operatorname{O}\nolimits\left(\frac{1}{\log\log x}\right)\right)\\
&& +\operatorname{O}\nolimits\left((\log\log x)^3\sum_{d\in D_3} |r_d|\right),
\end{eqnarray}
where in the last sum  ``$d\in D_4$'' is short for sum over  sum over all square-free products  $d=q_1\cdots q_4$ of at most four primes $q_i$ each $q_i<z$.

The sum $\sump x (\omega_z(\np{p})-\log\log x)^2$ is estimated by Theorem~\ref{th:liu} and is $\operatorname{O}\nolimits(\pi(x)\log\log x)$ and so by Cauchy-Schwartz inequality one has
$$|\sump x (\omega_z(\np{p})-\log\log x)|\ll \pi(x)\log\log x.$$
Let me indicate how to estimate the term (the estimation of other terms is similar to this one)
\begin{equation}\label{eq:error-erdos-kac}
\operatorname{O}\nolimits\left((\log\log x)^4\sum_{d}|r_d|\right)
\end{equation}
where the sum over $d$ runs over all square-free integers which are products of at most four primes $q<z$.  Each $d$ which contributes to the sum is at most a product of four primes $q<z$ so
\begin{equation}
d=q_1\cdots q_4 <z^4,
\end{equation}
and hence $d^{3/2}\log(d)\leq d^{5/2}<(z^4)^{(5/2)}=z^{10}$. Further since there are $\pi(z)$ primes less that $z$, so the number of such $d$ is $O(\pi(z)^4)$.  Thus one has
\begin{equation}
\sum|r_d| = O\left(x^{1/2}\log x z^{10}\pi(z)^4 \right)=O\left(x^{1/2}\log x \frac{z^{14}}{(\log(z))^4}\right)
\end{equation}
which is
\begin{equation}
O\left(\frac{x^{1/2} x^{14\beta}}{(\log x)^3}\right)
\end{equation}
Now choose $\beta$ so that $14\beta+\frac{1}{2}<1$ which gives $\beta<\frac{1}{28}$. So if we choose $\beta=\frac{1}{29}$ one  gets
\begin{equation}
\sum_d|r_d|=O\left(\frac{x^{57/58}}{(\log x)^3}\right),
\end{equation}
and finally with this estimate the last term in \eqref{eq:error-erdos-kac} is
\begin{equation}
\leq c_{\arabic{constantCounter}\stepcounter{constantCounter}} \frac{x^{57/58}(\log\log x)^4}{(\log x)^3} = o(\pi(x)(\log\log x)^2).
\end{equation}
The sum
\begin{equation}\label{eq:error-erdos-kac2}
\operatorname{O}\nolimits\left((\log\log x)^3\sum_{d\in D_3}|r_d|\right)
\end{equation}
can also be estimated similarly and is $o(\pi(x)(\log\log x)^2)$. The sum $\operatorname{O}\nolimits\left(\log\log x\sum_{d\in D_1}|r_d|\right)$ can also be estimated similarly and is certainly $o(\pi(x)(\log\log x)^2)$. This completes the proof of Lemma~\ref{le:fourth-moment}.

\end{proof}

The following theorem was proved in \citep{erdos36} for $k=3$.
\begin{theorem}\label{th:erdos-product}
Let $k\geq 3$ be an integer. Let $0<\vare < \frac{2}{(20(k^2+k))}$ be a  real number. Let
\begin{equation}
\setx{D}=\left\{n=n_1\cdots n_k\leq x: \omega(n_i)>(1-\vare)\llx \ {\rm for\ } 1\leq i\leq k \right\}.
\end{equation}
Then
\begin{equation}
\csetx{D} \leq c_{\arabic{constantCounter}}(\vare)\stepcounter{constantCounter} \frac{x}{(\log x)^{1+\delta}},
\end{equation}
with  $\delta=\log(2)(k-\vare\frac{k^2+k}{2})-2>\log(2)(k-\frac{1}{10})-2>0$.
\end{theorem}
\begin{proof}
Let $n\in \setx{D}$. Let us write $n=m_1\cdot m_2$ where $m_1$ is square-free and for every prime $q|m_2$ one has $q^2|m_2$. Let ${\mathcal{D}}_1(x)=\left\{n\in \setx{D}: \omega(m_2)>\vare\log\log x \right\}$ and let ${\mathcal{D}}_2(x)=\left\{n\in \setx{D}: \omega(m_2)\leq\vare\log\log x \right\}$. Then clearly $\setx{D}={\mathcal{D}}_1(x)\cup{\mathcal{D}}_2(x)$ and ${\mathcal{D}}_1(x)\cap{\mathcal{D}}_2(x)=\emptyset$. The cardinalities of each of ${\mathcal{D}}_1(x)$ and ${\mathcal{D}}_2(x)$ will be estimated separately.

If $n\in{\mathcal{D}}_1(x)$ then there is a divisor $d|n$ such that $\omega(d)=j$ where $j=[\vare\log\log x]$ and for each prime $q|d$ one has $q^2|n$. Hence  $|\dxo|$ is less than the number of such numbers less than $x$ i.e.,
\begin{equation}
|\dxo| <c_{\arabic{constantCounter}\stepcounter{constantCounter}} x\left(\sum_{q,\alpha\geq 2}\frac{1}{q^\alpha}\right)^j\frac{1}{j!},
\end{equation}
where the sum is over all primes $q$ and all integers $\alpha\geq 2$ and clearly the sum converges so let
$K=\sum_{q,\alpha\geq 2}\frac{1}{q^\alpha}$. Then
\begin{equation}
|\dxo| <c_{\arabic{constantCounter}\stepcounter{constantCounter}} \frac{x K^j}{j!}<c_{\arabic{constantCounter}\stepcounter{constantCounter}} x\left(\frac{K}{j}\right)^j.
\end{equation}
Now as $j=[\vare\log\log x]$ note that one has the trivial bounds (for all $x$ sufficiently large):
\begin{eqnarray}
 j & \leq & \vare\log\log x\\
 j & \geq & \frac{1}{2}\vare\log\log x.
\end{eqnarray}
Hence one has
\begin{eqnarray}
K^j & = & e^{j\log(K)}\leq e^{\vare\log\log x\log(K)},\\
    & = & e^{\log((\log x)^{\vare\log(K)})} \\
    & = & (\log x)^{\vare\log(K)},
\end{eqnarray}
and
\begin{equation}
\frac{1}{j}\leq \frac{2}{\vare\log\log x},
\end{equation}
so
\begin{equation}
\frac{1}{j^j} \leq \left(\frac{2}{\vare}\right)^j\frac{1}{(\log\log x)^j},
\end{equation}
and this gives
\begin{equation}
\frac{1}{j^j} \leq \frac{(\log x)^{\vare\log(2/\vare)}}{(\log\log x)^j}.
\end{equation}
Hence
\begin{equation}
\frac{K^j}{j^j}<\frac{(\log x)^{\vare\log(K)+\vare\log(2/\vare)}}{(\log\log x)^j}.
\end{equation}
Now as $ j  \geq  \frac{1}{2}\vare\log\log x$ one has
\begin{eqnarray}
\frac{1}{(\log\log x)^j} & \leq & \frac{1}{(\log\log x)^{(\vare/2)\log\log x}}\\
                   & \leq & \frac{1}{(\log x)^{(\vare/2)\log\llx}}.
\end{eqnarray}
So finally one sees that this gives
\begin{equation}
\frac{K^j}{j^j} \leq  \frac{(\log x)^{\vare\log(K)+\vare\log(2/\vare)}}{(\log x)^{(\vare/2)\log\llx}},
\end{equation}
which is certainly $\leq \frac{1}{(\log x)^{1+\delta}}$ for any $\delta>0$ for sufficiently large $x$.
Hence one deduces that for all $x$ sufficiently large one has
\begin{equation}\label{eq:cdxo-estmiate}
|\dxo| \leq c_{\arabic{constantCounter}\stepcounter{constantCounter}} x \frac{K^j}{j^j} \leq c_{\arabic{constantCounter}\stepcounter{constantCounter}} \frac{x}{(\log x)^{1+\delta}}
\end{equation}
for any positive $\delta$.

Now let me estimate $|\dxt|$. Each $n\in {\mathcal{D}}_2(x)$ is of the form $n=n_1\cdots n_k$ with $\omega(n_i)\geq (1-\vare)\llx$ for all $1\leq i\leq k$. Hence one has
\begin{equation}
\Omega(n) \geq (k-k\vare)\log\log x.
\end{equation}
As $n\in {\mathcal{D}}_2(x)$, one has $n=m_1\cdot m_2$ with $\omega(m_2)\leq \vare\log\log x$. Therefore each pair $n_i,n_j$ can have at most $\vare\log\log x$ common prime factors (as such factors  contribute to $m_2$). Therefore it follows that each $n\in{\mathcal{D}}_2(x)$ has at least
\begin{equation}
(k-k\vare-\vare k(k-1)/2)\log\log x = (k-\vare k(k+1)/2)\log\log x
\end{equation}
distinct prime factors. In other words
\begin{equation}
\omega(n)\geq (k-\vare k(k+1)/2)\log\log x.
\end{equation}
If $\vare <\frac{2}{20{(k^2+k)}}$ then the number on the right is positive. Thus each such $n$ has at least $2^{(k-\vare k(k+1)/2)\log\log x}$  distinct divisors.

On the other hand as
\begin{equation}
\sum_{n\leq x}d(n) \leqc_{\arabic{constantCounter}\stepcounter{constantCounter}} x\log x
\end{equation}
one has
\begin{equation}
2^{(k-\vare k(k+1)/2)\log\log x}\sum_{n\in{\mathcal{D}}_2(x)} 1 \leq \sum_{n\leq x}d(n) \leqc_{\arabic{constantCounter}\stepcounter{constantCounter}} x\log x
\end{equation}
which gives
\begin{eqnarray}
|\dxt| & \leq & c_{\arabic{constantCounter}\stepcounter{constantCounter}} \frac{x\log x}{2^{(k-\vare k(k+1)/2)\log\log x}} \\
      & \leq & c_{\arabic{constantCounter}\stepcounter{constantCounter}} \frac{x\log x}{e^{\log(2)(k-\vare k(k+1)/2)\log\log x}} \\
      & \leq & c_{\arabic{constantCounter}\stepcounter{constantCounter}} \frac{x\log x}{(\log x)^{\log(2)(k-\vare k(k+1)/2)}} \\
      & \leq & c_{\arabic{constantCounter}\stepcounter{constantCounter}} \frac{x}{(\log x)^{\log(2)(k-\vare k(k+1)/2)-1}}.
\end{eqnarray}
This will be $\ll \frac{x}{(\log x)^{1+\delta}}$ with $\delta>0$ if
\begin{equation}
\log(2)(k-\vare k(k+1)/2)-1 > 1.
\end{equation}
Now one choose $\vare$ as follows:
\begin{equation}
\vare < \frac{2}{20(k^2+k)}
\end{equation}
then it follows that
\begin{equation}
\log(2)\left(k-\vare k(k+1)/2\right)-1 > \log(2)\left(k-\frac{1}{20}\right)-1>1
\end{equation}
holds for all $k\geq 3$.  Hence in this case one has
\begin{equation}\label{eq:cdxt-estimate}
|\dxt|\leqc_{\arabic{constantCounter}}(\vare)\stepcounter{constantCounter} \frac{x}{(\log x)^{1+\delta(k)}}
\end{equation}
with $\delta(k)=\log(2)(k-\vare k(k+1)/2)-2>0$ for all $k\geq 3$. This is the only place where $k\geq 3$ is used in this proof. Unfortunately as $2\log(2)=1.38\cdots$ no choice of $\vare>0$ however small can give $\delta(2)>0$. Now one has $\csetx{D}=|\dxo|+|\dxt|$. Putting the estimates \eqref{eq:cdxo-estmiate} and \eqref{eq:cdxt-estimate} together one obtains that there is a $\delta>0$ (for $k\geq 3$) such that for all sufficiently large $x$ one has
\begin{equation}
\csetx{D} \leq c_{\arabic{constantCounter}}(\vare)\stepcounter{constantCounter}\frac{x}{(\log x)^{1+\delta}}.
\end{equation}
This proves the theorem.
\end{proof}

\section{Proof of Theorem~\ref{th:main}}
\begin{lemma}\label{le:basic-count1}
Let $E/\Q$ be an elliptic curve. If $E$ does not have complex multiplication assume GRH holds for $E$.\ Let $\vare$ be a fixed, sufficiently small positive real number.
Let $\setx{A}$ denote the number of primes $p\leq x$ such that
\begin{equation}
\omega(\np{p})<\lo
\end{equation} \
Then there is a positive constant $c_i(\vare)$ such that for all $x\geq x_0$ one has
$$
\csetx{A}\leqc_{\arabic{constantCounter}}(\vare)\stepcounter{constantCounter} \frac{\pi(x)}{(\log\log x)^2}.
$$
\end{lemma}
\begin{proof}
It is sufficient to prove that the number of $p\leq x$ such that
\begin{equation} \omega(\np{p})<(1-\vare)\llx\  {\rm or\ } \omega(\np{p})>(1+\vare)\llx\end{equation}
is
$\ll\frac{\pi(x)}{(\log\log x)^2}$ as the number of $p\leq x$ counted in the assertion of this Lemma is certainly less than the number of primes $p\leq x$   with this  property.
For $p\leq x$  satisfying the above property one has either $$\omega(\np{p})-\log\log x < -\vare\log\log x$$ i.e.
$$-(\omega(\np{p})-\log\log x)>\vare\log\log x$$ or $$\omega(\np{p})-\log\log x > \vare\log\log x$$ and hence at any rate
$$(\omega(\np{p})-\log\log x)^4 > \vare(\log\log x)^4.$$ Hence
\begin{equation}
\sum_{p\in \setx{A}}\vare^4(\log\log x)^4<\sump x(\omega(\np{p})-\log\log x)^4 \leq c_{\arabic{constantCounter}\stepcounter{constantCounter}}\pi(x)(\log\log x)^2,
\end{equation}
by Lemma~\ref{le:fourth-moment}. Thus the assertion follows.
\end{proof}

\begin{lemma}\label{le:convergence}
Let $E/\Q$ be an elliptic curve. If $E$ does not have complex multiplication assume GRH holds for $E$.\ The series $\sum_{p\in\setx{A}}\frac{1}{p}$ converges.
\end{lemma}
\begin{proof}
This follows from the previous lemma by Abel summation formula and the tautological bound ${\mathcal A}(t)\leq \pi(t)$ for all $t\geq 1$ as the integral $$\int_5^x\frac{dt}{t\log(t)(\log\log(t))^2}$$ converges as $x\to\infty$.
\end{proof}

\begin{lemma}\label{le:basic-count2}
Let $E/\Q$ be an elliptic curve. If $E$ does not have complex multiplication assume GRH holds for $E$.\ Let $\vare$ be a fixed, sufficiently small positive real number.
Let $\setx{C}$ denote the number of primes $p\leq x$ such that
\begin{equation}
 \omega(\np{p}) \geq (1-\vare)\llx.
\end{equation} \
Then there is a positive constant $c_i(\vare)$ such that for all $x\geq x_0$ one has
$$
\csetx{C}\geq \pi(x)-c_{\arabic{constantCounter}}(\vare)\stepcounter{constantCounter}\frac{\pi(x)}{(\log\log x)^2}>c_{\arabic{constantCounter}\stepcounter{constantCounter}}(\vare)\pi(x).
$$
\end{lemma}
\begin{proof}
By Lemma~\ref{le:basic-count1} the number of primes $p\leq x$ such that $ \omega(\np{p}) \geq (1-\vare)\llx$ does not hold is $o(\pi(x))$ hence the assertion follows.
\end{proof}

\begin{lemma}\label{le:series-lemma}
Let $E/\Q$ be an elliptic curve. If $E$ does not have complex multiplication assume GRH holds for $E$.\ Suppose $\vare$ is sufficiently small positive real number. Let $k\geq 1$ be an integer and $0<a<b<1$ be fixed positive real numbers such that $\log(b)-\log(a)>1$. Let
$$
\sum\nolimits_1= \sum_{x^a\leq p_1,\ldots,p_{k} <x^b} \frac{1}{p_1\cdots p_{k}},
$$
where the sum is over all the pairwise distinct primes $p_1,\ldots, p_k$ such that $x^a<p_i<x^b$ and  $\omega(\np{p_i})>(1-\vare)\llx$ for all $1\leq i\leq k$. Then there is a positive constant (depending on $a,b$) such that
$$\sum\nolimits_1 = c_{\arabic{constantCounter}\stepcounter{constantCounter}}(a,b,k)-\epsilon(x)>1,$$
where the right hand side is greater than one for all sufficiently large $x$ and $\epsilon(x)\to0$ as $x\to\infty$.
\end{lemma}
\begin{proof} For the purpose of this paper, I say that a \textit{prime $p$ is admissible if $\omega(\np{p})\geq(1-\vare)\llx$}, otherwise I say \textit{$p$ is inadmissible}.
Let me first remark that the proof shows that one may take $a,b$ satisfying
$0<b<\frac{1}{k-1}$ and $a<\frac{1}{e(k-1)}$ then one has $0<a<b<1$ and $\log(b/a)>1$. For $k=3$, in \citep{erdos36} $b=\frac{1}{4}<\frac{1}{2}$ and $a=\frac{1}{8}<\frac{1}{2e}$ so that $\log(b/a)=\log(4)>1$. Consider (the sum over all primes)
$$\sum_{x^a\leq p\leq x^b}\frac{1}{p}=\log\log(x^b)-\log\log(x^a)+o(1)=\log(b)-\log(a)+o(1)$$
by \citep[Mertens Theorem]{hardy-wright} and in particular, by the choice of $a,b$, this sum is greater than a positive constant for all $x$ sufficiently large. Now write this sum as
\begin{equation}
\sum_{x^a\leq p\leq x^b}\frac{1}{p}=\sum\nolimits_2\frac{1}{p}+\sum\nolimits_3\frac{1}{p}
\end{equation}
where $\sum\nolimits_2$ is over primes $p$ for which $\omega(\np{p})>(1-\vare)\llx$, in other words it is a sum over admissible primes, and
\begin{equation}
\sum\nolimits_3=\sum_{x^a\leq p\leq x^b,\omega(\np{p})\leq (1-\vare)\llx}\frac{1}{p}
\end{equation}
is a sum over inadmissible primes. By the convergence of the sum over inadmissible primes given by Lemma~\ref{le:convergence} one sees that if we write
\begin{equation}
\sum\nolimits_3=\sum_{x^a\leq p\leq x^b,\omega(\np{p})\leq (1-\vare)\llx}\frac{1}{p}=\epsilon(x)
\end{equation}
then $\epsilon(x)=o(1)$ as $x\to\infty$ and in particular for all sufficiently large $x$ one has $|\epsilon(x)|<\frac{1}{100}$. The choice of $\frac{1}{100}$ is arbitrary, the point being that $\epsilon(x)$ can be made smaller than any given positive real number by choosing $x$ sufficiently large. So one gets that the sum of the reciprocals of admissible primes (on the left)
\begin{equation}
	\sum\nolimits_2=\sum_{x^a\leq p\leq x^b}\frac{1}{p}-\epsilon(x)=\log(b/a)-\epsilon(x)
\end{equation}
and hence for sufficiently large $x$ the right hand side is greater than one as $\epsilon(x)$ can be made arbitrarily small by choosing $x$ sufficiently large. Thus the result is true for $k=1$.

Now before I give the general case let me illustrate the method of proof with the case $k=2$ as this also will allow me to set up notational conventions needed to do the general case. \textit{The sums over primes which occur below are sums over admissible primes.}
Since the result is true for $k=1$, one gets by squaring
\begin{equation}
\left(\sum\nolimits_2\right)^2=
(\log(b/a)-\epsilon(x))^2
\end{equation}
on the other hand by squaring the sum over admissible primes on the left one gets a sum over admissible primes
\begin{equation}
\left(\sum\nolimits_2\right)^2=\sum \frac{1}{p_1p_2}+2\sum \frac{1}{p_1^2}=(\log(b/a))^2-2\log(b/a)\epsilon(x)+\epsilon^2(x).
\end{equation}
Now the sum $\sum_p\frac{1}{p^2}$, taken over all primes (and not just admissible ones) is convergent and hence the sum $\sum_{x^a<p_1<x^b}\frac{1}{p_1^2}$ can be made arbitrarily small by choosing $x$ sufficiently large, in other words the sum is $\epsilon(x)$ for some function $\epsilon(x)$.  Since $\log(b/a)$ is fixed one may also write middle term in the above expression as $\epsilon(x)$ for yet another $\epsilon(x)$. Since there are finitely many such $\epsilon(x)$ it is possible to arrange $x$ so large that all of them be less than a given positive real number, similarly $\epsilon^2(x)$ is also $\epsilon(x)$. \textit{Note that I am systematically conflating all the $\epsilon(x)$--a diligent reader may easily workout the full notational genuflections required to separate them.}
Thus one gets from all this that
\begin{equation}
\sum_{x^a<p_1,p_2<x^b}\frac{1}{p_1p_2}=c_{\arabic{constantCounter}\stepcounter{constantCounter}}-\epsilon(x)>1
\end{equation}
for all sufficiently large $x$.

Now let me return to the general case which is proved by induction on $k$. Suppose the result is true for such a sum in the statement of the lemma and all natural numbers $\leq k$. Then I show that the result is also true for $k+1$. Since the result is true for $k$ one has
\begin{equation}
\sum\nolimits_1= \sum_{x^a\leq p_1,\ldots,p_{k} <x^b} \frac{1}{p_1\cdots p_{k}}=c_{\arabic{constantCounter}\stepcounter{constantCounter}}(a,b,k)-\epsilon(x),
\end{equation}
Then consider the product
\begin{equation} \label{eq:product-sum}
\left(\sum_{x^a\leq p_1,\ldots,p_{k} <x^b} \frac{1}{p_1\cdots p_{k}}\right)
\left(\sum_{x^a\leq p_1,\ldots,p_{k} <x^b} \frac{1}{p}\right)
\end{equation}
where both the sums are over all (pairwise) distinct admissible primes in the asserted range. Then by induction hypothesis (applied to both the sums) one has
\begin{equation}
\left(\sum_{x^a\leq p_1,\ldots,p_{k} <x^b} \frac{1}{p_1\cdots p_{k}}\right)
\left(\sum_{x^a\leq p_1,\ldots,p_{k} <x^b} \frac{1}{p}\right) = (c_{\arabic{constantCounter}\stepcounter{constantCounter}}-\epsilon(x))(c_{\arabic{constantCounter}\stepcounter{constantCounter}}-\epsilon(x))=c_{\arabic{constantCounter}\stepcounter{constantCounter}}-\epsilon(x),
\end{equation}
for some $\epsilon(x)$  and some constant depending on $a,b$ and in particular the right hand side is greater than one.

Now multiplying out the product one gets one term of the form
\begin{equation}
\sum\nolimits_1=\left(\sum_{x^a\leq p_1,\ldots,p_{k+1} <x^b} \frac{1}{p_1\cdots p_{k+1}}\right)
\end{equation}
a sum over pairwise distinct admissible primes in the asserted range
and the remaining terms are of the form
\begin{equation}
\sum\nolimits_1=\left(\sum_{x^a\leq p_1,\ldots,p_{k+1} <x^b} \frac{1}{p_{i_1}^{\ell_1}\cdots p_{i_j}^{\ell_j}}\right),
\end{equation}
where  the primes are admissible and $\ell_1+\cdots +\ell_j=k+1$ and $j\leq k$ and at least one of the exponents $\ell_1,\ldots,\ell_j$ is $\geq 2$, and which arise out of coincidences amongst the primes when one takes the product in \eqref{eq:product-sum}.
So one has to deal with these sort of sums.

I claim that these second sort of sums are always $\epsilon(x)$ for some $\epsilon(x)\to0$ as $x\to\infty$. Suppose this assertion is true for the moment. Then one gets

\begin{eqnarray}
\left(\sum_{x^a\leq p_1,\ldots,p_{k} <x^b} \frac{1}{p_1\cdots p_{k}}\right)
\left(\sum_{x^a\leq p_1,\ldots,p_{k} <x^b} \frac{1}{p}\right) & = & c_{\arabic{constantCounter}\stepcounter{constantCounter}}-\epsilon(x)\\
&=& \left(\sum_{x^a\leq p_1,\ldots,p_{k+1} <x^b} \frac{1}{p_1\cdots p_{k+1}}\right) + \\
&& \sum_{j,\ell_1,\ldots,\ell_j} \left(\sum_{x^a\leq p_{i_1},\ldots,p_{i_j} <x^b} \frac{1}{p_{i_1}^{\ell_1}\cdots p_{i_j}^{\ell_j}}\right)\\
&=& \left(\sum_{x^a\leq p_1,\ldots,p_{k+1} <x^b} \frac{1}{p_1\cdots p_{k+1}}\right)+\epsilon(x),
\end{eqnarray}
and hence rearranging terms
\begin{equation}
\left(\sum_{x^a\leq p_1,\ldots,p_{k+1} <x^b} \frac{1}{p_1\cdots p_{k+1}}\right)=c_{\arabic{constantCounter}\stepcounter{constantCounter}} -\epsilon(x)>1.
\end{equation}

Now let me return to prove that the sums where one has coincidences between some of the primes are small. I claim that for any $k$ and $\ell_1+\cdots+\ell_j=k$ and at least one $\ell_i\geq 2$,
\begin{equation}\label{le:subclaim}
\sum\nolimits_1=\left(\sum_{x^a\leq p_1,\ldots,p_{k+1} <x^b} \frac{1}{p_{i_1}^{\ell_1}\cdots p_{i_j}^{\ell_j}}\right)=\epsilon(x),
\end{equation}
This is again proved by induction on $j$. Suppose $j=1$. Then one has $\ell_1\geq 2$ as at least one of the exponents is required to be greater than one. Then the assertion follows as $\sum\frac{1}{p^2}$ converges. Hence assume the result is true when the number of primes factors is less than or equal to some $j$ and any $\ell_1,\ldots,\ell_j$ with one of them being greater that two. Then I prove the result is true for $j+1$ and any $\ell_1,\ldots,\ell_{j+1}$ with at least one of them being greater than two.

If at least one of $\ell_i\geq 2$ for some $1\leq i\leq j$ then the assertion follows by induction as $\sum_{x^a<p\leq x^b}\frac{1}{p^{\ell_{j+1}}}$ is bounded if $\ell_{j+1}=1$ and $o(1)$ if $\ell_{j+1}\geq 2$.

So now assume that $\ell_{i}=1$ for all $1\leq i\leq j$ and $\ell_{j+1}\geq 2$. Then the sum in consideration  occurs in the product
\begin{equation}
\left(\sum_{x^a<p_i\leq x^b}\frac{1}{p_1\cdots p_j}\right)\left(\sum_{p}\frac{1}{p^{\ell_{j+1}}}\right)
\end{equation}
and the first is a summand of $\left(\sum_{x^a<p\leq x^b}\frac{1}{p}\right)^j$ and hence is bounded (again one uses induction hypothesis to deal with series which coincidences as the number of factors will be $<j+1$) and the second sum is $o(1)$ as $\sum_p\frac{1}{p^{\ell_{j+1}}}$ converges as $\ell_{j+1}\geq 2$. Thus in all cases the claim~\ref{le:subclaim} is proved and hence Lemma~\ref{le:series-lemma} is proved.
\end{proof}

\begin{proof}[{Proof of Theorem~\ref{th:main}}]
Let $E/\Q$ be an elliptic curve. If $E$ does not have complex multiplication assume GRH holds for $E$.\ Let $\vare$ be a fixed, sufficiently small positive real number. Let $k\geq 1$ be an integer.  Let
\begin{equation}
\bx{k} =\{  p_1\cdot p_2\cdots p_k\leq x:   p_i\neq p_j\ {\rm if}\ i\neq j\  \ {\rm and}\ p_i\ {\rm admissible\ for\ all}\    1\leq i\leq k\}.
\end{equation}
Then I claim that it is sufficient to prove that for any $k\geq 1$ one has
\begin{equation}\label{eq:erdos-bound}
\cbx{k}\geq c_{\arabic{constantCounter}}(\vare)\stepcounter{constantCounter}\frac{x}{\log x}.
\end{equation}
Note that this claim does not use the assumption in Theorem~\ref{th:main} that $k\geq 3$. This assumption enters the proof only through the use of Theorem~\ref{th:erdos-product} which needs $k\geq 3$.

Suppose for the moment that \eqref{eq:erdos-bound} has been established. Then one can complete the proof of the theorem as follows. On one hand one has $\cbx{k}\geq c_{\arabic{constantCounter}}(\vare)\stepcounter{constantCounter} \frac{x}{\log x}$ on the other hand by Theorem~\ref{th:erdos-product} one has $\csetx{D}\leq c_{\arabic{constantCounter}}(\vare)\stepcounter{constantCounter}\frac{x}{(\log x)^{1+\delta}}$. So if one considers the mapping $\bx{k}\to \setx{D}$ given by $p_1\cdots p_k\in \bx{k} \mapsto \np{p_1}\cdots \np{p_k}$, then this map cannot have all fibers of bounded cardinality and the number $G_k(n)$ is precisely the cardinality of the fiber over $n$ under this mapping.

So let me now prove the bound asserted in \eqref{eq:erdos-bound}. If $k=1$ then this follows from Lemma~\ref{le:basic-count1} and Lemma~\ref{le:basic-count2}. So assume $k\geq 2$. If $m=p_1\cdots p_k\leq x$ then $p_k=\frac{m}{p_1\cdots p_{k-1}}\leq \frac{x}{p_1\cdots p_{k-1}}$ and hence $p_k\in \bxx{1}{x/(p_1\cdots p_{k-1})}$. Conversely if $p_1,\ldots,p_{k-1}$ are admissible primes and $p_1\cdots p_{k-1}\leq x$ and $p_k\leq \frac{x}{p_1\cdots p_{k-1}}$ is an admissible prime then $p_1\cdots p_k\leq x$ and is a product of admisible primes and hence is in $\bx{k}$. Moreover any of the $k!$ permutations of the factors $p_1,\ldots, p_k$ give the same $m=p_1\cdots p_k$. So one has
\begin{equation}
k! \cbx{k} \geq \sum \cbxx{1}{x/(p_1\cdots p_{k-1})},
\end{equation}
where the sum is over all the pairwise distinct $p_1,\ldots,p_{k-1}$ and excludes the terms corresponding to the coincidences $p_k=p_i$ for any of the $1\leq i\leq k-1$. Since one wants a lower bound, it is in fact enough to replace this sum by one of the same sort except where the primes $p_i$ are in some fixed range $x^a\leq p_i <x^b$ for $0<a<b<1$, and which satisfy the requirements of Lemma~\ref{le:series-lemma}.
Thus
\begin{equation}
k! \cbx{k} \geq \sum_{x^a\leq p_1,\ldots,p_{k-1} <x^b} \cbxx{1}{x/(p_1\cdots p_{k-1})},
\end{equation}

If one chooses $b$  so that
\begin{equation}\label{eq:b-bound}
1-b(k-1)>0
\end{equation}
 then $x/(p_1\cdots p_{k-1})\geq x/x^{b(k-1)}\gg 1$
and one can apply Lemma~\ref{le:basic-count2} to get
\begin{equation}
k! \cbx{k} \geq \sum_{x^a\leq p_1,\ldots,p_{k-1} <x^b} \cbxx{1}{x/(p_1\cdots p_{k-1})}\geq c_{\arabic{constantCounter}}(\vare)\stepcounter{constantCounter} \sum_{x^a\leq p_1,\ldots,p_{k-1} <x^b} \pi(x/(p_1\cdots p_{k-1})).
\end{equation}
Hence
\begin{equation}
k! \cbx{k} \geq c_{\arabic{constantCounter}}(\vare)\stepcounter{constantCounter}\sum_{x^a\leq p_1,\ldots,p_{k-1} <x^b} \frac{x/(p_1\cdots p_{k-1})}{\log(x/(p_1\cdots p_{k-1}))}.
\end{equation}
Since $\frac{1}{\log(x/(p_1\cdots p_{k-1}))}\geq \frac{1}{\log x}$ this gives
\begin{equation}
k! \cbx{k} \geq c_{\arabic{constantCounter}}(\vare)\stepcounter{constantCounter}\frac{x}{\log x}\left(\sum_{x^a\leq p_1,\ldots,p_{k-1} <x^b} \frac{1}{p_1\cdots p_{k-1}}\right).
\end{equation}
Now the proposition follows from Lemma~\ref{le:series-lemma} provided that in addition to the condition \eqref{eq:b-bound} one choose $a$ so that the conditions of Lemma \ref{le:series-lemma} are met. Thus on choosing $0<b<\frac{1}{k-1}$ and $0<a<\frac{1}{e(k-1)}$ one gets the assertion of  Theorem~\ref{th:main}.
\end{proof}

\section{The case $k=1$}\label{re:conj}
Now let us consider the case $k=1$. This is a case which does not occur in the classical context considered in \citep{erdos36}. Let $G_1(E,n)$ be the number of solutions (in primes $p$) to the equation $N_p(E)=n$. For example for the curve $E:y^2+y=x^3-x$ one has $$\np{1009}=\np{1063}=1057.$$ So $G_1(E,1057)\geq 2$ and a simple calculation using \eqref{eq:hw-bound} and a computer search shows that equality holds so $G_1(1057)=2$.

Here is a sample computation for another elliptic curve (chosen for no particular reason)
$$E:y^2 + 3 y = x^{3} -  x + 2$$ which tabulates a small list of solutions to the equation
$G_1(n)=3$ (in other words $n$ for which there are there primes $p_1,p_2,p_3$ with $\np{p}=n$):
\begin{center}
\begin{tabular}{|c|c|c|c|}
	\hline \rule[-2ex]{0pt}{5.5ex} $n$ & $p_1$ & $p_2$ & $p_3$ \\

	\hline 624 & 593 & 661 & 619 \\
	\hline 6495 & 6337 & 6449 & 6389 \\
	\hline 7440 & 7369 & 7523 & 7487 \\
	\hline 8568 & 8563 & 8423 & 8527 \\
	\hline 11422 & 11299 & 11617 & 11519 \\
	\hline 12312 & 12161 & 12421 & 12391 \\
	\hline 12672 & 12721 & 12791 & 12619 \\
	\hline 32022 & 31873 & 31699 & 32213 \\
	\hline 34240 & 34603 & 34217 & 34327 \\
	\hline 37464 & 37693 & 37571 & 37517 \\
	\hline
\end{tabular}
\end{center}

For the same curve $E:y^2 + 3 y = x^{3} -  x + 2$  here is an extract from the table of values of $G_1(n)\geq 2$:
\begin{center}

\begin{tabular}{|c|c|}
	\hline $n$ & $G_1(n)$ \\
	\hline
	\hline 10262 & 2\\
	\hline 10494 & 2\\
	\hline 10630 & 2\\
	\hline 10697 & 2\\
	\hline 10704 & 2\\
	\hline 11072 & 2\\
	\hline 11100 & 2\\
	\hline 11168 & 2\\
	\hline 11276 & 2\\
	\hline 11422 & 3\\
	\hline 11441 & 2\\
	\hline
\end{tabular}
\end{center}
I conjecture that for any elliptic curve $E$, one has $\varlimsup_{n\to\infty} G_1(E,n)=\infty$. In other words I conjecture that for any elliptic curve $E/{\mathbb Q}$, there exists integers $n$ (going to infinity) and arbitrarily large strings of primes $p_1<p_2<\cdots <p_m$ such that $$\np{p_1}=\np{p_2}=\cdots=\np{p_m}=n.$$  But this seems substantially more difficult at this point (even on GRH), and of course it would even more astounding if this turns out to be false. It seems natural to call such primes an \textit{elliptic progression progression of primes given by $E$} or simply \textit{primes in an elliptic progression} if the curve is unambiguously defined. At the very least I expect that $G_1(E,n)\geq 2$ infinitely often \textit{i.e.,} there are infinitely many pairs of primes in an elliptic progression. Let me summarize these conjectures:

\begin{conj}\label{con:k-one-case}
Let $E/{\mathbb Q}$ be an elliptic curve and for any integer $n\geq1$ let $G_1(n)$ be the number of good primes $p$ such that $\np{p}=n$. Then
\begin{enumerate}
\item $\varlimsup_{n\to\infty}G_1(n)=\infty$.
\item $\varlimsup_{n\to\infty}G_1(n)\geq 2$
\item There exists infinitely many $n$ such that $G_1(n)=2$.
\end{enumerate}
\end{conj}
Clearly \ref{con:k-one-case}(1)$\implies$\ref{con:k-one-case}(2) and \ref{con:k-one-case}(3)$\implies$\ref{con:k-one-case}(2). Note that \ref{con:k-one-case}(3) says that for any elliptic curve $E/{\mathbb Q}$, there are infinitely many distinct good primes $p,q$ such that $\np{p}=\np{q}$.
\bibliographystyle{plain}