n=int(input())
p=[-1]+list(map(int,input().split()))
mod=998244353


#############################
#############
cnb_max=10**6
#############

kai=[1]*(cnb_max+1)
rkai=[1]*(cnb_max+1)
for i in range(cnb_max):
    kai[i+1]=kai[i]*(i+1)%mod

rkai[cnb_max]=pow(kai[cnb_max],mod-2,mod)
for i in range(cnb_max):
    rkai[cnb_max-1-i]=rkai[cnb_max-i]*(cnb_max-i)%mod

def cnb(x,y):
    if y>x:
        return 0
    if x<0:return 0
    if y<0:return 0
    return (kai[x]*rkai[y]%mod)*rkai[x-y]%mod


def inv(n):
    return kai[n-1]*rkai[n]%mod

##################################




ind=[10**10]*(n+4)
for i in range(1,n+1):
    ind[p[i]]=i

l=ind[0]
r=ind[0]
ans=1
num=0
for i in range(1,n):
    if l<=ind[i]<=r:
        ans*=num
        ans%=mod
        num-=1
    else:
        if r<ind[i]:
            num+=ind[i]-r-1
            r=ind[i]
        else:
            num+=l-ind[i]-1
            l=ind[i]
print(ans)