#pragma region
#define _USE_MATH_DEFINES
#include <iostream>
#include <string>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <vector>
#include <map>
#include <queue>
#include <stack>
#include <set>
#include <list>
#include <iomanip>
#include <cstdint>
#include <bitset>
#include <sstream>
#include <regex>
#include <fstream>
#include <array>
#include <atcoder/all>

using namespace atcoder;
using mint = modint1000000007;
//using mint = modint998244353;
//using mint = modint;

#pragma region using
using namespace std;
typedef long long ll;
using vi = vector<int>;
using vvi = vector<vi>;
using vl = vector<ll>;
using vvl = vector<vl>;
using vmint = vector<mint>;
using vvmint = vector<vmint>;
using pint = pair<int, int>;
using pll = pair<ll, ll>;
using vpint = vector<pint>;
using vvpint = vector<vpint>;
using vd = vector<double>;
using vvd = vector<vd>;
//#define rep(i, s, e) for (int(i) = (s); (i) < (e); ++(i))
#define rep(i, e) for (int(i) = 0; (i) < (e); ++(i))
#define rrep(i, s) for (int(i) = (s) - 1; (i) >= 0; --(i))
#define all(x) x.begin(),x.end()
#define rall(x) x.rbegin(),x.rend()
#pragma endregion
#pragma region UnionFind
struct UnionFind
{
	vector<int> par;

	UnionFind(int n) : par(n, -1) {}
	void init(int n) { par.assign(n, -1); }

	int root(int x)
	{
		if (par[x] < 0) return x;
		else return par[x] = root(par[x]);
	}

	bool issame(int x, int y)
	{
		return root(x) == root(y);
	}

	bool merge(int x, int y)
	{
		x = root(x); y = root(y);
		if (x == y) return false;
		if (par[x] > par[y]) swap(x, y);
		par[x] += par[y];
		par[y] = x;
		return true;
	}

	int size(int x)
	{
		return -par[root(x)];
	}
};
#pragma endregion
#pragma region GCD
ll gcd(ll a, ll b)
{
	if (b == 0)return a;
	return gcd(b, a%b);
}
#pragma endregion
#pragma region LCM
ll lcm(ll a, ll b)
{
	return a / gcd(a, b) * b;
}
#pragma endregion
#pragma region chmin
template<class T> inline bool chmin(T& a, T b)
{
	if (a > b)
	{
		a = b;
		return true;
	}
	return false;
}
#pragma endregion
#pragma region chmax
template<class T> inline bool chmax(T& a, T b)
{
	if (a < b)
	{
		a = b;
		return true;
	}
	return false;
}
#pragma endregion
#pragma region Dijkstra
vl dijkstra(vector<vector<pair<int, ll>>> v, int s)
{
	ll INF = 1e18;
	int MAX = 1e6;
	vl res(MAX, INF);
	priority_queue<pair<ll, int>, vector<pair<ll, int>>, greater<pair<ll, int>>> q;
	q.push({ 0,s });
	while (!q.empty())
	{
		int now;
		ll d;
		tie(d, now) = q.top();
		q.pop();
		if (!chmin(res[now], d))continue;
		for (auto p : v[now])
		{
			int next;
			ll c;
			tie(next, c) = p;
			if (res[next] <= res[now] + c)continue;
			q.push({ res[now] + c,next });
		}
	}
	return res;
}
#pragma endregion
#pragma region degからrad
double dtor(double deg) { return deg * M_PI / 180; }
#pragma endregion
#pragma region radからdeg
double rtod(double r) { return r * 180 / M_PI; }
#pragma endregion
#pragma region グリッド内チェック
bool out(int x, int y, int h, int w)
{
	if (x < 0 || h <= x || y < 0 || w <= y)return true;
	else return false;
}
#pragma endregion
#pragma endregion

#pragma region 行列累乗
typedef vector<int> vec;
typedef vector<vector<int>> mat;
int m; //遷移行列のサイズ
//dpの更新
vec matmul(vec &dp, mat &mt)
{
	vec ret(m);
	rep(i, m)rep(j, m)
	{
		ret[i] += (mt[i][j] * dp[j]);
		ret[i] %= 10;
	}
	return ret;
}
//遷移行列の更新
mat update(mat &mt)
{
	mat ret(m, vec(m));
	rep(i, m)rep(j, m)rep(k, m)ret[i][j] += (mt[i][k] * mt[k][j]) % 10;
	return ret;
}

void matpow(vec &dp, mat &mt, ll k)
{
	while (k)
	{
		if (k & 1)dp = matmul(dp, mt);
		//正規化
		//double sum = dp[0] + dp[1];
		//dp[0] /= sum;
		//dp[1] /= sum;
		mt = update(mt);
		k /= 2;
	}
}
#pragma endregion

int main()
{
	int p, q, r; ll k;
	cin >> p >> q >> r >> k;
	p %= 10, q %= 10, r %= 10;
	m = 3;
	--k;
	vvi v = { {0,1,0},{0,0,1},{1,1,1} };
	vi a = { p,q,r };
	matpow(a, v, k);
	cout << a[0] % 10 << endl;
}