from numpy.polynomial import Polynomial

def c(n) :
    P,Q = [],[]
    for k in range(n+1) :
        if k%3 == 0 :
            P.append(1)
        else :
            P.append(0)
        if k%5 == 0 :
            Q.append(1)
        else :
            Q.append(0)
    s = 0
    for k in range(n+1) :
        s += P[k]*Q[n-k]
    return s
    
# for n in range(200) :
#     print(c(n))
for n in range(185) :
    print(c(n),c(n+15))
    
P = Polynomial([1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-1])*Polynomial([1,-1])
Q = Polynomial([1,0,0,-1])*Polynomial([1,0,0,0,0,-1])
print(P%Q)

def d(n) :
    s = 0
    for u in range(n) :
        for v in range(n) :
            if 3*u+5*v == n : s+=1
    return s
    
for n in range(200) :
    print(c(n),d(n))