#include using namespace std; using ll = long long; using ld = long double; using vl = vector; template using vc = vector; template using vvc = vector>; #define eb emplace_back #define all(x) (x).begin(), (x).end() #define rep(i, n) for (ll i = 0; i < (n); i++) #define repr(i, n) for (ll i = (n)-1; i >= 0; i--) #define repe(i, l, r) for (ll i = (l); i < (r); i++) #define reper(i, l, r) for (ll i = (r)-1; i >= (l); i--) #define repa(i,n) for (auto& i: n) template inline bool chmax(T1 &a, const T2 &b) {if (a inline bool chmin(T1 &a, const T2 &b) {if (b void verr(const set &st) { repa(a, st) cerr << a << " "; cerr << endl; } template void verr(const map &mp) { repa(a, mp) cerr << "{" << a.first << ", " << a.second << "} "; cerr << endl; } template void verr(const vector>& a, const N& n) { rep(i, n) cerr << "{" << a[i].first << ", " << a[i].second << "} "; cerr << endl; } template void verr(const vector& a, const N& n) { rep(i, n) cerr << a[i] << " "; cerr << endl; } ll dbgt = 1; void err() { cerr << "passed " << dbgt++ << endl; } template void err(H&& h,T&&... t){ cerr<< h << (sizeof...(t)?" ":"\n") << flush; if(sizeof...(t)>0) err(forward(t)...); } #endif const ll INF = 4e18; const ld EPS = 1e-11; const ld PI = acos(-1.0L); const ll MOD = 1e9 + 7; // const ll MOD = 998244353; //--------------------------------------------------------------------------------// template struct Matrix { vector> mat; int row, col; Matrix(int r, int c, T v = 0) : mat(r, vector(c, v)), row(r), col(c){}; inline const vector &operator[](int k) const { return mat[k]; } inline vector &operator[](int k) { return mat[k]; } Matrix &operator+=(const Matrix &A) { for (int i = 0; i < row; i++) for (int j = 0; j < col; j++) (*this)[i][j] += A[i][j]; return (*this); } Matrix &operator-=(const Matrix &A) { for (int i = 0; i < row; i++) for (int j = 0; j < col; j++) (*this)[i][j] -= A[i][j]; return (*this); } Matrix &operator*=(const Matrix &A) { vector> res(row, vector(A.col, 0)); for (int i = 0; i < row; i++) for (int j = 0; j < A.col; j++) for (int k = 0; k < col; k++) (res[i][j] += mat[i][k] * A[k][j]) %= 10; mat.swap(res); (*this).col = A.col; return (*this); } Matrix &operator^=(long long p) { Matrix res = Matrix(row, col); for (int i = 0; i < row; i++) res[i][i] = 1; while(p > 0){ if (p & 1) res *= *this; *this *= *this; p >>= 1; } mat.swap(res.mat); return (*this); } Matrix operator+(const Matrix &B) const { return (Matrix(*this) += B); } Matrix operator-(const Matrix &B) const { return (Matrix(*this) -= B); } Matrix operator*(const Matrix &B) const { return (Matrix(*this) *= B); } Matrix operator^(const long long k) const { return (Matrix(*this) ^= k); } friend ostream &operator<<(ostream &os, Matrix &p) { for (int i = 0; i < p.row; i++) { os << "["; for (int j = 0; j < p.col; j++) { os << p[i][j] << (j + 1 == p.col ? "]\n" : ","); } } return (os); } }; int main() { ll p, q, r, K; cin >> p >> q >> r >> K; Matrix A(3, 3); rep(i, 3) A[0][i] = 1; A[1][0] = A[2][1] = 1; Matrix X(3, 1); X[0][0] = r % 10, X[1][0] = q % 10, X[2][0] = p % 10; A ^= K - 3; X = A * X; cout << X[0][0] << endl; }