#include using namespace std; #define rep(i, a, n) for(int i = a; i < n; i++) #define rrep(i, a, n) for(int i = a; i >= n; i--) #define ll long long #define pii pair #define pll pair constexpr ll MOD = 1000000007; constexpr int IINF = 1001001001; constexpr ll INF = 1LL<<60; template void chmax(t&a,u b){if(a void chmin(t&a,u b){if(b(chrono::system_clock::now().time_since_epoch()).count(); int n; cin >> n; vector> card(n), init(n); rep(i, 0, n){ ll a, b; cin >> a >> b; card[i] = {a, b, i}; init[i] = {a, b, i}; } sort(card.begin(), card.end()); vector> ope; rep(_, 0, 2){ rep(__, 0, 25){ // 1~25回目は表、26~50回目は裏が最大, 最小となるカードについて操作を行う auto [a, c, i] = card[0]; auto [b, d, j] = card[n-1]; card[0] = {(a+b)/2, (c+d)/2, i}; card[n-1] = {(a+b)/2, (c+d)/2, j}; ope.push_back({i, j}); sort(card.begin(), card.end()); } rep(i, 0, n){ auto [a, c, j] = card[i]; card[i] = {c, a, j}; } } // 1番のカード auto get_card1 = [&](vector> card_) -> pair { rep(i, 0, n){ auto [a, c, j] = card_[i]; if(!j){ return {a, c}; } } return {0, 0}; }; // スコア auto score = [&](pair card_) -> ll { ll base = 5e17; return (ll)(2e6-1e5*log10(max(abs(card_.first-base), abs(card_.second-base))+1)); }; // 操作 auto simulation = [&](vector> ope_) -> ll { vector> card_ = init; for(auto [u, v]: ope_){ auto [a, c, i] = card_[u]; auto [b, d, j] = card_[v]; card_[u] = {(a+b)/2, (c+d)/2, i}; card_[v] = {(a+b)/2, (c+d)/2, j}; } return score(get_card1(card_)); }; ll best_score = simulation(ope); // 山登り int cnt = 0; while(true){ auto now_time = chrono::duration_cast(chrono::system_clock::now().time_since_epoch()).count(); if(now_time-start_time >= TIME_LIMIT) break; vector> new_ope = ope; if(cnt&1){ int num = rand()%(ope.size()); int u = rand()%n, v = rand()%n; while(u == v) v = rand()%n; new_ope[num] = {u, v}; ll new_score = simulation(new_ope); if(best_score < new_score){ swap(ope, new_ope); swap(best_score, new_score); } }else{ int u = rand()%(ope.size()), v = rand()%(ope.size()); while(u == v) v = rand()%(ope.size()); swap(new_ope[u], new_ope[v]); ll new_score = simulation(new_ope); // 温度関数 double temp = start_temp+(end_temp-start_temp)*((double)(now_time-start_time))/TIME_LIMIT; // 遷移確率関数(最大化の場合) double prob = exp((double)(new_score-best_score)/temp); // 確率probで遷移する if (prob > (rand()%IINF)/(double)IINF) { swap(ope, new_ope); swap(best_score, new_score); } } cnt++; } //cout << simulation(ope) << endl; // 出力 cout << ope.size() << endl; for(auto [u, v]: ope){ cout << u+1 << ' ' << v+1 << endl; } return 0; }