#include <cstdio>
#include <algorithm>
#include <numeric>
#include <cstring>
#include <queue>
typedef long long LL;
typedef unsigned long long ULL;
typedef std::pair<int, int> pii;
#define MP std::make_pair
#define fi first
#define se second
#define clr(f, n) memset(f, 0, (n) << 2)
#define cpy(f, g, n) memcpy(f, g, (n) << 2)
int read()
{
	int s = 0, f = 1;
	char c = getchar();
	while (!(c >= '0' && c <= '9'))
		f ^= (c == '-'), c = getchar();
	while (c >= '0' && c <= '9')
		s = s * 10 + (c ^ 48), c = getchar();
	return f ? s : -s;
}
const int MAXN = 200005, MAXD = 30;
int Q, K, m, n = 1;
auto fplus = [](int x, int y){ return x + y >= K ? x + y - K : x + y; };
auto fminus = [](int x, int y){ return x >= y ? x - y : x + K - y; };
auto Fplus = [](int &x, int y){ return x = fplus(x, y); };
auto Fminus = [](int &x, int y){ return x = fminus(x, y); };
void swap(int a[], int b[], int n)
{
	for (int i = 0; i < n; i++)
		std::swap(a[i], b[i]);
}
int exgcd(int x, int y, LL &a, LL &b)
{
	if (!y) return a = 1, b = 0, x;
	int d = exgcd(y, x % y, b, a);
	b -= a * (x / y);
	return d;
}
int gcd(int x, int y)
{
	for (; y; std::swap(x, y)) x %= y;
	return x;
}
struct Linear_Basis
{
	int p[MAXD][MAXD];
	void insert(int v[])
	{
		for (int i = m - 1; i >= 0; i--)
			if (v[i])
			{
				// printf("!%d\n", i);
				// if (!p[i][i])
				// {
				// 	cpy(p[i], v, m);
				// 	int d = gcd(K, v[i]);
				// 	for (int j = 0; j < m; j++)
				// 		v[j] = (LL)v[j] * (K / d) % K;
				// }
				// else
				{
					while (true)
					{
						for (int j = 0; j < m; j++)
							Fminus(v[j], (LL)(v[i] / p[i][i]) * p[i][j] % K);
						if (!v[i]) break;
						swap(p[i], v, m);
					}

				}
			}
	}
	void query(int v[])
	{
		// for (int i = 0; i < m; i++)
		// 	for (int j = 0; j < m; j++)
		// 		printf("%d%c", p[i][j], " \n"[j == m - 1]);
		for (int i = m - 1; i >= 0; i--)
			if (v[i] && p[i][i])
			{
				LL a, b;
				int g = gcd(p[i][i], K);
				// 0/0;
				exgcd(p[i][i] / g, K / g, a, b);//gcd(p[i][i], K);
				// printf("%d %lld %lld\n", g, a, b);
				// int g = exgcd(p[i][i], K, a, b);//gcd(p[i][i], K);
				// printf("%d %lld %lld\n", g, a, b);
				int c = (K - v[i] + g - 1) / g;
				a = ((K / g) + c * a % (K / g)) % (K / g);
				for (int j = 0; j < m; j++)
					v[j] = (a * p[i][j] + v[j]) % K;
			}
	}
} lb;
int val[MAXN][MAXD];
int tmp[MAXD], tm[MAXD];
void calc(int f[], int v)
{
	clr(f, m);
	for (int i = 0; v; v /= K, i++)
		f[i] = v % K;
}
void pplus(int f[], int g[])
{
	for (int i = 0; i < m; i++)
		Fplus(f[i], g[i]);
}
int calc(int f[])
{
	int v = 0;
	for (int i = 0, j = 1; i < m; i++, j *= K)
		v += f[i] * j;
	return v;
}
int main()
{
	Q = read(), K = read(), m = read();
	for (int _ = 1; _ <= Q; _++)
	{
		int o = read(), x = read(), y = read();
		if (o == 1)
		{
			calc(tmp, y), cpy(val[++n], tmp, m), pplus(val[n], val[x]);
			for (int i = 0; i < m; i++) Fplus(tmp[i], tmp[i]);
			// printf("%d: ", n);
			// for (int i = 0; i < m; i++) printf("%d ", val[n][i]);
			// putchar('\n');
			lb.insert(tmp);
		}
		else if (o == 2)
		{
			calc(tmp, read());
			for (int i = 0; i < m; i++) tm[i] = fplus(tmp[i], tmp[i]);
			lb.insert(tm);
			for (int i = 0; i < m; i++)
				Fplus(tmp[i], fplus(val[x][i], val[y][i]));
			lb.insert(tmp);
		}
		else
		{
			for (int i = 0; i < m; i++)
				tmp[i] = fplus(val[x][i], val[y][i]);
			lb.query(tmp);
			printf("%d\n", calc(tmp));
		}
	}
	return 0;
}