|
|
A246261
|
|
Numbers n such that A003961(n) is of the form 4k+1.
|
|
17
|
|
|
1, 3, 4, 9, 10, 11, 12, 13, 14, 16, 23, 25, 27, 30, 31, 33, 34, 35, 36, 37, 38, 39, 40, 42, 44, 47, 48, 49, 52, 56, 58, 59, 64, 69, 71, 75, 81, 82, 83, 85, 86, 89, 90, 92, 93, 95, 97, 99, 100, 102, 105, 106, 107, 108, 109, 110, 111, 114, 117, 119, 120, 121, 122, 124, 126, 130, 131, 132, 133, 134, 136, 139, 140, 141, 143, 144
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Equally: numbers n for which A048673(n) is odd.
Also, after 1, numbers n for which A243062(n) is even.
|
|
LINKS
|
|
|
FORMULA
|
|
|
MATHEMATICA
|
{1}~Join~Select[Range[144], Mod[Times @@ Map[NextPrime[#1]^#2 & @@ # &, FactorInteger[#]], 4] == 1 &] (* Michael De Vlieger, Mar 12 2021 *)
|
|
PROG
|
(define A246261 (MATCHING-POS 1 1 (lambda (n) (= 1 (modulo (A003961 n) 4)))))
(Python)
from sympy import factorint, prime, primepi
from operator import mul
def a003961(n):
f=factorint(n)
return 1 if n==1 else reduce(mul, [prime(primepi(i) + 1)**f[i] for i in f])
print([n for n in range(1, 201) if a003961(n)%4==1]) # Indranil Ghosh, Jun 12 2017
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|