#include #include #include using namespace std; using namespace atcoder; using mint=modint998244353; int mem=0; mint Q_data[62][7]; unordered_map mp; mint Bostan_Mori(long long N){ if(mp.count(N)){ return mp[N]; } if(N==0){ mint ans=1; mp[N]=ans; return ans; } if(N<=6){ mint ans=mint(2).pow(N-1); mp[N]=ans; return ans; } long long M=N; mint Q[7]; if(mem>0){ copy(begin(Q_data[0]),end(Q_data[0]),begin(Q)); } if(mem==0){ Q[0]=1; for(int i=1;i<=6;i++){ Q[i]=-1; } copy(begin(Q),end(Q),begin(Q_data[0])); mem++; } mint P[6]; P[0]=1; for(int i=1;i<6;i++){ P[i]=0; } int count=0; while(N>0){ count++; if(mem<=count){ mint nQ[7]; for(int i=0;i<=6;i++){ nQ[i]=0; for(int j=max(0,2*i-6);j<=min(2*i,6);j++){ if(j%2==0){ nQ[i]+=Q[j]*Q[2*i-j]; continue; } nQ[i]-=Q[j]*Q[2*i-j]; } } copy(begin(nQ),end(nQ),begin(Q_data[mem])); mem++; } mint nP[6]; if(N%2==0){ for(int i=0;i<12;i+=2){ nP[i/2]=0; for(int j=max(0,i-5);j<=min(6,i);j++){ if(j%2==0){ nP[i/2]+=Q[j]*P[i-j]; continue; } nP[i/2]-=Q[j]*P[i-j]; } } copy(begin(nP),end(nP),begin(P)); } if(N%2==1){ for(int i=1;i<12;i+=2){ nP[i/2]=0; for(int j=max(0,i-5);j<=min(6,i);j++){ if(j%2==0){ nP[i/2]+=Q[j]*P[i-j]; continue; } nP[i/2]-=Q[j]*P[i-j]; } } copy(begin(nP),end(nP),begin(P)); } copy(begin(Q_data[count]),end(Q_data[count]),begin(Q)); N/=2; continue; } mp[M]=P[0]; return P[0]; } template struct matrix{ int H,W; T table[6][6]; matrix(int h,int w) : H(h),W(w){ for(int i=0;i=0;i--){ W+=Bostan_Mori(N*M+i)*(i-4); } for(;k--;){ int C; scanf("%d",&C); matrix P(6,6); matrix Q(6,1); for(int i=0;i<6;i++){ if(i==C){ continue; } if(C=0;i--){ if(i==C){ continue; } if(j==C){ continue; } if(i=0;i--){ if(P.table[i][0]==0){ continue; } if(i=0;i--){ ans+=R[i]*(i-4); } ans=W-ans; //printf("%d\n",ans.val()); } printf("%d\n",(int)(mp.size())); }