From 8426cd2ce320586f8af64e0c4d33e29b72275eea Mon Sep 17 00:00:00 2001
From: drbaggy <js5@sanger.ac.uk>
Date: Mon, 4 Oct 2021 10:56:58 +0100
Subject: smith numbers with opt!

---
 challenge-133/james-smith/perl/ch-2.pl | 49 ++++++++++++++++++++++++++++++++++
 1 file changed, 49 insertions(+)
 create mode 100644 challenge-133/james-smith/perl/ch-2.pl

diff --git a/challenge-133/james-smith/perl/ch-2.pl b/challenge-133/james-smith/perl/ch-2.pl
new file mode 100644
index 0000000000..5f7bed5d67
--- /dev/null
+++ b/challenge-133/james-smith/perl/ch-2.pl
@@ -0,0 +1,49 @@
+#!/usr/local/bin/perl
+
+use strict;
+
+use warnings;
+use feature qw(say);
+use Data::Dumper qw(Dumper);
+
+my @primes    = (2,3);
+my %ds;
+my %comp;
+my @t =  smith_numbers(@ARGV?$ARGV[0]:100);
+print "@t\n";
+
+sub smith_numbers {
+  my $C = shift;
+  my @sn = ();
+  my $n = 4;
+  while(1) {
+    my @q = prime_factors($n);
+    if( @q ) {
+      my $t = sum_digits($n);
+      $t-= sum_digits($_) foreach @q;
+      unless($t) {
+        push @sn, $n;
+        return @sn if @sn==$C;
+      }
+    }
+    $n++;
+  }
+}
+
+sub sum_digits {
+  return $ds{$_[0]} if exists $ds{$_[0]};
+  my $t = 0;
+  $t+= $_ foreach split //, $_[0];
+  return $ds{$_[0]} = $t;
+}
+
+sub prime_factors {
+  my $N = shift;
+  my @factors;
+  foreach (@primes) {
+    return @{ $comp{$N} = [ $_, exists $comp{$N/$_}? @{$comp{$N/$_}} : $N/$_ ] } unless $N % $_;
+  }
+  push @primes,$N;
+  return;
+}
+
-- 
cgit