%I #18 Dec 05 2023 08:25:08
%S 2,3,6,7,10,11,12,14,15,18,19,22,23,24,27,30,31,32,34,35,38,39,40,44,
%T 46,47,48,51,52,55,56,58,59,60,63,64,66,67,70,71,72,74,75,76,78,79,82,
%U 83,86,87,88,91,92,94,95,96,98,99,100,102,104,106,108,110,112,114,115,116,118,119,120,122
%N Numbers k such that A014285(k) and A007504(k) are coprime.
%C Numbers k such that A306834(k) = A014285(k).
%C No terms == 1 (mod 4).
%C Numbers k such that A309036(k)=1. - _Robert Israel_, Jul 09 2019
%H Robert Israel, <a href="/A307414/b307414.txt">Table of n, a(n) for n = 1..10000</a>
%e a(3) = 6 is a term because A007504(6) = 41 and A014285(6) = 184 are coprime.
%p N:= 1000: # for terms <= N
%p Primes:= map(ithprime, [$1..N]):
%p S1:= ListTools:-PartialSums(Primes):
%p S2:= ListTools:-PartialSums(zip(`*`,Primes, [$1..N])):
%p select(t -> igcd(S1[t],S2[t])=1, [$1..N]);
%t okQ[n_] := With[{pp = Prime[Range[n]]}, CoprimeQ[Total[pp], Total[pp.Range[n]]]];
%t Select[Range[200], okQ] (* _Jean-François Alcover_, Dec 05 2023 *)
%o (PARI) isok(k) = my(vp=primes(k)); gcd(sum(i=1, k, vp[i]), sum(i=1, k, i*vp[i])) == 1; \\ _Michel Marcus_, Apr 07 2019
%Y Cf. A007504, A014285, A306834, A309036.
%K nonn
%O 1,1
%A _Robert Israel_, Apr 07 2019
|