#1189557 (C++23) No.3744 XY Tiling

• 問題: No.3744 XY Tiling
• ユーザー: noya2
• 提出日時: 2026-09-19 17:59:55
• 言語: C++23 (gcc 15.3.0 + boost 1.92.0 + ACL)
• 結果: RE
• 実行時間: -
• コード長: 12,769 bytes
• コンパイル時間: 3,224 ms
• コンパイル使用メモリ: 361,952 KB
• 実行使用メモリ: 10,052 KB
• 最終ジャッジ日時: 2026-09-19 18:01:09
• 合計ジャッジ時間: 6,518 ms
• ジャッジサーバーID （参考情報）: judge4_0 / judge2_0

サブタスク	配点	結果
部分点	60 %	AC * 19
満点	40 %	AC * 39 WA * 16 RE * 5
合計	5 * 60% = 300 点

ソースコード

#line 2 "/Users/noya2/Desktop/Noya2_library/template/template.hpp"
using namespace std;

#include<bits/stdc++.h>
#line 1 "/Users/noya2/Desktop/Noya2_library/template/inout_old.hpp"
namespace noya2 {

template <typename T, typename U>
ostream &operator<<(ostream &os, const pair<T, U> &p){
    os << p.first << " " << p.second;
    return os;
}
template <typename T, typename U>
istream &operator>>(istream &is, pair<T, U> &p){
    is >> p.first >> p.second;
    return is;
}

template <typename T>
ostream &operator<<(ostream &os, const vector<T> &v){
    int s = (int)v.size();
    for (int i = 0; i < s; i++) os << (i ? " " : "") << v[i];
    return os;
}
template <typename T>
istream &operator>>(istream &is, vector<T> &v){
    for (auto &x : v) is >> x;
    return is;
}

void in() {}
template <typename T, class... U>
void in(T &t, U &...u){
    cin >> t;
    in(u...);
}

void out() { cout << "\n"; }
template <typename T, class... U, char sep = ' '>
void out(const T &t, const U &...u){
    cout << t;
    if (sizeof...(u)) cout << sep;
    out(u...);
}

template<typename T>
void out(const vector<vector<T>> &vv){
    int s = (int)vv.size();
    for (int i = 0; i < s; i++) out(vv[i]);
}

struct IoSetup {
    IoSetup(){
        cin.tie(nullptr);
        ios::sync_with_stdio(false);
        cout << fixed << setprecision(15);
        cerr << fixed << setprecision(7);
    }
} iosetup_noya2;

} // namespace noya2
#line 1 "/Users/noya2/Desktop/Noya2_library/template/const.hpp"
namespace noya2{

const int iinf = 1'000'000'007;
const long long linf = 2'000'000'000'000'000'000LL;
const long long mod998 =  998244353;
const long long mod107 = 1000000007;
const long double pi = 3.14159265358979323;
const vector<int> dx = {0,1,0,-1,1,1,-1,-1};
const vector<int> dy = {1,0,-1,0,1,-1,-1,1};
const string ALP = "ABCDEFGHIJKLMNOPQRSTUVWXYZ";
const string alp = "abcdefghijklmnopqrstuvwxyz";
const string NUM = "0123456789";

void yes(){ cout << "Yes\n"; }
void no(){ cout << "No\n"; }
void YES(){ cout << "YES\n"; }
void NO(){ cout << "NO\n"; }
void yn(bool t){ t ? yes() : no(); }
void YN(bool t){ t ? YES() : NO(); }

} // namespace noya2
#line 2 "/Users/noya2/Desktop/Noya2_library/template/utils.hpp"

#line 6 "/Users/noya2/Desktop/Noya2_library/template/utils.hpp"

