結果
| 問題 | No.3736 Purely Bool Hell |
| コンテスト | |
| ユーザー |
 uruzunyaa
|
| 提出日時 | 2026-09-19 02:59:01 |
| 言語 | C++23 (gcc 15.3.0 + boost 1.92.0 + ACL) |
| 結果 |
AC
不安定
|
| 実行時間 | 938 ms / 3,000 ms |
| + 118µs | |
| コード長 | 6,432 bytes |
| 記録 | |
| コンパイル時間 | 2,853 ms |
| コンパイル使用メモリ | 364,196 KB |
| 実行使用メモリ | 7,552 KB |
| 最終ジャッジ日時 | 2026-09-19 13:26:06 |
| 合計ジャッジ時間 | 12,200 ms |
|
ジャッジサーバーID （参考情報） |
judge1_0 / judge3_0 |
（要ログイン）
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 |
| other | AC * 39 |
ソースコード
//#pragma GCC optimize("O3")
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define rep(i,n) for (ll i=0;i<(ll)n;i++)
#define rrep(i,n) for (ll i=(n)-1;i>=(ll)0;i--)
#define loop(i,m,n) for(ll i=m;i<=(ll)n;i++)
#define rloop(i,m,n) for(ll i=m;i>=(ll)n;i--)
#define vl vector<ll>
#define vvl vector<vl>
#define vvvl vector<vvl>
#define vdbg(a) rep(ii,a.size()){cout<<a[ii]<<" ";}cout<<endl;
#define vpdbg(a) rep(ii,a.size()){cout<<"{"<<a[ii].first<<","<<a[ii].second<<"} ";}cout<<endl;
#define vvdbg(a) rep(ii,a.size()){rep(jj,a[ii].size()){cout<<a[ii][jj]<<" ";}cout<<endl;}
#define setdbg(a) for(const auto & ii:a){cout<<ii<<" ";}cout<<endl;
#define inf 4000000000000000000LL
#define mod 998244353LL
//#define mod 1000000007LL
#define eps 0.000000001
#define circlepi 3.14159265358979323846
random_device rnd;// 非決定的な乱数生成器
mt19937 mt(rnd());// メルセンヌ・ツイスタの32ビット版、引数は初期シード
//#include<boost/multiprecision/cpp_int.hpp>
//#define bbi boost::multiprecision::cpp_int
//#include<atcoder/lazysegtree>
//整数同士の累乗の計算をする。
ll power(ll A, ll B) {
ll result = 1;
for (ll i=0;i<B;i++){
result *= A;
}
return result;
}
// nのk乗をmodで割った余りを計算
ll power_mod(ll n, ll k){
n%=mod;
ll ans = 1;
while (k > 0){
if ((k&1) ==1)ans=(ans*n)%mod;
n=n*n%mod;
k >>= 1;
}
return ans;
}
//受け取った2次元文字の外側に、文字pをコーティングする。
vector<string> pad(vector<string> &s,char p){
ll h=s.size();
ll w=s[0].size();
vector<string> res(h+2,string(w+2,p));
rep(i,h)rep(j,w)res[i+1][j+1]=s[i][j];
return res;
}
// Union-Find
struct UnionFind {
vector<int> par, siz;
UnionFind(int n) : par(n, -1) , siz(n, 1) { }
// 根を求める
int root(int x) {
if (par[x] == -1) return x;
else return par[x] = root(par[x]);
}
// x と y が同じグループに属するかどうか (根が一致するかどうか)
bool issame(int x, int y) {
return root(x) == root(y);
}
// x を含むグループと y を含むグループとを併合する
bool unite(int x, int y) {
x = root(x), y = root(y);
if (x == y) return false;
if (siz[x] < siz[y]) swap(x, y);
par[y] = x;
siz[x] += siz[y];
return true;
}
// x を含むグループのサイズ
int size(int x) {
return siz[root(x)];
}
};
//グリッド問題等用
vl dx={1,0,-1,0};
vl dy={0,1,0,-1};
void solve(){
ll n;
cin>>n;
vl x(n),y(n),z(n+n-1);
//横ANDがX,縦ORがY,左下→右上向き斜めXORがZ
rep(i,n)cin>>x[i];
rep(i,n)cin>>y[i];
rep(i,n+n-1)cin>>z[i];
vvl ans(n,vl(n,0));
rep(b,30){
vl xx=x,yy=y,zz=z;
vl xcnt=x,ycnt=y,zcnt=z;
rep(i,n){
xx[i]=(xx[i]>>b)%2;
yy[i]=(yy[i]>>b)%2;
xcnt[i]=n;
ycnt[i]=n;
}
rep(i,n+n-1){
zz[i]=(zz[i]>>b)%2;
zcnt[i]=min(i+1,n+n-i-1);
}
vvl tmp(n,vl(n,-1));
ll type=2;
rep(i,n){
if(xx[i]==1){
if(type==1){
cout<<-1<<endl;
return;
}
type=0;
rep(j,n){
tmp[i][j]=1;
xcnt[i]=0;
ycnt[j]=0;
zcnt[i+j]--;
zz[i+j]^=1;
}
}
if(yy[i]==0){
if(type==0){
cout<<-1<<endl;
return;
}
type=1;
rep(j,n){
tmp[j][i]=0;
xcnt[j]=0;
ycnt[i]=0;
zcnt[i+j]--;
}
}
}
if(type==2){
if(n<=2){
//3以下は全探索
rep(i,1LL<<(n*n)){
vvl testing(n,vl(n));
rep(j,n*n){
if(i&(1LL<<j))testing[j/n][j%n]=1;
else testing[j/n][j%n]=0;
}
vl ansx(n,1),ansy(n,0),ansz(n+n-1,0);
rep(i,n){
rep(j,n){
ansx[i]&=testing[i][j];
ansy[j]|=testing[i][j];
ansz[i+j]^=testing[i][j];
}
}
if(ansx==xx&&ansy==yy&&ansz==zz){
tmp=testing;
break;
}
}
if(tmp[0][0]==-1){
cout<<-1<<endl;
return;
}
rep(i,n)rep(j,n)ans[i][j]+=tmp[i][j]*(1LL<<b);
continue;
}
if(n==3){
//3以下は全探索
rep(i,1LL<<((n-1)*(n-1))){
vvl testing(n,vl(n));
rep(j,(n-1)*(n-1)){
if(i&(1LL<<j))testing[j/(n-1)][j%(n-1)+1]=1;
else testing[j/(n-1)][j%(n-1)+1]=0;
}
testing[0][0]=zz[0];
testing[1][0]=zz[1]^testing[0][1];
testing[2][0]=zz[2]^testing[1][1]^testing[0][2];
testing[2][1]=zz[3]^testing[1][2];
testing[2][2]=zz[4];
vl ansx(n,1),ansy(n,0);
rep(i,n){
rep(j,n){
ansx[i]&=testing[i][j];
ansy[j]|=testing[i][j];
}
}
if(ansx==xx&&ansy==yy){
tmp=testing;
break;
}
}
if(tmp[0][0]==-1){
cout<<-1<<endl;
return;
}
rep(i,n)rep(j,n)ans[i][j]+=tmp[i][j]*(1LL<<b);
continue;
}
//n>=4のみ解く。
tmp[0][1]=0;
tmp[1][1]=0;
tmp[n-2][1]=1;
tmp[n-2][0]=1;
zz[n-2]^=1;
zz[n-1]^=1;
zcnt[1]--;
zcnt[2]--;
zcnt[n-2]--;
zcnt[n-1]--;
loop(i,2,n-1){
tmp[1][i]=1;
zz[1+i]^=1;
zcnt[1+i]--;
tmp[i][n-2]=0;
zcnt[n-2+i]--;
}
//XORを調整する。
rep(i,n)rep(j,n){
if(tmp[i][j]==-1){
tmp[i][j]=zz[i+j];
zz[i+j]=0;
}
}
}else{
if(type==0){
rep(i,n)rep(j,n){
if(tmp[i][j]==-1&&zcnt[i+j]>=2){
tmp[i][j]=0;
xcnt[i]--;
zcnt[i+j]--;
}
}
}else{
rep(j,n)rep(i,n){
if(tmp[i][j]==-1&&zcnt[i+j]>=2){
tmp[i][j]=1;
zz[i+j]^=1;
ycnt[j]--;
zcnt[i+j]--;
}
}
}
//zcnt埋めが終わった奴をチェック
rep(i,n+n-1){
if(zcnt[i]==0){
if(zz[i]==1){
cout<<-1<<endl;
return;
}
}
}
//zcntが1な場所を埋める
rep(i,n)rep(j,n){
if(zcnt[i+j]==1&&tmp[i][j]==-1){
tmp[i][j]=zz[i+j];
zcnt[i+j]=0;
xcnt[i]=max(0LL,xcnt[i]-1);
ycnt[j]=max(0LL,ycnt[j]-1);
}
}
}
vl ansx(n,1),ansy(n,0);
rep(i,n){
rep(j,n){
ansx[i]&=tmp[i][j];
ansy[j]|=tmp[i][j];
}
}
ll hokan=0;
rep(i,n){
if(tmp[i][n-i-1]==type){
hokan=i;
break;
}
}
rep(i,n){
if(ansx[i]!=xx[i]){
swap(tmp[hokan][n-hokan-1],tmp[i][n-i-1]);
}
if(ansy[i]!=yy[i]){
swap(tmp[hokan][n-hokan-1],tmp[n-i-1][i]);
}
}
rep(i,n)rep(j,n)ans[i][j]+=tmp[i][j]*(1LL<<b);
}
vl ansx(n,(1LL<<30)-1),ansy(n,0);
rep(i,n){
rep(j,n){
ansx[i]&=ans[i][j];
ansy[j]|=ans[i][j];
}
}
if(ansx!=x||ansy!=y){
cout<<-1<<endl;
return;
}
vvdbg(ans);
return;
}
//メイン
int main(){
ll t;
cin>>t;
rep(i,t)solve();
}