Submission #42968021


Source Code Expand

/**
 *    author:  otera
**/
#include<bits/stdc++.h>

namespace otera {
    class random_gen {
        public:
            random_gen() : mt64(std::chrono::steady_clock::now().time_since_epoch().count()) {}

            // [a, b)
            long long operator()(long long a, long long b) {
                std::uniform_int_distribution<long long> dist(a, b - 1);
                return dist(mt64);
            }
        private:
            std::mt19937_64 mt64;
    };
} // namespace otera

namespace otera {
    struct timer {
        public:
        timer() {
            start();
        }
        void start() {
            time = std::clock();
        }   
        double elapsed() {
            return 1.0 * (std::clock() - time) / CLOCKS_PER_SEC;
        }
        private:
        int time;
    };
} // namespace otera

namespace otera {}
using namespace std;
using namespace otera;

struct io_setup {
    io_setup(int precision = 20) {
        ios::sync_with_stdio(false);
        cin.tie(0);
        cout << fixed << setprecision(precision);
    }
} io_setup_ {};

#define int long long

using ll = long long;
using ld = long double;
using ull = unsigned long long;
using i128 = __int128_t;
using u128 = __uint128_t;
#define repa(i, n) for(int i = 0; i < (n); ++ i)
#define repb(i, a, b) for(int i = (a); i < (b); ++ i)
#define repc(i, a, b, c) for(int i = (a); i < (b); i += (c))
#define overload4(a, b, c, d, e, ...) e
#define overload3(a, b, c, d, ...) d
#define rep(...) overload4(__VA_ARGS__, repc, repb, repa)(__VA_ARGS__)
#define rep1a(i, n) for(int i = 0; i <= (n); ++ i)
#define rep1b(i, a, b) for(int i = (a); i <= (b); ++ i)
#define rep1c(i, a, b, c) for(int i = (a); i <= (b); i += (c))
#define rep1(...) overload4(__VA_ARGS__, rep1c, rep1b, rep1a)(__VA_ARGS__)
#define rev_repa(i, n) for(int i=(n)-1;i>=0;i--)
#define rev_repb(i, a, b) assert((a) > (b));for(int i=(a);i>(b);i--)
#define rev_rep(...) overload3(__VA_ARGS__, rev_repb, rev_repa)(__VA_ARGS__)
#define rev_rep1a(i, n) for(int i=(n);i>=1;i--)
#define rev_rep1b(i, a, b) assert((a) >= (b));for(int i=(a);i>=(b);i--)
#define rev_rep1(...) overload3(__VA_ARGS__, rev_rep1b, rev_rep1a)(__VA_ARGS__)
#define for_subset(t, s) for(ll t = (s); t >= 0; t = (t == 0 ? -1 : (t - 1) & (s)))
typedef pair<int, int> P;
typedef pair<ll, ll> LP;
#define pb push_back
#define pf push_front
#define ppb pop_back
#define ppf pop_front
#define eb emplace_back
#define mkp make_pair
#define mkt make_tuple
#define fr first
#define sc second
#define all(c) std::begin(c), std::end(c)
#define rall(c) std::rbegin(c), std::rend(c)
#define lb(c, x) distance(std::begin(c), lower_bound(all(c), (x)))
#define ub(c, x) distance(std::begin(c), upper_bound(all(c), (x)))
#define Sort(a) sort(all(a))
#define Rev(a) reverse(all(a))
#define Uniq(a) sort(all(a));(a).erase(unique(all(a)),std::end(a))
#define si(c) (int)(c).size()
inline ll popcnt(ull a){ return __builtin_popcountll(a); }
#define kth_bit(x, k) ((x>>k)&1)
#define unless(A) if(!(A))
#define elif else if
ll modulo(ll n, ll d){ return (n % d + d) % d; };
ll intpow(ll a, ll b){ ll ans = 1; while(b){ if(b & 1) ans *= a; a *= a; b /= 2; } return ans; }
ll intpow(ll a, ll b, ll m) {ll ans = 1; while(b){ if(b & 1) (ans *= a) %= m; (a *= a) %= m; b /= 2; } return ans; }
template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; }
template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; }
#define INT(...) int __VA_ARGS__;in(__VA_ARGS__)
#define LL(...) ll __VA_ARGS__;in(__VA_ARGS__)
#define ULL(...) ull __VA_ARGS__;in(__VA_ARGS__)
#define STR(...) string __VA_ARGS__;in(__VA_ARGS__)
#define CHR(...) char __VA_ARGS__;in(__VA_ARGS__)
#define DBL(...) double __VA_ARGS__;in(__VA_ARGS__)
#define LD(...) ld __VA_ARGS__;in(__VA_ARGS__)
#define vec(type,name,...) vector<type>name(__VA_ARGS__)
#define VEC(type,name,size) vector<type>name(size);in(name)
#define vv(type,name,h,...) vector<vector<type>>name(h,vector<type>(__VA_ARGS__))
#define VV(type,name,h,w) vector<vector<type>>name(h,vector<type>(w));in(name)
#define vvv(type,name,h,w,...) vector<vector<vector<type>>>name(h,vector<vector<type>>(w,vector<type>(__VA_ARGS__)))
template <class T> using vc = vector<T>;
template <class T> using vvc = vector<vc<T>>;
template <class T> using vvvc = vector<vvc<T>>;
template <class T> using vvvvc = vector<vvvc<T>>;
template <class T> using pq = priority_queue<T>;
template <class T> using pqg = priority_queue<T, vector<T>, greater<T>>;
template <class T, class U> using umap = unordered_map<T, U>;
template<class T> void scan(T& a){ cin >> a; }
template<class T> void scan(vector<T>& a){ for(auto&& i : a) scan(i); }
void in(){}
template <class Head, class... Tail> void in(Head& head, Tail&... tail){ scan(head); in(tail...); }
void print(){ cout << ' '; }
template<class T> void print(const T& a){ cout << a; }
template<class T> void print(const vector<T>& a){ if(a.empty()) return; print(a[0]); for(auto i = a.begin(); ++i != a.end(); ){ cout << ' '; print(*i); } }
int out(){ cout << '\n'; return 0; }
template<class T> int out(const T& t){ print(t); cout << '\n'; return 0; }
template<class Head, class... Tail> int out(const Head& head, const Tail&... tail){ print(head); cout << ' '; out(tail...); return 0; }
#define CHOOSE(a) CHOOSE2 a
#define CHOOSE2(a0,a1,a2,a3,a4,x,...) x
#define debug_1(x1) cout<<#x1<<": "<<x1<<endl
#define debug_2(x1,x2) cout<<#x1<<": "<<x1<<", "#x2<<": "<<x2<<endl
#define debug_3(x1,x2,x3) cout<<#x1<<": "<<x1<<", "#x2<<": "<<x2<<", "#x3<<": "<<x3<<endl
#define debug_4(x1,x2,x3,x4) cout<<#x1<<": "<<x1<<", "#x2<<": "<<x2<<", "#x3<<": "<<x3<<", "#x4<<": "<<x4<<endl
#define debug_5(x1,x2,x3,x4,x5) cout<<#x1<<": "<<x1<<", "#x2<<": "<<x2<<", "#x3<<": "<<x3<<", "#x4<<": "<<x4<<", "#x5<<": "<<x5<<endl
#ifdef DEBUG
#define debug(...) CHOOSE((__VA_ARGS__,debug_5,debug_4,debug_3,debug_2,debug_1,~))(__VA_ARGS__)
#define dump(...) { print(#__VA_ARGS__); print(":"); out(__VA_ARGS__); }
#else
#define debug(...)
#define dump(...)
#endif

