OFFSET
1,1
COMMENTS
The corresponding n are a subsequence of A001333; see example.
LINKS
Paolo Xausa, Table of n, a(n) for n = 1..19
Joerg Arndt, Arctan relations for Pi (the author's starting point for this sequence).
EXAMPLE
The first few terms correspond to
7^2 + 1 = 2 * 5^2 = 2 * 25,
41^2 + 1 = 2 * 29^2 = 2 * 841,
239^2 + 1 = 2 * 13^4 = 2 * 28561,
8119^2 + 1 = 2 * 5741^2 = 2 * 32959081,
47321^2 + 1 = 2 * 33461^2 = 2 * 1119638521,
63018038201^2+1 = 2 * 44560482149^2 = 2 * 1985636569351347658201.
MATHEMATICA
With[{r=Range[100]}, Select[((ChebyshevT[r, I]/I^r)^2+1)/2, !PrimeQ[#]&&PrimePowerQ[#]&]] (* Paolo Xausa, Nov 13 2023, after Joerg Arndt *)
PROG
(PARI) forstep(n=1, 10^9, 2, t=(n^2+1)/2; if( !isprime(t) && isprimepower(t), print1(t, ", ")));
(PARI) /* much more efficient: */
{b(n) = polchebyshev(n, 1, I) / I^n}
for(n=1, 10^3, t=(b(n)^2+1)/2; if(!isprime(t)&&isprimepower(t), print1(t, ", ")));
CROSSREFS
KEYWORD
nonn
AUTHOR
Joerg Arndt, Dec 11 2014
STATUS
approved