Let $f(n)$ be an integer function such that $$ f(n) = \sum\limits_{i=1}^{n} \sum\limits_{j=1}^{i} \sum\limits_{k=1}^{j} [(3ijk - (ij + ik + jk)) = n]. $$
Here square bracket denotes Iverson bracket.
I conjecture that if $f(n) = 0$, then $2n+1$ is always prime.
Note that these primes are just a subsequence of all primes and the sequence is not in the OEIS.
Is there a way to prove it?