random_gen rng;
timer ti;

const int n = 30;
const int dx[6] = {-1, -1, 0, 0, 1, 1};
const int dy[6] = {-1, 0, -1, 1, 0, 1};

const int down_dx[2] = {1, 1};
const int down_dy[2] = {0, 1};

const int up_dx[2] = {-1, -1};
const int up_dy[2] = {0, -1};

const int K = 10000;

const int M = 465;

struct S {
    int px, py;
    int nx, ny;
};

bool inside(int x, int y) {
    if(x < 0 or x >= n) return false;
    if(y < 0 or y >= x + 1) return false;
    return true;
}

void solve() {
    vvc<int> b(n);
    vc<S> ans;
    vc<P> pos(M);
    // input
    {
        rep(i, n) {
            b[i].resize(i + 1);
            rep(j, i + 1) {
                in(b[i][j]);
                pos[b[i][j]] = mkp(i, j);
            }
        }
    }
    auto greedy2 = [&]() -> void {
        auto simulate = [&](const int val) -> int {
            auto b_ = b;
            auto pos_ = pos;
            int ret = 0;
            for(int i = val; i < M; ++ i) {
                while(1) {
                    auto [x, y] = pos_[i];
                    bool found = 0;
                    int nxt_x = -1, nxt_y = -1;
                    if(inside(x - 1, y) and inside(x - 1, y - 1) and b_[x - 1][y] > b_[x][y] and b_[x - 1][y - 1] > b_[x][y]) {
                        if(b_[x - 1][y] > b_[x - 1][y - 1]) {
                            nxt_x = x - 1, nxt_y = y;
                        } else {
                            nxt_x = x - 1, nxt_y = y - 1;
                        }
                        ++ ret;
                        swap(pos_[b_[x][y]], pos_[b_[nxt_x][nxt_y]]);
                        swap(b_[nxt_x][nxt_y], b_[x][y]);
                        found = 1;
                    }
                    rep(k, 2) {
                        if(found) break;
                        int nx = x + up_dx[k], ny = y + up_dy[k];
                        if(inside(nx, ny) and b_[nx][ny] > b_[x][y]) {
                            ++ ret;
                            swap(pos_[b_[x][y]], pos_[b_[nx][ny]]);
                            swap(b_[nx][ny], b_[x][y]);
                            found = 1;
                        }
                    }
                    if(!found) break;
                }
            }
            return ret;
        };
        rep(i, M) {
            if(si(ans) >= K) break;
            while(1) {
                if(si(ans) >= K) break;
                auto [x, y] = pos[i];
                bool found = 0;
                int nxt_x = -1, nxt_y = -1;
                if(inside(x - 1, y) and inside(x - 1, y - 1) and b[x - 1][y] > b[x][y] and b[x - 1][y - 1] > b[x][y]) {
                    swap(pos[b[x][y]], pos[b[x - 1][y]]);
                    swap(b[x - 1][y], b[x][y]);
                    int cost = simulate(i);
                    swap(pos[b[x][y]], pos[b[x - 1][y]]);
                    swap(b[x - 1][y], b[x][y]);
                    // 
                    swap(pos[b[x][y]], pos[b[x - 1][y - 1]]);
                    swap(b[x - 1][y - 1], b[x][y]);
                    int ncost = simulate(i);
                    debug(cost, ncost);
                    if(cost < ncost) {
                        swap(pos[b[x][y]], pos[b[x - 1][y - 1]]);
                        swap(b[x - 1][y - 1], b[x][y]);
                        // 
                        swap(pos[b[x][y]], pos[b[x - 1][y]]);
                        swap(b[x - 1][y], b[x][y]);
                        ans.pb(S{x - 1, y, x, y});
                    } else {
                        ans.pb(S{x - 1, y - 1, x, y});
                    }
                    found = 1;
                }
                rep(k, 2) {
                    if(found) break;
                    int nx = x + up_dx[k], ny = y + up_dy[k];
                    if(inside(nx, ny) and b[nx][ny] > b[x][y]) {
                        ans.eb(S{nx, ny, x, y});
                        swap(pos[b[x][y]], pos[b[nx][ny]]);
                        swap(b[nx][ny], b[x][y]);
                        found = 1;
                    }
                }
                if(!found) break;
            }
        }
    };
    // solution
    {
        greedy2();
    }
    // output
    {
        int k = si(ans);
        out(k);
        for(auto [px, py, nx, ny]: ans) {
            out(px, py, nx, ny);
        }
    }
}

