#include // clang-format off using namespace std; using ll=long long; using ull=unsigned long long; using pll=pair; const ll INF=4e18; void print0(){}; template void print0(H h,T... t){cout<void print(H h,T... t){print0(h);if(sizeof...(T)>0)print0(" ");print(t...);} void perr0(){}; template void perr0(H h,T... t){cerr<void perr(H h,T... t){perr0(h);if(sizeof...(T)>0)perr0(" ");perr(t...);} void ioinit() { cout<; mt19937 engine(0); clock_t start_time; double now() { return 1000.0 * (clock() - start_time) / CLOCKS_PER_SEC; } void marathon_init() { start_time = clock(); random_device seed_gen; engine.seed(seed_gen()); } struct point { int i; int j; }; const int PLANET = 1; const int STATION = 2; const int N = 100; const int M = 8; const int ALPHA = 5; vector planets(N); int distance2(point a, point b) { return (a.i - b.i) * (a.i - b.i) + (a.j - b.j) * (a.j - b.j); } int calc_score(vector ops, vector stations) { // round(10^9/(1000+sqrt(S))) double S = 0; int m = ops.size(); for (int i = 1; i < m; i++) { pii pre = ops[i - 1]; pii cur = ops[i]; double ratio = 1; point pp; point cp; if (pre.first == STATION) { pp = stations[pre.second]; } else { pp = planets[pre.second]; ratio *= ALPHA; } if (cur.first == STATION) { cp = stations[cur.second]; } else { cp = planets[cur.second]; ratio *= ALPHA; } S += ratio * distance2(cp, pp); } return 1e9 / (1000.0 + sqrt(S)) + 0.5; } void solve() { // k-means // 初期値は固定(TODO) vector stations = { {166, 166}, {333, 500}, {166, 833}, {500, 666}, {833, 833}, {666, 500}, {833, 166}, {500, 333}, }; assert(stations.size() == M); vector pre_assign(N, -1); for (int iter = 0; iter < 10000; iter++) { vector assign(N); for (int pid = 0; pid < N; pid++) { auto pl = planets[pid]; int mindist = 1e9; int minstation = -1; for (int sid = 0; sid < M; sid++) { auto st = stations[sid]; int d = distance2(pl, st); if (mindist > d) { mindist = d; minstation = sid; } } assign[pid] = minstation; } if (pre_assign == assign) { break; } vector sums(M, {0, 0}); vector counts(M); for (int pid = 0; pid < N; pid++) { auto pl = planets[pid]; auto sid = assign[pid]; sums[sid].i += pl.i; sums[sid].j += pl.j; counts[sid]++; } for (int sid = 0; sid < M; sid++) { int cnt = counts[sid]; auto su = sums[sid]; if (cnt == 0) { // TODO 初期値の選び方によってこれを起こさないことはできるが... // stationの位置に変更なし continue; } stations[sid] = {su.i / cnt, su.j / cnt}; } pre_assign = assign; } // クラスタの移動は適当に決め打ち。 // TODO bitdpでできる // station1->planet1.1->station1->planet1.2->station1->station2 のように移動 // だが、planet1.1->planet1.2 のようにやるほうが短いならそれを選んでもよい vector ops; { int sid = pre_assign[0]; int sid_init = sid; ops.push_back({PLANET, 0}); while (true) { // 座圧? vector i2p; map p2i; for (int pid = 0; pid < N; pid++) { if (pre_assign[pid] == sid) { p2i[pid] = i2p.size(); i2p.push_back(pid); } } if (i2p.size() > 0) { // 距離に偏りのあるTSP. まずは貪欲 / TODO 焼きなましとかしたほうがよさそう int m = i2p.size(); point st = stations[sid]; vector done(m); int donenum = 0; while (true) { pii cur = ops.back(); if (cur.first == STATION) { // 適当に選ぶ for (int i = 0; i < m; i++) { if (!done[i]) { ops.push_back({PLANET, i2p[i]}); break; } } continue; } int pid = cur.second; int i = p2i[pid]; done[i] = true; donenum++; if (donenum == m) { ops.push_back({STATION, sid}); break; } int d_station = ALPHA * distance2(planets[pid], st); int mind = d_station; int minj = -1; for (int j = 0; j < m; j++) { if (done[j]) continue; int dj = ALPHA * ALPHA * distance2(planets[pid], planets[i2p[j]]); if (mind > dj) { mind = dj; minj = j; } } if (minj == -1) { ops.push_back({STATION, sid}); continue; } else { ops.push_back({PLANET, i2p[minj]}); continue; } } } sid = (sid + 1) % M; if (sid == sid_init) break; ops.push_back({STATION, sid}); } ops.push_back({STATION, sid_init}); ops.push_back({PLANET, 0}); } // 出力 for (auto st : stations) { print(st.i, st.j); } print(ops.size()); for (auto op : ops) { print(op.first, op.second + 1); } { perr("time=", now()); perr("score=", calc_score(ops, stations)); } } int main() { marathon_init(); ioinit(); int n, m; cin >> n >> m; for (int i = 0; i < n; i++) { int a, b; cin >> a >> b; planets[i] = {i : a, j : b}; } solve(); return 0; }