#x,y=map(int, input().split())
#A=list(map(int, input().split()))
#D=[[]for i in range(N)]
S=input()
mod=998244353
C=[0]*10
for s in S:
  C[int(s)]+=1

mod=998244353;n=10**6
inv_t=[0]+[1]
for i in range(2,n):
  inv_t+=[inv_t[mod%i]*(mod-int(mod/i))%mod]

kai=[1,1]
rev_kai=[1,inv_t[1]]
for i in range(2,n):
  kai.append(kai[-1]*i%mod)
  rev_kai.append(rev_kai[-1]*inv_t[i]%mod)

def cmb(n,r):
  return kai[n]*rev_kai[r]*rev_kai[n-r]%mod

ansa=1
al=sum(C)-C[0]
ansa*=cmb(al-1+C[0],al-1)
for i in range(1,10):
  ansa*=cmb(al,C[i])
  ansa%=mod 
  al-=C[i]
print(ansa)