Yeah! I'm reaching a similar conclusion..now if only I could work out to use Algorithm::Loops, it'd be done in a trice :)

Here's a recursive solution to find the closest set of factors to a target (with wrapper script). I think it's solid.

use strict; 
use warnings;
use List::Util qw[ reduce ];

my @pfs = ( 2,3,3,3,7, 997 );
my $num = reduce { $a * $b } @pfs;
my $root = sqrt $num;
print "N=$num, R=$root\n";

sub closest {
  # Args are target and factor-list
  my ($target, $car, @cdr) = @_;
  @cdr or return $car;
  reduce { abs($target - $a) < abs($target - $b) ? $a : $b }
         ($car*closest($target/$car, @cdr), closest($target, @cdr), $c
+ar);
}

print closest($root, @pfs), "\n";
[download]

---

Caution: Contents may have been coded under pressure.

It seems to work well as far as I can verify it, with one exception. Like one of my other attempts, if N is prime, it produces N as the nearest factor rather than 1.

However, being recursive, it suffers in performance one you really start making it work.

Try 8585937356. I get a nearest factor of 4 with my dumb crunch code. I gave up waiting for your's :)

---

Examine what is said, not who speaks.

Silence betokens consent.

Love the truth but pardon error.

BrowserUk,
The following code demonstrates that this is not actually true:
Read more... (2 kB)
• If a subset has been seen before, it is skipped prior to determining the product of the subset
• If a subset is the start of a progression of subsets previously seen, the entire progression is skipped
• If a factor is determined to be larger than the sqrt, determining if it is the closest is skipped
Unfortunately, I do not know if this is the fastest implementation. I have several variations ranging from bitstrings to avoiding push/slice. The trouble is that I have found an apparent bug in Math::Pari so I am not sure which solution should be fastest. The bug only seems to manifest itself when I add the sqrt() code. To see the bug for yourself using 24777695232, add print "$_\n" for @factor; right after the array is created. To see it go away, comment out the $sqrt definition and change the return statement to return $f.

Cheers - L~R

I got an entirely different error when I did a sqrt($N) before factoring:

  DB<1> n
main::(factsqrt.pl:5):    my $N = $_;
  DB<1> n
main::(factsqrt.pl:6):    my $sqrt2 = sqrt($N);
  DB<1> n
main::(factsqrt.pl:7):    my $a = factorint($N);
  DB<1> n
  DB<1> p $a
Mat([-1,1])
[download]

However, this is not a bug in Math::Pari, but a side effect of perl's conversion of string to numbers.

The number "24777695232" first comes in as a string from the command line. Applying the function "sqrt" to it converts it to a double. When you pass it to Math::Pari, it goes in as a double. Naturally, the module fails to do the right thing. If you call a Math::Pari function first, like factorint, it automatically converts the string into a long Pari integer and operates on that.

The bug only seems to manifest itself when I add the sqrt() code. .... To see it go away, comment out the $sqrt definition and change the return statement to return $f.

Sorry Limbic~Region, I do not follow these instructions?

The only definition of $sqrt I can see is     my $sqrt = sqrt( $N );

But I do not understand what you mean by "change the return statement to return $f."?

---

Examine what is said, not who speaks.

Silence betokens consent.

Love the truth but pardon error.

Here's an Algorithm::Loops solution, then, based on lidden's nifty solution:

use strict; 
use warnings;
use List::Util qw[ reduce ];
use Algorithm::Loops 'NextPermuteNum';

sub closest{
    my $exact = sqrt(shift);
    my $best = 1;
    do {
        my $guess = 1;
        for (@_) {
        $guess *= $_;
            $best = $guess if (abs($exact - $guess) < abs($exact - $be
+st));
            last if $guess > $exact;
        }
    } while NextPermuteNum(@_);
    $best;
}

my @pfs = (2,2,3,3,5,5,7);
# Compute product of factors
our ($a, $b);
my $num = reduce { $a * $b } @pfs;
my $root = sqrt $num;
print "N=$num, R=$root\n";

print closest($num, @pfs), "\n";
[download]

---

Caution: Contents may have been coded under pressure.

• Prior to determining the product of each combination, you can abort if you have already seen this combination
• When determining the product of each combination, you can abort if the value is greater than the sqrt
• If combinations are gathered left to right, once you have determined that the product of a given combination is greater than the sqrt() all future combinations that include this combination can be skipped