login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A226020 Composite squarefree numbers n such that the ratio (n + 1/2)/(p(i) + 1/2) is an integer, where p(i) are the prime factors of n. 10

%I #18 Jun 04 2013 00:16:00

%S 13702,42997,1004062,1684462,38447662,40243549,70801087,107728582,

%T 409055062,594021862,760767262,1045475437,1104435202,1471700587,

%U 1686747562,1920806662,3136180162,3469071937,5291041297,7239716347,7903353667,12738885862,22711489762

%N Composite squarefree numbers n such that the ratio (n + 1/2)/(p(i) + 1/2) is an integer, where p(i) are the prime factors of n.

%C Also composite squarefree numbers n such that (2*p(i)+1) | (2*n+1).

%H Giovanni Resta, <a href="/A226020/b226020.txt">Table of n, a(n) for n = 1..65</a> (terms < 3*10^12)

%e The prime factors of 13702 are 2, 13, 17 and 31. We see that (13702 + 1)/(2 + 1/2) = 5481, (13702 + 1/2)/(13 + 1/2) = 1015, (13702 + 1)/(17 + 1/2) = 783 and ( 13702 + 1/2)/(31 + 1/2) = 435. Hence 13702 is in the sequence.

%e The prime factors of 1123545 are 3, 5 and 74903. We see that

%e (1123545 + 1/2)/(3 + 1/2) = 321013, (1123545 + 1/2)/(5 + 1/2) = 204281 but (1123545 + 1/2)/(74903+ 1/2) = 321013/21401. Hence 1123545 is not in the sequence.

%p with(numtheory); A226020:=proc(i, j) local c, d, n, ok, p;

%p for n from 2 to i do if not isprime(n) then p:=ifactors(n)[2]; ok:=1;

%p for d from 1 to nops(p) do if p[d][2]>1 or not type((n+j)/(p[d][1]+j),integer) then ok:=0; break; fi; od;

%p if ok=1 then print(n); fi; fi; od; end: A226020(10^9,1/2);

%Y Cf. A208728, A225702-A225720, A226111-A226114.

%K nonn,hard

%O 1,1

%A _Paolo P. Lava_, May 23 2013

%E a(9)-a(23) from _Giovanni Resta_, Jun 02 2013

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 07:20 EDT 2024. Contains 371235 sequences. (Running on oeis4.)