#include<bits/stdc++.h>
using namespace std;
constexpr int maxn=5005;
int n,a[maxn],k,b[maxn];
typedef long long ll;
constexpr ll mod=998244353;
ll dp[maxn][maxn],ans;
inline void solve(){
	int t=a[k+1];
	memset(dp,0,sizeof dp);
	dp[0][0]=1;
	for(int i=1;i<=k;i++)
	for(int j=0;j<=t;j++) dp[i][j]=(dp[i-1][j]*(a[i]-j)+(j?dp[i-1][j-1]*(t-j+1):0))%mod,cout<<i<<' '<<j<<'\n';
	ans+=dp[k][t];
}
int main(){
	cin>>n;for(int i=1;i<=n;i++) cin>>a[i],k=max(k,min(a[i],n-i+1)),b[i]=a[i];
	reverse(a+1,a+n+1);solve();
	for(int i=1;i<=n;i++) a[i]=0;
	for(int i=1;i<=n;i++) for(int j=1;j<=b[i];j++) a[j]++;
	solve();
	ll w=1;
	for(int i=1;i<=k;i++) w=w*i%mod;
	cout<<k<<' '<<(ans+mod-w)%mod<<'\n';
}