Построить дороги из единиц между городами из "-1"

#include <iostream>
#include <cmath>
 
using namespace std;
 
double s(int x1, int y1, int x2, int y2, int x, int y){
 
    return sqrt((x1 - x)*(x1 - x) + (y1 - y)*(y1 - y))
         + sqrt((x2 - x)*(x2 - x) + (y2 - y)*(y2 - y));
}
 
int main(){
    
    int x1, y1, x2, y2;
    const int H = 34;
    const int W = 55;
    
    do{
        cout << "x1(1 - " << W << "): y1(1 - " << H << "): ";
        cin >> x1 >> y1;
    }while (x1 < 1 || W < x1 || y1 < 1 || H < y1);  
    do{
        cout << "x2(1 - " << W << "): y2(1 - " << H << "): ";
        cin >> x2 >> y2;
        if (x1 == x2 && y1 == y2)
            cout << "Ââåäèòå òî÷êó îòëè÷íóþ îò ïåðâîé:" << endl;            
    }while (x2 < 1 || W < x2 || y2 < 1 || H < y2 || (x1 == x2 && y1 == y2));
    
    int map[H * W] = {};
    
    --x1; --x2;// --y1; --y2;
    y1 = H - y1;
    y2 = H - y2;
    
    map[y1 * W + x1] = 2;
    map[y2 * W + x2] = 2;
    
    if (x1 > x2){
        int buf = x1; x1 = x2; x2 = buf;
            buf = y1; y1 = y2; y2 = buf;
    }
 
    int x = x1;
    int y = y1;
    if (y2 > y1)
        while (x != x2 && y != y2){
                    
            double min;
            if (s(x1, y1, x2, y2, x + 1, y) < s(x1, y1, x2, y2, x, y + 1))
                min = s(x1, y1, x2, y2, x + 1, y);
            else
                min = s(x1, y1, x2, y2, x, y + 1);
                
            if (s(x1, y1, x2, y2, x + 1, y + 1) < min){
                ++x;
                ++y;
            }else if (s(x1, y1, x2, y2, x + 1, y) < s(x1, y1, x2, y2, x, y + 1))
                ++x;
            else
                ++y;
            
            if (x != x2 && y != y2)
                map[y * W + x] = 1;
        }
    else
        while (x != x2 && y != y2){
            
            double min;
            if (s(x1, y1, x2, y2, x + 1, y) < s(x1, y1, x2, y2, x, y - 1))
                min = s(x1, y1, x2, y2, x + 1, y);
            else
                min = s(x1, y1, x2, y2, x, y - 1);
                
            if (s(x1, y1, x2, y2, x + 1, y - 1) < min){
                ++x;
                --y;
            }else if (s(x1, y1, x2, y2, x + 1, y) < s(x1, y1, x2, y2, x, y - 1))
                ++x;
            else
                --y;
            
            if (x != x2 && y != y2)
                map[y * W + x] = 1;
        }
        
    for (int i = 0; i < H; ++i){
        for (int j = 0; j < W; ++j)
            cout << map[i * W + j];
        cout << endl;
    }   
}