function [ a, more, m, m2 ] = ksub_next ( n, k, a, more, m, m2 )
%*****************************************************************************80
%
%% ksub_next() generates the subsets of size K from a set of size N, one at a time.
%
% Licensing:
%
% This code is distributed under the MIT license.
%
% Modified:
%
% 28 May 2015
%
% Author:
%
% Original FORTRAN77 version by Albert Nijenhuis, Herbert Wilf.
% This version by John Burkardt.
%
% Reference:
%
% Albert Nijenhuis, Herbert Wilf,
% Combinatorial Algorithms,
% Academic Press, 1978, second edition,
% ISBN 0-12-519260-6.
%
% Input:
%
% integer N, the size of the set from which subsets are drawn.
%
% integer K, the desired size of the subsets. K must
% be between 0 and N.
%
% integer A(K). A(I) is the I-th element of the
% subset. Thus A(I) will be an integer between 1 and N.
% Note that the routine will return the values in A
% in sorted order: 1 <= A(1) < A(2) < ... < A(K) <= N
%
% logical MORE. Set MORE = FALSE before first call
% for a new sequence of subsets. It then is set and remains
% TRUE as long as the subset computed on this call is not the
% final one. When the final subset is computed, MORE is set to
% FALSE as a signal that the computation is done.
%
% integer M, M2, two variables used by this
% procedure for bookkeeping. The user must declare these variables,
% and the output values from one call must be used as the input values
% on the next. The user should not change these values.
%
% Output:
%
% integer A(K). A(I) is the I-th element of the
% subset. Thus A(I) will be an integer between 1 and N.
% Note that the routine will return the values in A
% in sorted order: 1 <= A(1) < A(2) < ... < A(K) <= N
%
% logical MORE. Set MORE = FALSE before first call
% for a new sequence of subsets. It then is set and remains
% TRUE as long as the subset computed on this call is not the
% final one. When the final subset is computed, MORE is set to
% FALSE as a signal that the computation is done.
%
% integer M, M2, two variables used by this
% procedure for bookkeeping. The user must declare these variables,
% and the output values from one call must be used as the input values
% on the next. The user should not change these values.
%
if ( k < 0 )
fprintf ( 1, '\n' );
fprintf ( 1, 'KSUB_NEXT - Fatal error!\n' );
fprintf ( 1, ' K = %d\n', k );
fprintf ( 1, ' but 0 <= K is required!\n' );
error ( 'KSUB_NEXT - Fatal error!' );
end
if ( n < k )
fprintf ( 1, '\n' );
fprintf ( 1, 'KSUB_NEXT - Fatal error!\n' );
fprintf ( 1, ' N = %d\n', n );
fprintf ( 1, ' K = %d\n', k );
fprintf ( 1, ' but K <= N is required!\n' );
error ( 'KSUB_NEXT - Fatal error!' );
end
if ( ~ more )
m2 = 0;
m = k;
else
if ( m2 < n - m )
m = 0;
end
m = m + 1;
m2 = a(k+1-m);
end
for j = 1 : m
a(k+j-m) = m2 + j;
end
more = ( a(1) ~= ( n - k + 1 ) );
return
end