A055131 - OEIS

A055131

Those composite s for which A055095[s] = 2.

15, 39, 51, 87, 111, 123, 159, 183, 219, 267, 291, 303, 327, 339, 411, 447, 471, 519, 543, 579, 591, 687, 699, 723, 771, 807, 831, 843, 879, 939, 951, 1011, 1047, 1059, 1119, 1167, 1191, 1203, 1227, 1263, 1299, 1347, 1371, 1383, 1527, 1563, 1623, 1671

(list; graph; refs; listen; history; text; internal format)

OFFSET

0,1

LINKS

Table of n, a(n) for n=0..47.

FORMULA

a(n) = 3*((4*A005098[n])+1) = 3*A002144[n] ??? (Conjecture, not yet proved)

MAPLE

find_A055095_is_2_composites := proc(upto_n) local j, a; a := []; for j from 1 to upto_n do if(-1 = (j - wt(GrayCode(qrs2bincode((2*j)+1))))) then if(not isprime((2*j)+1)) then a := [op(a), ((2*j)+1)]; fi; fi; od; RETURN(a); end;

MATHEMATICA

A005811[n_] := Length[Length /@ Split[IntegerDigits[n, 2]]];

A055094[n_] := With[{rr = Table[Mod[k^2, n], {k, 1, n-1}] // Union}, Boole[MemberQ[rr, #]] & /@ Range[n-1]] // FromDigits[#, 2]&;

A055095[1] = 0; A055095[n_] := 2*A005811[A055094[n]] - (n-1);

A055131 = Position[Array[A055095, 2000], 2] // Flatten // Select[#, CompositeQ]& (* Jean-François Alcover, Mar 06 2016 *)

CROSSREFS

Sequence in context: A272189 A336560 A176257 * A121051 A139042 A020142

Adjacent sequences: A055128 A055129 A055130 * A055132 A055133 A055134

KEYWORD

nonn

AUTHOR

Antti Karttunen, Apr 04 2000

EXTENSIONS

More terms from James A. Sellers, Apr 21 2000

STATUS

approved