function [h,hh] = genstrict1(P,m,n)
%P = P + 1e-12;
%P = P/sum(P);
hh = [];
for i=2:2^n-1
  setstr = dec2bin(i-1,n);
  ii = find(setstr=='1');
  hh(end+1) = strictii(P,m,n,ii);
end
h = max(hh);
return

function h = strictii(P,m,n,ii)
jj = setdiff(1:n,ii);
PP = [];

lii = length(ii);

dims = repmat([m],[1 lii]);
for xi = 1:m^lii
  X = ind2coord(xi,dims);
  [pX,Xinds] = getPX(P,m,n,X,ii);
  Pii = P(Xinds)/(pX+realmin);
  PP = [PP Pii];
end
h = strict(PP);
return


function h = strict(A)
[n,n] = size(A);
h = 0;
for i=1:n
  for j=i+1:n
    x = A(:,i);
    y = A(:,j);
    d = .5*sum(abs(x-y));
    if d>h
      h=d;
    end
  end
end
return