# Liber Abaci XI.7.13 a

# SOLUZIONI INTERE POSITIVE DI
#       x1 + x2 + x3 + x4 = 100 
#       3x1 + 2x2+ x3/2 + x4 = 100     

print("\tLiber Abaci XI.7.13\n\tSOLUZIONI INTERE POSITIVE DI\n\t\tx1 + x2 + x3 + x4 = 100\n\t\t3x1 + 2x2+ x3/2 + x4 = 100\n")

n=1
for x4 in range (1,100):
    T=100-x4
    u=int(T/3)+1
    while u<2*T/5:
        a=3*u-T
        b=2*T-5*u
        x1=a
        x2=b
        x3=4*a+2*b
        if (n-1)%4 == 0: print()	 # a capo dopo 4 risultati		
        print('  sol.',n,':  (',x1,x2,x3,x4,')',end="\t")
#        print("verifica:..",3*x1+2*x2+x3/2+x4,"== 100")
        u=u+1
        n=n+1