#953285 (C++23) No.5020 Averaging

• 問題: No.5020 Averaging
• ユーザー: butsurizuki
• 提出日時: 2024-02-20 14:39:47
• 言語: C++23 (gcc 13.3.0 + boost 1.87.0)
• 結果: AC
• 実行時間: 612 ms / 1,000 ms
• コード長: 2,295 bytes
• コンパイル時間: 5,333 ms
• コンパイル使用メモリ: 298,012 KB
• 実行使用メモリ: 6,676 KB
• スコア: 24,278,744
• 最終ジャッジ日時: 2024-02-25 12:43:59
• 合計ジャッジ時間: 37,689 ms
• ジャッジサーバーID （参考情報）: judge11 / judge15

ファイルパターン	結果
other	AC * 50

ソースコード

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
#pragma GCC target("avx2")
#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")
#include <bits/stdc++.h>
#include <sys/time.h>
using namespace std;
long long start_time;
void start_clock(){
  struct timeval tv;
  gettimeofday(&tv, NULL);
  start_time=(tv.tv_sec*1000000+tv.tv_usec);
}
long long current_clock(){
  struct timeval tv;
  gettimeofday(&tv, NULL);
  long long current_time=(tv.tv_sec*1000000+tv.tv_usec);
  // cout << current_time-start_time << "(us)\n";
  return current_time-start_time;
}
using pl=pair<long long,long long>;
long long target=500000000000000000;
double score(pl a){
  long long mg=max(abs(a.first-target),abs(a.second-target));
  return 1.0-0.05*log10(mg+1);
}
pl f(pl a,pl b){
  return {(a.first+b.first)/2,(a.second+b.second)/2};
}
typedef struct{
  vector<pl> card;
  vector<int> hand;
}Game;
double eval(Game &g){
  double x=0.0;
  x=score(g.card[0]);
  return x;
}
#define BSIZE 500
using pdi=pair<double,int>;
int main(){
  start_clock();
  int n;
  cin >> n;
  Game inig;
  inig.card.resize(n);
  for(auto &nx : inig.card){
    cin >> nx.first >> nx.second;
  }
  inig.hand.clear();
  vector<Game> gs={inig};
  double rr=-100.0;
  Game res=inig;
  for(int tr=0;tr<50;tr++){
    vector<Game> ng;
    priority_queue<pdi,vector<pdi>,greater<pdi>> pq;
    for(auto &nx : gs){
      for(int i=0;i<n;i++){
        for(int j=i+1;j<n;j++){
          pl ci=nx.card[i];
          pl cj=nx.card[j];
          pl ck=f(ci,cj);
          nx.card[i]=ck;
          nx.card[j]=ck;
          nx.hand.push_back((i<<6)+j);
          double ce=eval(nx);
          if(rr<ce){
            rr=ce;
            res=nx;
          }
          if(ng.size()<BSIZE){
            pq.push({ce,ng.size()});
            ng.push_back(nx);
          }
          else if(pq.top().first < ce){
            pdi od=pq.top();pq.pop();
            od.first=ce;
            pq.push(od);
            ng[od.second]=nx;
          }
          nx.card[i]=ci;
          nx.card[j]=cj;
          nx.hand.pop_back();
        }
      }
    }
    gs=ng;
  }
  cout << res.hand.size() << "\n";
  for(auto &nx : res.hand){
    cout << (nx/64)+1 << " " << (nx%64)+1 << "\n";
  }
  //for(auto &nx : res.card){
  //  cout << nx.first << " " << nx.second << "\n";
  //}
  return 0;
}
 
 
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX