#include <bits/stdc++.h>

using namespace std;

int n, v[5005], a[5005], dp1[5005][5005], dp2[5005][5005], t, k;

long long fact (int n)
{
    int p=1;
    for (int i=2; i<=n; i++)
    {
        p*=i;
        p%=998244353;
    }
    return p%998244353;
}

void pd (int v[], int dp[][5005])
{
    int x=v[k+1];
    dp[0][x]=1;
    for (int i=0; i<=k; i++)
    {
        for (int j=0; j<=v[i+1] && j<=x; j++)
        {
            dp[i+1][j]+=(1LL*(v[i+1]-j)*dp[i][j]%998244353);
            dp[i+1][j]%=998244353;
            if (j>0)
            {
                dp[i+1][j-1]+=(1LL*j*dp[i][j]%998244353);
                dp[i+1][j-1]%=998244353;
            }
        }
    }
}

void rotire ()
{
    for (int i=1; i<=n; i++)
    {
        for (int j=1; j<=n; j++)
        {
            if (v[j]>=i)
                a[i]++;
            else
                break;
        }
    }
}

int main()
{
    //9 1 1 2 5 6 6 6 7 9
    cin >> n;
    for (int i=1; i<=n; i++)
        cin >> v[n-i+1];
    for (int i=1; i<=n; i++)
    {
        bool ok=true;
        for (int j=1; j<=i; j++)
        {
            if (v[j]<i)
                ok=false;
        }
        if (!ok)
            break;
        else
            k=i;
    }
    rotire();
    pd(v,dp1);
    pd(a,dp2);
    //cout << dp1[k][0] << " " << dp2[k][0];
    long long rez=(dp1[k][0]%998244353+dp2[k][0]%998244353-1LL*fact(k)%998244353)%998244353;
    cout << k << " " << 1LL*rez%998244353;
    return 0;
}