Powerful digit counts

AUTHOR

polettix

https://projecteuler.net/problem=63

The 5-digit number, 16807=7^5, is also a fifth power. Similarly, the 9-digit number, 134217728=8^9, is a ninth power.

How many n-digit positive integers exist which are also an nth power?

use v6;



# As of August 24th, 2009 we don't have big integers, so we'll have
# to conjure up something. We'll represent each number with an
# array of digits, base 10-exp for ease of length computation. The most
# significant part is at the end of the array, i.e. the array should
# be read in reverse.
# Setting '1' for the number of digits means representing the base-10
# system with one digit in each array position.
my $digits = 5;
my $limit = 10 ** $digits;

my $count = 0;

# 9 is the maximum possible base for this problem. 9**22 has 21 digits
sub MAIN(Bool :$verbose = False) {
    for 1 .. 9 -> $x {
        my @x = (1);
        for 1 .. * -> $y {
            @x = multby(@x, $x);
            my $px = printable(@x);
            if ($px.encode('utf-8').bytes == $y) {
                say "$x ** $y = $px (", $px.encode('utf-8').bytes, ')'
                    if $verbose;
                $count++;
            }
            elsif ($px.encode('utf-8').bytes < $y) {
                last;
            }
        }
    }
    say $count;
}

sub printable (@x is copy) {
    my $msb = pop @x;
    return $msb ~ @x.reverse.map({sprintf '%0'~$digits~'d', $_ }).join('');
}

# Add a "number" to another, modifies first parameter in place.
# This assumes that length(@y) <= length(@x), which will be true in
# our program because @y is lower than @x
sub add (@x is copy, @y) {
    my $rest = 0;
    return add(@y, @x) if +@x < +@y;
    for @x Z (@y, 0, *) -> $x is rw, $y {
        $x += $y + $rest;
        $rest = int($x / $limit);
        $x %= $limit;
    }
    push @x, $rest if $rest;
    return @x;
}

sub multby (@x is copy, $y) {
    my $rest = 0;
    for @x -> $x is rw {
        $x = $x * $y + $rest;
        $rest = $x div $limit;
        $x %= $limit;
    }
    push @x, $rest if $rest;
    return @x;
}

# Not really needed...
sub mult (@x is copy, @y) {
    my @result = (0);
    for @y -> $y {
        my @addend = multby(@x, $y);
        @result = add(@result, @addend);
        @x.unshift(0);
    }
    return @result;
}

# vim: expandtab shiftwidth=4 ft=perl6

See Also

prob001-cspencer.raku

Multiples of 3 and 5

prob001-eric256.raku

Multiples of 3 and 5

prob001-grondilu.raku

Multiples of 3 and 5

prob001-hexmode.raku

Multiples of 3 and 5

prob001-unobe.raku

Multiples of 3 and 5

prob002-eric256.raku

Even Fibonacci numbers

prob002-gerdr.raku

Even Fibonacci numbers

prob002-hexmode.raku

Even Fibonacci numbers

prob003-eric256.raku

Largest prime factor

prob003-gerdr.raku

Largest prime factor

prob003-hexmode.raku

Largest prime factor

prob003-lanny.raku

Largest prime factor

prob004-unobe.raku

Largest palindrome product

prob004-xfix.raku

Largest palindrome product

prob005-unobe.raku

Smallest multiple

prob005-xfix.raku

Smallest multiple

prob006-polettix.raku

Sum square difference

prob007-polettix.raku

10001st prime

prob008-duff.raku

Largest product in a series

prob008-duff2.raku

Largest product in a series

prob009-gerdr-feeds.raku

Special Pythagorean triplet

prob009-gerdr.raku

Special Pythagorean triplet

prob009-polettix.raku

Special Pythagorean triplet

prob010-polettix.raku

Summation of primes

prob011-moritz.raku

Largest product in a grid

prob012-polettix.raku

Highly divisible triangular number

prob013-grondilu.raku

Large sum

prob014-felher.raku

Longest Collatz sequence

prob015-felher.raku

Lattice paths

prob016-grondilu.raku

Power digit sum

prob017-duff.raku

Number letter counts

prob018-felher.raku

Maximum path sum I

prob019-grondilu.raku

Counting Sundays

prob020-grondilu.raku

Factorial digit sum

prob021-gerdr.raku

Amicable numbers

prob022-grondilu.raku

Names scores

prob023-shlomif.raku

Non-abundant sums

prob024-moritz.raku

Lexicographic permutations

prob025-polettix.raku

1000-digit Fibonacci number

prob026-shlomif.raku

Reciprocal cycles

prob027-shlomif.raku

Quadratic primes

prob028-shlomif.raku

Number spiral diagonals

prob029-gerdr.raku

Distinct powers

prob029-polettix.raku

Distinct powers

prob031-shlomif.raku

Coin sums

prob033-andreoss.raku

Digit cancelling fractions

prob034-quinny.raku

Digit factorials

prob036-xenu.raku

Double-base palindromes

prob038-andreoss.raku

Pandigital multiples

prob039-quinny.raku

Integer right triangles

prob041-heyajulia-alternative.raku

Pandigital Prime

prob041-heyajulia.raku

Pandigital Prime

prob042-shlomif.raku

Coded triangle numbers

prob047-gerdr.raku

Distinct primes factors

prob052-duff.raku

Permuted multiples

prob053-duff.raku

Combinatoric selections

prob053-gerdr.raku

Combinatoric selections

prob054-andreoss.raku

Poker hands

prob055-shlomif.raku

Lychrel numbers

prob056-shlomif.raku

prob059-andreoss.raku

XOR decryption

prob063-moritz.raku

Powerful digit counts

prob065-andreoss.raku

Convergents of e

prob065-grondilu.raku

prob066-andreoss.raku

Diophantine equation

prob067-felher.raku

Maximum path sum II

prob080-andreoss.raku

Square root digital expansion

prob081-moritz.raku

Path sum: two ways

prob089-andreoss.raku

Roman numerals

prob092-moritz.raku

Square digit chains

prob097-andreoss.raku

Large non-Mersenne prime

prob098-andreoss.raku

Anagramic squares

prob099-andreoss.raku

Largest exponential

README.md

The Camelia image is copyright 2009 by Larry Wall. "Raku" is trademark of the Yet Another Society. All rights reserved.

Built with Podlite — the markup and publishing tools behind this site.