while($0 < H * W)
  $1 = $0 / W + 1
  $2 = $0 % W + 1
  $3 = $3 + (((A[$1, $2] - B[$1, $2]) * (A[$1, $2] - B[$1, $2]) + H * W - 1) / (H * W)) * $0
  $0 = $0 + 1
end

return($3 / W + 1, $3 % W + 1)