function hh = spcstrict(P,m,n)
% conditioned on X1=0

hh = [];
for i=2:n
  hh(end+1) = stricti(P,m,n,i);
end
return

function h = stricti(P,m,n,i)
h = 0;

P0 = getcondP(P,m,n,[i:n],[1],[1]);
P1 = getcondP(P,m,n,[i:n],[ ],[ ]);
h = .5*sum(abs(P0-P1));

return