login

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”).

A325968
a(n) is the sum k of such a subset of divisors of n with the largest sum and for which n-k and n-(sigma(n)-k) are relatively prime.
6
1, 3, 4, 7, 6, 7, 8, 15, 13, 17, 12, 27, 14, 23, 20, 31, 18, 38, 20, 41, 32, 35, 24, 59, 31, 39, 40, 29, 30, 71, 32, 63, 44, 53, 48, 91, 38, 59, 56, 89, 42, 95, 44, 83, 74, 69, 48, 123, 57, 93, 68, 95, 54, 119, 72, 119, 80, 89, 60, 167, 62, 95, 104, 127, 84, 143, 68, 125, 92, 143, 72, 194, 74, 113, 124, 137, 96, 167, 80, 185, 121, 125
OFFSET
1,2
FORMULA
a(n) = A000203(n) - A325967(n).
a(n) = n + A325969(n).
For all n:
a(A000040(n)) = A000040(n)+1.
a(A000396(n)) = A000396(n)+1.
a(n) <= A325818(n).
EXAMPLE
For n=15, its divisors are [1, 3, 5, 15]. If we take the full set [1, 3, 5, 15] and its complement [], their sums are 24 and 0, but gcd(15-0, 24-15) = gcd(15, 9) = 3 > 1. If we take subsets [1] and [3, 5, 15], then their sums are 1 and 23, but gcd(15-1, 23-15) = gcd(14,8) = 2 > 1. If we take subsets [3] and [1, 5, 15], their sums are 3 and 21, but gcd(15-3, 21-15) = gcd(12, 6) = 6 > 1. Only when we take the subset with the four smallest sum, [1, 3] and its complement [5, 15], we get such sums 4 and 20 for which gcd(15-4, 20-15) = gcd(11, 5) = 1. Thus a(15) = 20, the size of the subset with larger sum.
PROG
(PARI)
A325968(n) = { my(divs=divisors(n), s=sigma(n), r, ms=0); for(b=0, (2^(length(divs)))-1, r=sumbybits(divs, b); if(1==gcd(n-(s-r), n-r), ms=max(r, ms))); (ms); };
sumbybits(v, b) = { my(s=0, i=1); while(b>0, s += (b%2)*v[i]; i++; b >>= 1); (s); };
KEYWORD
nonn
AUTHOR
Antti Karttunen, May 29 2019
STATUS
approved