#!/usr/bin/env python3
from sympy import *
def L(n, p):
    N = prod([n - k for k in range(n) if k != p])
    D = prod([p - k for k in range(n) if k != p])
    return N//D
for n in range(2, 10):
    sgn = "" if sign(L(n, 0)) >=0 else "-"
    print(f"u{n} = {sgn}u0", end="")
    for p in range(1, n):
        a = L(n, p)
        sgn = " +" if sign(a) >= 0 else " -"
        print(f"{sgn} {abs(a)}u{p}", end="")
    print()