namespace noya2{

unsigned long long inner_binary_gcd(unsigned long long a, unsigned long long b){
    if (a == 0 || b == 0) return a + b;
    int n = __builtin_ctzll(a); a >>= n;
    int m = __builtin_ctzll(b); b >>= m;
    while (a != b) {
        int mm = __builtin_ctzll(a - b);
        bool f = a > b;
        unsigned long long c = f ? a : b;
        b = f ? b : a;
        a = (c - b) >> mm;
    }
    return a << std::min(n, m);
}

template<typename T> T gcd_fast(T a, T b){ return static_cast<T>(inner_binary_gcd(std::abs(a),std::abs(b))); }

long long sqrt_fast(long long n) {
    if (n <= 0) return 0;
    long long x = sqrt(n);
    while ((x + 1) * (x + 1) <= n) x++;
    while (x * x > n) x--;
    return x;
}

template<typename T> T floor_div(const T n, const T d) {
    assert(d != 0);
    return n / d - static_cast<T>((n ^ d) < 0 && n % d != 0);
}

template<typename T> T ceil_div(const T n, const T d) {
    assert(d != 0);
    return n / d + static_cast<T>((n ^ d) >= 0 && n % d != 0);
}

template<typename T> void uniq(std::vector<T> &v){
    std::sort(v.begin(),v.end());
    v.erase(unique(v.begin(),v.end()),v.end());
}

template <typename T, typename U> inline bool chmin(T &x, U y) { return (y < x) ? (x = y, true) : false; }

template <typename T, typename U> inline bool chmax(T &x, U y) { return (x < y) ? (x = y, true) : false; }

template<typename T> inline bool range(T l, T x, T r){ return l <= x && x < r; }

} // namespace noya2
#line 8 "/Users/noya2/Desktop/Noya2_library/template/template.hpp"

#define rep(i,n) for (int i = 0; i < (int)(n); i++)
#define repp(i,m,n) for (int i = (m); i < (int)(n); i++)
#define reb(i,n) for (int i = (int)(n-1); i >= 0; i--)
#define all(v) (v).begin(),(v).end()

using ll = long long;
using ld = long double;
using uint = unsigned int;
using ull = unsigned long long;
using pii = pair<int,int>;
using pll = pair<ll,ll>;
using pil = pair<int,ll>;
using pli = pair<ll,int>;

namespace noya2{

/*　~ (. _________ . /)　*/

}

using namespace noya2;


#line 2 "c.cpp"

#line 2 "/Users/noya2/Desktop/Noya2_library/graph/grid.hpp"

#line 5 "/Users/noya2/Desktop/Noya2_library/graph/grid.hpp"
#include <ranges>

namespace noya2 {

struct grid {
    int h, w;
    grid (int _h = 0, int _w = 0) : h(_h), w(_w) {}
    int idx(int x, int y) const {
        return x * w + y;
    }
    std::pair<int,int> pos(int id) const {
        return {id / w, id % w};
    }
    bool operator()(int x, int y) const {
        return 0 <= x && x < h && 0 <= y && y < w;
    }
    bool operator()(int id) const {
        return operator()(id / w, id % w);
    }
    static constexpr std::array<std::pair<int,int>, 4> dxys = {std::pair<int,int>{0, 1}, {-1, 0}, {0, -1}, {1, 0}};
    auto adj4(int x, int y) const {
        return dxys | std::views::transform([x, y](std::pair<int,int> dxy){
            return std::pair<int,int>(dxy.first + x, dxy.second + y);
        }) | std::views::filter([this](std::pair<int,int> xy){
            return operator()(xy.first, xy.second);
        });
    }
    auto adj4(int id) const {
        return adj4(id / w, id % w) | std::ranges::views::transform([this](std::pair<int,int> xy){
            return idx(xy.first, xy.second);
        });
    }
    int size() const {
        return h * w;
    }
};

} // namespace noya2
#line 4 "c.cpp"

struct dom {
    int x1, y1, c1;
    int x2, y2, c2;
    friend ostream &operator<<(ostream &os, dom d){
        return os << d.x1+1 << ' ' << d.y1+1 << ' ' << d.c1 << ' ' << d.x2+1 << ' ' << d.y2+1 << ' ' << d.c2;
    }
};

