#include using namespace std; #define all(a) a.begin(),a.end() #define pb push_back #define sz(a) ((int)a.size()) using ll=long long; using u32=unsigned int; using u64=unsigned long long; using i128=__int128; using u128=unsigned __int128; using f128=__float128; using pii=pair; using pll=pair; template using vc=vector; template using vvc=vc>; template using vvvc=vc>; using vi=vc; using vll=vc; using vvi=vc; using vvll=vc; #define vv(type,name,n,...) \ vector> name(n,vector(__VA_ARGS__)) #define vvv(type,name,n,m,...) \ vector>> name(n,vector>(m,vector(__VA_ARGS__))) template using min_heap=priority_queue,greater>; template using max_heap=priority_queue; // https://trap.jp/post/1224/ #define rep1(n) for(ll i=0; i<(ll)(n); ++i) #define rep2(i,n) for(ll i=0; i<(ll)(n); ++i) #define rep3(i,a,b) for(ll i=(ll)(a); i<(ll)(b); ++i) #define rep4(i,a,b,c) for(ll i=(ll)(a); i<(ll)(b); i+=(c)) #define cut4(a,b,c,d,e,...) e #define rep(...) cut4(__VA_ARGS__,rep4,rep3,rep2,rep1)(__VA_ARGS__) #define per1(n) for(ll i=((ll)n)-1; i>=0; --i) #define per2(i,n) for(ll i=((ll)n)-1; i>=0; --i) #define per3(i,a,b) for(ll i=((ll)a)-1; i>=(ll)(b); --i) #define per4(i,a,b,c) for(ll i=((ll)a)-1; i>=(ll)(b); i-=(c)) #define per(...) cut4(__VA_ARGS__,per4,per3,per2,per1)(__VA_ARGS__) #define rep_subset(i,s) for(ll i=(s); i>=0; i=(i==0?-1:(i-1)&(s))) template constexpr T ifloor(const T a, const S b){return a/b-(a%b&&(a^b)<0);} template constexpr T iceil(const T a, const S b){return ifloor(a+b-1,b);} template void sort_unique(vector &vec){ sort(vec.begin(),vec.end()); vec.resize(unique(vec.begin(),vec.end())-vec.begin()); } template constexpr bool chmin(T &a, const S b){if(a>b) return a=b,true; return false;} template constexpr bool chmax(T &a, const S b){if(a istream& operator >> (istream& i, pair &p){return i >> p.first >> p.second;} template ostream& operator << (ostream& o, const pair &p){return o << p.first << ' ' << p.second;} #ifdef i_am_noob #define bug(...) cerr << "#" << __LINE__ << ' ' << #__VA_ARGS__ << "- ", _do(__VA_ARGS__) template void _do(vector x){for(auto i: x) cerr << i << ' ';cerr << "\n";} template void _do(set x){for(auto i: x) cerr << i << ' ';cerr << "\n";} template void _do(unordered_set x){for(auto i: x) cerr << i << ' ';cerr << "\n";} template void _do(T && x) {cerr << x << endl;} template void _do(T && x, S&&...y) {cerr << x << ", "; _do(y...);} #else #define bug(...) 777771449 #endif template void print(vector x){for(auto i: x) cout << i << ' ';cout << "\n";} template void print(set x){for(auto i: x) cout << i << ' ';cout << "\n";} template void print(unordered_set x){for(auto i: x) cout << i << ' ';cout << "\n";} template void print(T && x) {cout << x << "\n";} template void print(T && x, S&&... y) {cout << x << ' ';print(y...);} template istream& operator >> (istream& i, vector &vec){for(auto &x: vec) i >> x; return i;} vvi read_graph(int n, int m, int base=1){ vvi adj(n); for(int i=0,u,v; i> u >> v,u-=base,v-=base; adj[u].pb(v),adj[v].pb(u); } return adj; } vvi read_tree(int n, int base=1){return read_graph(n,n-1,base);} template pair operator + (const pair &a, const pair &b){return {a.first+b.first,a.second+b.second};} template constexpr T inf=0; template<> constexpr int inf = 0x3f3f3f3f; template<> constexpr ll inf = 0x3f3f3f3f3f3f3f3f; template vector operator += (vector &a, int val){for(auto &i: a) i+=val; return a;} template T isqrt(const T &x){T y=sqrt(x+2); while(y*y>x) y--; return y;} #define ykh mt19937 rng(chrono::steady_clock::now().time_since_epoch().count()) #include using namespace atcoder; using mint=modint998244353; //using mint=modint1000000007; template using poly=vector>; template poly add(poly a, poly b){ poly res=a; if(b.size()>a.size()) res.resize(b.size()); for(int i=0; i poly sub(poly a, poly b){ poly res=a; if(b.size()>a.size()) res.resize(b.size()); for(int i=0; i poly inverse(poly a){ int n=a.size(); assert(a[0]!=0); poly res(1,a[0].inv()); for(int m=1; m v1(a.begin(),a.begin()+m*2),v2=res; v1=convolution(convolution(v1,v2),v2); res.resize(m*2); for(int i=0; i pair,poly> divide(poly a, poly b){ int n=a.size(),m=b.size(),k=n-m+1; if(n ra=a,rb=b; reverse(all(ra)),reverse(all(rb)); ra.resize(k),rb.resize(k); poly Q=convolution(ra,inverse(rb)); Q.resize(k); reverse(all(Q)); poly res=convolution(b,Q),R(m-1); for(int i=0; i poly Derivative(poly a){ int n=a.size(); poly res(n-1); for(int i=0; i poly Integral(poly a){ int n=a.size(); poly res(n+1); for(int i=0; i poly Ln(poly a){ // a[0] = 1 int n=a.size(); if(n==1) return {}; poly d=Derivative(a); a.pop_back(); poly res=convolution(d,inverse(a)); res.resize(n-1); return Integral(res); } template poly Exp(poly a){ // a[0] = 0 int n=a.size(); poly q(1,1); a[0]+=1; for(int m=1; m g(a.begin(),a.begin()+m*2),h(all(q)); h.resize(m*2),h=Ln(h); for(int i=0; i poly pow(poly a, ll k){ int n=a.size(),m=0; poly res(n,0); while(m=n||k*m>=n)) return res; if(m==n) return res[0]=1,res; int lead=k*m; poly b(a.begin()+m,a.end()); static_modint base=b[0].pow(k),inv=b[0].inv(); for(int i=0; i> n >> m; vi a(n); cin >> a; string s; cin >> s; reverse(all(s)); s+='a'; vi tmp1,tmp2; int cur=0; for(auto c: s){ if(c=='s') cur++; else{ tmp1.pb(cur),tmp2.pb(cur+1); } } auto solve=[&](auto &self, int l, int r, bool flag) -> vc{ if(l+1==r){ if(flag) return {1,-tmp2[l]}; else return {1,-tmp1[l]}; } int mid=l+r>>1; return convolution(self(self,l,mid,flag),self(self,mid,r,flag)); }; vc vec1=solve(solve,0,sz(tmp1),0),vec2=solve(solve,0,sz(tmp2),1); vec1.resize(n+1),vec2.resize(n+1); // 1 / (1 - cur * x) vc vec3(n+1); vec3[0]=1; rep(n) vec3[i+1]=vec3[i]*cur; vc tmp=convolution(vec1,vec3); tmp.resize(n+1); vc P=convolution(tmp,inverse(vec2)); mint res; rep(n) res+=a[i]*P[i]*P[n-1-i]; print(res.val()); } signed main(){ ios_base::sync_with_stdio(0),cin.tie(0); cout << fixed << setprecision(20); int t=1; //cin >> t; while(t--) ahcorz(); }