def S(p,a): s=0 for k in range(0,a+1): s=s+p**k return s