#include<bits/stdc++.h>
using namespace std;
namespace P {
const int MOD=998244353,N=1<<19,G=3;
int rev[N],inv[N],fac[N],ifac[N],w[N<<1];
int ksm(int a,int b=MOD-2) {
	int ret=1;
	for(;b;a=1ll*a*a%MOD,b=b>>1) if(b&1) ret=1ll*ret*a%MOD;
	return ret;
}
void poly_init() {
	inv[1]=1;
	for(int i=2;i<N;++i) inv[i]=1ll*(MOD-MOD/i)*inv[MOD%i]%MOD;
	fac[0]=ifac[0]=1;
	for(int i=1;i<N;++i) fac[i]=1ll*fac[i-1]*i%MOD,ifac[i]=1ll*ifac[i-1]*inv[i]%MOD;
	for(int k=1;k<=N;k<<=1) {
		int x=ksm(G,(MOD-1)/k); w[k]=1;
		for(int i=1;i<k;++i) w[i+k]=1ll*x*w[i+k-1]%MOD;
	}
}
int plen(int x) { int y=1; for(;y<x;y<<=1); return y;  }
void ntt(int *f,bool idft,int n) {
	for(int i=0;i<n;++i) {
		rev[i]=(rev[i>>1]>>1);
		if(i&1) rev[i]|=n>>1;
	}
	for(int i=0;i<n;++i) if(rev[i]<i) swap(f[i],f[rev[i]]);
	for(int k=2,x,y;k<=n;k<<=1) {
		for(int i=0;i<n;i+=k) {
			for(int j=i;j<i+k/2;++j) {
				x=f[j],y=1ll*f[j+k/2]*w[k+j-i]%MOD;
				f[j]=(x+y>=MOD)?x+y-MOD:x+y,f[j+k/2]=(x>=y)?x-y:x+MOD-y;
			}
		}
	}
	if(idft) {
		reverse(f+1,f+n);
		for(int i=0,x=ksm(n);i<n;++i) f[i]=1ll*f[i]*x%MOD;
	}
}
void poly_mul(const int *f,const int *g,int *h,int n,int m) {
	static int a[N],b[N];
	for(int i=0;i<n;++i) a[i]=f[i];
	for(int i=0;i<m;++i) b[i]=g[i];
	int len=plen(n+m-1);
	ntt(a,0,len),ntt(b,0,len);
	for(int i=0;i<len;++i) h[i]=1ll*a[i]*b[i]%MOD;
	ntt(h,1,len);
	memset(a,0,sizeof(int)*len);
	memset(b,0,sizeof(int)*len);
}
}
const int MAXN=4e5+5,N=1<<19,MOD=998244353;
int n,a[MAXN],st[MAXN],tp,L[MAXN],R[MAXN],nxt[MAXN];
int m,b[MAXN],X[N],Y[N],Z[N],f[MAXN];
void calc() {
	for(int i=0;i<=m;++i) {
		X[i]=b[i+1]-b[i]-(i==m);
		Y[i]=b[m+1-i]-b[m-i]-(i==m);
	}
	P::poly_mul(X,Y,Z,m+1,m+1);
	for(int i=0;i<m;++i) f[m-i]=(f[m-i]+Z[i])%MOD;
}
void solve() {
	cin>>n;
	map <int,int> pos;
	for(int i=1;i<=n;++i) cin>>a[i],f[i]=0;
	for(int i=n;i>=1;--i) {
		if(!pos.count(a[i])) nxt[i]=n+2;
		else nxt[i]=pos[a[i]];
		pos[a[i]]=i;
	}
	for(int i=1;i<=n;++i) {
		while(tp&&a[st[tp]]<a[i]) R[st[tp--]]=i;
		st[++tp]=i;
	}
	while(tp) R[st[tp--]]=n+1;
	for(int i=n;i>=1;--i) {
		while(tp&&a[st[tp]]<a[i]) L[st[tp--]]=i;
		st[++tp]=i;
	}
	while(tp) L[st[tp--]]=0;
	for(auto it:pos) {
		int u=it.second;
		while(u<=n) {
			int v=u; b[0]=L[u],m=0;
			while(v<R[u]) b[++m]=v,v=nxt[v];
			b[m+1]=R[u],calc(),u=v;
		}
	}
	int ans=0;
	for(int i=1;i<=n;++i) ans=(ans+1ll*f[i]*f[i]%MOD*i)%MOD;
	cout<<ans<<"\n";
}
signed main() {
	P::poly_init();
	ios::sync_with_stdio(false);
	int T; cin>>T;
	while(T--) solve();
	return 0;
}