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

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))

根据题目中的代码，这个程序的主要功能是计算所有满足以下条件的数的和：该数是质数且该数的形式为形如3n+1的数。我们可以按照以下步骤分析这个程序：

首先，函数f(n)初始化一个变量s为0，然后通过一个循环遍历从2到n+1的所有整数（包括2和n+1）。对于循环中的每一个数i，代码进入一个内层循环，该内层循环用于检查i是否可被其小于等于其平方根的所有正整数整除。如果存在这样的数，那么i就不是质数，变量p被设置为False并跳出内层循环。如果内层循环完成后，变量p仍为True，那么我们知道i是一个质数。然后，代码检查这个质数是否满足形式为3n+1的条件。如果满足这个条件，那么该质数就被加入到变量s中。最后，函数返回变量s的值。

在给定的例子中，我们调用函数f(13)。我们需要找出所有在2到14之间（包括2和14）的质数，且这些数要满足形式为3n+1的条件。这些数是：2（虽然它是最小的质数，但它不满足形式为3n+1的条件），7和13（它们都是质数并且满足形式为3n+1的条件）。因此，执行这个程序的结果是：7 + 13 = 20。所以正确答案是B，即输出结果为20。