#include <bits/stdc++.h>

using namespace std;

const long long mod=998244353;
long long n,a[5001],k,x,dp1[5001][5001],dp2[5001][5001],v[5001];
long long factorial(long long n)
{
    long long nr=1;
    for(long long i=1;i<=n;i++)
        nr=(nr*i)%mod;
    return nr;
}
void rotire()
{
    for(long long i=1;i<=n;i++)
    {
        long long j=1;
        while(j<=n and a[j]>=i)
            j++;
        v[i]=j-1;
    }
    for(int i=1;i<=n;i++)
        a[i]=v[i];
}
int main()
{
    cin>>n;
    for(long long i=1;i<=n;i++)
        cin>>a[i];
    reverse(a+1,a+n+1);
    for(long long i=1;i<=n;i++)
    {
        if(a[i]>=i)
            k=i;
        else
            break;
    }
    x=a[k+1];
    dp1[0][x]=1;
    for(long long i=0;i<=k;i++)
    {
        for(long long j=0;j<=x;j++)
        {
            if(j>0)
                dp1[i+1][j-1]+=dp1[i][j]*j;
            dp1[i+1][j]+=dp1[i][j]*(a[i+1]-j);
        }
    }
    rotire();
    x=a[k+1];
    dp2[0][x]=1;
    for(long long i=0;i<=k;i++)
    {
        for(long long j=0;j<=x;j++)
        {
            if(j>0)
                dp2[i+1][j-1]+=dp2[i][j]*j;
            dp2[i+1][j]+=dp2[i][j]*(a[i+1]-j);
        }
    }
    long long total=0;
    total=dp1[k][0]+dp2[k][0]-factorial(k);
    total%=mod;
    while(total<0)
        total+=mod;
    cout<<k<<" "<<total;
    return 0;
}