function y = myobjm(x)
global nstate lenseq
global A C
A = reshape(x,[nstate nstate lenseq-1]);
C = zeros(nstate,1); C(2)=1;
Ef1 = E(1);
Ef0 = E([]);
Y_ac = Ef1 - Ef0;
%y = abs(Y_ac);
y = -Y_ac^2;
return

function r = F(x)
r = length(find(x~=1));
return

function Ef = E(Xi)
% E[f(X)|Xi]
global ALLSEQ
%global A
Ef = 0;
i = length(Xi);
n = size(ALLSEQ,2);
for k=1:size(ALLSEQ,1)
  x = ALLSEQ(k,:);
  p = P(x);
  if (i & (i<n))
    %pxi = A(x(i+1),x(i));
    %pXi = A(x(i+1),Xi(i));
    pxi = p_n(x(i+1),x (i),i);
    pXi = p_n(x(i+1),Xi(i),i);
    p = p / pxi * pXi;
  end
  f = F([Xi x(i+1:end)]);
  Ef = Ef + f*p;
end
return

function p = P(x)
%global A
%global C
if isempty(x)
  p = 1;
else
  %p = C(x(1));
  p = p_0(x(1));
end
for j=2:length(x)
  %p = p * A(x(j),x(j-1)); % markov case
  p = p * p_n(x(j),x(j-1),j-1); % markov case
  %p = p * C(x(j)); % independent case
end
return

function p = p_0(x)
global C
p = C(x);
return

function p = p_n(i,j,n)
% primitive cond prob
global A
if i
  p = A(i,j,n);
else
  p = 1;
end
return