OFFSET
1,1
COMMENTS
Integers whose Jacobi-vector does not form a valid Motzkin-path.
LINKS
Antti Karttunen and J. Moyer, C-program for computing the initial terms of this sequence
FORMULA
a(n) = 4*A095275(n) + 3.
PROG
(Sage)
def is_Motzkin(n, k):
s = 0
for i in (1..k) :
s += jacobi_symbol(i, n)
if s < 0 : return False
return True
def A095101_list(n):
return [m for m in range(3, n+1, 4) if not is_Motzkin(m, m//2)]
A095101_list(467) # Peter Luschny, Aug 08 2012
(PARI) isok(m) = {my(s=0); if(m%4==3, for(i=1, m-1, if((s+=kronecker(i, m))<0, return(1)))); 0; } \\ Jinyuan Wang, Jul 20 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jun 01 2004
STATUS
approved