#include<bits/stdc++.h>
#define inl inline
#define CL(s) memset(s,0,sizeof(s))
using namespace std;
const int N=2005,P=998244353;
int n,x,y,an,w[N],s0[N],s1[N],s2[N],s3[N];
struct F {int a[N],b[N],c[N]; inl void G() {CL(a); CL(b); CL(c); } }f[2];
inl void Az(int &x,int y) {(x+=y)>=P&&(x-=P); }
int main()
{
	for(;~scanf("%d",&n);printf("%d\n",an),an=0)
	{
		f[x=0].G(); y=1; for(int i=1;i<=n;++i)
		{
			scanf("%d",&w[i]);
			f[0].a[i]=w[i], f[0].b[i]=(i-1ll)*w[i]%P, f[0].c[i]=(i-1ll)*(i-1)*w[i]%P;
		}
		for(int i=1;i<n-2;++i,x^=1,y^=1)
		{
			f[y].G(); CL(s0); CL(s1); CL(s2); CL(s3);
			for(int j=1;j<=n;++j)
			{
				int t0=f[x].a[j],t1=f[x].b[j],t2=f[x].c[j];
				Az(s0[j],t0); Az(s1[j],t0); Az(s2[1],t0); Az(s2[j+1],P-t0); s3[j+1]=(1ll*j*j*t0+s3[j+1])%P;
				Az(f[y].a[j],t1); s3[j+1]=(1ll*j*t1+s3[j+1])%P; Az(f[y].b[j],t2); Az(s3[j+1],t2);
			}
			for(int j=n;j;--j) Az(s0[j],s0[j+1]), Az(s1[j],s1[j+1]);
			for(int j=1;j<=n;++j)
			{
				Az(s2[j],s2[j-1]); Az(s3[j],s3[j-1]);
				f[y].a[j]=(1ll*s0[j]+f[y].a[j])*w[j]%P;
				f[y].b[j]=(s1[j]*(j-1ll)+f[y].b[j])%P*w[j]%P;
				f[y].c[j]=(s2[j]*(j-1ll)*(j-1)+s3[j])%P*w[j]%P;
			}
		}
		for(int i=1;i<=n;++i) an=(1ll*f[x].a[i]*i*i+1ll*f[x].b[i]*i+f[x].c[i]+an)%P;
	}
	return 0;
}