in reply to algorithm for 'best subsets'

Ugly, off-the-top-of-my-head idea. This depends on having no space characters in the keywords - I suppose you can use \0 as a separator if you need to, but I found ' ' easier to type:
use Math::Combinatorics; my $nwordsatonce = 4; my ($k,$v); my %totals = (); local $" = ' '; # just in case while (($k,$v) = each %items) { next unless $nwordsatonce <= @$v; my @words = sort @$v; do {$totals{"@$_"} += 1;} for combine($nwordsatonce,@words); } my @comb = sort {$totals{$b} <=> $totals{$a}} keys %totals; print "Top $nwordsatonce - word combinations:\n"; do {print "$_\n";} for @comb[0..4];
I'm a bit curious where this data is coming from.
-- @/=map{[/./g]}qw/.h_nJ Xapou cets krht ele_ r_ra/; map{y/X_/\n /;print}map{pop@$_}@/for@/

Replies are listed 'Best First'.
Re^2: algorithm for 'best subsets'
by fizbin (Chaplain) on Mar 03, 2005 at 19:50 UTC
    And here's a sub version that I've actually been using for testing, using the test data code you posted. It finds 3-at-once in under 6 seconds and 4-at-once in under 11.5 seconds on my laptop, using the supplied test data.
    sub countcomb { my $nwordsatonce = shift; my ($k,$v); my %totals = (); local $" = ' '; # just in case while (($k,$v) = each %Items) { next unless $nwordsatonce <= @$v; do {$totals{"@$_"} += 1;} for combine($nwordsatonce,sort @$v); } my @comb = sort {$totals{$b} <=> $totals{$a} or $a cmp $b} keys %tot +als; my $topn = 5; my $toptot = $totals{$comb[$topn-1]}; while ($toptot <= $totals{$comb[$topn]}) {$topn++;} print "Top $nwordsatonce - word combinations: (cuttoff $toptot)\n"; do {print "$_ ($totals{$_})\n";} for @comb[0..$topn-1]; }
    Unfortunately, it does have a tendency to die of out-of-memory errors if you up either the number of keywords or the average length of a set, and tieing %totals to a db file doesn't seem to prevent it, so there must be some other sort of memory leak going on. (I suppose that it could also be the sort exploding, but working around that should be relatively straightforward)
    -- @/=map{[/./g]}qw/.h_nJ Xapou cets krht ele_ r_ra/; map{y/X_/\n /;print}map{pop@$_}@/for@/