//* #pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") //*/ #include // #include using namespace std; // using namespace atcoder; // #define _GLIBCXX_DEBUG #define DEBUG(x) cerr << #x << ": " << x << endl; #define DEBUG_VEC(v) \ cerr << #v << ":"; \ for (int iiiiiiii = 0; iiiiiiii < v.size(); iiiiiiii++) \ cerr << " " << v[iiiiiiii]; \ cerr << endl; #define DEBUG_MAT(v) \ cerr << #v << endl; \ for (int iv = 0; iv < v.size(); iv++) { \ for (int jv = 0; jv < v[iv].size(); jv++) { \ cerr << v[iv][jv] << " "; \ } \ cerr << endl; \ } typedef long long ll; // #define int ll #define vi vector #define vl vector #define vii vector> #define vll vector> #define vs vector #define pii pair #define pis pair #define psi pair #define pll pair template pair operator+(const pair &s, const pair &t) { return pair(s.first + t.first, s.second + t.second); } template pair operator-(const pair &s, const pair &t) { return pair(s.first - t.first, s.second - t.second); } template ostream &operator<<(ostream &os, pair p) { os << "(" << p.first << ", " << p.second << ")"; return os; } #define X first #define Y second #define rep(i, n) for (int i = 0; i < (int)(n); i++) #define rep1(i, n) for (int i = 1; i <= (int)(n); i++) #define rrep(i, n) for (int i = (int)(n)-1; i >= 0; i--) #define rrep1(i, n) for (int i = (int)(n); i > 0; i--) #define REP(i, a, b) for (int i = a; i < b; i++) #define in(x, a, b) (a <= x && x < b) #define all(c) c.begin(), c.end() void YES(bool t = true) { cout << (t ? "YES" : "NO") << endl; } void Yes(bool t = true) { cout << (t ? "Yes" : "No") << endl; } void yes(bool t = true) { cout << (t ? "yes" : "no") << endl; } void NO(bool t = true) { cout << (t ? "NO" : "YES") << endl; } void No(bool t = true) { cout << (t ? "No" : "Yes") << endl; } void no(bool t = true) { cout << (t ? "no" : "yes") << endl; } template bool chmax(T &a, const T &b) { if (a < b) { a = b; return 1; } return 0; } template bool chmin(T &a, const T &b) { if (a > b) { a = b; return 1; } return 0; } #define UNIQUE(v) v.erase(std::unique(v.begin(), v.end()), v.end()); const ll inf = 1000000001; const ll INF = (ll)1e18 + 1; const long double pi = 3.1415926535897932384626433832795028841971L; int popcount(ll t) { return __builtin_popcountll(t); } // int dx[4] = {1, 0, -1, 0}, dy[4] = {0, 1, 0, -1}; // int dx2[8] = { 1,1,0,-1,-1,-1,0,1 }, dy2[8] = { 0,1,1,1,0,-1,-1,-1 }; vi dx = {0, 0, -1, 1}, dy = {-1, 1, 0, 0}; vi dx2 = {1, 1, 0, -1, -1, -1, 0, 1}, dy2 = {0, 1, 1, 1, 0, -1, -1, -1}; struct Setup_io { Setup_io() { ios_base::sync_with_stdio(0), cin.tie(0), cout.tie(0); cout << fixed << setprecision(25); } } setup_io; const ll MOD = 1000000007; // const ll MOD = 998244353; // #define mp make_pair //#define endl '\n' struct dice { mt19937 mt; dice() : mt(chrono::steady_clock::now().time_since_epoch().count()) {} // [0, x)の一様乱数 ll operator()(ll x) { return this->operator()(0, x); } // [x, y)の一様乱数 ll operator()(ll x, ll y) { uniform_int_distribution dist(x, y - 1); return dist(mt); } vl operator()(int n, ll x, ll y) { vl res(n); for (int i = 0; i < n; i++) res[i] = this->operator()(x, y); return res; } } rnd; constexpr int n = 100, m = 8; constexpr int alpha = 5; vector ab(n); vector cd(m); void input() { int nn, mm; cin >> nn >> mm; rep(i, n) { cin >> ab[i].first >> ab[i].second; } rep(i, m) { int y = i / 3; int x = i % 3; cd[i] = pii(250 * (y + 1), 250 * (x + 1)); } } int calc_dist(pii xy1, pii xy2) { pll dxy = xy1 - xy2; return dxy.first * dxy.first + dxy.second * dxy.second; } int ori_hyou[n][n]; int hyou[n][n]; pii ori_inter[n][n]; pii inter[n][n]; void fill_hyou(vector &jun) { rep(i, n) { rep(j, n) { ori_hyou[i][j] = hyou[i][j]; ori_inter[i][j] = inter[i][j]; } } rep(i, n) { REP(j, i, n) { hyou[i][j] = calc_dist(ab[i], ab[j]) * alpha * alpha; inter[i][j] = pii(-1, -1); rep(k, m) { if (chmin(hyou[i][j], (calc_dist(ab[i], cd[k]) + calc_dist(cd[k], ab[j])) * alpha)) { inter[i][j] = pii(k, -1); } } rep(k1, m) { int temp = calc_dist(ab[i], cd[k1]) * alpha; if (temp >= hyou[i][j]) { continue; } rep(k2, m) { int add = calc_dist(cd[k1], cd[k2]) + calc_dist(cd[k2], ab[j]) * alpha; if (chmin(hyou[i][j], temp + add)) { inter[i][j] = pii(k1, k2); } } } hyou[j][i] = hyou[i][j]; inter[j][i] = inter[i][j]; if (inter[j][i].second != -1) { swap(inter[j][i].first, inter[j][i].second); } } } } int calc_score(const vi &dist) { ll sum = 0; rep(i, n) sum += dist[i]; return round(1e9 / (1000 + sqrt(sum))); } int two_opt(vector &jun) { fill_hyou(jun); static vector dist(n); rep(i, n) { int s = jun[i], t = jun[(i + 1) % n]; dist[i] = hyou[s][t]; } bool updated = true; while (updated) { updated = false; for (int i = 1; i < n; i++) { for (int j = i + 1; j < n; j++) { int sub = 0, add = 0; sub += hyou[jun[i - 1]][jun[i]]; add += hyou[jun[i - 1]][jun[j]]; sub += hyou[jun[j]][jun[(j + 1) % n]]; add += hyou[jun[i]][jun[(j + 1) % n]]; if (sub > add) { reverse(jun.begin() + i, jun.begin() + j + 1); reverse(dist.begin() + i, dist.begin() + j); dist[i - 1] = hyou[jun[i - 1]][jun[i]]; dist[j] = hyou[jun[j]][jun[(j + 1) % n]]; updated = true; } } } } return calc_score(dist); } int double_bridge(vector &jun) { static vi erase_idx; while (true) { erase_idx.clear(); rep(i, 4) { erase_idx.push_back(rnd(n)); } sort(all(erase_idx)); UNIQUE(erase_idx); if (erase_idx.size() == 4) { break; } } static vi new_jun; new_jun.clear(); rep(i, erase_idx[0] + 1) { new_jun.emplace_back(i); } REP(i, erase_idx[2] + 1, erase_idx[3] + 1) { new_jun.emplace_back(i); } REP(i, erase_idx[1] + 1, erase_idx[2] + 1) { new_jun.emplace_back(i); } REP(i, erase_idx[0] + 1, erase_idx[1] + 1) { new_jun.emplace_back(i); } REP(i, erase_idx[3] + 1, n) { new_jun.emplace_back(i); } swap(jun, new_jun); return two_opt(jun); } void output(vector &jun) { rep(i, m) { cout << cd[i].first << " " << cd[i].second << endl; } vector ans; ans.emplace_back(1, 0); rep(i, n) { int s = jun[i], t = jun[(i + 1) % n]; pii keiyu = inter[s][t]; if (keiyu.first != -1) { ans.emplace_back(2, keiyu.first); } if (keiyu.second != -1) { ans.emplace_back(2, keiyu.second); } ans.emplace_back(1, t); } cout << ans.size() << endl; for (auto [t, u] : ans) { cout << t << " " << u + 1 << endl; } } signed main() { clock_t start_c = clock(); input(); vector jun(n); rep(i, n) { jun[i] = i; } int score = two_opt(jun); int best_score = score; vi best_jun = jun; vector best_cd = cd; rep(try_num, inf) { if (try_num % 256 == 0) { double sec = (double)(clock() - start_c) / CLOCKS_PER_SEC; if (sec > 0.9) { DEBUG(sec); DEBUG(try_num); break; } } if (try_num % 64 == 0) { DEBUG(score); score = double_bridge(jun); DEBUG(score); continue; } int idx = rnd(m); pii ori_cd = cd[idx]; cd[idx].first += rnd(-30, 30); cd[idx].second += rnd(-30, 30); cd[idx].first = min(max(cd[idx].first, 0), 1000); cd[idx].second = min(max(cd[idx].second, 0), 1000); static vi ori_jun; ori_jun = jun; int new_score = two_opt(jun); if (new_score < score) { cd[idx] = ori_cd; swap(jun, ori_jun); rep(i, n) { rep(j, n) { hyou[i][j] = ori_hyou[i][j]; inter[i][j] = ori_inter[i][j]; } } } else { score = new_score; if (score > best_score) { best_score = score; best_jun = jun; best_cd = cd; } } } cd = best_cd; jun = best_jun; DEBUG(score); DEBUG(best_score); DEBUG(two_opt(jun)); output(jun); }