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.

A339567 Numbers k such that A007088(k) == 1 (mod k).

This page as a plain text file.
%I A339567 #22 Feb 10 2021 13:53:30
%S A339567 1,5,15,25,55,91,137,525,625,925,3967,5995,7625,10767,25087,57225,
%T A339567 68817,565027,591415,2515825,2757625,4162019,5276309,96689255,
%U A339567 115686005,133890625,242899421,492029715,588620625,1839399055,7786281065,11231388063,17251448809,71050380625
%N A339567 Numbers k such that A007088(k) == 1 (mod k).
%C A339567 All terms are odd.
%e A339567 a(3) = 15 is a term because 15 = 1111_2 and 1111 == 1 (mod 15).
%p A339567 filter:= t -> convert(t,binary) mod t = 1: filter(1):= true:
%p A339567 select(filter, [seq(i,i=1..10^7,2)]);
%t A339567 Block[{a = {1}, k}, Do[If[Mod[FromDigits@ IntegerDigits[i, 2], i] == 1, AppendTo[a, i]], {i, 2, 10^7}]; a] (* _Michael De Vlieger_, Dec 12 2020 *)
%o A339567 (PARI) isok(n) = Mod(fromdigits(binary(n)), n) == 1;
%o A339567 forstep(k=1, 10^7, 2, if(isok(k), print1(k, ", "))); \\ _Daniel Suteu_, Dec 12 2020
%Y A339567 Cf. A007088, A339566.
%K A339567 nonn,base
%O A339567 1,2
%A A339567 _Robert Israel_, Dec 09 2020
%E A339567 a(30)-a(34) from _Daniel Suteu_, Dec 12 2020