signed main() {
    int testcase = 1;
    // in(testcase);
    while(testcase--) solve();
    return 0;
}

Submission Info

Submission Time
Task A - Pyramid Sorting
User otera
Language C++ (GCC 9.2.1)
Score 13516495
Code Size 10718 Byte
Status AC
Exec Time 55 ms
Memory 3840 KiB

Compile Error

./Main.cpp: In lambda function:
./Main.cpp:230:21: warning: unused variable ‘nxt_x’ [-Wunused-variable]
  230 |                 int nxt_x = -1, nxt_y = -1;
      |                     ^~~~~
./Main.cpp:230:33: warning: unused variable ‘nxt_y’ [-Wunused-variable]
  230 |                 int nxt_x = -1, nxt_y = -1;
      |                                 ^~~~~

Judge Result

Set Name test_ALL
Score / Max Score 13516495 / 15000000
Status
AC × 150
Set Name Test Cases
test_ALL test_0000.txt, test_0001.txt, test_0002.txt, test_0003.txt, test_0004.txt, test_0005.txt, test_0006.txt, test_0007.txt, test_0008.txt, test_0009.txt, test_0010.txt, test_0011.txt, test_0012.txt, test_0013.txt, test_0014.txt, test_0015.txt, test_0016.txt, test_0017.txt, test_0018.txt, test_0019.txt, test_0020.txt, test_0021.txt, test_0022.txt, test_0023.txt, test_0024.txt, test_0025.txt, test_0026.txt, test_0027.txt, test_0028.txt, test_0029.txt, test_0030.txt, test_0031.txt, test_0032.txt, test_0033.txt, test_0034.txt, test_0035.txt, test_0036.txt, test_0037.txt, test_0038.txt, test_0039.txt, test_0040.txt, test_0041.txt, test_0042.txt, test_0043.txt, test_0044.txt, test_0045.txt, test_0046.txt, test_0047.txt, test_0048.txt, test_0049.txt, test_0050.txt, test_0051.txt, test_0052.txt, test_0053.txt, test_0054.txt, test_0055.txt, test_0056.txt, test_0057.txt, test_0058.txt, test_0059.txt, test_0060.txt, test_0061.txt, test_0062.txt, test_0063.txt, test_0064.txt, test_0065.txt, test_0066.txt, test_0067.txt, test_0068.txt, test_0069.txt, test_0070.txt, test_0071.txt, test_0072.txt, test_0073.txt, test_0074.txt, test_0075.txt, test_0076.txt, test_0077.txt, test_0078.txt, test_0079.txt, test_0080.txt, test_0081.txt, test_0082.txt, test_0083.txt, test_0084.txt, test_0085.txt, test_0086.txt, test_0087.txt, test_0088.txt, test_0089.txt, test_0090.txt, test_0091.txt, test_0092.txt, test_0093.txt, test_0094.txt, test_0095.txt, test_0096.txt, test_0097.txt, test_0098.txt, test_0099.txt, test_0100.txt, test_0101.txt, test_0102.txt, test_0103.txt, test_0104.txt, test_0105.txt, test_0106.txt, test_0107.txt, test_0108.txt, test_0109.txt, test_0110.txt, test_0111.txt, test_0112.txt, test_0113.txt, test_0114.txt, test_0115.txt, test_0116.txt, test_0117.txt, test_0118.txt, test_0119.txt, test_0120.txt, test_0121.txt, test_0122.txt, test_0123.txt, test_0124.txt, test_0125.txt, test_0126.txt, test_0127.txt, test_0128.txt, test_0129.txt, test_0130.txt, test_0131.txt, test_0132.txt, test_0133.txt, test_0134.txt, test_0135.txt, test_0136.txt, test_0137.txt, test_0138.txt, test_0139.txt, test_0140.txt, test_0141.txt, test_0142.txt, test_0143.txt, test_0144.txt, test_0145.txt, test_0146.txt, test_0147.txt, test_0148.txt, test_0149.txt
Case Name Status Exec Time Memory
test_0000.txt AC 54 ms 3776 KiB
test_0001.txt AC 45 ms 3624 KiB
test_0002.txt AC 47 ms 3660 KiB
test_0003.txt AC 45 ms 3732 KiB
test_0004.txt AC 44 ms 3676 KiB
test_0005.txt AC 50 ms 3776 KiB
test_0006.txt AC 46 ms 3680 KiB
test_0007.txt AC 45 ms 3660 KiB
test_0008.txt AC 46 ms 3704 KiB
test_0009.txt AC 45 ms 3812 KiB
test_0010.txt AC 45 ms 3732 KiB
test_0011.txt AC 46 ms 3628 KiB
test_0012.txt AC 51 ms 3664 KiB
test_0013.txt AC 47 ms 3660 KiB
test_0014.txt AC 46 ms 3812 KiB
test_0015.txt AC 43 ms 3776 KiB
test_0016.txt AC 48 ms 3764 KiB
test_0017.txt AC 48 ms 3760 KiB
test_0018.txt AC 49 ms 3776 KiB
test_0019.txt AC 48 ms 3732 KiB
test_0020.txt AC 41 ms 3704 KiB
test_0021.txt AC 46 ms 3840 KiB
test_0022.txt AC 48 ms 3780 KiB
test_0023.txt AC 52 ms 3676 KiB
test_0024.txt AC 46 ms 3780 KiB
test_0025.txt AC 43 ms 3624 KiB
test_0026.txt AC 50 ms 3836 KiB
test_0027.txt AC 46 ms 3840 KiB
test_0028.txt AC 50 ms 3732 KiB
test_0029.txt AC 47 ms 3648 KiB
test_0030.txt AC 51 ms 3632 KiB
test_0031.txt AC 51 ms 3840 KiB
test_0032.txt AC 49 ms 3764 KiB
test_0033.txt AC 44 ms 3656 KiB
test_0034.txt AC 48 ms 3732 KiB
test_0035.txt AC 48 ms 3732 KiB
test_0036.txt AC 48 ms 3680 KiB
test_0037.txt AC 40 ms 3760 KiB
test_0038.txt AC 52 ms 3732 KiB
test_0039.txt AC 41 ms 3760 KiB
test_0040.txt AC 46 ms 3836 KiB
test_0041.txt AC 54 ms 3764 KiB
test_0042.txt AC 51 ms 3840 KiB
test_0043.txt AC 46 ms 3756 KiB
test_0044.txt AC 48 ms 3704 KiB
test_0045.txt AC 46 ms 3780 KiB
test_0046.txt AC 43 ms 3812 KiB
test_0047.txt AC 50 ms 3760 KiB
test_0048.txt AC 44 ms 3624 KiB
test_0049.txt AC 45 ms 3708 KiB
test_0050.txt AC 47 ms 3732 KiB
test_0051.txt AC 46 ms 3724 KiB
test_0052.txt AC 45 ms 3780 KiB
test_0053.txt AC 42 ms 3724 KiB
test_0054.txt AC 47 ms 3780 KiB
test_0055.txt AC 49 ms 3676 KiB
test_0056.txt AC 45 ms 3732 KiB
test_0057.txt AC 50 ms 3632 KiB
test_0058.txt AC 49 ms 3660 KiB
test_0059.txt AC 48 ms 3772 KiB
test_0060.txt AC 47 ms 3664 KiB
test_0061.txt AC 50 ms 3808 KiB
test_0062.txt AC 45 ms 3764 KiB
test_0063.txt AC 44 ms 3652 KiB
test_0064.txt AC 52 ms 3680 KiB
test_0065.txt AC 46 ms 3668 KiB
test_0066.txt AC 48 ms 3676 KiB
test_0067.txt AC 47 ms 3676 KiB
test_0068.txt AC 54 ms 3776 KiB
test_0069.txt AC 47 ms 3780 KiB
test_0070.txt AC 49 ms 3676 KiB
test_0071.txt AC 45 ms 3836 KiB
test_0072.txt AC 45 ms 3652 KiB
test_0073.txt AC 47 ms 3728 KiB
test_0074.txt AC 44 ms 3736 KiB
test_0075.txt AC 43 ms 3676 KiB
test_0076.txt AC 49 ms 3760 KiB
test_0077.txt AC 49 ms 3836 KiB
test_0078.txt AC 47 ms 3840 KiB
test_0079.txt AC 44 ms 3736 KiB
test_0080.txt AC 50 ms 3768 KiB
test_0081.txt AC 43 ms 3676 KiB
test_0082.txt AC 43 ms 3676 KiB
test_0083.txt AC 50 ms 3708 KiB
test_0084.txt AC 46 ms 3732 KiB
test_0085.txt AC 44 ms 3660 KiB
test_0086.txt AC 47 ms 3664 KiB
test_0087.txt AC 48 ms 3708 KiB
test_0088.txt AC 50 ms 3816 KiB
test_0089.txt AC 50 ms 3652 KiB
test_0090.txt AC 46 ms 3764 KiB
test_0091.txt AC 46 ms 3836 KiB
test_0092.txt AC 46 ms 3708 KiB
test_0093.txt AC 51 ms 3676 KiB
test_0094.txt AC 48 ms 3760 KiB
test_0095.txt AC 53 ms 3624 KiB
test_0096.txt AC 50 ms 3812 KiB
test_0097.txt AC 43 ms 3728 KiB
test_0098.txt AC 45 ms 3680 KiB
test_0099.txt AC 47 ms 3772 KiB
test_0100.txt AC 42 ms 3728 KiB
test_0101.txt AC 49 ms 3656 KiB
test_0102.txt AC 49 ms 3732 KiB
test_0103.txt AC 51 ms 3664 KiB
test_0104.txt AC 45 ms 3736 KiB
test_0105.txt AC 52 ms 3840 KiB
test_0106.txt AC 48 ms 3760 KiB
test_0107.txt AC 49 ms 3732 KiB
test_0108.txt AC 50 ms 3764 KiB
test_0109.txt AC 47 ms 3676 KiB
test_0110.txt AC 44 ms 3664 KiB
test_0111.txt AC 49 ms 3656 KiB
test_0112.txt AC 49 ms 3816 KiB
test_0113.txt AC 44 ms 3736 KiB
test_0114.txt AC 53 ms 3680 KiB
test_0115.txt AC 39 ms 3840 KiB
test_0116.txt AC 48 ms 3776 KiB
test_0117.txt AC 46 ms 3772 KiB
test_0118.txt AC 45 ms 3680 KiB
test_0119.txt AC 46 ms 3628 KiB
test_0120.txt AC 55 ms 3816 KiB
test_0121.txt AC 48 ms 3676 KiB
test_0122.txt AC 48 ms 3780 KiB
test_0123.txt AC 48 ms 3732 KiB
test_0124.txt AC 45 ms 3732 KiB
test_0125.txt AC 54 ms 3772 KiB
test_0126.txt AC 47 ms 3756 KiB
test_0127.txt AC 54 ms 3768 KiB
test_0128.txt AC 47 ms 3840 KiB
test_0129.txt AC 52 ms 3764 KiB
test_0130.txt AC 49 ms 3804 KiB
test_0131.txt AC 45 ms 3732 KiB
test_0132.txt AC 51 ms 3728 KiB
test_0133.txt AC 47 ms 3736 KiB
test_0134.txt AC 48 ms 3676 KiB
test_0135.txt AC 48 ms 3656 KiB
test_0136.txt AC 46 ms 3816 KiB
test_0137.txt AC 48 ms 3736 KiB
test_0138.txt AC 53 ms 3624 KiB
test_0139.txt AC 46 ms 3816 KiB
test_0140.txt AC 47 ms 3672 KiB
test_0141.txt AC 51 ms 3676 KiB
test_0142.txt AC 45 ms 3772 KiB
test_0143.txt AC 48 ms 3812 KiB
test_0144.txt AC 46 ms 3680 KiB
test_0145.txt AC 46 ms 3732 KiB
test_0146.txt AC 46 ms 3840 KiB
test_0147.txt AC 49 ms 3664 KiB
test_0148.txt AC 50 ms 3732 KiB
test_0149.txt AC 51 ms 3676 KiB