#include using namespace std; typedef long long ll; typedef __int128_t lll; typedef long double ld; typedef pair pl; typedef tuple tt; typedef vector vl; typedef vector vvl; typedef vector vvvl; typedef vector vd; typedef vector vvd; typedef vector vvvd; typedef vector vb; typedef vector vvb; typedef vector vvvb; typedef vector vc; typedef vector vvc; typedef vector vvvc; typedef vector vp; typedef vector vt; typedef vector vs; const ll INF = 1000000010; const ll INFL = INF * INF; const double MYPI = acos(-1); const ld epsilon = 1e-10; #define ovl(a, b, c, d, e, ...) e #define pb push_back #define eb emplace_back #define MP make_pair #define DEBUG(...) DEBUG_(#__VA_ARGS__, __VA_ARGS__) #define REP1(i, n) for (int i = 0; i < (n); i++) #define REP2(i, l, r) for (int i = (l); i < (r); i++) #define REP3(i, l, r, d) for (int i = (l); i < (r); i += (d)) #define REP(...) ovl(__VA_ARGS__, REP3, REP2, REP1)(__VA_ARGS__) #define RREP1(i, n) for (int i = (n)-1; i >= 0; i--) #define RREP2(i, l, r) for (int i = (r)-1; i >= (l); i--) #define RREP3(i, l, r, d) for (int i = (r)-1; i >= (l); i -= (d)) #define RREP(...) ovl(__VA_ARGS__, RREP3, RREP2, RREP1)(__VA_ARGS__) #define EACH1(e, a) for (auto &e : a) #define EACH2(x, y, a) for (auto &[x, y] : a) #define EACH3(x, y, z, a) for (auto &[x, y, z] : a) #define EACH(...) ovl(__VA_ARGS__, EACH3, EACH2, EACH1)(__VA_ARGS__) #define ALL(x) begin(x), end(x) #define RALL(x) (x).rbegin(), (x).rend() #define sz(x) (ll) x.size() #define LB(a, x) (lower_bound(ALL(a), x) - a.begin()) #define UB(a, x) (upper_bound(ALL(a), x) - a.begin()) #define FLG(x, i) (((x) >> (i)) & 1) #define CNTBIT(x) __builtin_popcountll(x) #define TOPBIT(t) (t == 0 ? -1 : 63 - __builtin_clzll(t)) #define IN(x, a, b) ((a) <= (x) && (x) < (b)) template void OUT(const T &a, const Ts &...b) { cout << a; (cout << ... << (cout << " ", b)); cout << endl; } template void DEBUG_(string_view name, const T &a) { cout << name << ": " << a << endl; } template void DEBUG_(string_view name, const vector &a) { cout << name << ": "; REP(i, sz(a)) cout << a[i] << " "; cout << endl; } template bool chmax(T1 &x, const T2 &y) { return x < y ? x = y, 1 : 0; } template bool chmin(T1 &x, const T2 &y) { return x > y ? x = y, 1 : 0; } // struct UnionFind { // explicit UnionFind(int n) : _n(n), par_or_sz(n, -1) {} // int merge(int a, int b) { // assert(IN(a, 0, _n) && IN(b, 0, _n)); // int x = find(a), y = find(b); // if (x == y) return x; // if (-par_or_sz[x] < -par_or_sz[y]) swap(x, y); // par_or_sz[x] += par_or_sz[y]; // return par_or_sz[y] = x; // } // int find(int a) { // assert(IN(a, 0, _n)); // return par_or_sz[a] < 0 ? a : par_or_sz[a] = find(par_or_sz[a]); // } // bool same(int a, int b) { // assert(IN(a, 0, _n) && IN(b, 0, _n)); // return find(a) == find(b); // } // protected: // int _n; // vl par_or_sz; // }; int main() { cin.tie(0); ios::sync_with_stdio(false); int N; cin >> N; vector> XYs(N); REP(i, N) { int x, y; cin >> x >> y; XYs[i] = {x, y}; } vector> edges; REP(i, N) { auto [xi, yi] = XYs[i]; REP2(j, i + 1, N) { auto [xj, yj] = XYs[j]; edges.eb((xi - xj) * (xi - xj) + (yi - yj) * (yi - yj), i, j); } } sort(ALL(edges)); set seen; vp ans; EACH3(w, i, j, edges) { if (seen.find(i) == seen.end() && seen.find(j) == seen.end()) { ans.eb(i, j); seen.insert(i); seen.insert(j); } } OUT(ans.size()); EACH(e, ans) { OUT(e.first + 1, e.second + 1); } return 0; }