Perl – 페이지 3 – blog.update.sh

프로젝트 오일러 40번

2013/09/17
Perl

소수점 뒤에 양의 정수를 차례대로 붙여 나가면 아래와 같은 무리수를 만들 수 있습니다.

0.123456789101112131415161718192021…

이 무리수의 소수점 아래 12번째 자리에는 1이 옵니다 (위에서 붉게 표시된 숫자).

소수점 아래 n번째 숫자를 d_n이라고 했을 때, 아래 식의 값은 얼마입니까?

d₁ × d₁₀ × d₁₀₀ × d₁₀₀₀ × d₁₀₀₀₀ × d₁₀₀₀₀₀ × d_1000000

perl의 주특기인 문자열 처리!

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

my $s = join("", (0..1000000));
say substr($s, 1, 1)
	*substr($s, 10, 1)
	*substr($s, 100, 1)
	*substr($s, 1000, 1)
	*substr($s, 10000, 1)
	*substr($s, 100000, 1)
	*substr($s, 1000000, 1);

#!/usr/bin/env perl

use 5.010;

use strict;

use warnings;

my $s = join("", (0..1000000));

say substr($s, 1, 1)

*substr($s, 10, 1)

*substr($s, 100, 1)

*substr($s, 1000, 1)

*substr($s, 10000, 1)

*substr($s, 100000, 1)

*substr($s, 1000000, 1);

프로젝트 오일러 39번

2013/09/17
Perl

세 변의 길이가 모두 자연수 {a, b, c}인 직각삼각형의 둘레를 p 로 둘 때, p = 120 을 만족하는 직각삼각형은 아래와 같이 세 개가 있습니다.

{20, 48, 52}, {24, 45, 51}, {30, 40, 50}

1000 이하의 둘레 p에 대해서, 직각삼각형이 가장 많이 만들어지는 p의 값은 얼마입니까?

30초가 넘게 걸린다.
식을 좀 예쁘게 정리하면 좀 빠르게 구할 수 있겠다.

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

my $r = 0;
my $c = 0;

for my $g (3..1000) {
	my $t = 0;
	for my $i (1..$g/3) {
		for my $j (1..($g-$i)/2) {
			my @a = sort {$a <=> $b} ($i ,$j, $g-$i-$j);
			$t++ if ($a[2] <= $a[1] + $a[0] and $a[2]**2 == $a[1]**2 + $a[0]**2);
		}
	}
	if ($t > $c) {
		$c = $t;
		$r = $g;
	}
}

say $r;

#!/usr/bin/env perl

use 5.010;

use strict;

use warnings;

my $r = 0;

my $c = 0;

