bindi:=proc(n) if n<1 then [] else [op(bindi(floor(n/2))), n mod 2] fi: end: idnib:=proc(a) local i:if nops(a)=0 then 0 else a[-1]+2*idnib([seq(a[i],i=1..nops(a)-1)]) fi:end: oy:=proc(n,m) local S,T,s,t,U,i: S:=bindi(n):T:=bindi(m):s:=nops(S):t:=nops(T): if s