#if !__INCLUDE_LEVEL__ #include __FILE__ ll Euclid2d(ll x1,ll x2,ll y1,ll y2){ return (x1-x2)*(x1-x2)+(y1-y2)*(y1-y2); }; int main() { int N,M;cin >> N >> M; vector Planet(N,vector(2)); rep(i,N){ cin >> Planet[i][0] >> Planet[i][1]; } vector Cluster(N); random_device seed_gen; mt19937_64 rnd(seed_gen()); uniform_int_distribution Clu_N(0, M-1); vector Center(M,vector(2)); rep(i,N){ Cluster[i]=Clu_N(rnd); } rep(k,1000){ vector ClusterNum(M,0); vector tmp_center(M,vector(2)); rep(i,N){ tmp_center[Cluster[i]][0] += Planet[i][0]; tmp_center[Cluster[i]][1] += Planet[i][1]; ClusterNum[Cluster[i]]+=1; } rep(m,M){ if(ClusterNum[m]==0)continue; tmp_center[m][0]/=ClusterNum[m]; tmp_center[m][1]/=ClusterNum[m]; } bool change = false; rep(i,N){ int c = Cluster[i]; ll D = Euclid2d(Planet[i][0],tmp_center[c][0],Planet[i][1],tmp_center[c][1]); rep(m,M){ ll dist = Euclid2d(Planet[i][0],tmp_center[m][0],Planet[i][1],tmp_center[m][1]); if(dist> Out; int start = 0; vector Cplanet(8,vector(0)); rep(i,N){ Cplanet[Cluster[i]].emplace_back(i); if(i==0)start = Cluster[i]; } Out.emplace_back(1,1); Out.emplace_back(2,start+1); rep(m,M){ int station = (m+start)%M; if(station==start){ range(j,1,Cplanet[station].size()){ Out.emplace_back(1,Cplanet[station][j]+1); Out.emplace_back(2,station+1); } } else{ rep(j,Cplanet[station].size()){ Out.emplace_back(1,Cplanet[station][j]+1); Out.emplace_back(2,station+1); } } Out.emplace_back(2,(station+1)%M+1); } Out.emplace_back(1,1); ll V = Out.size(); cout << V << "\n"; rep(v,V){ cout << Out[v].first <<" " << Out[v].second << "\n"; } } #else //tie(a, b, c) = t //repをllにしているので時間ギリギリならintに変換 #include #include using namespace std; using namespace atcoder; #define rep(i, n) for(long long int i = 0; i < n; i++) #define rrep(i, n) for(long long int i = n-1; i >= 0; i--) #define range(i, m, n) for(long long int i = m; i < n; i++) #define fore(i,a) for(auto &i:a) #define all(v) v.begin(), v.end() #define rall(v) v.rbegin(), v.rend() #define Sum(v) accumulate(all(v),0LL) #define minv(v) *min_element(all(v)) #define maxv(v) *max_element(all(v)) typedef long long ll; typedef vector vl; typedef vector> vvl; const ll INF = 1e16; const ll MOD1 = 1000000007; const ll MOD2 = 998244353; template inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } ll SN(char s){return ll(s-'0');} ll SN(string s){return stoll(s);} int alpN(char s){return int(s-'a');} int AlpN(char s){return int(s-'A');} int Nalp(int n){return char(n+97);} int NAlp(int n){return char(n+65);} using Graph = vector>; using GraphCost = vector>>; using mint = modint; using mint1 = modint1000000007; using mint2 = modint998244353; using pll = pair; template ostream &operator<<(ostream &o,const vector&v){for(int i=0;i<(int)v.size();i++)o<<(i>0?" ":"")< bool contain(const std::string& s, const T& v) { return s.find(v) != std::string::npos; } ll max(int x,ll y){return max((ll)x,y);} ll max(ll x,int y){return max(x,(ll)y);} ll min(int x,ll y){return min((ll)x,y);} ll min(ll x,int y){return min(x,(ll)y);} template struct edge { int src, to; T cost; edge(int to, T cost) : src(-1), to(to), cost(cost) {} edge(int src, int to, T cost) : src(src), to(to), cost(cost) {} edge& operator=(const int& x) { to = x; return *this; } operator int() const { return to; } }; template using Edges = vector >; template using WeightedGraph = vector >; using UnWeightedGraph = vector >; template using Matrix = vector >; //unorderd_mapの拡張…https://qiita.com/hamamu/items/4d081751b69aa3bb3557 template size_t HashCombine(const size_t seed,const T &v){ return seed^(std::hash()(v)+0x9e3779b9+(seed<<6)+(seed>>2)); } /* pair用 */ template struct std::hash>{ size_t operator()(const std::pair &keyval) const noexcept { return HashCombine(std::hash()(keyval.first), keyval.second); } }; /* vector用 */ template struct std::hash>{ size_t operator()(const std::vector &keyval) const noexcept { size_t s=0; for (auto&& v: keyval) s=HashCombine(s,v); return s; } }; /* tuple用 */ template struct HashTupleCore{ template size_t operator()(const Tuple &keyval) const noexcept{ size_t s=HashTupleCore()(keyval); return HashCombine(s,std::get(keyval)); } }; template <> struct HashTupleCore<0>{ template size_t operator()(const Tuple &keyval) const noexcept{ return 0; } }; template struct std::hash>{ size_t operator()(const tuple &keyval) const noexcept { return HashTupleCore>::value>()(keyval); } }; ll ceil(ll a, ll b){ return (a + b - 1) / b; } #endif