vector<dom> make_sml_even(int h, int w){
    assert(h % 2 == 0);
    assert(h <= w);
    vector<dom> ans;
    int lx = 0, rx = h, ly = 0, ry = w;
    int p = 1;
    while (true){
        if (rx - lx == 2){
            int u = lx, d = lx+1;
            repp(j,ly,ry){
                ans.push_back({u,j,p,d,j,2});
            }
            break;
        }
        if (rx - lx == 4){
            int uc = p, dc = 4-p;
            // ans.push_back({lx+2,ly,dc,lx+2,ly+1,2});
            // ans.push_back({lx+3,ly,dc,lx+3,ly+1,2});
            // ans.push_back({lx,ly,uc,lx+1,ly,2});
            // ans.push_back({lx,ly+1,uc,lx+1,ly+1,2});
            repp(j,ly,ry){
                ans.push_back({lx,j,uc,lx+1,j,2});
                ans.push_back({lx+2,j,2,lx+3,j,dc});
            }
            break;
        }
        repp(i,lx,rx-2){
            ans.push_back({i,ly,p,i,ly+1,2});
        }
        repp(j,ly,ry-2){
            ans.push_back({rx-2,j,p,rx-1,j,2});
        }
        repp(i,lx+2,rx){
            ans.push_back({i,ry-2,2,i,ry-1,4-p});
        }
        repp(j,ly+2,ry){
            ans.push_back({lx,j,2,lx+1,j,4-p});
        }
        lx += 2, rx -= 2;
        ly += 2, ry -= 2;
        p = 4-p;
    }
    return ans;
}

vector<dom> make1(int h){
    vector<dom> ans;
    int lx = 0, rx = h, ly = 0, ry = h+1;
    int p = 1;
    while (true){
        if (rx - lx == 1){
            assert(ry - ly == 2);
            ans.push_back({lx,ly,4-p,lx,ly+1,2});
            break;
        }
        repp(i,lx,rx-2){
            ans.push_back({i,ly,2,i,ly+1,p});
        }
        repp(j,ly,ry-2){
            ans.push_back({rx-2,j,2,rx-1,j,4-p});
        }
        repp(i,lx+2,rx){
            ans.push_back({i,ry-2,4-p,i,ry-1,2});
        }
        repp(j,ly+2,ry){
            ans.push_back({lx,j,p,lx+1,j,2});
        }
        lx += 2, rx -= 2;
        ly += 2, ry -= 2;
        p = 4-p;
    }
    return ans;
}

vector<dom> make3(int h){
    vector<dom> ans;
    int lx = 0, rx = h, ly = 0, ry = h+1;
    int p = 1;
    while (true){
        if (rx - lx == 3){
            assert(ry - ly == 4);
            ans.push_back({lx,ly,2,lx,ly+1,p});
            ans.push_back({lx,ly+2,p,lx,ly+3,2});
            repp(j,ly,ry){
                ans.push_back({lx+1,j,p,lx+2,j,2});
            }
            break;
        }
        repp(i,lx,rx-2){
            ans.push_back({i,ly,p,i,ly+1,2});
        }
        repp(j,ly,ry-2){
            ans.push_back({rx-2,j,p,rx-1,j,2});
        }
        repp(i,lx+2,rx){
            ans.push_back({i,ry-2,2,i,ry-1,4-p});
        }
        repp(j,ly+2,ry){
            ans.push_back({lx,j,2,lx+1,j,4-p});
        }
        lx += 2, rx -= 2;
        ly += 2, ry -= 2;
        p = 4-p;
    }
    return ans;
}

vector<dom> make(int h, int w){
    if (h > w){
        auto ans = make(w,h);
        for (auto &[x1, y1, c1, x2, y2, c2] : ans){
            swap(x1,y1);
            swap(x2,y2);
        }
        return ans;
    }
    if (h % 2 == 0){
        return make_sml_even(h,w);
    }
    assert(w % 2 == 0);
    vector<dom> ans = [&]{
        if (h % 4 == 1){
            return make1(h);
        }
        else {
            return make3(h);
        }
    }();
    int nw = w-(h+1);
    if (nw == 0){
        return ans;
    }
    auto adds = [&]{
        if (nw != 4){
            return make(h,nw);
        }
        vector<dom> ret;
        ret.push_back({0,0,2,1,0,1});
        ret.push_back({0,1,2,1,1,1});
        ret.push_back({0,2,2,0,3,3});
        ret.push_back({1,2,2,1,3,3});
        repp(i,2,h){
            ret.push_back({i,0,1,i,1,2});
            ret.push_back({i,2,2,i,3,3});
        }
        return ret;
    }();
    for (auto [x1, y1, c1, x2, y2, c2] : adds){
        y1 = nw-1-y1;
        y2 = nw-1-y2;
        y1 += h+1;
        y2 += h+1;
        ans.push_back({x1,y1,c1,x2,y2,c2});
    }
    return ans;
}

