/*
#include <boost/functional/hash.hpp>
#include <boost/multiprecision/cpp_int.hpp>
typedef boost::multiprecision::cpp_int mlint;
*/
#include <bits/stdc++.h>
using namespace std;

#define esc(x) cout << (x) << endl,exit(0)
#define inf(T) (numeric_limits<T>::max() / 2)
#define each(i,v) for(auto i = begin(v),end_copy = end(v); i != end_copy; ++i)
#define reach(i,v) for(auto i = rbegin(v),rend_copy = rend(v); i != rend_copy; ++i)
#define urep(i,s,t) for(i64 i = (s),t_copy = (t); i <= t_copy; ++i)
#define drep(i,s,t) for(i64 i = (s),t_copy = (t); i >= t_copy; --i)
#define rep(i,n) urep(i,0,(n) - 1)
#define rep1(i,n) urep(i,1,(n))
#define rrep(i,n) drep(i,(n) - 1,0)
#define all(v) begin(v),end(v)
#define rall(v) rbegin(v),rend(v)
#define vct vector
#define prique priority_queue
#define umap unordered_map
#define uset unordered_set
#define l_bnd lower_bound
#define u_bnd upper_bound
#define rsz resize
#define ers erase
#define emp emplace
#define emf emplace_front
#define emb emplace_back
#define pof pop_front
#define pob pop_back
#define mkp make_pair
#define mkt make_tuple
#define fir first
#define sec second

typedef int64_t i64;
typedef double db;
typedef vct<vct<i64>> mat;
typedef pair<int,int> pii;
typedef pair<i64,int> pli;
typedef pair<db,int> pdi;
typedef map<int,int> mpii;


#ifdef Local
#define debug 1
#define print(x) cout << #x << " = " << (x) << " (L" << __LINE__ << ")" << endl
#include <sys/time.h>
struct timeval start_time;
#endif

void setup(short prec = 10) {
    cin.tie(0);
    ios::sync_with_stdio(false);
    cout << fixed << setprecision(prec);

    #ifdef Local
    gettimeofday(&start_time,NULL);
    cout << "-- Execution At Local ---\n" << "\n<Standard Output>\n\n";
    auto print_atexit = []() {
        struct timeval end_time;
        gettimeofday(&end_time,NULL);
        double runtime_millis = 1000.0 * (end_time.tv_sec - start_time.tv_sec) + 0.001 * (end_time.tv_usec - start_time.tv_usec);
        cout << "\n----------------\n";
        cout << "\nExec time : " << runtime_millis << " ms\n\n";
    };
    atexit(print_atexit);
    #endif
}

template <class T, class U> ostream& operator << (ostream& s, pair<T,U> p) { return s << p.fir << " " << p.sec; }
template <class T> ostream& operator << (ostream& s, vector<T>& v) { for(auto i = v.begin(); i != v.end(); i++) { if(i != begin(v)) s << " "; s << *i; } return s; }

template <class T> void hash_combine(size_t &seed, T const &key) {
    hash<T> hasher;
    seed ^= hasher(key) + 0x9e3779b9 + (seed << 6) + (seed >> 2);
}

namespace std {
    template <class T, class U> struct hash<pair<T,U>> {
        size_t operator()(pair<T,U> const &p) const {
            size_t seed(0);
            hash_combine(seed,p.first);
            hash_combine(seed,p.second);
            return seed;
        }
    };
}

bool odd(i64 x) { return x & 1; }
bool even(i64 x) { return !odd(x); }
bool parity(i64 x, i64 y) { return odd(x) ^ even(y); }
bool bit(i64 n, uint8_t e) { return n >> e & 1; }
int ilog(i64 x, uint64_t b) { if(x) return ilog(x / b,b) + 1; return -1; }
i64 qceil(i64 x, i64 y) { return x > 0 ? (x - 1) / y + 1 :  x / y; }
i64 gcd(i64 a, i64 b) { if(!b) return abs(a); if(a % b) return gcd(b, a % b); return b; }
i64 lcm(i64 a, i64 b) { if(a & b) return a / gcd(a, b) * b; return 0; }
template <class T, class U> bool chmax(T& m, U x) { if(m < x) { m = x; return 1; } return 0; }
template <class T, class U> bool chmin(T& m, U x) { if(m > x) { m = x; return 1; } return 0; }
template <class T> T cmprs(T& v) {
    T tmp = v;
    sort(all(tmp));
    tmp.erase(unique(all(tmp)),end(tmp));
    each(i,v) *i = l_bnd(all(tmp),*i) - begin(tmp) + 1;
    return v;
}

const int dir[8][2] = { {1,0},{0,1},{-1,0},{0,-1},{1,1},{-1,1},{-1,-1},{1,-1} };
const i64 mod = 1e9 + 7;
const db eps = 1e-15;

int n,x,y,z,w,mask;
i64 dp[41][1 << 17] = {};
i64 pw[41];
i64 ans=0;

void erabit(int &s) {
	s &= (1 << w) - 1;
}

bool check(int s) {
	return (s & mask) == mask;
}

int main() {
	setup();
	cin>>n>>x>>y>>z;
	w=x+y+z;
	mask = (1 << z-1) + (1 << y+z-1) + (1 << w-1);
	pw[0] = 1;
	rep(i,40) pw[i + 1] = pw[i] * 10 % mod;
	dp[0][0] = 1;
	rep(i,n) {
		i64 *crt = dp[i],*nxt = dp[i + 1];
		rep(s,1 << w) {
			int cp = (s << 1) + 1;
			erabit(cp);
			rep(k,10) {
				if(check(cp)) (ans+=pw[n-1-i]*crt[s])%=mod;
				else (nxt[cp]+=crt[s])%=mod;
				erabit(cp<<=1);
			}
		}
	}
	esc(ans);
}

Measured memory_bytes: 41216