view release on metacpan or  search on metacpan

eg/permutations  view on Meta::CPAN

#!/usr/bin/env perl
use strict;
use warnings;

# WARNING: Running this could take a very long time...

use Acme::Affinity;
use Algorithm::Combinatorics qw/ variations_with_repetition /;

my $questions = [
    { 'how messy are you' => [ 'very messy', 'average', 'very organized' ] },
    { 'do you like to be the center of attention' => [ 'no', 'maybe', 'yes' ] },
];
my $importance = {
    'irrelevant'         => 0,
    'a little important' => 1,
    'somewhat important' => 10,
    'very important'     => 50,
    'mandatory'          => 250,
};

my @responses1 = get_responses( $questions->[0], $importance ); # 45 responses
my @responses2 = get_responses( $questions->[1], $importance ); # "

my @response_pairs;
for my $r1 ( @responses1 ) {
    for my $r2 ( @responses2 ) {
        push @response_pairs, [ $r1, $r2 ];
    }
}

my $iter = variations_with_repetition( \@response_pairs, 2 ); # 4.1M variations

my %scores;

my $n = 0;

while ( my $v = $iter->next ) {
    $n++;

    my $affinity = Acme::Affinity->new(
        questions  => $questions,
        importance => $importance,
        me         => $v->[0],
        you        => $v->[1],
    );

#    print $n, '. ', $affinity->score, "\n";

    $scores{ $affinity->score }++;
}

$n = 0;

for my $score ( sort { $a <=> $b } keys %scores ) {
    print ++$n, '. ', $score, ' => ', $scores{$score}, "\n";
}

sub get_responses {
    my ( $question, $importance ) = @_;

    my @responses;

    for my $v ( values %$question ) {
        my $iter = variations_with_repetition( $v, 2 );

        while ( my $answer = $iter->next ) {
            for my $import ( sort { $importance->{$a} <=> $importance->{$b} } keys %$importance ) {
                push @responses, [ @$answer, $import ];
            }
        }
    }

    return @responses;
}

__END__
1. 0 => 3224529
2. 0.398406374501992 => 1296
3. 0.883848951807801 => 2592
4. 1.23787406856394 => 2592
5. 1.90312273167978 => 2592
6. 1.96078431372549 => 1296
7. 2.57684036015037 => 5184
8. 2.74617518190545 => 2592
9. 3.84615384615385 => 1296
10. 4.22200330920749 => 2592
11. 4.46321842677452 => 10368
12. 5.7166195047503 => 5184
13. 5.76199021246126 => 5184
14. 5.91312395989083 => 2592
15. 6.01820249894957 => 2592
16. 6.18937034281971 => 2592
17. 6.24975587367918 => 2592
18. 6.29935788878163 => 2592
19. 6.31194403097803 => 25920
20. 8.00640769025437 => 5184
21. 9.09090909090909 => 1296
22. 9.90147542976674 => 10368
23. 12.3091490979333 => 5184
24. 12.7827498141228 => 5184
25. 13.3511467458638 => 2592
26. 13.7308759095272 => 2592
27. 13.864838846795 => 2592
28. 13.8675049056307 => 10368
29. 13.9748789763252 => 2592
30. 14.0028008402801 => 25920
31. 16.6666666666667 => 5184
32. 17.9028718509858 => 5184
33. 18.6989398001691 => 2592
34. 19.2307692307692 => 2592
35. 19.4183909345154 => 2592
36. 19.5725075656075 => 2592
37. 19.6116135138184 => 25920
38. 21.320071635561 => 10368
39. 27.524094128159 => 5184
40. 28.7479787288034 => 2592
41. 28.8675134594813 => 20736
42. 29.5656197994541 => 2592