結果

問題 No.3739 Stronger Network
コンテスト
ユーザー 👑 Nachia
提出日時 2026-09-19 15:33:21
言語 C++17
(gcc 15.3.0 + boost 1.92.0 + ACL)
コンパイル:
g++-15 -O2 -lm -std=c++17 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 34 ms / 2,000 ms
+ 468µs
コード長 2,116 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 632 ms
コンパイル使用メモリ 107,760 KB
実行使用メモリ 9,952 KB
最終ジャッジ日時 2026-09-19 15:33:29
合計ジャッジ時間 4,692 ms
ジャッジサーバーID
（参考情報） 		judge4_0 / judge5_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 48
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#ifdef NACHIA
#define _GLIBCXX_DEBUG
#else
// disable assert
#define NDEBUG
#endif
#include <iostream>
#include <string>
#include <vector>
#include <algorithm>
using namespace std;
using ll = long long;
const ll INF = 1ll << 60;
#define REP(i,n) for(ll i=0; i<ll(n); i++)
template <class T> using V = vector<T>;
template <class A, class B> void chmax(A& l, const B& r){ if(l < r) l = r; }
template <class A, class B> void chmin(A& l, const B& r){ if(r < l) l = r; }

V<ll> seq(ll n){
  if(n == 1) return {0};
  auto A = seq(n / 2);
  V<ll> ans;
  REP(i,A.size()){
    ans.push_back(A[i] * 2 + (0 ^ (i%2)));
    ans.push_back(A[i] * 2 + (1 ^ (i%2)));
  }
  return ans;
}

V<ll> row(ll N){
  V<ll> ans;
  ll n = 1; while(n*2 <= N) n *= 2;
  auto buf = seq(n);
  REP(i,buf.size()) if(i%2 == 0){
    ans.push_back(buf[i]);
    if(buf[i] + n < N){
      ans.push_back(buf[i] + n);
      ans.push_back(buf[i+1] + n);
    }
    ans.push_back(buf[i+1]);
  }
  return ans;
}

struct Node {
  ll val = INF;
  V<ll> A, B;
};

bool operator<(const Node& l, const Node& r){ return l.val < r.val; }

Node do_dp(ll A, ll B){
  ll a = 0; while((1ll << a) <= A) a++;
  ll b = 0; while((1ll << b) <= B) b++;
  V<V<Node>> dp(a+1, V<Node>(b+1));
  dp[0][0].val = 0;
  REP(i,a+1) REP(j,b+1){
    if(i){
      Node nx = dp[i-1][j];
      nx.A.push_back(1ll << (i+j-1));
      if(A & (1ll << (i-1))) nx.val += (1ll << (i+j-1));
      chmin(dp[i][j], nx);
    }
    if(j){
      Node nx = dp[i][j-1];
      nx.B.push_back(1ll << (i+j-1));
      if(B & (1ll << (j-1))) nx.val += (1ll << (i+j-1));
      chmin(dp[i][j], nx);
    }
  }
  return dp[a][b];
}

void testcase(){
  ll H, W; cin >> H >> W;
  if(H%2 != 0 || W%2 != 0){ cout << "-1\n"; return; }

  auto rows = row(H);
  auto cols = row(W);
  auto mp = do_dp(H-1, W-1);

  REP(y,H){
    REP(x,W){
      if(x) cout << " ";
      ll val = 0;
      REP(a,mp.A.size()) val += (rows[y] >> a) % 2 * mp.A[a];
      REP(a,mp.B.size()) val += (cols[x] >> a) % 2 * mp.B[a];
      cout << val;
    } cout << "\n";
  }
}

int main(){
  cin.tie(0)->sync_with_stdio(0);
  testcase();
  return 0;
}
0