#include<cstdio>

#define N 500050

using namespace std;

struct lhy{
    int x,y,next;
}edge[N];

int tot,n,k,A,B,Ans,ans;
int son[N],f[N][11],g[N][11];

inline void add(int x,int y)
{
    edge[++tot].x=x;edge[tot].y=y;edge[tot].next=son[x];son[x]=tot;
}

void dfs(int x)
{
    f[x][0]=1;
    for(int i=1;i<=k;i++)f[x][i]=0;
    for(int i=son[x];i;i=edge[i].next)
    {
        dfs(edge[i].y);
        for(int j=0;j<k;j++)
            f[x][j+1]+=f[edge[i].y][j];
    }
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d%d%d",&n,&k,&A,&B);
        tot=0;
        for(int i=1;i<=n;i++)son[i]=0;
        for(int i=2;i<=n;i++)
            add((1ll*A*i%(i-1)+B)%(i-1)+1,i);
        dfs(1);
        Ans=0;
        for(int x=1;x<=n;x++)
        {
            for(int i=son[x];i;i=edge[i].next)
            {
                int y=edge[i].y;
                for(int j=0;j<=k;j++)
                    g[y][j]=0;
                g[y][1]=f[x][0];
                for(int j=1;j<k;j++)
                    g[y][j+1]+=(f[x][j]-f[y][j-1])+g[x][j];
            }
            ans=0;
            for(int j=0;j<=k;j++)
                ans+=g[x][j]+f[x][j];
            Ans^=ans;
        }
        printf("%d\n",Ans); 
    }
}