Search: id:A072037 Results 1-1 of 1 results found. %I A072037 %S A072037 4,8,9,121,343,1331,10201,14641,94249,1030301,104060401, %T A072037 900075181570009,10022212521222001,12124434743442121,12323244744232321, %U A072037 12341234943214321 %N A072037 Palindromic powers of a prime but not a prime (A025475). %e A072037 E.g. 94249=307*307 %t A072037 a = {}; Do[pp = Prime[n]^i; d = IntegerDigits[pp]; If[d == Reverse[d], a = Append[a, pp]], {n, 1, PrimePi[ Sqrt[10^21]]}, {i, 2, Floor[ Log[ Prime[n], 10^21]]}]; Sort[a] (Robert G. Wilson v) %o A072037 (PARI) {a=10^15; v=[]; m=sqrt(a); forprime(p=2,m,q=p; while((q=q*p)0,d=divrem(n,10); n=d[1]; rev=10*rev+d[2]); if(q==rev, v=concat(v,q)))); v=vecsort(v); for(j=1,matsize(v)[2],print1(v[j], ","))} %Y A072037 Cf. A025475, A002385. %Y A072037 Sequence in context: A075786 A046450 A077271 this_sequence A084093 A076703 A141220 %Y A072037 Adjacent sequences: A072034 A072035 A072036 this_sequence A072038 A072039 A072040 %K A072037 base,easy,nonn %O A072037 1,1 %A A072037 Labos E. (labos(AT)ana.sote.hu), Jun 07 2002 %E A072037 Two more term from Klaus Brockhaus (klaus-brockhaus(AT)t-online.de), Jun 07 2002 %E A072037 Four more terms from Robert G. Wilson v (rgwv(AT)rgwv.com), Oct 31 2002 Search completed in 0.001 seconds