#include #include using namespace std; using namespace atcoder; using ll=long long; #pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #define speedup std::cin.tie(0)->sync_with_stdio(0); #define _GLIBCXX_DEBUG //優先度付きキュー使って時間制限厳しい場合コメントアウト #define int ll #define int_max 2147483647 #define int_min -2147483647 #define uint_max 4294967295 #define ll_max 9223372036854775807 #define ll_min -9223372036854775807 #define ull_max 18446744073709551615 #define rep(i,n) for(ll i=0;i<(n);i++) #define reps(i,n) for(ll i=1;i<=(n);i++) #define REP(i,j,n) for(ll i=(j);i<(n);i++) #define all(a) (a).begin(), (a).end() #define veci vector #define veci2 vector> //veci2 A(a,veci(b)) #define vecll vector #define vecll2 vector> #define vecst vector #define repc(i,n,A) rep(i,n)cin>>A[i] #define repc2(i,n,A,B) rep(i,n)cin>>A[i]>>B[i] #define repc2vec(i,j,a,b,A) rep(i,a)rep(j,b)cin>>A[i][j] #define repair(i,n,A) rep(i,n)cin>>A[i].F>>A[i].S>> #define ST(A) sort(all(A)) #define RV(A) reverse(all(A)); #define pb push_back #define mp make_pair #define Endl endl #define F first #define S second #define yes(b) ((b)?"yes":"no") #define Yes(b) ((b)?"Yes":"No") #define YES(b) ((b)?"YES":"NO") #define TA(b) ((b)?"Takahashi":"Aoki") #define AB(b) ((b)?"Alice":"Bob") template inline T gcd(T a,T b) {return (b==0)?a:gcd(b,a%b);}//最大公約数 template inline T lcm(T a, T b) {return (a*b)/gcd(a,b);}//最小公倍数 veci dx={-1,0,1,-1,1,-1,0,1}; veci dy={-1,-1,-1,0,0,1,1,1}; veci DX={-1,1,0,0}; veci DY={0,0,-1,1}; ll mod=998244353; //mod=1000000007; //ここからコード入力(関数使用時用) //関数 /*ll nibutan(ll K){ //二分探索 a=要素数 [m]の前配列名 ll ng=-1; ll ok=a; while(ok-ng>1){ ll m=(ng+ok)/2; if([m]>K)ok=m; else ng=m; } return ok; }*/ vector>soinsuubunkai(int N){ vector>res; for(int a=2;a*a<=N;a++){ if(N%a!=0)continue; int ex=0; // 指数 // 割れる限り割り続ける while(N%a==0){ ex++; N/=a; } // その結果を push res.pb({a, ex}); } // 最後に残った数について if(N!=1)res.pb({N,1}); return res; //N=6 {2,1} {3,1}(2^1+3^1)のようになる autoで受け取ろう } /*//コンビネーション ll C_MAX=510000; vecll fac(C_MAX); vecll finv(C_MAX); vecll inv(C_MAX); // テーブルを作る前処理 void COMinit(){ fac[0]=fac[1]=1; finv[0]=finv[1]=1; inv[1]=1; REP(i,2,C_MAX){ fac[i]=fac[i-1]*i%mod; inv[i]=mod-inv[mod%i]*(mod/i)%mod; finv[i]=finv[i-1]*inv[i]%mod; } } //二項係数計算 ll COM(ll n,ll k){ if(n>A>>B; ll now1=1,now2=1,now3=A,now4=B,end=0; vecll2 C(A+2,vecll(B+2)); rep(i,A+2)rep(j,B+2)C[i][j]=-1; reps(i,A)reps(j,B)cin>>C[i][j]; priority_queue,vector>,greater>>D; setD2; priority_queue,vector>,greater>>E; setE2; vecll2 color(A+2,vecll(B+2)); D.push(mp(0,1*(A+2)+1)); E.push(mp(0,A*(A+2)+B)); D2.insert(1*(A+2)+1); E2.insert(A*(A+2)+B); rep(i,A*B){ if(i%2==0){ now1=D.top().S/(A+2); now2=D.top().S%(A+2); D.pop(); color[now1][now2]=1; rep(j,4){ if(color[now1+DX[j]][now2+DY[j]]==2){ end++; break; } if(!(D2.count((now1+DX[j])*(A+2)+now2+DY[j]))&&C[now1+DX[j]][now2+DY[j]]!=-1){ D2.insert((now1+DX[j])*(A+2)+now2+DY[j]); D.push(mp(C[now1+DX[j]][now2+DY[j]],(now1+DX[j])*(A+2)+now2+DY[j])); } } } if(i%2==1){ now3=E.top().S/(A+2); now4=E.top().S%(A+2); E.pop(); color[now3][now4]=2; rep(j,4){ if(color[now3+DX[j]][now4+DY[j]]==1){ end++; //cout<<"!!!"<