結果
| 問題 |
No.2159 Filling 4x4 array
|
| コンテスト | |
| ユーザー |
hotman78
|
| 提出日時 | 2022-12-02 15:41:45 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 1,677 ms / 5,000 ms |
| コード長 | 1,591 bytes |
| コンパイル時間 | 2,041 ms |
| コンパイル使用メモリ | 200,412 KB |
| 最終ジャッジ日時 | 2025-02-09 03:12:59 |
|
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 5 |
| other | AC * 45 |
ソースコード
#include<bits/stdc++.h>
#include<atcoder/modint>
using namespace std;
using lint=long long;
using mint=atcoder::static_modint<998244353>;
#define rep(i,n) for(int i=0;i<int(n);++i)
int main(){
array<lint,8> v;
rep(i,8)cin>>v[i];
rep(i,8)v[i]-=4;
if(v[0]+v[1]+v[2]+v[3]!=v[4]+v[5]+v[6]+v[7]){
cout<<0<<endl;
return 0;
}
swap(v[3],v[6]);
swap(v[6],v[7]);
vector<vector<mint>> dp(31,vector<mint>(1<<14));
dp[0][0]=1;
array<lint,(1<<9)>a;
a.fill(0);
rep(i,1<<9){
rep(j,3)rep(k,3){
if(i>>(j*3+k)&1){
a[i]+=1<<(j*2);
a[i]+=1<<((k+3)*2);
}
}
}
array<lint,30>v2;
v2.fill(0);
rep(d,30){
rep(j,7){
v2[d]+=(v[j]>>d&1)<<(j*2);
}
}
// 今考えている桁
rep(d,30){
//flag : 解説におけるh',w' を2bitずつに割り当てたもの
rep(flag,1<<14){
// 左上9マスの選び方
rep(i,1<<9){
// tmp : 前半3制約の内遇奇が合わない物
// 21 = 1+4+4^2
const lint tmp=((flag^a[i]^v2[d])&21LL);
// cnt : 7つ目の制約の変化bit
const lint cnt=__builtin_popcountll(tmp)<<12;
// tmp2 : 遷移の値を加算したもの
// 5461 = 1+4+...+4^16
const lint tmp2=(flag+a[i]+cnt)+(((flag^(a[i]+cnt))^v2[d])&5461LL)-v2[d];
dp[d+1][tmp2/2]+=dp[d][flag];
}
}
}
cout<<dp[30][0].val()<<endl;
}
hotman78