We may assume m≥n. It is well known that
gcd(2p−1,2q−1)=2gcd(p,q)−1,
so that
gcd(2m−2n,2m2+mn+n2−1)=gcd(2m−n−1,2m2+mn+n2−1)=2gcd(m−n,m2+mn+n2)−1.
Next, consider a divisor d∣m−n. We must have gcd(m,d)=1, since m and n are relatively prime. It follows that 0≡m2+mn+n2≡3m2(modd) is equivalent to d∣3, and we infer that
gcd(m−n,m2+mn+n2)=gcd(m−n,3),
which is 1 or 3.
Hence gcd(2m−2n,2m2+mn+n2−1) may only assume the values 1 and 7. Both values are possible, since m=2,n=1 gives
gcd(22−21,222+2⋅1+12−1)=gcd(2,27−1)=1,
and m=1,n=1 gives
gcd(21−21,212+1⋅1+12−1)=gcd(0,23−1)=7.