cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A277075 Numbers n for which A276711(n) = 1.

This page as a plain text file.
%I A277075 #13 Jun 02 2025 12:19:39
%S A277075 2,4,5,8,9,10,16,22,24,25,26,36,58,64,76,82,85,120,170,196,202,214,
%T A277075 324,328,370,412,505,562,676,706,730,799,841,1024,1204,1243,1549,1681,
%U A277075 1849,2146,2986,10404,46656,52900,112896,122500
%N A277075 Numbers n for which A276711(n) = 1.
%C A277075 Is the sequence finite?
%C A277075 a(47) > 10^7 if it exists.
%e A277075 122500 is in the sequence because 122500 = 121169 + 11^3 is the only way to write 122500 as the sum of a prime and a nonnegative perfect power.
%p A277075 N:= 10^6: # to get all terms <= N
%p A277075 Primes:= select(isprime, [2,seq(i,i=3..N,2)]):
%p A277075 Pows:= {0,1,seq(seq(b^k,k=2..floor(log[b](N))),b=2..floor(sqrt(N)))}:
%p A277075 G:= expand(add(x^p,p=Primes)*add(x^r,r=Pows)):
%p A277075 B:= [seq(coeff(G,x,i),i=1..N)]:
%p A277075 select(t -> B[t]=1, [$1..N]);
%Y A277075 Cf. A276711.
%K A277075 nonn,more
%O A277075 1,1
%A A277075 _Robert Israel_, Sep 27 2016