def nPr(n, r, mod): if ( r<0 or r>n ): return 0 return g1[n] * g2[n-r] % mod def nCr(n, r, mod): if ( r<0 or r>n ): return 0 r = min(r, n-r) return (g1[n] * g2[r] % mod) * g2[n-r] % mod import sys input = sys.stdin.readline MOD = 998244353 N = int(input()) A = list(map(int, input().split())) g1 = [1, 1] # 元テーブル g2 = [1, 1] #逆元テーブル inverse = [0, 1] #逆元テーブル計算用テーブル fact = [1, 1] for i in range( 2, 2*10**5 + 1 ): g1.append( ( g1[-1] * i ) % MOD ) inverse.append( ( -inverse[MOD % i] * (MOD//i) ) % MOD ) g2.append( (g2[-1] * inverse[-1]) % MOD ) fact.append( (fact[-1] * i) % MOD ) A.append(0) ans = dict() for i in range(N): a = A[i] if i not in ans: ans[i] = nCr(i+a, a, MOD) a2 = A[i+1] for aa in range(a2, a): ans[i+1] = nCr(i+1+aa, aa, MOD) #print(ans) ans2 = 0 bk = -1 for i in range(N, -1, -1): if bk==A[i]: continue ans2 += ans[i] ans2 %= MOD bk = A[i] print(ans2)