Each table row shows performance measurements for this PyPy 3 program with a particular command-line input value N.
| N | CPU secs | Elapsed secs | Memory KB | Code B | ≈ CPU Load |
|---|---|---|---|---|---|
| 10 | 2.78 | 0.68 | 80,052 | 894 | 49% 93% 52% 52% 48% 52% 53% 46% |
Read the ↓ make, command line, and program output logs to see how this program was run.
Read fannkuch-redux benchmark to see what this program should do.
# 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]))
Thu, 08 Sep 2022 07:54:40 GMT COMMAND LINE: /usr/bin/pypy3 fannkuchredux.pypy3-3.pypy3 10 PROGRAM OUTPUT: 73196 Pfannkuchen(10) = 38