comment on

I had an easier time thinking of this as manipulating the characteristic sequence of the subsets (i.e, a string of 0's and 1's, where the first bit is a 1 if the first factor is included in the subset, etc). If you look at the characteristic sequences given in your example, you can see that there are three rules being applied (you can print out $str to see what's going on under the covers). Anyway, this code does its manipulations on the characteristic sequence and converts it to the appropriate subset:

sub iter {
    my @factors = @_;
    my $str;

    return sub {
        if (not defined $str) {
            $str =  "1" . ("0" x $#factors);
            return map { substr($str, $_, 1) ? $factors[$_] : () }
                   0 .. $#factors;
        }
        for ($str) {
            s/0(0*)$/1$1/
              or s/11$/01/
              or s/^(.*)10(.*)$/"${1}01" . "0" x length $2/e
              or return;
        }

        return map { substr($str, $_, 1) ? $factors[$_] : () }
               0 .. $#factors;
    };
}   

my $i = iter( 2, 3, 5, 7 );
while (my @s = $i->()) {
    print "@s\n";
}

__END__
2
2 3
2 3 5
2 3 5 7
2 3 7
2 5
2 5 7
2 7
3
3 5
3 5 7
3 7
5
5 7
7
[download]

I'd be interested to see this simplified a bit. I know you mentioned in the CB that tye had an idea for a solution, and I wonder if he's able to implement it with manipulations to the subsets themselves instead of their characteristic sequences.

You also mentioned wanting to know when the "next phase" of iterations began (the horizontal lines in your example). You can figure this out by checking which substitution rule was actually applied.

Update: Here is the code with the "next phase" markers:

sub iter {
    my @factors = @_;
    my $str;
    my $break = 0;

    return sub {
        if (not defined $str) {
            $str =  "1" . ("0" x $#factors);
            return map { substr($str, $_, 1) ? $factors[$_] : () }
                   0 .. $#factors;
        }
        for ($str) {
            s/0(0*)$/1$1/ and last;

            return "BREAK" if $break = !$break;


            s/11$/01/
              or s/^(.*)10(.*)$/"${1}01" . "0" x length $2/e
              or return;
        }

        return map { substr($str, $_, 1) ? $factors[$_] : () }
               0 .. $#factors;
    };
}


my $i = iter( 2, 3, 5, 7 );

while (my @s = $i->()) {
    print "@s\n";
}

__END__
2
2 3
2 3 5
2 3 5 7
BREAK
2 3 7
BREAK
2 5
2 5 7
BREAK
2 7
BREAK
3
3 5
3 5 7
BREAK
3 7
BREAK
5
5 7
BREAK
7
BREAK
[download]

blokhead

In reply to Re: Generating powerset with progressive ordering by blokhead
in thread Generating powerset with progressive ordering by Limbic~Region

Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!

Titles consisting of a single word are discouraged, and in most cases are disallowed outright.

Read Where should I post X? if you're not absolutely sure you're posting in the right place.

Please read these before you post! —

Posts may use any of the Perl Monks Approved HTML tags:

a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr

You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)

	For:		Use:
	&		`&`
	<		`<`
	>		`>`
	[		`[`
	]		`]`

Link using PerlMonks shortcuts! What shortcuts can I use for linking?

See Writeup Formatting Tips and other pages linked from there for more info.