#include <bits/stdc++.h>

using namespace std;
int q;
struct point
{
    double x;
    double y;
};
point vt(point a,point b)
{
    point vt;
    vt.x=b.x-a.x;
    vt.y=b.y-a.y;
    return vt;
}
double scalar(point a,point b)
{
    double scalar=0;
    scalar=sqrt((b.x-a.x)*(b.x-a.x)+(b.y-a.y)*(b.y-a.y));
    return scalar;
}
double mul1(point a,point b)
{
    double mul=a.x*b.x+a.y*b.y;
    return mul;
}
point mul2(point a,double b)
{
    point c;
    c.x=a.x*b;
    c.y=a.y*b;
    if (c.x==-0) c.x=0;
    if (c.y==-0) c.y=0;
    return c;
}
int main()
{
    point p1,p2;
    cin>>p1.x>>p1.y;
    cin>>p2.x>>p2.y;
    cin>>q;
    while (q--)
    {
        point p;
        cin>>p.x>>p.y;
        point v=vt(p1,p);
        point u=vt(p1,p2);
        //cout<<v.x<< " "<<v.y<<endl;
        //cout<<u.x<< " "<<u.y<<endl;
        double su=scalar(p1,p2);
        //cout<<su<<endl;
        point ans;
        ans=mul2(u,mul1(u,v)/(su*su));
        cout<<setprecision(10)<<fixed<<ans.x+p1.x<<" "<<ans.y+p1.y<<endl;
    }
    return 0;
}