for my $g (3..1000) {

my $t = 0;

for my $i (1..$g/3) {

for my $j (1..($g-$i)/2) {

my @a = sort {$a <=> $b} ($i ,$j, $g-$i-$j);

$t++ if ($a[2] <= $a[1] + $a[0] and $a[2]**2 == $a[1]**2 + $a[0]**2);

}

if ($t > $c) {

$c = $t;

$r = $g;

}

say $r;

프로젝트 오일러 38번

2013/09/16
Perl

숫자 192에 1, 2, 3을 각각 곱합니다.

192 × 1 = 192
192 × 2 = 384
192 × 3 = 576

곱한 결과를 모두 이어보면 192384576 이고, 이것은 1 ~ 9 팬디지털(pandigital)인 숫자입니다. 이런 과정을 편의상 ‘곱해서 이어붙이기’라고 부르기로 합니다.

같은 식으로 9와 (1, 2, 3, 4, 5)를 곱해서 이어붙이면 918273645 라는 1 ~ 9 팬디지털 숫자를 얻습니다.

어떤 정수와 (1, 2, … , n)을 곱해서 이어붙였을 때 얻을 수 있는 가장 큰 아홉자리의 1 ~ 9 팬디지털 숫자는 무엇입니까? (단 n > 1)

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

my $r = 0;
for my $i (1..9999) {
	my $s = "";
	for (my $j = 1; 1; $j++) {
		if (length($s) < 9) {
			$s .= $i*$j;
		} elsif (length($s) == 9) {
			if ("123456789" eq join("", sort split(//, $s))) { 
				if ($r < $s) {
					$r = $s;
				}
			}
			last;
		} else {
			last;
		}
	}
}

say $r;

#!/usr/bin/env perl

use 5.010;

use strict;

use warnings;

my $r = 0;

for my $i (1..9999) {

my $s = "";

for (my $j = 1; 1; $j++) {

if (length($s) < 9) {

$s .= $i*$j;

} elsif (length($s) == 9) {

if ("123456789" eq join("", sort split(//, $s))) {

if ($r < $s) {

$r = $s;

}

last;

} else {

last;

}

say $r;

프로젝트 오일러 37번

2013/09/16
Perl

소수 3797에는 왼쪽부터 자리수를 하나씩 없애거나 (3797, 797, 97, 7) 오른쪽부터 없애도 (3797, 379, 37, 3) 모두 소수가 되는 성질이 있습니다.
이런 성질을 가진 소수는 단 11개만이 존재합니다. 이것을 모두 찾아서 합을 구하세요.
(참고: 2, 3, 5, 7은 제외합니다)

그냥 하자.
안느리다.

#!/usr/bin/env perl
use 5.010;
use strict;
use warnings;
use Math::Prime::Util ':all';

my $c = 11;
my $r = 0;
my @a;
my $p = 7;
while (1) {
	$p = next_prime($p);
	my @a = my @b = split (//, $p);
	my $b = 1;
	for (0..length($p)-2) {
		shift(@a);
		pop(@b);
		if (!is_prime(join("", @a)) or !is_prime(join("", @b))) {
			$b = 0;
		}
	}
	if ($b) {
		$r += $p;
		last if (--$c == 0);
	}
}

say $r;

#!/usr/bin/env perl

use 5.010;

use strict;

use warnings;

use Math::Prime::Util ':all';

my $c = 11;

my $r = 0;

my @a;

my $p = 7;

while (1) {

$p = next_prime($p);

my @a = my @b = split (//, $p);

my $b = 1;

for (0..length($p)-2) {

shift(@a);

pop(@b);

if (!is_prime(join("", @a)) or !is_prime(join("", @b))) {

$b = 0;

}

if ($b) {

$r += $p;

last if (--$c == 0);

}

say $r;

프로젝트 오일러 36번

2013/09/16
Perl

대칭수(palindrome)인 585는 2진수로 나타내도 1001001001₂가 되어 여전히 대칭수입니다.
10진법과 2진법으로 모두 대칭수인 1,000,000 이하 숫자의 합을 구하세요.
(주의: 첫번째 자리가 0이면 대칭수가 아님)

sprintf를 이용해 2진수로 변환한다.

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

my $r = 0;
for (1..1000000) {
	if ($_ == reverse $_) {
		my $b = sprintf ("%b", $_);
		$r += $_ if ($b == reverse $b);
	}
}

say $r;

#!/usr/bin/env perl

use 5.010;

use strict;

use warnings;

my $r = 0;

for (1..1000000) {

if ($_ == reverse $_) {

my $b = sprintf ("%b", $_);

$r += $_ if ($b == reverse $b);

}

say $r;

프로젝트 오일러 35번

2013/09/14
Perl

소수 중에서 각 자리의 숫자들을 순환시켜도 여전히 소수인 것을 circular prime이라고 합니다. 예를 들어 197은 971, 719가 모두 소수이므로 여기에 해당합니다.

이런 소수는 100 밑으로 2, 3, 5, 7, 11, 13, 17, 31, 37, 71, 73, 79, 97 처럼 13개가 있습니다.

그러면 1,000,000 밑으로는 모두 몇 개나 있을까요?

뭐 별 수 있겠는가? 그냥 짜야지…
역시 소수를 가져오는데는 모듈을 사용했다.

#!/usr/bin/env perl
use 5.010;
use strict;
use warnings;
use Math::Prime::Util ':all';

my @a;
for (@{primes(1000000)}) {
	my @t = split (//, $_);
	my $b = 1;
	for (0..length($_)-2) {
		push (@t, shift(@t));
		next if ($t[0] == 0);
		my $s = join("", @t);
		if (!is_prime($s)) {
			$b = 0;
		}
	}
	push (@a, join("", @t)) if ($b);
}

say $#a+1;

#!/usr/bin/env perl

use 5.010;

use strict;

use warnings;

use Math::Prime::Util ':all';

my @a;

for (@{primes(1000000)}) {

my @t = split (//, $_);

my $b = 1;

for (0..length($_)-2) {

push (@t, shift(@t));

next if ($t[0] == 0);

my $s = join("", @t);

if (!is_prime($s)) {

$b = 0;

}

push (@a, join("", @t)) if ($b);

}

say $#a+1;

프로젝트 오일러 34번

2013/09/14
Perl

숫자 145에는 신기한 성질이 있습니다. 각 자릿수의 팩토리얼(계승)을 더하면 1! + 4! + 5! = 1 + 24 + 120 = 145 처럼 자기 자신이 됩니다.

이렇게 각 자릿수의 팩토리얼을 더하면 자기 자신이 되는 모든 수의 합을 구하세요.

단, 1! = 1 과 2! = 2 의 경우는 덧셈이 아니므로 제외합니다.

$l 에 한계값 넣고 돌린다.
오래걸린다. 약 28초…

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

my $r = 0;
my $l = 9;
my @f;
push @f, f($_) for (0..9); 
$l .= 9 while ($l < length($l)*($f[9]));

for my $i (3..$l) {
	my $a = 0;
	$a += $f[$_] for (split (//, $i));
	$r += $i if ($i == $a);
}

say $r;

sub f {
	my $r = 1;
	$r *= $_ for (1..$_[0]);
	return $r;
}

#!/usr/bin/env perl

use 5.010;

use strict;

use warnings;

my $r = 0;

my $l = 9;

my @f;

push @f, f($_) for (0..9);

$l .= 9 while ($l < length($l)*($f[9]));

for my $i (3..$l) {

my $a = 0;

$a += $f[$_] for (split (//, $i));

$r += $i if ($i == $a);

}

say $r;

sub f {

my $r = 1;

$r *= $_ for (1..$_[0]);

return $r;

}

프로젝트 오일러 33번

2013/09/14
Perl

분수 ⁴⁹/₉₈에는 재미있는 성질이 있습니다. 수학을 잘 모르는 사람이 분모와 분자에서 9를 각각 지워서 간단히 하려고 ⁴⁹/₉₈ = ⁴/₈ 처럼 계산해도 올바른 결과가 됩니다.

이에 비해 ³⁰/₅₀ = ³/₅ 같은 경우는 다소 진부한 예라고 볼 수 있습니다.

위와 같은 성질을 가지면서 ‘진부하지 않은’ 분수는, 값이 1보다 작고 분자와 분모가 2자리 정수인 경우 모두 4개가 있습니다.

이 4개의 분수를 곱해서 약분했을 때 분모는 얼마입니까?

… 그냥 풀자.
사실 손으로 풀어도 된다.

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

my $c = 1;
my $p = 1;
for my $i (1..9) {
	for my $j (0..9) {
		my $n = (9*$i*$j)/(10*$i-$j);
		my $a = (10*$i+$n);
		my $b = (10*$n+$j);
		if ($i != $n and $a > 9 and $a < 100 and $b > 9 and $b < 100 and isint($n)) {
			$c *= $a;
			$p *= $b;
		}
	}
}

say $p/gcd($c, $p);

sub isint{
  my $val = shift;
  return ($val =~ m/^d+$/);
}

sub gcd {
  my $a = $_[0];
  my $b = $_[1];

  while ($b != 0) {
    my $temp = $a % $b;
    $a = $b;
    $b = $temp;
  }

  return abs($a);
}

#!/usr/bin/env perl

use 5.010;

use strict;

use warnings;

my $c = 1;

my $p = 1;

for my $i (1..9) {

for my $j (0..9) {

my $n = (9*$i*$j)/(10*$i-$j);

my $a = (10*$i+$n);

my $b = (10*$n+$j);

if ($i != $n and $a > 9 and $a < 100 and $b > 9 and $b < 100 and isint($n)) {

$c *= $a;

$p *= $b;

}

say $p/gcd($c, $p);

sub isint{

my $val = shift;

return ($val =~ m/^d+$/);

}

sub gcd {

my $a = $_[0];

my $b = $_[1];

while ($b != 0) {

my $temp = $a % $b;

$a = $b;

$b = $temp;

}

return abs($a);

}

프로젝트 오일러 32번

2013/09/14
Perl

1부터 n까지의 각 숫자를 한번씩만 써서 만들 수 있는 숫자를 팬디지털(pandigital)이라고 합니다.
예를 들면 15234는 1부터 5의 숫자가 한번씩만 쓰였으므로 1 ~ 5 팬디지털입니다.

7254라는 숫자는 그런 면에서 특이한데, 39 × 186 = 7254 라는 곱셈식을 만들 때 이것이 1 ~ 9 팬디지털이 되기 때문입니다.

이런 식으로 a × b = c 가 1 ~ 9 팬디지털이 되는 모든 c의 합은 얼마입니까?

(참고: 어떤 c는 두 개 이상의 (a, b)쌍에 대응될 수도 있는데, 이런 경우는 하나로 칩니다)

…오래걸린다.
약 21초 걸렸다.

#!/usr/bin/env perl
use 5.010;
use strict;
use warnings;
use Math::Permute::Array;

my @a = (1..9);
my $r = 0;
my %h;
my $p = new Math::Permute::Array(@a);
for my $i (0..$p->cardinal()-1) {
	my @t = @{$p->permutation($i)};
	for my $i (0..7) {
		for my $j (0..7-$i) {
			no warnings;
			if (join("", @t[0..$i])*join("", @t[$i+1..$j]) == join("", @t[$j+1..8])) {
				$h{join("", @t[$j+1..8])} = 1;
			}
		}
	}
}

$r += $_ foreach (keys %h);
say $r;

#!/usr/bin/env perl

use 5.010;

use strict;

use warnings;

use Math::Permute::Array;

my @a = (1..9);

my $r = 0;

my %h;

my $p = new Math::Permute::Array(@a);

for my $i (0..$p->cardinal()-1) {

my @t = @{$p->permutation($i)};

for my $i (0..7) {

for my $j (0..7-$i) {

no warnings;

if (join("", @t[0..$i])*join("", @t[$i+1..$j]) == join("", @t[$j+1..8])) {

$h{join("", @t[$j+1..8])} = 1;

}

$r += $_ foreach (keys %h);

say $r;

프로젝트 오일러 31번

2013/09/14
Perl

영국의 화폐 단위는 파운드(£)와 펜스(p)이고, 동전의 종류는 아래와 같습니다.

1p, 2p, 5p, 10p, 20p, 50p, £1 (100p), £2 (200p)
이 동전들을 가지고 2파운드를 만드는 방법은 다양할 것입니다. 예를 하나 들면 이렇습니다.

1×£1 + 1×50p + 2×20p + 1×5p + 1×2p + 3×1p
2파운드를 만드는 서로 다른 방법은 모두 몇가지나 있습니까?

사실 슬슬 귀찮았다.
동적계획법을 쓰면 좋았겠지만, 귀찮으니까 그냥 바보같이 짰다.
실행시간은 약 2초

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

my $r = 0;
for my $h (-1..200) {
	next if ($h == -1);
for my $g (-1..max(0, 200 - (200*$h))) {
	next if ($g == -1);
for my $f (-1..max(0, 200 - (100*$g + 200*$h))) {
	next if ($f == -1);
for my $e (-1..max(0, 200 - (40*$f + 100*$g + 200*$h))) {
	next if ($e == -1);
for my $d (-1..max(0, 200 - (20*$e + 40*$f + 100*$g + 200*$h))) {
	next if ($d == -1);
for my $c (-1..max(0, 200 - (10*$d + 20*$e + 40*$f + 100*$g + 200*$h))) {
	next if ($c == -1);
for my $b (-1..max(0, 200 - (4*$c + 10*$d + 20*$e + 40*$f + 100*$g + 200*$h))) {
	next if ($b == -1);
for my $a (-1..max(0, 200 - (2*$b + 4*$c + 10*$d + 20*$e + 40*$f + 100*$g + 200*$h))) {
	next if ($a == -1);
	$r++ if (200 == 1*$a + 2*$b + 5*$c + 10*$d + 20*$e + 50*$f + 100*$g + 200*$h);
}}}}}}}} # sorry...!

say $r;

sub max {
	return $_[0] > $_[1] ? $_[0] : $_[1];
}

#!/usr/bin/env perl

use 5.010;

use strict;

use warnings;

my $r = 0;

for my $h (-1..200) {

next if ($h == -1);

for my $g (-1..max(0, 200 - (200*$h))) {

next if ($g == -1);

for my $f (-1..max(0, 200 - (100*$g + 200*$h))) {

next if ($f == -1);

for my $e (-1..max(0, 200 - (40*$f + 100*$g + 200*$h))) {

next if ($e == -1);

for my $d (-1..max(0, 200 - (20*$e + 40*$f + 100*$g + 200*$h))) {

next if ($d == -1);

for my $c (-1..max(0, 200 - (10*$d + 20*$e + 40*$f + 100*$g + 200*$h))) {

next if ($c == -1);

for my $b (-1..max(0, 200 - (4*$c + 10*$d + 20*$e + 40*$f + 100*$g + 200*$h))) {

next if ($b == -1);

for my $a (-1..max(0, 200 - (2*$b + 4*$c + 10*$d + 20*$e + 40*$f + 100*$g + 200*$h))) {

next if ($a == -1);

$r++ if (200 == 1*$a + 2*$b + 5*$c + 10*$d + 20*$e + 50*$f + 100*$g + 200*$h);

}}}}}}}} # sorry...!

say $r;

sub max {

return $_[0] > $_[1] ? $_[0] : $_[1];

}

blog.update.sh

이 블로그는 더 이상 업데이트되지 않습니다.

Perl (63)

프로젝트 오일러 40번

프로젝트 오일러 39번

프로젝트 오일러 38번

프로젝트 오일러 37번

프로젝트 오일러 36번

프로젝트 오일러 35번

프로젝트 오일러 34번

프로젝트 오일러 33번

프로젝트 오일러 32번

프로젝트 오일러 31번

카테고리

보관함

최신 글

최신 댓글

메타