#include <iostream>
#include <fstream>
#include <algorithm>

using namespace std;
#define MaxN 5000
#define MOD 998244353

int a[MaxN], b[MaxN];
long long int dp[MaxN][MaxN];


long long int dinamic(int a[], int k, int x)
{
    dp[0][x]=1;
    int i, j;
    for(i=0; i<k; i++)
    {
        for(j=x; j>=0; j--)
        {
            dp[i+1][j]+=(a[i+1]-j)*dp[i][j];
            dp[i+1][j]%=MOD;
            if(j>0)
            {
                dp[i+1][j-1]+=j*dp[i][j];
                dp[i+1][j-1]%=MOD;
            }

        }
    }
    return dp[k][0];
}

long long int fact(int k)
{
    long long int prod=1;
    while(k>0)
    {
        prod=prod*k;
        prod%=MOD;
        k--;
    }
    return prod;
}
int main()
{
    int n, i, j, k, x=0;
    long long int res;
    cin>>n;
    for(i=1; i<=n; i++)
    {
        cin>>a[i];
    }
    reverse(a+1, a+n+1);

    k=0;
    while(k<=n && a[k]>=k)
    {
        k++;
    }
    k--;
    cout<<k<<" ";

    x=a[k+1];
    res+=dinamic(a, k, x);
    res%=MOD;

    for(i=n; i>0; i--)
    {
        j=1;
        while(a[j]>=i)
        {
            j++;
        }
        b[i]=j-1;
    }

    x=b[k+1];

    for(i=0; i<=n; i++)
    {
        for(j=0; j<=n; j++)
        {
            dp[i][j]=0;
        }
    }
    res+=dinamic(b, k, x);
    res%=MOD;
    res-=fact(k);
    while(res<0)
        res+=MOD;
    cout<<res;
    return 0;
}