#include #include #include using namespace std; using namespace atcoder; using mint=modint998244353; int mem=0; vector> Q_data(60); mint Bostan_Mori(long long N){ int d=6; vector Q(d+1); if(mem>0){ Q=Q_data[0]; } if(mem==0){ Q[0]=1; for(int i=1;i<=d;i++){ Q[i]=-1; } Q_data[0]=Q; mem++; } vector P(d); P[0]=1; int count=0; while(N>0){ count++; if(mem<=count){ vector nQ(d+1); for(int i=0;i<=d;i++){ for(int j=max(0,2*i-d);j<=min(2*i,d);j++){ if(j%2==0){ nQ[i]+=Q[j]*Q[2*i-j]; continue; } nQ[i]-=Q[j]*Q[2*i-j]; } } Q_data[mem]=nQ; mem++; } vector nP(d); if(N%2==0){ for(int i=0;i<2*d;i+=2){ for(int j=max(0,i-d+1);j<=min(d,i);j++){ if(j%2==0){ nP[i/2]+=Q[j]*P[i-j]; continue; } nP[i/2]-=Q[j]*P[i-j]; } } P=nP; } if(N%2==1){ for(int i=1;i<2*d;i+=2){ for(int j=max(0,i-d+1);j<=min(d,i);j++){ if(j%2==0){ nP[i/2]+=Q[j]*P[i-j]; continue; } nP[i/2]-=Q[j]*P[i-j]; } } P=nP; } Q=Q_data[count]; N/=2; continue; } return P[0]; } template struct matrix{ int H,W; vector> table; matrix(int h,int w) : H(h),W(w){ table.resize(h,vector(w)); } int size(){ return H; } matrix pow(long long N){ assert(H==W); assert(0<=N); matrix x=*this; matrix r(H,H); for(int i=0;i>=1; } return r; } vector& operator[](int index){ assert(index A(6); A[0]=1; A[1]=1; A[2]=2; A[3]=4; A[4]=8; A[5]=16; vector c(6,mint(1)); for(;k--;){ int C; scanf("%d",&C); vector> P(6,vector(6)); //P[i][j]:X=iからスタートして,X=N+j(ピッタリ)で終わるサイコロの振り方であって,常にX \not \equiv C (mod N)である場合 for(int i=0;i<6;i++){ for(int j=0;j<6;j++){ if(i==C){ P[i][j]=0; continue; } if(j==C){ P[i][j]=0; continue; } if(i Q(6,6); for(int i=0;i<6;i++){ for(int j=0;j<6;j++){ Q[i][j]=P[j][i]; } } Q=Q.pow(M-2);//to do:M=1の場合に場合分けして書く matrix R(6,1); for(int i=0;i<6;i++){ R[i][0]=P[0][i]; } Q*=R; vector S(6); for(int i=0;i<6;i++){ if(Q[i][0]==0){ continue; } for(int j=0;j<6;j++){ //iからN+jに動く //N~N+jの間のCについては無視 if(i=0;i--){ for(int j=i-1;j>=0;j--){ S[i]-=S[j]; } ans+=S[i]; } printf("%d\n",ans.val()); } }