OFFSET
1,5
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..1000
FORMULA
EXAMPLE
a(1) = 1 mod 1 = 0;
a(2) = 9 mod 3 = 0;
a(3) = 15 mod 5 = 0;
a(4) = 21 mod 7 = 0;
a(5) = 25 mod 9 = 7.
MAPLE
A014076 := proc(n) option remember; local a; if n = 1 then 1; else for a from procname(n-1)+2 by 2 do if not isprime(a) then return a; end if; end do: end if; end proc:
A163505 := proc(n) A014076(n) mod (2*n-1); end: seq(A163505(n), n=1..80) ; # R. J. Mathar, Oct 10 2009
MATHEMATICA
A014076 := Select[Range[1, 10000, 2], PrimeOmega[#] != 1 &]; Table[Mod[A014076[[n]], 2*n - 1], {n, 1, 50}] (* G. C. Greubel, Jul 27 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Jul 29 2009
EXTENSIONS
Missing term between a(53) and a(54) inserted by G. C. Greubel, Jul 27 2017
STATUS
approved