//#ifdef ONLINE_JUDGE #pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include #include //#endif #include #define bgn(a) begin(a) #define rbg(a) rbegin(a) #define fin(...) exit(pp(__VA_ARGS__)) #define ep(...) emplace(__VA_ARGS__) #define ef(...) emplace_front(__VA_ARGS__) #define eb(...) emplace_back(__VA_ARGS__) #define pf(...) pop_front(__VA_ARGS__) #define pb(...) pop_back(__VA_ARGS__) #define ei(...) else if(__VA_ARGS__) #define el else #define wh(...) while(__VA_ARGS__) #define lb(...) lower_bound(__VA_ARGS__) #define ub(...) upper_bound(__VA_ARGS__) #define itmax(...) max_element(al(__VA_ARGS__)) #define itmin(...) min_element(al(__VA_ARGS__)) #define srt(...) sort(al(__VA_ARGS__)) #define rv(...) reverse(al(__VA_ARGS__)) #define uq(a) srt(a);a.erase(unique(al(a)),end(a)) #define sw(a,b) swap(a,b) #define ov3(a,b,c,d,...) d #define ov4(a,b,c,d,e,...) e #define ov5(a,b,c,d,e,f,...) f #define ov6(a,b,c,d,e,f,g,...) g #define al1(v) bgn(v),end(v) #define al2(v,b) bgn(v),bgn(v)+b #define al3(v,a,b) bgn(v)+a,bgn(v)+b #define al(...) ov3(__VA_ARGS__,al3,al2,al1)(__VA_ARGS__) #define fo1(b) for(ll ii=0;ii<(ll)(b);ii++) #define fo2(i,b) for(ll i=0;i<(ll)(b);i++) #define fo3(i,a,b) for(ll i=(ll)(a);i<(ll)(b);i++) #define fo4(i,a,b,c) for(ll i=(ll)(a);i<(ll)(b);i+=(ll)(c)) #define fo(...) ov4(__VA_ARGS__,fo4,fo3,fo2,fo1)(__VA_ARGS__) #define of1(a) fo(a) #define of2(i,a) for(ll i=(ll)(a)-1;i>=0;i--) #define of3(i,a,b) for(ll i=(ll)(a)-1;i>=(ll)(b);i--) #define of4(i,a,b,c) for(ll i=(ll)(a)-1;i>=(ll)(b);i-=(ll)(c)) #define of(...) ov4(__VA_ARGS__,of4,of3,of2,of1)(__VA_ARGS__) #define fe1(a) fo(len(a)) #define fe2(a,v) for(au&&a:v) #define fe3(a,b,v) for(au&&[a,b]:v) #define fe4(a,b,c,v) for(au&&[a,b,c]:v) #define fe5(a,b,c,d,v) for(au&&[a,b,c,d]:v) #define fe(...) ov5(__VA_ARGS__,fe5,fe4,fe3,fe2,fe1)(__VA_ARGS__) #define LL(...) ll __VA_ARGS__;li(__VA_ARGS__) #define DD(...) dd __VA_ARGS__;li(__VA_ARGS__) #define ST(...) str __VA_ARGS__;li(__VA_ARGS__) #define UL2(n,a) u1 a(n);li(a) #define UL3(n,a,b) u1 a(n),b(n);li(a,b) #define UL4(n,a,b,c) u1 a(n),b(n),c(n);li(a,b,c) #define UL5(n,a,b,c,d) u1 a(n),b(n),c(n),d(n);li(a,b,c,d) #define UL6(n,a,b,c,d,e) u1 a(n),b(n),c(n),d(n),e(n);li(a,b,c,d,e) #define UL(...) ov6(__VA_ARGS__,UL6,UL5,UL4,UL3,UL2)(__VA_ARGS__) #define UV2(n,a) u1 a(n);vi(a) #define UV3(n,a,b) u1 a(n),b(n);vi(a,b) #define UV4(n,a,b,c) u1 a(n),b(n),c(n);vi(a,b,c) #define UV5(n,a,b,c,d) u1 a(n),b(n),c(n),d(n);vi(a,b,c,d) #define UV6(n,a,b,c,d,e) u1 a(n),b(n),c(n),d(n),e(n);vi(a,b,c,d,e) #define UV(...) ov6(__VA_ARGS__,UV6,UV5,UV4,UV3,UV2)(__VA_ARGS__) #define DL2(n,a) v1
a(n);li(a) #define DL3(n,a,b) v1
a(n),b(n);li(a,b) #define DL4(n,a,b,c) v1
a(n),b(n),c(n);li(a,b,c) #define DL5(n,a,b,c,d) v1
a(n),b(n),c(n),d(n);li(a,b,c,d) #define DL(...) ov5(__VA_ARGS__,DL5,DL4,DL3,DL2)(__VA_ARGS__) #define DV2(n,a) v1
a(n);vi(a) #define DV3(n,a,b) v1
a(n),b(n);vi(a,b) #define DV4(n,a,b,c) v1
a(n),b(n),c(n);vi(a,b,c) #define DV5(n,a,b,c,d) v1
a(n),b(n),c(n),d(n);vi(a,b,c,d) #define DV(...) ov5(__VA_ARGS__,DV5,DV4,DV3,DV2)(__VA_ARGS__) #define U23(h,w,a) u2 a;rs(a,h,w);li(a) #define U24(h,w,a,b) u2 a,b;rs(a,h,w),rs(b,h,w);li(a,b) #define U25(h,w,a,b,c) u2 a,b,c;rs(a,h,w),rs(b,h,w),rs(c,h,w);li(a,b,c) #define U2(...) ov5(__VA_ARGS__,U25,U24,U23)(__VA_ARGS__) #define S22(h,a) strs a(h);li(a) #define S23(h,a,b) strs a(h),b(h);li(a,b) #define S2(...) ov3(__VA_ARGS__,S23,S22)(__VA_ARGS__) #define aa2(a,b) ((a)&&(b)) #define aa3(a,b,c) ((a)&&(b)&&(c)) #define aa4(a,b,c,d) ((a)&&(b)&&(c)&&(d)) #define aa5(a,b,c,d,e) ((a)&&(b)&&(c)&&(d)&&(e)) #define aa6(a,b,c,d,e,f) ((a)&&(b)&&(c)&&(d)&&(e)&&(f)) #define aa(...) ov6(__VA_ARGS__,aa6,aa5,aa4,aa3,aa2)(__VA_ARGS__) #define oo2(a,b) ((a)||(b)) #define oo3(a,b,c) ((a)||(b)||(c)) #define oo4(a,b,c,d) ((a)||(b)||(c)||(d)) #define oo5(a,b,c,d,e) ((a)||(b)||(c)||(d)||(e)) #define oo6(a,b,c,d,e,f) ((a)||(b)||(c)||(d)||(e)||(f)) #define oo(...) ov6(__VA_ARGS__,oo6,oo5,oo4,oo3,oo2)(__VA_ARGS__) #define I template #define J class #define rr return #define ss struct #define uu using #define as assert #define au auto #define bk break #define cs const #define ct continue //#define endl "\n" #define fd friend #define op operator #define sc static #define th this namespace my{ //#ifdef ONLINE_JUDGE uu namespace atcoder; uu namespace boost::multiprecision; uu mint=modint998244353; uu cint=cpp_int; //endif uu namespace std; uu ll=long long; uu dd=long double; uu ui=unsigned int; uu cha=char; uu str=string; uu bo=bool; uu vo=void; uu is=istream; uu os=ostream; str sp=" ",nc="",nl="\n"; au np=nullptr; Iss cp{ A a;B b; cp():a(A()),b(B()){} cp(A a,B b):a(a),b(b){} cp op+()cs{rr*th;} cp op-()cs{rr cp(-a,-b);} cp&op++(){a++,b++;rr*th;}cp op++(int){cp r=*th;++*th;rr r;} cp&op--(){a--,b--;rr*th;}cp op--(int){cp r=*th;--*th;rr r;} cp&op+=(cs cp&c){a+=c.a;b+=c.b;rr*th;}cp op+(cs cp&c)cs{rr cp(*th)+=c;} cp&op-=(cs cp&c){a-=c.a;b-=c.b;rr*th;}cp op-(cs cp&c)cs{rr cp(*th)-=c;} cp&op*=(cs cp&c){a*=c.a;b*=c.b;rr*th;}cp op*(cs cp&c)cs{rr cp(*th)*=c;} cp&op+=(cs A&d){a+=d;b+=d;rr*th;}cp op+(cs A&d)cs{rr cp(*th)+=d;} cp&op-=(cs A&d){a-=d;b-=d;rr*th;}cp op-(cs A&d)cs{rr cp(*th)-=d;} cp&op*=(cs A&d){a*=d;b*=d;rr*th;}cp op*(cs A&d)cs{rr cp(*th)*=d;} bo op==(cs cp&c)cs{rr aa(a==c.a,b==c.b);}bo op!=(cs cp&c)cs{rr!(*th==c);} bo op<(cs cp&c)cs{rr a!=c.a?a=(cs cp&c)cs{rr!(*th(cs cp&c)cs{rr a!=c.a?a>c.a:b>c.b;}bo op<=(cs cp&c)cs{rr!(*th>c);} fd is&op>>(is&i,cp&c){rr i>>c.a>>c.b;} fd os&op<<(os&o,cs cp&c){rr o<ss tr{ A a;B b;C c; tr():a(A()),b(B()),c(C()){} tr(A a,B b,C c):a(a),b(b),c(c){} bo op==(cs tr&t)cs{rr aa(a==t.a,b==t.b,c==t.c);} bo op!=(cs tr&t)cs{rr!(*th==t);} bo op<(cs tr&t)cs{rr a!=t.a?a(cs tr&t)cs{rr a!=t.a?a>t.a:b!=t.b?b>t.b:c>t.c;} fd is&op>>(is&i,tr&t){rr i>>t.a>>t.b>>t.c;} fd os&op<<(os&o,cs tr&t){rr o<ss qu{ A a;B b;C c;D d; qu():a(A()),b(B()),c(C()),d(D()){} qu(A a,B b,C c,D d):a(a),b(b),c(c),d(d){} bo op==(cs qu&q)cs{rr aa(a==q.a,b==q.b,c==q.c,d==q.d);} bo op!=(cs qu&q)cs{rr!(*th==q);} bo op<(cs qu&q)cs{rr a!=q.a?a(cs qu&q)cs{rr a!=q.a?a>q.a:b!=q.b?b>q.b:c!=q.c?c>q.c:d>q.d;} fd is&op>>(is&i,qu&q){rr i>>q.a>>q.b>>q.c>>q.d;} fd os&op<<(os&o,cs qu&q){rr o<uu v1=vector;uu u1=v1; Iuu v2=v1>;uu u2=v1; Iuu v3=v1>;uu u3=v1; Iuu v4=v1>;uu u4=v1; Iuu mset=multiset; Iuu fn=function; Iuu pqmax=priority_queue; Iuu pqmin=priority_queue,greater>; uu strs=v1; uu cl=cp; uu tl=tr; uu cd=cp; uu ql=deque; uu sl=set; uu msl=mset; uu ml=map; ss edg{ll t,w;edg(){}edg(ll t,ll w=1):t(t),w(w){}}; uu graph=v2; ll inf=3e18; dd e=1e-10,pi=acosl(-1); u1 dx={-1,0,1,0,-1,1,1,-1},dy={0,-1,0,1,-1,-1,1,1}; str first(bo t=1){rr t?"first":"second";} str First(bo t=1){rr t?"First":"Second";} str yes(bo t=1){rr t?"yes":"no";}str no(){rr yes(0);} str Yes(bo t=1){rr t?"Yes":"No";}str No(){rr Yes(0);} str YES(bo t=1){rr t?"YES":"NO";}str NO(){rr YES(0);} str possible(bo t=1){rr t?"possible":"impossible";} str Possible(bo t=1){rr t?"Possible":"Impossible";} str POSSIBLE(bo t=1){rr t?"POSSIBLE":"IMPOSSIBLE";} ui rnd(){ sc ui x=123456789,y=362436069,z=521288629,w=88675123; ui t=(x^(x<<11));x=y,y=z,z=w; rr(w=(w^(w>>19))^(t^(t>>8))); } dd rnd_prb(){rr (dd)rnd()/UINT_MAX;} ll rnd_rng(ll r){rr rnd()%r;} ll rnd_rng(ll l,ll r){rr rnd()%(r-l)+l;} bo od(ll x){rr x&1;} bo ev(ll x){rr !od(x);} ll bit(ll a){rr 1LL<ll len(cs T&a){rr a.size();} Ibo yu(cs T&a){rr len(a);} Ibo mu(cs T&a){rr !yu(a);} IT sq(T a){rr a*a;} ll at(ll S,ll i){rr S>>i&1;} IT at(cs v1&v,ll i){ll n=len(v);rr v[(i%n+n)%n];} Ill sgn(cs T&a){rr(a>e)-(a<-e);} Ill sgn(cs T&a,cs U&b){rr sgn(a-b);} IT cei(T x,U y){as(y);rr(y<0?cei(-x,-y):(x>0?(x+y-1)/y:x/y));} IT flo(T x,U y){as(y);rr(y<0?flo(-x,-y):(x>0?x/y:x/y-!(x%y)));} str sb(cs str&s,ll i,ll n){rr s.substr(i,n);} str tos(cs cha&a){rr a+nc;} str tos(cs cha*a){rr a;} str tos(cs str&a){rr a;} Istr tos(cs T&a){rr to_string(a);} Istr tos(cs T&a,cs A&...b){rr tos(a)+tos(b...);} u1 tou(cs str&s,cha b='a'){u1 a;fe(c,s)a.eb(c-b);rr a;} u2 tou(cs strs&s,cha b='a'){u2 a;fe(t,s)a.eb(tou(t,b));rr a;} u1 tou(cs str&s,cs str&t){ ll n=len(s); u1 a(n); fo(i,n)fo(j,len(t))if(s[i]==t[j]){a[i]=j;bk;} rr a; } u2 tou(cs strs&s,cs str&t){ll h=len(s);u2 a(h);fo(i,h)a[i]=tou(s[i],t);rr a;} Iau op+(v1v){rr v;} Iau op-(v1v){fe(x,v)x=-x;rr v;} Iau&op++(v1&v){fe(x,v)x++;rr v;} Iau&op--(v1&v){fe(x,v)x--;rr v;} Iau op++(v1&v,int){au r=v;++v;rr r;} Iau op--(v1&v,int){au r=v;--v;rr r;} Iau&op+=(v1&v,cs v1&u){fo(i,len(v))v[i]+=u[i];rr v;} Iau&op-=(v1&v,cs v1&u){fo(i,len(v))v[i]-=u[i];rr v;} Iau&op*=(v1&v,cs v1&u){fo(i,len(v))v[i]*=u[i];rr v;} Iau op+(v1v,cs v1&u){rr v+=u;} Iau op-(v1v,cs v1&u){rr v-=u;} Iau op*(v1v,cs v1&u){rr v*=u;} Iau&op+=(v1&v,cs U&a){fe(x,v)x+=a;rr v;} Iau&op-=(v1&v,cs U&a){fe(x,v)x-=a;rr v;} Iau&op*=(v1&v,cs U&a){fe(x,v)x*=a;rr v;} Iau&op/=(v1&v,cs U&a){fe(x,v)x/=a;rr v;} Iau&op%=(v1&v,cs U&a){fe(x,v)x%=a;rr v;} Iau op+(v1v,cs U&a){rr v+=a;} Iau op-(v1v,cs U&a){rr v-=a;} Iau op*(v1v,cs U&a){rr v*=a;} Iau op/(v1v,cs U&a){rr v/=a;} Iau op%(v1v,cs U&a){rr v%=a;} Iau jo(v1&a,cs v1&b){copy(al(b),back_inserter(a));} Ibo amax(T&a,cs U&b){rr abo amin(T&a,cs U&b){rr a>b?a=b,1:0;} IT max(cs v1&a){rr *itmax(a);} IT min(cs v1&a){rr *itmin(a);} Iau max(cs v2&a){rr max(*itmax(a));} Iau min(cs v2&a){rr min(*itmin(a));} Iau max(T...a){rr max(initializer_list>{a...});} Iau min(T...a){rr min(initializer_list>{a...});} Ill argmax(cs v1&a){rr itmax(a)-bgn(a);} Ill argmin(cs v1&a){rr itmin(a)-bgn(a);} IT gcd(cs T&a,cs A&...b){rr gcd(a,gcd(b...));} IT sum(cs v1&a){rr accumulate(al(a),T());} Ill lbs(cs v1&a,cs U&b){rr lb(al(a),b)-bgn(a);} Ill ubs(cs v1&a,cs U&b){rr ub(al(a),b)-bgn(a);} Iau minmax(cs T&a,cs U&b){rr cp(min(a,b),max(a,b));} Iau minmax(cs cp&p){rr minmax(p.a,p.b);} Iau mp(cs fn&f,v1x){v1y;fe(a,x)y.eb(f(a));rr y;} Iau keys(cs map&m){v1r;fe(k,v,m)r.eb(k);rr r;} Iau vals(cs map&m){v1r;fe(k,v,m)r.eb(v);rr r;} u1 io(ll n,ll a=0){u1 i(n);iota(al(i),a);rr i;} Ivo emplace_front(v1&v,U x=U()){v.ep(bgn(v),x);} Ivo emplace_front(v1&v,T x=T()){v.ep(bgn(v),x);} Ivo pop_front(v1&v){v.erase(bgn(v));} Iau&bg(T&a){rr*bgn(a);} Iau&rb(T&a){rr*rbg(a);} IT bg(pqmax&q){rr q.top();} IT bg(pqmin&q){rr q.top();} IT pop(v1&v){T r=rb(v);v.pb();rr r;} IT pop_front(deque&q){T r=bg(q);q.pf();rr r;} IT pop_back(deque&q){T r=rb(q);q.pb();rr r;} IT pop(deque&q){rr pop_front(q);} IT pop(pqmax&q){T r=bg(q);q.pop();rr r;} IT pop(pqmin&q){T r=bg(q);q.pop();rr r;} ss z{z(){cin.tie(0);ios::sync_with_stdio(0);cout<>(is&i,mint&x){ll t;i>>t;x=t;rr i;} os&op<<(os&o,cs mint&x){rr o<is&op>>(is&i,pair&p){rr i>>p.first>>p.second;} Ios&op<<(os&o,cs pair&p){rr o<is&op>>(is&i,v1&v){fe(x,v)i>>x;rr i;} Ios&op<<(os&o,cs v1&v){fe(x,v)o<os&op<<(os&o,cs v2&v){fe(x,v)o<os&op<<(os&o,cs map&m){fe(p,m)o<os&op<<(os&o,cs set&s){fe(x,s)o<os&op<<(os&o,cs mset&s){fe(x,s)o<os&op<<(os&o,dequeq){wh(yu(q))o<os&op<<(os&o,pqmaxq){wh(yu(q))o<os&op<<(os&o,pqminq){wh(yu(q))o<ll pp(cs T&a,cs A&...b){cout<ll li(T&...a){(cin>>...>>a);rr 0;} Ivo vi(ll i,T&...a){(cin>>...>>a[i]);} Ivo vi(v1&a,A&...b){fo(i,len(a))vi(i,a,b...);} Ivo rs(v1&v,ll n){v.resize(n);} Ivo rs(v1&V,ll n,cs A&...a){rs(V,n);fe(v,V)rs(v,a...);} Ivo fl(T&x,cs U&a){x=a;} Ivo fl(v1&V,cs U&a){fe(v,V)fl(v,a);} vo gi(u2&g,ll m,bo d=1){fo(m){LL(a,b);a--,b--;g[a].eb(b);if(d)g[b].eb(a);}} vo gi(graph&g,ll m,bo w=1,bo d=1){ fo(m){LL(a,b);a--,b--;ll c=1;if(w)li(c);g[a].eb(b,c);if(d)g[b].eb(a,c);} } vo ti(u2&g,bo d=1){gi(g,len(g)-1,d);} vo ti(graph&g,bo w=1,bo d=1){gi(g,len(g)-1,w,d);} Iss rcl:private F{ explicit rcl(F&&f):F(forward(f)){} Idecltype(au)op()(T&&...a)cs{rr F::op()(*th,forward(a)...);} }; vo pz(cs u1&v,u1&h){fe(x,h)x=lbs(v,x);} Ivo pz(cs u1&v,u1&h,T&...t){pz(v,h);pz(v,t...);} Iu1 zp(T&...a){u1 v;fe(e,{a...})jo(v,e);uq(v);pz(v,a...);rr v;} Ivo sv(cs u1&o,v1&a){au c=a;fo(i,len(a))a[i]=c[o[i]];} Ivo sv(cs u1&o,v1&a,A&...b){sv(o,a);sv(o,b...);} Iau vs(cs fn&f,T&...a){ ll n=0;fe(e,{a...})n=len(e); u1 o=io(n),p(n); sort(al(o),f); fo(i,n)p[o[i]]=i; sv(o,a...); rr cp(o,p); } ll bsl(cs fn&j,ll o,ll n){wh(abs(o-n)>1)(j((o+n)/2)?o:n)=(o+n)/2;rr o;} dd bsd(cs fn&j,dd o,dd n){wh(abs(o-n)>e)(j((o+n)/2)?o:n)=(o+n)/2;rr o;} Iau zt(v1a){fo(i,1,len(a))a[i]+=a[i-1];rr a;} Iau mb(v1a){of(i,len(a),1)a[i]-=a[i-1];rr a;} Iau zt(v2a){fe(v,a)v=zt(v);fo(i,1,len(a))a[i]+=a[i-1];rr a;} Iau mb(v2a){fe(v,a)v=mb(v);of(i,len(a),1)a[i]-=a[i-1];rr a;} au kth(cs u1&a,ll k){u1 s=zt(a);rr cp(lbs(s,k),s);} au kth2(cs u1&a,ll&k){au[i,s]=kth(a,k);k-=s[i-1];rr i;} IT esq(T x,ll n,T e){wh(n){if(n&1)e*=x;x*=x;n/=2;}rr e;} vo wafl(u2&d){ll n=len(d);fo(k,n)fo(i,n)fo(j,n)amin(d[i][j],d[i][k]+d[k][j]);} u1 divs(ll n){ml m;fo(i,1,n/i+1)if(n%i==0)m[i]++,m[n/i]++;rr keys(m);} ml fact(ll n){ml m;fo(i,2,n/i+1)wh(n%i==0)m[i]++,n/=i;if(n>1)m[n]++;rr m;} ll pow(ll x,ll n,ll m=inf){ll r=1;wh(n){if(n&1)r*=x,r%=m;x*=x,x%=m,n/=2;}rr r;} ll inv(ll a,ll m){ ll b=m,u=1,v=0; wh(b)u-=a/b*v,sw(u,v),a-=a/b*b,sw(a,b); rr(u%m+m)%m; } ll mrg(u1&a){ ll n=len(a); if(n<=1)rr 0; u1 b(al(a,n/2)),c(al(a,n/2,n)); ll cnt=mrg(b)+mrg(c); for(ll i=0,j=0,k=0;ill jobsch(v1>s){ vs([&](ll i,ll j){rr s[i].bdp(n+1),su(n+1);dp[0]=su[0]=1; u1 p(len(zp(a))); fo(i,1,n+1){ ll&k=p[a[i]]; dp[i]=su[i-1]-(k?su[k-1]:0); su[i]=su[i-1]+dp[i]; k=i; } rr su[n]-1; } u2 bfs2(cs u2&g,ll sx=0,ll sy=0){ ll h=len(g),w=len(g[0]); u2 d;rs(d,h,w);fl(d,inf);d[sx][sy]=0; dequeq;q.eb(sx,sy); wh(yu(q)){ au[x,y]=pop(q); fo(k,4){ ll nx=x+dx[k],ny=y+dy[k]; if(oo(nx<0,ny<0,nx>=h,ny>=w))ct; if(g[nx][ny])ct; if(d[nx][ny]==inf)d[nx][ny]=d[x][y]+1,q.eb(nx,ny); } } rr d; } u1 bfs(cs u2&g,ll s=0){ ll n=len(g); u1 d(n,inf);d[s]=0; ql q;q.eb(s); wh(yu(q)){ ll u=pop(q); fe(v,g[u])if(d[v]==inf)d[v]=d[u]+1,q.eb(v); } rr d; } u1 djks(cs graph&g,ll s=0){ ll n=len(g); u1 d(n,inf);d[s]=0; pqminq;q.ep(0,s); wh(yu(q)){ au[du,u]=pop(q); if(d[u]ss com{ v1fa,rf,bs; v2m; v2u; com(ll n){ n++; rs(fa,n);rs(rf,n);rs(bs,n); fa[0]=rf[0]=bs[0]=1; fo(i,1,n)fa[i]=fa[i-1]*i,rf[i]=rf[i-1]/i; fo(i,1,n)bs[i]=bs[i-1]+rf[i]*pow(-1,i&1); if(n<=5e3)rs(m,n,n),rs(u,n,n); } T c(ll n,ll k){rr oo(n=0)rr u[n][k]=1,m[n][k]=p(n,k-1)+p(n-k,k); rr u[n][k]=1,m[n][k]=p(n,k-1); } }; ss ufnd{ ll n; u1 d; ufnd():n(0){} ufnd(ll n):n(n),d(n,-1){} ll ldr(ll a){as(aa(0<=a,a&f,cs fn&e>ss sgt{ ll n,A,log; v1d; sgt(ll n):sgt(v1(n,e())){} sgt(cs v1&v):n(len(v)){ log=lg2(n*2-1)-1; A=bit(log); d=v1(A*2,e()); fo(i,n)d[i+A]=v[i]; of(i,A,1)upd(i); } ll size()cs{rr n;} vo upd(ll k){d[k]=f(d[k*2],d[k*2+1]);} vo set(ll i,T x){as(aa(0<=i,i>j);} vo add(ll i,T x){as(aa(0<=i,ida;rs(da,N+1,A+1,B+1,C+1); v4
db;rs(db,N+1,A+1,B+1,C+1); v4
dc;rs(dc,N+1,A+1,B+1,C+1); u4 used;rs(used,N+1,A+1,B+1,C+1); au f=rcl([&](au f,ll i,ll a,ll b,ll c){ if(i==0)rr; if(used[i][a][b][c])rr; used[i][a][b][c]=1; f(i-1,a,b,c); ll s=a+b+c; dd pa=(s>=2?(dd)a*(a-1)/(s*(s-1)):0); dd pb=(s>=2?(dd)b*(b-1)/(s*(s-1)):0); dd pc=(s>=2?(dd)c*(c-1)/(s*(s-1)):0); if(a>=2)f(i-1,a-2,b,c); if(b>=2)f(i-1,a,b-2,c); if(c>=2)f(i-1,a,b,c-2); da[i][a][b][c]+=da[i-1][a][b][c]*(1-pa-pb-pc); da[i][a][b][c]+=pa; if(a>=2)da[i][a][b][c]+=da[i-1][a-2][b][c]*pa; if(b>=2)da[i][a][b][c]+=da[i-1][a][b-2][c]*pb; if(c>=2)da[i][a][b][c]+=da[i-1][a][b][c-2]*pc; db[i][a][b][c]+=db[i-1][a][b][c]*(1-pa-pb-pc); db[i][a][b][c]+=pb; if(a>=2)db[i][a][b][c]+=db[i-1][a-2][b][c]*pa; if(b>=2)db[i][a][b][c]+=db[i-1][a][b-2][c]*pb; if(c>=2)db[i][a][b][c]+=db[i-1][a][b][c-2]*pc; dc[i][a][b][c]+=dc[i-1][a][b][c]*(1-pa-pb-pc); dc[i][a][b][c]+=pc; if(a>=2)dc[i][a][b][c]+=dc[i-1][a-2][b][c]*pa; if(b>=2)dc[i][a][b][c]+=dc[i-1][a][b-2][c]*pb; if(c>=2)dc[i][a][b][c]+=dc[i-1][a][b][c-2]*pc; }); f(N,A,B,C); pp(da[N][A][B][C],db[N][A][B][C],dc[N][A][B][C]); } }int main(){uu namespace my;ll T=1;(0)&&li(T);fo(T)solve();}