pascal.py

import subprocess

def pascalMod2(row):
    rv = [1]
    for i, x in enumerate(row, start=1):
        if i == len(row):
            break
        rv.append((x + row[i]) % 2)
    rv.append(1)
    assert len(rv) == len(row) + 1
    return rv

def rowToBin(row):
    rv = 0
    for x in row:
        rv <<= 1
        if x:
            rv |= 1
    return rv

def factor(x):
    # XXX hax
    return subprocess.check_output(["factor", str(x)]).strip()

row = [1]
for i in range(40):
    print row
    b = rowToBin(row)
    print b, '=', factor(b)
    row = pascalMod2(row)