A080113 Positions of A080115 in A000040.
1, 7, 8, 10, 13, 14, 16, 18, 19, 21, 24, 25, 26, 28, 29, 30, 31, 33, 34, 35, 37, 38, 40, 41, 42, 44, 45, 47, 48, 49, 50, 51, 53, 55, 57, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 77, 78, 79, 80, 82, 84, 86, 87, 88, 89, 90, 91, 93, 94, 95, 97, 98, 99, 100, 101
Offset: 1
Keywords
Programs
-
Maple
A080113 := proc(n) option remember; local i; if(1 = n) then RETURN(1); fi; i := A080113(n-1)+1; while(i > 0) do if(A080111(i) > 0) then RETURN(i); fi; i := i+1; od; end;
-
Mathematica
s[p_, u_] := Sum[JacobiSymbol[j, p], {j, 1, u}]; Join[{1}, Select[Range[ 100], !(p = Prime[#]; AllTrue[Range[(p - 1)/2], s[p, #] >= 0 &]) &]] (* Jean-François Alcover, Mar 07 2016 *)