#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=1e6+6;
int n,K,Q,T,dfn[N];
bool cmp(const int x,const int y)
{
    return dfn[x]<dfn[y];
}
struct Tree{
    int number;
    int p[N],rt;
    int fa[N][20];
    struct Edge{
        int v,nxt;
    }e[N];
    int h[N],tot1;
    void addEdge(int u,int v)
    {
        tot1++;
        e[tot1].v=v;e[tot1].nxt=h[u];
        h[u]=tot1;
    }
    void input()
    {
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&p[i]);
            if(p[i]>=1)addEdge(p[i],i),fa[i][0]=p[i];
            else rt=i,fa[i][0]=0;
        }
    }
    int dfn[N],idx,dep[N],rnk[N];
    void dfs(int u)
    {
        dfn[u]=++idx;rnk[idx]=u;
        for(int i=h[u];i;i=e[i].nxt)
        {
            int v=e[i].v;dep[v]=dep[u]+1;
            dfs(v);
        }
    }
    int a[N];
    LL d[N];
    int lca(int u,int v)
    {
        if(dep[u]<dep[v])
        {
            swap(u,v);
        }
        for(int j=19;j>=0;j--)
        {
            if(dep[u]-dep[v]>>j&1)
            {
                u=fa[u][j];
            }
        }
        if(u==v)
        {
            return u;
        }
        for(int j=19;j>=0;j--)
        {
            if(fa[u][j]!=fa[v][j])
            {
                u=fa[u][j];
                v=fa[v][j];
            }
        }
        return fa[u][0];
    }
    void init()
    {
        dfs(rt);
        for(int j=1;j<=19;j++)
        {
            for(int i=1;i<=n;i++)
            {
                fa[i][j]=fa[fa[i][j-1]][j-1];
            }
        }
        for(int i=1;i<=K;i++)
        {
            a[i]=i;
        }
        int tp=0;
        for(int i=1;i<=n;i++)
        {
            if(rnk[i]<=K)
            {
                a[++tp]=rnk[i];
            }
        }
        for(int i=1;i<K;i++)
        {
            cout<<"? "<<a[i]<<' '<<a[i+1]<<endl;
        }
    }
    void work()
    {
        d[a[1]]=0;
        for(int i=1;i<K;i++)
        {
            LL X,Y,Z,W;
            scanf("%lld%lld%lld%lld",&X,&Y,&Z,&W);
            int x=lca(a[i],a[i+1]);
            if(number==1)
            {
                d[x]=d[a[i]]-X;
                d[a[i+1]]=d[x]+Y;
            }else{
                d[x]=d[a[i]]-Z;
                d[a[i+1]]=d[x]+W;
            }
        //    cout<<a[i]<<' '<<x<<' '<<a[i+1]<<endl;
        //    cout<<d[a[i]]<<' '<<d[x]<<' '<<d[a[i+1]]<<endl;
        }
    /*    for(int i=1;i<=n;i++)
        {
            cout<<d[i]<<' ';
        }
        cout<<endl;*/
    }
    LL query(int u,int v)
    {
        return d[u]+d[v]-2*d[lca(u,v)];
    }
}t1,t2;
int main()
{
    int u,v;
    scanf("%d%d%d%d",&n,&K,&Q,&T);
    if(Q!=K-1||n>5e5)return 0;
    t1.number=1;t2.number=2;
    t1.input();t2.input();
    for(int i=1;i<=K;i++)
    {
        cout<<i<<' ';
    }
    cout<<endl;
    t1.init();
    cout<<"!\n";
    fflush(stdout);
    t1.work();
    while(T--)
    {
        scanf("%d%d",&u,&v);
        printf("%lld %lld\n",t1.query(u,v),t1.query(u,v));
    }
    fflush(stdout);
    return 0;
}