执行以下程序后所输出的结果是？（    ）

def f(n):

   s=0

   for i in range(2, n+1):

       p = True

       for j in range(2, int(i**0.5)+1):

           if i % j == 0:

               p = False

               break    

       if p==True and i%3==1:

           s=s+i

   return s

print(f(13))