martijn 4

d = [ ];
for i = 1 : length( a )
    d = [ d ones( 1, b( i ) ) * a( i ) ];
end
e = perms( d );
f = find( all( diff( e' ) ) );
c = e( f( 1 ), : );