login
Numbers > 1 in A000695 which are not in A152021.
6

%I #6 Sep 01 2019 02:36:28

%S 16,20,64,65,68,80,84,85,256,260,272,273,276,320,324,325,336,340,1024,

%T 1025,1028,1029,1040,1044,1045,1088,1092,1104,1105,1108,1280,1281,

%U 1284,1285,1296,1300,1344,1345,1348,1360,1364,1365,4096,4097,4100

%N Numbers > 1 in A000695 which are not in A152021.

%t f[n_] := FromDigits[IntegerDigits[n, 2], 4]; s = Array[f, 100, 2]; div[a_, b_] := Divisible[a, b] && a > b; n = 1; s0 = s; While[Length[s] > n, s = Select[s, ! div[#, s[[n]]] &]; n++]; Complement[s0, s] (* _Amiram Eldar_, Aug 31 2019 *)

%Y Cf. A000695, A152021.

%K nonn

%O 1,1

%A _Vladimir Shevelev_, Nov 20 2008