For a positive integer n, let us denote by S(n) the sum of all the positive factors of n whose remainder when divided by 4 is not equal to 2. Let us first determine S(n).
Suppose the prime factorization of n is given by
2mp1m1⋯pkmk(p1,…,pk are distinct odd primes, m≥0,m1,…,mk≥1).
Since the fact that an integer has a remainder 2 when divided by 4 is equivalent to the fact that it is divisible by 2 only once, we can see that S(n) is the sum of all numbers of the form
2lp1l1⋯pklk(where 0≤l≤m,1=l=1,0≤l1≤m1,…,0≤lk≤mk).
Consequently, S(n) equals
l=0,l=1∑m2ll1=0∑m1p1l1⋯lk=0∑mkpklk.
(Note that because of the distributive law the number of terms in each sum corresponds to the number of possible values for each of the exponents.)
For the sake of simplicity, let for each non-negative integer m,
f(2,m)=l=0,l=1∑m2l;f(p,m)=l=0∑mpl(when p is a prime =2).
Then, if n=2mp1m1⋯pkmk, we have S(n)=f(2,m)f(p1,m1)⋯f(pk,mk). In order to determine positive integers n for which S(n)=1000, let us first determine the pairs (p,m), where p is a prime and m is a positive integer for which f(p,m) is a factor of 1000.
When p=2, we get that if m≥9, then f(2,m)≥f(2,9)=1021. So, it is sufficient to consider the cases for m≤8, and we can conclude that f(2,1)=1, f(2,2)=5, f(2,6)=125 are the only cases which give a factor of 1000 for f(2,m).
When 3≤p≤31, we can similarly check that f(3,1)=4, f(3,3)=40, f(7,1)=8, f(19,1)=20 are the only cases for this range of primes p for which f(p,m) is a factor of 1000.
When p≥32, we get if m≥2 f(p,m)≥f(m,2)=1+p+p2≥1+32+322>1000. So, it is enough to check the cases for m=1 only for this range of p, and we get f(199,1)=200, f(499)=500 as the only possibilities for a factor of 1000.
Finally, we search for combinations of these values of f(p,m)'s which yield the product 1000, and we find that the desired answer is given by 26×71=448 and 22×1991=796.