#include <bits/stdc++.h>
using namespace std;
const int dx[]={0,1,0,-1,1,-1,-1,1};
const int dy[]={-1,0,1,0,1,1,-1,-1};
const int INF = 1<<30;
const double EPS = 1e-15;
#define PB push_back
#define mk make_pair
#define fi first
#define se second
#define ll long long
#define reps(i,j,k) for(int i = (j); i < (k); i++)
#define rep(i,j) reps(i,0,j)
#define MOD 1000000007
typedef pair<int,int> Pii;
typedef vector<int> vi;
typedef vector<vi> vvi;
struct Point{
    double x,y;
    Point(){}
    Point(double _x,double _y){
        x = _x;
        y = _y;
    }
    Point operator + (Point p){
        return Point(x+p.x,y+p.y);
    }
  
    Point operator - (Point p){
        return Point(x-p.x,y-p.y);
    }
  
    Point operator * (double k){
        return Point(k*x,k*y);
    }
     
    Point operator * (Point p){
        return Point(p.x*x-p.y*y,p.x*y+p.y*x);
    }

    Point operator / (double k){
        return Point(x/k,y/k);
    }
    
    Point operator = (Point p){
        x = p.x;
        y = p.y;
        return *this;
    }

    bool operator<(const Point &p)const{
        return x!=p.x?x<p.x:y<p.y;
    }

};
double dot(Point p, Point q){
    return p.x*q.x+p.y*q.y;
}
  
double cross(Point p,Point q){
    return p.x*q.y-q.x*p.y;
}
  
double norm(Point p){
    //?????????????????§?????????2???
    return p.x*p.x+p.y*p.y;
}
  
double pabs(Point p){
    //?????????????????§??????
    return sqrt(norm(p));
}
Point UnitVector(Point a){
    //??????????????????
    return a/pabs(a);
}
typedef vector<Point> Polygon;
struct Line : Polygon{
    Line(){}
    Line(Point p,Point q){
        PB(p);
        PB(q);  
    }
};
struct Circle{
    Point p; //?????????
    double r; //??????
    Circle(){}
    Circle(Point _p,double _r){
        p = _p;
        r = _r;
    }
};
int ccw(Point a,Point b,Point c){
    b = b-a;
    c = c-a;
    if(cross(b,c) > 0)return +1; //a???b ??§???????¨?????????????c ccw
    if(cross(b,c) < 0)return -1; //a???b ??§????¨?????????????c ccw
    if(dot(b,c) < 0)return +2; // c--a--b;
    if(pabs(b)<pabs(c))return -2; //a--b--c
    return 0; //a--c--b
}

Point projection(Line l, Point p){
    double t = dot(p-l[0],l[0]-l[1]) / norm(l[0]-l[1]);
    return l[0] + (l[0]-l[1])*t;
}

bool intersectSS(Line l,Line m){
    //l??¨m??????????????????????????????
    return (ccw(l[0],l[1],m[0])*ccw(l[0],l[1],m[1])<EPS)&&
            (ccw(m[0],m[1],l[0])*ccw(m[0],m[1],l[1])<EPS);
}

bool intersectSP(Line l, Point p){
    //??????l??¨???p?????????????????????
    return !ccw(l[0],l[1],p);
}

double distanceLP(Line l,Point p){
    //??????l??¨???p????????¢
    p = p - l[0];
    l[1] = l[1] - l[0];
    if(intersectSP(l,p))return 0.0;
    return fabs(cross(p,l[1])/pabs(l[1]));
}

bool intersectLL(Line l,Line m){
    return fabs(cross(l[1]-l[0],m[1]-m[0])) > EPS || //????????§??????
            fabs(cross(l[1]-l[0],m[0]-l[0])) < EPS; //????????????
}

bool intersectCP(Circle c,Point p){
    //???c??¨???p???????????????
    //true->???c????????¨??????p????????¨????????????.
    return fabs(pabs(p-c.p)-c.r)<EPS;
}

bool intersectCL(Circle c,Line l,Point *p,Point *q){
    
}

int intersectCC(Circle a,Circle b,Point *p,Point *q){
    double dist = pabs(a.p-b.p);
    if(dist < fabs(a.r-b.r) || dist > a.r+b.r){
        //????????????????????´??????????????¨???
        //????????????????????¨???????????¨???
        *p = Point(1e15,1e15);
        *q = Point(1e15,1e15);
        return 0;
    }
    else if(fabs(dist-a.r-b.r) < EPS){
        //??????1?????§????????£???????????¨???
        Point way;
        way = b.p-a.p;
        *p = a.p+UnitVector(way)*a.r;
        *q = *p;
        return 1;
    }
    else{
        //??????2?????§????????£???????????¨???
        if(a.r < b.r){
            swap(a.r,b.r);
            swap(a.p,b.p);
        }
        double s = (a.r + b.r + dist) / 2.0 ;
        double area = sqrt(s * (s - a.r) * (s - b.r) * (s - dist));
        double h = 2.0 * area / dist;

        Point  v = b.p-a.p;
        v = UnitVector(v);
        Point  m = a.p + (v*sqrt(a.r*a.r-h*h));
        Point  n = v * Point(0.0,1.0);
        
        *p = m + (n * h);
        *q = m - (n * h);
        return 2;
    }
}

double TriangleArea(Point a,Point b,Point c){
    b = b - a;
    c = c - a;
    return cross(b,c)*0.5; //?????????????????¨??????????????§??¨???
}
 
double VertexArea(Polygon v){
    double subArea = 0.0;
    rep(i,v.size()-1){
        subArea += TriangleArea(v[0],v[i],v[i+1]);
    }
    return subArea; //?????????????????¨??????????????§??¨???
}

Point crosspoint(Line l,Line m){
    //??????l??¨??????m?????????
    if(intersectSS(l,m) == false)return Point(INF,INF);
    return m[0]+(m[1]-m[0])*cross(l[1]-l[0],l[1]-m[0])/cross(l[1]-l[0],m[1]-m[0]);
}
 
Polygon ConvexHull(Polygon v){
    //??????????????????
    int n = v.size();
    int k = 0;
    sort(v.begin(),v.end());
    Polygon e(2*n);
    for(int i=0;i<n;e[k++]=v[i++])
        while(k>=2&&ccw(e[k-2],e[k-1],v[i])<=0)k--;
    for(int i=n-2,t=k+1;i>=0;e[k++]=v[i--])
        while(k>=t&&ccw(e[k-2],e[k-1],v[i])<=0)k--;
    e.resize(k-1);
    return e;
}

int main(){
    double x,y;
    Line l;
    cin >> x >> y;
    l.PB(Point(x,y));
    cin >> x >> y;
    l.PB(Point(x,y));
    int Q;
    cin >> Q;
    while(Q--){
        cin >> x >> y;
        Point ret = projection(l,Point(x,y));
        printf("%.8lf %.8lf\n",ret.x,ret.y);
    }
    return 0;
}

Measured memory_bytes: 1256