結果

問題 No.621 3 x N グリッド上のドミノの置き方の数
ユーザー rickythetarickytheta
提出日時 2017-12-21 00:22:44
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 12 ms / 3,000 ms
コード長 2,007 bytes
コンパイル時間 1,356 ms
コンパイル使用メモリ 145,184 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-08-22 11:07:31
合計ジャッジ時間 4,381 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 2 ms
4,380 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 2 ms
4,380 KB
testcase_04 AC 2 ms
4,380 KB
testcase_05 AC 2 ms
4,376 KB
testcase_06 AC 2 ms
4,376 KB
testcase_07 AC 2 ms
4,376 KB
testcase_08 AC 2 ms
4,380 KB
testcase_09 AC 2 ms
4,380 KB
testcase_10 AC 2 ms
4,376 KB
testcase_11 AC 2 ms
4,380 KB
testcase_12 AC 2 ms
4,376 KB
testcase_13 AC 3 ms
4,376 KB
testcase_14 AC 3 ms
4,376 KB
testcase_15 AC 4 ms
4,376 KB
testcase_16 AC 4 ms
4,376 KB
testcase_17 AC 5 ms
4,376 KB
testcase_18 AC 5 ms
4,376 KB
testcase_19 AC 5 ms
4,376 KB
testcase_20 AC 6 ms
4,376 KB
testcase_21 AC 6 ms
4,380 KB
testcase_22 AC 7 ms
4,380 KB
testcase_23 AC 7 ms
4,376 KB
testcase_24 AC 7 ms
4,376 KB
testcase_25 AC 8 ms
4,376 KB
testcase_26 AC 8 ms
4,376 KB
testcase_27 AC 9 ms
4,376 KB
testcase_28 AC 9 ms
4,376 KB
testcase_29 AC 10 ms
4,376 KB
testcase_30 AC 9 ms
4,380 KB
testcase_31 AC 9 ms
4,376 KB
testcase_32 AC 9 ms
4,380 KB
testcase_33 AC 10 ms
4,380 KB
testcase_34 AC 9 ms
4,380 KB
testcase_35 AC 9 ms
4,376 KB
testcase_36 AC 9 ms
4,376 KB
testcase_37 AC 10 ms
4,376 KB
testcase_38 AC 10 ms
4,380 KB
testcase_39 AC 10 ms
4,376 KB
testcase_40 AC 10 ms
4,376 KB
testcase_41 AC 10 ms
4,376 KB
testcase_42 AC 10 ms
4,380 KB
testcase_43 AC 10 ms
4,376 KB
testcase_44 AC 10 ms
4,376 KB
testcase_45 AC 10 ms
4,376 KB
testcase_46 AC 10 ms
4,376 KB
testcase_47 AC 10 ms
4,380 KB
testcase_48 AC 10 ms
4,376 KB
testcase_49 AC 10 ms
4,376 KB
testcase_50 AC 10 ms
4,376 KB
testcase_51 AC 10 ms
4,380 KB
testcase_52 AC 10 ms
4,380 KB
testcase_53 AC 10 ms
4,376 KB
testcase_54 AC 10 ms
4,380 KB
testcase_55 AC 10 ms
4,376 KB
testcase_56 AC 10 ms
4,376 KB
testcase_57 AC 9 ms
4,376 KB
testcase_58 AC 9 ms
4,380 KB
testcase_59 AC 12 ms
4,380 KB
testcase_60 AC 12 ms
4,376 KB
testcase_61 AC 12 ms
4,380 KB
testcase_62 AC 7 ms
4,376 KB
testcase_63 AC 7 ms
4,376 KB
testcase_64 AC 7 ms
4,380 KB
testcase_65 AC 7 ms
4,380 KB
testcase_66 AC 7 ms
4,380 KB
testcase_67 AC 7 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;

typedef int _loop_int;
#define REP(i,n) for(_loop_int i=0;i<(_loop_int)(n);++i)
#define FOR(i,a,b) for(_loop_int i=(_loop_int)(a);i<(_loop_int)(b);++i)
#define FORR(i,a,b) for(_loop_int i=(_loop_int)(b)-1;i>=(_loop_int)(a);--i)

#define DEBUG(x) cout<<#x<<": "<<x<<endl
#define DEBUG_VEC(v) cout<<#v<<":";REP(i,v.size())cout<<" "<<v[i];cout<<endl
#define ALL(a) (a).begin(),(a).end()

#define CHMIN(a,b) a=min((a),(b))
#define CHMAX(a,b) a=max((a),(b))

// mod
const ll MOD = 1000000007ll;
#define FIX(a) ((a)%MOD+MOD)%MOD

// floating
typedef double Real;
const Real EPS = 1e-11;
#define EQ0(x) (abs(x)<EPS)
#define EQ(a,b) (abs(a-b)<EPS)
typedef complex<Real> P;

ll n;

ll A[27][27];
ll tmp[27][27];
ll r[27][27];

#define UNPACK(a,b,c,x) int a=(x)%3, b=((x)/3)%3, c=(x)/9

void matmul(ll r[27][27], ll a[27][27], ll b[27][27]){
  REP(i,27)REP(j,27)tmp[i][j] = 0;
  REP(i,27)REP(j,27)REP(k,27)(tmp[i][j] += a[i][k] * b[k][j] % MOD) %= MOD;
  REP(i,27)REP(j,27)r[i][j] = tmp[i][j];
}

int main(){
  REP(i,27){
    UNPACK(a,b,c,i);
    REP(j,27){
      UNPACK(x,y,z,j);
      // has space
      if(x==0 && y==0)continue;
      if(y==0 && z==0)continue;
      if(a==0 && x==0)continue;
      if(b==0 && y==0)continue;
      if(c==0 && z==0)continue;
      // not connect check
      if(a==1 && x!=2)continue;
      if(b==1 && y!=2)continue;
      if(c==1 && z!=2)continue;
      // 2 put check
      bool conx = a!=1 && x==2;
      bool cony = b!=1 && y==2;
      bool conz = c!=1 && z==2;
      if(conx && cony && conz)continue;
      if(conx && !cony)continue;
      if(cony && !conx && !conz)continue;
      if(conz && !cony)continue;
      // ok
      A[j][i] = 1;
    }
  }
  REP(i,27)r[i][i] = 1;
  cin>>n;
  n++;
  while(n){
    if(n&1)matmul(r,r,A);
    matmul(A,A,A);
    n>>=1;
  }
  printf("%lld\n", r[13][26]);
  return 0;
}
0