login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

a(n) = (number of x with largest prime factor of the form 4k+3 less than or equal to n) minus (number of x with largest prime factor of the form 4k+1 less than or equal to n).
1

%I #15 Nov 05 2024 21:15:02

%S 0,0,1,1,0,1,2,2,3,2,3,4,3,4,3,3,2,3,4,3,4,5,6,7,6,5,6,7,6,5,6,6,7,6,

%T 7,8,7,8,7,6,5,6,7,8,7,8,9,10,11,10,9,8,7,8,9,10,11,10,11,10,9,10,11,

%U 11,10,11,12,11,12,13,14,15,14,13,12,13,14,13,14,13,14,13,14,15,14,15,14,15

%N a(n) = (number of x with largest prime factor of the form 4k+3 less than or equal to n) minus (number of x with largest prime factor of the form 4k+1 less than or equal to n).

%C Weak conjecture: sequence gives positive values only; strong conjecture: for n large enough a(n) > 2*sqrt(n).

%H Robert Israel, <a href="/A071825/b071825.txt">Table of n, a(n) for n = 1..10000</a>

%H Robert Israel, <a href="/A071825/a071825.png">Log-log plot of (n, a(n)) for 1 <= n <= 10^6</a>

%F a(n) = Card( x<=n : A006530(x) == 3 mod 4) - Card( x<=n : A006530(x) == 1 mod 4)

%F a(n) = A071823(n) - A071824(n). - _Sean A. Irvine_, Aug 16 2024

%p f:= proc(n) local t; t:= max(numtheory:-factorset(n)) mod 4; if t = 3 then 1 elif t= 1 then -1 else 0 fi end proc:

%p ListTools:-PartialSums(map(f, [$1..100])); # _Robert Israel_, Nov 05 2024

%o (PARI) for(n=1,200,print1(sum(i=2,n,if((component(component(factor(i),1),omega(i))-3)%4,0,1))-sum(i=2,n,if((component(component(factor(i),1),omega(i))-1)%4,0,1)),","))

%Y Cf. A071823, A071824.

%K easy,nonn,look

%O 1,7

%A _Benoit Cloitre_, Jun 07 2002

%E Missing a(1)=0 inserted by _Sean A. Irvine_, Aug 16 2024