#define _CRT_SECURE_NO_WARNINGS #include "bits/stdc++.h" #include #include #include using namespace std; //呪文 #define DUMPOUT cerr #define dump(...) DUMPOUT<<" ";DUMPOUT<<#__VA_ARGS__<<" :["<<__LINE__<<":"<<__FUNCTION__<<"]"< pii; typedef pair pll; typedef pair pdd; typedef pair pss; template ostream& operator << (ostream& ostr, const pair<_KTy, _Ty>& m) { ostr << "{" << m.first << ", " << m.second << "}"; return ostr; } template ostream& operator << (ostream& ostr, const map<_KTy, _Ty>& m) { if (m.empty()) { ostr << "{ }"; return ostr; } ostr << "{" << *m.begin(); for (auto itr = ++m.begin(); itr != m.end(); itr++) { ostr << ", " << *itr; } ostr << "}"; return ostr; } template ostream& operator << (ostream& ostr, const vector<_Ty>& v) { if (v.empty()) { ostr << "{ }"; return ostr; } ostr << "{" << v.front(); for (auto itr = ++v.begin(); itr != v.end(); itr++) { ostr << ", " << *itr; } ostr << "}"; return ostr; } template ostream& operator << (ostream& ostr, const set<_Ty>& s) { if (s.empty()) { ostr << "{ }"; return ostr; } ostr << "{" << *(s.begin()); for (auto itr = ++s.begin(); itr != s.end(); itr++) { ostr << ", " << *itr; } ostr << "}"; return ostr; } template ostream& operator << (ostream& ostr, const stack<_Ty>& s) { if (s.empty()) { ostr << "{ }"; return ostr; } stack<_Ty> t(s); ostr << "{" << t.top(); t.pop(); while (!t.empty()) { ostr << ", " << t.top(); t.pop(); } ostr << "}"; return ostr; } template ostream& operator << (ostream& ostr, const list<_Ty>& l) { if (l.empty()) { ostr << "{ }"; return ostr; } ostr << "{" << l.front(); for (auto itr = ++l.begin(); itr != l.end(); ++itr) { ostr << ", " << *itr; } ostr << "}"; return ostr; } template istream& operator >> (istream& istr, pair<_KTy, _Ty>& m) { istr >> m.first >> m.second; return istr; } template istream& operator >> (istream& istr, vector<_Ty>& v) { for (size_t i = 0; i < v.size(); i++) istr >> v[i]; return istr; } namespace aux { // print tuple template struct tp { static void print(ostream& os, const Ty& v) { os << get(v) << ", "; tp::print(os, v); } }; template struct tp { static void print(ostream& os, const Ty& value) { os << get(value); } }; } template ostream& operator<<(ostream& os, const tuple& t) { os << "{"; aux::tp, 0, sizeof...(Tys)-1>::print(os, t); os << "}"; return os; } template inline void Fill(A(&array)[N], const T &val) { std::fill((T*)array, (T*)(array + N), val); } void dump_func() { DUMPOUT << endl; } template void dump_func(Head&& head, Tail&&... tail) { DUMPOUT << head; if (sizeof...(Tail) == 0) { DUMPOUT << " "; } else { DUMPOUT << ", "; } dump_func(std::move(tail)...); } #define PI 3.14159265358979323846 #define EPS 1e-11 #define FOR(i,a,b) for(int i=(a);i<(b);++i) #define REP(i,n) FOR(i,0,n) #define all(x) (x).begin(), (x).end() #define SZ(x) ((int)(x).size()) #define fake false struct SXor128 { unsigned x, y, z, w; SXor128() { x = 123456789; y = 362436069; z = 521288629; w = 88675123; } SXor128(int _w) { x = 123456789; y = 362436069; z = 521288629; w = _w; } void setSeed() { x = 123456789; y = 362436069; z = 521288629; w = 88675123; } void setSeed(int _w) { x = 123456789; y = 362436069; z = 521288629; w = _w; } unsigned nextUInt() { unsigned t = (x ^ (x << 11)); x = y; y = z; z = w; return (w = (w ^ (w >> 19)) ^ (t ^ (t >> 8))); } // [0, mod) unsigned nextUInt(unsigned mod) { unsigned t = (x ^ (x << 11)); x = y; y = z; z = w; w = (w ^ (w >> 19)) ^ (t ^ (t >> 8)); return w % mod; } // [l, r) unsigned nextUInt(unsigned l, unsigned r) { unsigned t = (x ^ (x << 11)); x = y; y = z; z = w; w = (w ^ (w >> 19)) ^ (t ^ (t >> 8)); return w % (r - l) + l; } }; int N, K; vector L; SXor128 rnd; struct State { vector> board; int score; vector> op; State(vector> init_board) : board(init_board), score(0) {} friend ostream& operator << (ostream& os, const State& s) { cout << "Score = " << s.score << endl; for (int i = 0; i < s.board.size(); i++) cout << s.board[i] << endl; return os; } void randomFlip(int index) { int length = L[index]; bool row = rnd.nextUInt() & 1; int sx, sy; if (row) { sy = rnd.nextUInt(N); sx = rnd.nextUInt(N - length + 1); } else { sy = rnd.nextUInt(N - length + 1); sx = rnd.nextUInt(N); } if (row) { for (int x = sx; x < sx + length; x++) { score += (board[sy][x] ? 1 : -1); board[sy][x] = !board[sy][x]; } } else { for (int y = sy; y < sy + length; y++) { score += (board[y][sx] ? 1 : -1); board[y][sx] = !board[y][sx]; } } if (row) op.push_back({ sy, sx, sy, sx + length - 1 }); else op.push_back({ sy, sx, sy + length - 1, sx }); } }; int main() { cin.tie(0); ios::sync_with_stdio(false); cin >> N >> K; L.resize(K); cin >> L; vector input(N); cin >> input; vector> init_board(N, vector(N)); // あとで long long で for (int i = 0; i < N; i++) for (int j = 0; j < N; j++) init_board[i][j] = input[i][j] - '0'; State init_state(init_board); State bestState = init_state; while ((double)clock() / CLOCKS_PER_SEC < 0.99) { State state = init_state; for (int i = 0; i < K; i++) { state.randomFlip(i); } if (state.score > bestState.score) { bestState = state; } } for (int i = 0; i < bestState.op.size(); i++) { cout << bestState.op[i][0] << " " << bestState.op[i][1] << " " << bestState.op[i][2] << " " << bestState.op[i][3] << endl; } return 0; }