Python Interpreters Benchmarks
x64 ArchLinux : AMD® Ryzen 7 4700U®

 performance measurements

Each table row shows performance measurements for this program with a particular command-line input value N.

 N  CPU secs Elapsed secs Memory KB Code B ≈ CPU Load
103.190.4696,672894  78% 92% 87% 89% 85% 89% 87% 91%

Read the ↓ make, command line, and program output logs to see how this program was run.

Read  benchmark to see what this program should do.

 notes

  source code

# The Computer Language Benchmarks Game
# http://benchmarksgame.alioth.debian.org/
#
# contributed by Joerg Baumann
# many thanks to Oleg Mazurov for his helpful description

from sys import argv
from math import factorial
from multiprocessing import cpu_count, Pool
from itertools import islice, starmap

def permutations(n, start, size):
    p = bytearray(range(n))

    remainder = start
    for v in range(n - 1, 0, -1):
        rotation_count, remainder = divmod(remainder, factorial(v))
        for _ in range(rotation_count):
            p[:v], p[v] = p[1:v + 1], p[0]

    if size < 2:
        yield p[:]
    else:
        rotations = [(v, factorial(v)) for v in range(3, n)]
        for i in range(start + 2, start + size + 2, 2):
            yield p[:]
            p[0], p[1] = p[1], p[0]
            yield p[:]
            p[1], p[2] = p[2], p[1]
            for v, modulo in rotations:
                if i % modulo != 0:
                    break
                p[:v], p[v] = p[1:v + 1], p[0]

def alternating_flips_generator(n, start, size):
    maximum_flips = 0
    alternating_factor = 1
    for permutation in permutations(n, start, size):
        first = permutation[0]
        if first:
            flips_count = 1
            while True:
                permutation[:first + 1] = permutation[first::-1]
                first = permutation[0]
                if not first: break
                flips_count += 1
            if maximum_flips < flips_count:
                maximum_flips = flips_count
            yield flips_count * alternating_factor
        else:
            yield 0
        alternating_factor = -alternating_factor
    yield maximum_flips

def task(n, start, size):
    alternating_flips = alternating_flips_generator(n, start, size)
    return sum(islice(alternating_flips, size)), next(alternating_flips)

def fannkuch(n):
    assert(n > 0)

    task_count = cpu_count()
    total = factorial(n)
    task_size = (total + task_count - 1) // task_count

    if task_size < 20000:
        task_size = total
        task_count = 1

    assert(task_size % 2 == 0)

    task_args = [(n, i * task_size, task_size) for i in range(task_count)]

    if task_count > 1:
        with Pool() as pool:
            checksums, maximums = zip(*pool.starmap(task, task_args))
    else:
        checksums, maximums = zip(*starmap(task, task_args))

    checksum, maximum = sum(checksums), max(maximums)
    print("{0}\nPfannkuchen({1}) = {2}".format(checksum, n, maximum))

if __name__ == "__main__":
    fannkuch(int(argv[1]))

 make, command-line, and program output logs

 Wed, 28 Sep 2022 08:57:52 GMT

COMMAND LINE:
 /usr/bin/pyston fannkuchredux.pyston-3.pyston 10

PROGRAM OUTPUT:
73196
Pfannkuchen(10) = 38

Revised BSD license

  Home   Conclusions   License   Play