#line 2 "/Users/noya2/Desktop/Noya2_library/data_structure/dsu.hpp"

#line 6 "/Users/noya2/Desktop/Noya2_library/data_structure/dsu.hpp"

namespace noya2{

struct dsu {
  public:
    dsu() : _n(0), _cc(0) {}
    dsu(int n) : _n(n), _cc(n), parent_or_size(n, -1) {}

    int merge(int a, int b) {
        assert(0 <= a && a < _n);
        assert(0 <= b && b < _n);
        int x = leader(a), y = leader(b);
        if (x == y) return x;
        if (-parent_or_size[x] < -parent_or_size[y]) std::swap(x, y);
        parent_or_size[x] += parent_or_size[y];
        parent_or_size[y] = x;
        _cc--;
        return x;
    }

    bool same(int a, int b) {
        assert(0 <= a && a < _n);
        assert(0 <= b && b < _n);
        return leader(a) == leader(b);
    }

    int leader(int a) {
        assert(0 <= a && a < _n);
        if (parent_or_size[a] < 0) return a;
        return parent_or_size[a] = leader(parent_or_size[a]);
    }

    int size(int a) {
        assert(0 <= a && a < _n);
        return -parent_or_size[leader(a)];
    }

    std::vector<std::vector<int>> groups() {
        std::vector<int> leader_buf(_n), group_size(_n);
        for (int i = 0; i < _n; i++) {
            leader_buf[i] = leader(i);
            group_size[leader_buf[i]]++;
        }
        std::vector<std::vector<int>> result(_n);
        for (int i = 0; i < _n; i++) {
            result[i].reserve(group_size[i]);
        }
        for (int i = 0; i < _n; i++) {
            result[leader_buf[i]].push_back(i);
        }
        result.erase(
            std::remove_if(result.begin(), result.end(),
                           [&](const std::vector<int>& v) { return v.empty(); }),
            result.end());
        return result;
    }

    int group_count() const {
        return _cc;
    }

  private:
    int _n, _cc;
    // root node: -1 * component size
    // otherwise: parent
    std::vector<int> parent_or_size;
};

} // namespace noya2
#line 171 "c.cpp"

vector<dom> recolor(int h, int w, vector<dom> a){
    dsu d(h*w);
    vector c(h,vector<int>(w));
    for (auto [x1, y1, c1, x2, y2, c2] : a){
        c[x1][y1] = c1;
        c[x2][y2] = c2;
    }
    grid gr(h,w);
    rep(i,h) rep(j,w){
        if (i > 0 && c[i-1][j] == c[i][j]){
            d.merge(gr.idx(i-1,j),gr.idx(i,j));
        }
        if (j > 0 && c[i][j-1] == c[i][j]){
            d.merge(gr.idx(i,j-1),gr.idx(i,j));
        }
    }
    auto grs = d.groups();
    int sz = grs.size();
    rep(col,sz){
        for (auto id : grs[col]){
            auto [x, y] = gr.pos(id);
            c[x][y] = col+1;
        }
    }
    for (auto &[x1, y1, c1, x2, y2, c2] : a){
        c1 = c[x1][y1];
        c2 = c[x2][y2];
    }
    return a;
}

void solve(){
    int h, w; in(h,w);
    auto ans = make(h,w);
    // out("ok"); cout << flush;
    ans = recolor(h,w,ans);
    int ma = 0;
    for (auto d : ans){
        chmax(ma,d.c1);
        chmax(ma,d.c2);
    }
    out(ma);
    for (auto d : ans){
        out(d);
    }
}

int main(){
    int t = 1; //in(t);
    while (t--) { solve(); }
}