function [phi1,K1,maxs] = boundsumdefmu(K,m,n)

maxs = 0;
dims = repmat([m],[1 n]);

K1 = getK1(K,m,n);
n1 = n-1;
[Ainf,Binf] = lipab(m,n1);
LB =    zeros(m^n1,1);
UB = n1+zeros(m^n1,1);

[x,fval] = linprog(-K1,Ainf,Binf,[],[],LB,UB);

phi1 = round(x-min(x))';


%K0 = getK0(K,m,n);
%[Ainf,Binf] = lipab(m,1);
%[x,fval] = linprog(-K0,Ainf,Binf,[],[],[],[]);
%
%mu = round(x-min(x))';

mu = (K>0);

s2 = sum(K.*mu);

maxs = phi1*K1 + s2;

return