%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Copyright: 2010 Integrated Sytems Laboratory, ETH Zurich
%%            http://www.iis.ee.ethz.ch/~sha3
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [s] = expand1(j)       % j in {16,17}
global qnt32 Qb;

    s = mod(...
        bin2num(qnt32,sprintf('%-1d', s1(Qb(j-16+1,:)))) + ...
        bin2num(qnt32,sprintf('%-1d', s2(Qb(j-15+1,:)))) + ...
        bin2num(qnt32,sprintf('%-1d', s3(Qb(j-14+1,:)))) + ...
        bin2num(qnt32,sprintf('%-1d', s0(Qb(j-13+1,:)))) + ...
        bin2num(qnt32,sprintf('%-1d', s1(Qb(j-12+1,:)))) + ...
        bin2num(qnt32,sprintf('%-1d', s2(Qb(j-11+1,:)))) + ...
        bin2num(qnt32,sprintf('%-1d', s3(Qb(j-10+1,:)))) + ...
        bin2num(qnt32,sprintf('%-1d', s0(Qb(j-9+1,:)))) + ...
        bin2num(qnt32,sprintf('%-1d', s1(Qb(j-8+1,:)))) + ...
        bin2num(qnt32,sprintf('%-1d', s2(Qb(j-7+1,:)))) + ...
        bin2num(qnt32,sprintf('%-1d', s3(Qb(j-6+1,:)))) + ...
        bin2num(qnt32,sprintf('%-1d', s0(Qb(j-5+1,:)))) + ...
        bin2num(qnt32,sprintf('%-1d', s1(Qb(j-4+1,:)))) + ...
        bin2num(qnt32,sprintf('%-1d', s2(Qb(j-3+1,:)))) + ...
        bin2num(qnt32,sprintf('%-1d', s3(Qb(j-2+1,:)))) + ...
        bin2num(qnt32,sprintf('%-1d', s0(Qb(j-1+1,:)))) + ...
        AddEl(j-16), 2^32);