fortran中的combinations_with_replacement函数？

Question

我用Python语言编写了一个小模块，用于获取适合y总成本的x产品的所有可能性。该模块运行良好，但速度较慢。计算6个产品到每个产品的30次迭代大约需要6个小时。所以，我在考虑用FORTRAN重写脚本，看看能不能挤出更快的速度。不幸的是，我是第一次接触FORTRAN，现在还不了解大多数的库等等。

在FORTRAN语言中，有没有类似于python的itertools.combinations_with_replacement(pool, r)的模块/函数，或者实现相同功能的模块？

Answer 1

别干那事。当你需要进行算法优化时，你正在尝试进行微优化。这就是为什么我认为should not use the exponential-time solution是基于itertools而不是better, recursive solution的原因。

Answer 2

一言以蔽之，itertools documentation对combinations_with_replacement()有一个纯python等效项。它很简洁，翻译成Fortran应该不难

def combinations_with_replacement(iterable, r):
    # combinations_with_replacement('ABC', 2) --> AA AB AC BB BC CC
    pool = tuple(iterable)
    n = len(pool)
    if not n and r:
        return
    indices = [0] * r
    yield tuple(pool[i] for i in indices)
    while True:
        for i in reversed(range(r)):
            if indices[i] != n - 1:
                break
        else:
            return
        indices[i:] = [indices[i] + 1] * (r - i)
        yield tuple(pool[i] for i in indices)