Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%I #39 Apr 24 2018 06:00:02
%S 0,2,4,6,6,6,6,6,10,18,18,22,22,24,24,24,24,24,24,24,24,26,28,34,44,
%T 46,46,46,46,46,57,58,61,61,61,61,61,61,61,61,61,61,61,61,61,61,61,62,
%U 62,62,62,62,62,70,74,78,82,82,82,82,82,90,110,110,110,110,126,130,136,138,138,142,142,142,142
%N Maximum remainder of prime(p) + prime(q) divided by p + q with p <= q <= n.
%C Odd numbers k which are terms of this sequence are 57, 61, 353, 2113, ...
%C Approximate self-similar growing patterns appear at different scales which suggest a fractal-like structure, see plots in Links section.
%H Altug Alkan, <a href="/A302754/b302754.txt">Table of n, a(n) for n = 1..10000</a>
%H Andres Cicuttin, <a href="/A302754/a302754.pdf">Several plots showing similar stair-like patterns</a>
%e a(1) = 0 because only option is p = q = 1.
%e a(4) = a(8) = 6 because (prime(4) + prime(4)) mod 8 = (prime(8) + prime(7)) mod 15 = 6 is the largest remainder for both.
%e a(31) = 57 because (prime(28) + prime(31)) mod 59 = 57 is the largest remainder.
%t a[n_]:=Table[Table[Mod[Prime[j]+Prime[i],i+j],{i,1,j}],{j,1,n}]//Flatten//Max;
%t Table[a[n],{n,1,100}]
%o (PARI) a(n) = vecmax(vector(n, q, vecmax(vector(q, p, (prime(p)+prime(q)) % (p+q)))));
%Y Cf. A247824, A302245, A302446.
%K nonn,look
%O 1,2
%A _Andres Cicuttin_ and _Altug Alkan_, Apr 12 2018