S=input()
v=[0]*10
mod=998244353
for i in range(len(S)):
    v[int(S[i])]+=1
result=0
u=[1]*10**6
u2=[1]*10**6
for i in range(1,10**6):
    u[i]=u[i-1]*i
    u[i]%=mod
    u2[i]=pow(u[i],-1,mod)
def ncm(x,y):
    ans=u[x]*u2[y]
    ans%=mod
    ans*=u2[x-y]
    ans%=mod
    return ans
for x in range(1,10):
    if v[x]==0:
        continue
    h=v[:]
    h[x]-=1
    rest=sum(h)
    w=1
    for y in range(10):
        n=h[y]
        w*=ncm(rest,n)
        w%=mod
        rest-=n
    result+=w
    result%=mod
print(result)