cp's OEIS Frontend

A337121 a(n) is the number of ways the n-th prime number prime(n) can be represented as sum of two smaller odd prime numbers p1, p2 with prime(n) > p1 > (p2 minus the maximum odd prime factor of (p1-p2)).

%I A337121 #20 Jan 09 2021 02:39:23
%S A337121 1,1,2,3,2,2,4,3,3,5,3,6,4,4,4,7,4,4,5,6,6,6,9,7,8,8,7,7,6,11,4,11,9,
%T A337121 7,8,9,7,13,12,6,10,15,10,9,7,13,13,11,13,10,15,10,13,14,11,13,13,12,
%U A337121 14,17,13,13,19,9,14,19,12,8,14,22,17,14,13,16,9,15
%N A337121 a(n) is the number of ways the n-th prime number prime(n) can be represented as sum of two smaller odd prime numbers p1, p2 with prime(n) > p1 > (p2 minus the maximum odd prime factor of (p1-p2)).
%C A337121 This sequence counts the cases such that prime(n) = p1 + p2 - MaxOddPrimeFactor(p1-p2), where MaxOddPrimeFactor(m) is defined as the maximum odd prime factor of the positive integer m. If there is no odd prime factor of m, MaxOddPrimeFactor(m) is defined as 1.
%C A337121 Conjecture: a(n) > 0 when n >= 4.
%C A337121 Some nonprime odd numbers, like 27, cannot be partitioned into the form of p1 + p2 - MaxOddPrimeFactor(p1-p2).
%e A337121 When n=4, prime(4)=7, MaxOddPrimeFactor(5-3)=1, 7=5+3-1. This is the only case, so a(4)=1.
%e A337121 When n=5, prime(5)=11, MaxOddPrimeFactor(7-5)=1, 11=7+5-1. This is the only case, so a(5)=1.
%e A337121 When n=6, prime(6)=13, MaxOddPrimeFactor(11-3)=1, 13=11+3-1; and MaxOddPrimeFactor(11-5)=3, 13=11+5-3. Two cases found, so a(6)=2.
%t A337121 MaxOddPrimeFactor[m_] :=
%t A337121 Module[{factors, l, res}, factors = FactorInteger[m];
%t A337121   l = Length[factors]; res = factors[[l, 1]]; If[res == 2, res = 1];
%t A337121   res]
%t A337121 Table[p = Prime[n]; p1 = NextPrime[p/2, -1]; ct = 0;
%t A337121 While[p1 = NextPrime[p1]; p1 < p, p2 = NextPrime[p - p1, -1];
%t A337121   While[p2 = NextPrime[p2]; p2 < p1,
%t A337121    If[p == (p1 + p2 - MaxOddPrimeFactor[p1 - p2]), ct++]]]; ct, {n, 4,
%t A337121    79}]
%K A337121 nonn
%O A337121 4,3
%A A337121 _Lei Zhou_, Aug 17 2020