%I #59 Feb 13 2024 06:58:42
%S 1,3,8,13,23,38,39,41,43,48,50,53,56,57,58,66,68,70,73,77,84,90,94,98,
%T 126,128,134,140,143,145,149,151,153,157,160,164,167,168,172,174,176,
%U 182,191,194,196,200,210,212,215,217,218,219,222,225,228,229
%N Numbers k such that sum of first k primes divides product of first k primes.
%H Amiram Eldar, <a href="/A051838/b051838.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1000 from T. D. Noe)
%F From _Reinhard Zumkeller_, Oct 03 2011: (Start)
%F A002110(a(n)) mod A007504(a(n)) = 0.
%F A116536(n) = A002110(a(n)) / A007504(a(n)). (End)
%e Sum of first 8 primes is 77 and product of first 8 primes is 9699690. 77 divides 9699690 therefore a(3)=8.
%t p = Prime@ Range@ 250; Flatten@ Position[ Mod[ First@#, Last@#] & /@ Partition[ Riffle[ Rest[ FoldList[ Times, 1, p]], Accumulate@ p], 2], 0] (* _Harvey P. Dale_, Dec 19 2010 *)
%o (Haskell)
%o import Data.List (elemIndices)
%o a051838 n = a051838_list !! (n-1)
%o a051838_list =
%o map (+ 1) $ elemIndices 0 $ zipWith mod a002110_list a007504_list
%o -- _Reinhard Zumkeller_, Oct 03 2011
%o (PARI) for(n=1,100,P=prod(i=1,n,prime(i));S=sum(i=1,n,prime(i));if(!(P%S),print1(n,", "))) \\ _Derek Orr_, Jul 19 2015
%o (PARI) isok(n) = my(p = primes(n)); (vecprod(p) % vecsum(p)) == 0; \\ _Michel Marcus_, Dec 20 2018
%o (GAP) P:=Filtered([1..2000],IsPrime);;
%o Filtered([1..Length(P)],n->Product([1..n],i->P[i]) mod Sum([1..n],i->P[i])=0); # _Muniru A Asiru_, Dec 20 2018
%Y Cf. A007504, A002110, A159639, A196415.
%Y A116536 gives the quotients, A140763 the divisors and A159578 the dividends.
%K nonn
%O 1,2
%A _G. L. Honaker, Jr._, Dec 12 1999