-3

我正在尝试制作一个可以打印 x,y 的 n 阶多项式的函数

即poly(x,y,1)会输出c[0] + c[1]*x + c[2]*y

即poly(x,y,2)会输出c[0] + c[1]*x + c[2]*y + c[3]*x**2 + c[4]*y**2 + c[5]*x*y

你能给我一些想法吗？也许是迭代工具？

4

3 回答 3

3

你可以尝试从类似的东西开始

def poly(x,y,n):
    counter = 0
    for nc in range(n+1):
        for i in range(nc+1):
            print "c[", counter, "]",
            print " * ", x, "**", i,
            print " * ", y, "**", nc-i,
            print " + ",
            counter += 1

例如

poly("x", "y", 2)

会产生

c[ 0 ]  *  x ** 0  *  y ** 0  +  c[ 1 ]  *  x ** 0  *  y ** 1  +  c[ 2 ]  *  x ** 1  *  y ** 0  +  c[ 3 ]  *  x ** 0  *  y ** 2  +  c[ 4 ]  *  x ** 1  *  y ** 1  +  c[ 5 ]  *  x ** 2  *  y ** 0  +

内置ifs，如果你想抑制不需要的输出。

于 2012-07-02T19:24:14.013 回答
2

既然您想要一个带有 itertools 的功能性解决方案，这里有一个单行：

import itertools as itt
from collections import Counter
n = 3
xy = ("x", "y") # list of variables may be extended indefinitely
poly = '+'.join(itt.starmap(lambda u, t: u+"*"+t if t else u, zip(map(lambda v: "C["+str(v)+"]", itt.count()),map(lambda z: "*".join(z), map(lambda x: tuple(map(lambda y: "**".join(map(str, filter(lambda w: w!=1, y))), x)), map(dict.items, (map(Counter, itt.chain.from_iterable(itt.combinations_with_replacement(xy, i) for i in range(n+1))))))))))

那会给你

C[0]+C[1]*x+C[2]*y+C[3]*x**2+C[4]*y*x+C[5]*y**2+C[6]*x**3+C[7]*y*x**2+C[8]*y**2*x+C[9]*y**3

请注意，系数的顺序略有不同。这不仅适用于任何 n，而且适用于任意数量的变量（x、y、z 等......）

只为一笑

于 2012-07-02T23:58:12.690 回答
1

稍微笼统一点：

from itertools import product

def make_clause(c, vars, pows):
    c = ['c[{}]'.format(c)]
    vp = (['', '{}', '({}**{})'][min(p,2)].format(v,p) for v,p in zip(vars,pows))
    return '*'.join(c + [s for s in vp if s])

def poly(vars, max_power):
    res = (make_clause(c, vars, pows) for c,pows in enumerate(product(*(range(max_power+1) for v in vars))))
    return ' + '.join(res)

然后poly(['x', 'y'], 2)返回

"c[0] + c[1]*y + c[2]*(y**2) + c[3]*x + c[4]*x*y + c[5]*x*(y**2) + c[6]*(x**2) + c[7]*(x**2)*y + c[8]*(x**2)*(y**2)"
于 2012-07-02T23:59:19.077 回答