#include <bits/stdc++.h>
using namespace std;
#define fir first
#define sec second

#define lson(i) (tr[i].ls)
#define rson(i) (tr[i].rs)

#define FIO(file) freopen(file".in","r",stdin), freopen(file".out","w",stdout)

typedef long long ll;
typedef double db;
typedef long double ldb;
typedef unsigned int uint;
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
typedef __int128_t i128;
typedef __uint128_t u128;
ll LIM,m,q,mod;
ll f[15],g[15];
vector<ll> val;
ll fd[1000000],gd[1000000];
inline ll getf(int x){int pos=lower_bound(val.begin(),val.end(),x)-val.begin();if(pos<0||pos>=val.size()||val[pos]!=x)return 0;else return fd[pos];}
inline ll getg(int x){int pos=lower_bound(val.begin(),val.end(),x)-val.begin();if(pos<0||pos>=val.size()||val[pos]!=x)return 0;else return gd[pos];}
int main()
{
    cin>>f[0]>>g[0]>>q>>mod;
    LIM=1e15;
    for(int i=1;i<=13;i++)
    {
        f[i]=max((ll)i,g[i/2]+g[i/3]+g[i/5]+g[i/7]);
        g[i]=max((ll)i,f[i/2]+f[i/3]+f[i/4]+f[i/5]);
    }
    for(int i=0;i<=13;i++)f[i]%=mod,g[i]%=mod;
    val.push_back(0);
    for(ll t=1;t<=min(13ll,LIM);t++)
        for(ll z=t;z<=LIM;z*=2)
            for(ll t=z;t<=LIM;t*=3)
                for(ll z=t;z<=LIM;z*=5)
                    for(ll t=z;t<=LIM;t*=7)
                        val.push_back(t);
    sort(val.begin(),val.end());
    val.resize(unique(val.begin(),val.end())-val.begin());
    for(int i=0;i<val.size();i++)
    {
        if(val[i]<=13){
            if(val[i]==0)fd[i]=f[0],gd[i]=g[0];
            else fd[i]=(f[val[i]]-f[val[i]-1]+mod)%mod,gd[i]=(g[val[i]]-g[val[i]-1]+mod)%mod;
        }
        else 
        {
            fd[i]=(val[i]%2==0?getg(val[i]/2):0)+(val[i]%3==0?getg(val[i]/3):0)+(val[i]%5==0?getg(val[i]/5):0)+(val[i]%7==0?getg(val[i]/7):0);fd[i]%=mod;
            gd[i]=(val[i]%2==0?getf(val[i]/2):0)+(val[i]%3==0?getf(val[i]/3):0)+(val[i]%4==0?getf(val[i]/4):0)+(val[i]%5==0?getf(val[i]/5):0);gd[i]%=mod;
        }
    }
    for(int i=1;i<val.size();i++)
        fd[i]+=fd[i-1],gd[i]+=gd[i-1],fd[i]%=mod,gd[i]%=mod;
    for(int i=1;i<=q;i++)
    {
        ll x;
        cin>>x;
        int pos=upper_bound(val.begin(),val.end(),x)-val.begin()-1;
        //cout<<"!"<<pos<<'\n';
        cout<<fd[pos]<<' '<<gd[pos]<<'\n';
    }
}