%I A133322
%S A133322 5,13,25,41,61,113,181,313,421,613,761,841,1013,1201,1301,1741,1861,
%T A133322 2113,2381,2521,3121,3613,4513,5101,7321,8581,9661,9941,10513,12641,
%U A133322 13613,14281,14621,15313,16381,19013,19801,20201,21013,21841,23981
%N A133322 Centered square numbers that are prime powers of the form (4n+1)^k.
%t A133322 a = {}; For[i = 1, i < 140, i++, k = 1; b = (2*i*(i + 1) + 1); While[b^(1/
               k) > 4, If[PrimeQ[b^(1/k)], If[Mod[b^(1/k), 4] == 1, AppendTo[a, 
               b]]]; k++ ]]; a - Stefan Steinerberger (stefan.steinerberger(AT)gmail.com), 
               Oct 21 2007
%Y A133322 Subsequence of A001844 and A120960.
%Y A133322 Sequence in context: A096891 A001844 A099776 this_sequence A146590 A098483 
               A147205
%Y A133322 Adjacent sequences: A133319 A133320 A133321 this_sequence A133323 A133324 
               A133325
%K A133322 nonn
%O A133322 1,1
%A A133322 Mats Granvik (mgranvik(AT)abo.fi), Oct 18 2007, Oct 19 2007
%E A133322 More terms from Stefan Steinerberger (stefan.steinerberger(AT)gmail.com), 
               Oct 21 2007