﻿from sympy import *

def f(x):
    y=x**2-x+41
    return y

def Test(n):
    c=0
    for k in range(0,n+1):
        if isprime(f(k)):
            c=c+1
    return c
