sum(dispterms(n)+n:n+?gcd(i,2^i),i,1+n:1,501);