結果

問題 No.5020 Averaging
ユーザー A8pfA8pf
提出日時 2024-02-25 16:26:24
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 940 ms / 1,000 ms
コード長 5,859 bytes
コンパイル時間 3,022 ms
コンパイル使用メモリ 234,252 KB
実行使用メモリ 7,552 KB
スコア 18,961,420
最終ジャッジ日時 2024-02-25 16:27:17
合計ジャッジ時間 51,163 ms
ジャッジサーバーID
(参考情報)
judge15 / judge10
純コード判定しない問題か言語
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 894 ms
7,364 KB
testcase_01 AC 905 ms
7,552 KB
testcase_02 AC 907 ms
7,400 KB
testcase_03 AC 903 ms
7,416 KB
testcase_04 AC 898 ms
7,292 KB
testcase_05 AC 891 ms
7,380 KB
testcase_06 AC 939 ms
7,416 KB
testcase_07 AC 920 ms
7,352 KB
testcase_08 AC 881 ms
7,316 KB
testcase_09 AC 889 ms
7,552 KB
testcase_10 AC 918 ms
7,552 KB
testcase_11 AC 918 ms
7,552 KB
testcase_12 AC 900 ms
7,276 KB
testcase_13 AC 874 ms
7,552 KB
testcase_14 AC 905 ms
7,552 KB
testcase_15 AC 923 ms
7,416 KB
testcase_16 AC 880 ms
7,552 KB
testcase_17 AC 888 ms
7,552 KB
testcase_18 AC 918 ms
7,392 KB
testcase_19 AC 896 ms
7,312 KB
testcase_20 AC 918 ms
7,552 KB
testcase_21 AC 905 ms
7,380 KB
testcase_22 AC 892 ms
7,552 KB
testcase_23 AC 883 ms
7,552 KB
testcase_24 AC 911 ms
7,424 KB
testcase_25 AC 899 ms
7,508 KB
testcase_26 AC 920 ms
7,344 KB
testcase_27 AC 881 ms
7,552 KB
testcase_28 AC 935 ms
7,552 KB
testcase_29 AC 878 ms
7,404 KB
testcase_30 AC 928 ms
7,408 KB
testcase_31 AC 886 ms
7,420 KB
testcase_32 AC 926 ms
7,392 KB
testcase_33 AC 901 ms
7,552 KB
testcase_34 AC 920 ms
7,340 KB
testcase_35 AC 883 ms
7,552 KB
testcase_36 AC 940 ms
7,552 KB
testcase_37 AC 899 ms
7,552 KB
testcase_38 AC 884 ms
7,360 KB
testcase_39 AC 921 ms
7,304 KB
testcase_40 AC 912 ms
7,296 KB
testcase_41 AC 924 ms
7,320 KB
testcase_42 AC 918 ms
7,404 KB
testcase_43 AC 929 ms
7,552 KB
testcase_44 AC 903 ms
7,368 KB
testcase_45 AC 892 ms
7,332 KB
testcase_46 AC 884 ms
7,288 KB
testcase_47 AC 897 ms
7,552 KB
testcase_48 AC 890 ms
7,392 KB
testcase_49 AC 885 ms
7,552 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include<bits/stdc++.h>
using namespace std;
using ull=unsigned long long;

std::mt19937 mt_for_action(0); // 行動選択用の乱数生成器を初期化

// 128bit整数
std::ostream &operator<<(std::ostream &dest, __int128_t value) {
  std::ostream::sentry s(dest);
  if (s) {
    __uint128_t tmp = value < 0 ? -value : value;
    char buffer[128];
    char *d = std::end(buffer);
    do {
      --d;
      *d = "0123456789"[tmp % 10];
      tmp /= 10;
    } while (tmp != 0);
    if (value < 0) {
      --d;
      *d = '-';
    }
    int len = std::end(buffer) - d;
    if (dest.rdbuf()->sputn(d, len) != len) {
      dest.setstate(std::ios_base::badbit);
    }
  }
  return dest;
}

__int128 parse(string &s) {
  __int128 ret = 0;
  for (int i = 0; i < s.length(); i++)
    if ('0' <= s[i] && s[i] <= '9')
      ret = 10 * ret + s[i] - '0';
  return ret;
}

// 128bit整数ここまで
constexpr int MAX_TURN = 50;
constexpr __int128 TARGET = 500000000000000000;

struct Card {
  __int128 omote, ura;
};

class GameState
{
  private:
  int _turn;
  vector<Card> _cards;
  vector<pair<int, int>> _commandHistory;

  public:
  __int128 nowScore;
  pair<int, int> first_action_;
  GameState(){};
  GameState(vector<unsigned long long> omotes, vector<unsigned long long> uras): _turn(0), nowScore(0) {
    for (int i = 0; i < omotes.size(); i++) {
      _cards.push_back({omotes[i], uras[i]});
    }
  }

  __int128 getEstimatedScore() {
    __int128 score = 0;
    // 1枚目のみ別に評価
    __int128 omotediff = TARGET - _cards[0].omote;
    if (omotediff > 0) {
      omotediff *= -1;
    }
    score += omotediff;
    __int128 uradiff = TARGET - _cards[0].ura;
    if (uradiff > 0) {
      uradiff *= -1;
    }
    score += uradiff;
    // 満点ならその時点で終了
    if (omotediff == 0 && uradiff == 0) {
      score += 9999999999999999;
    }
    // 2枚目以降
    const Card OTHER_TARGET = { TARGET * 2 - _cards[0].omote, TARGET * 2 - _cards[0].ura };
    for (int i = 1; i < _cards.size(); i++) {
      __int128 omotediff = OTHER_TARGET.omote - _cards[i].omote;
      if (omotediff > 0) {
        omotediff *= -1;
      }
      score += omotediff * _turn;
      __int128 uradiff = OTHER_TARGET.ura - _cards[i].ura;
      if (uradiff > 0) {
        uradiff *= -1;
      }
      score += uradiff * _turn;
      // 満点ならその時点で終了
      if (omotediff == 0 && uradiff == 0) {
        score += 99999999;
      }
    }
    return score;
  }

  // 1ターン進める
  void advance(int u, int v)
  {
    _commandHistory.push_back({u, v});
    _turn++;
    Card newCard = {(_cards[u].omote + _cards[v].omote) / 2, (_cards[u].ura + _cards[v].ura) / 2};
    _cards[u] = newCard;
    _cards[v] = newCard;
    nowScore = getEstimatedScore();
  }

  // 合法手を列挙する
  vector<pair<int, int>> legalActions() const
  {
    vector<pair<int, int>> ret;
    for (int i = 1; i < _cards.size(); i++) {
      ret.push_back({0, i});
    }
    return ret;
  }
  
  // 1番目に触らないランダムなアクションを返す
  pair<int, int> randomAction() const
  {
    uniform_int_distribution<int> dist(1, (int)_cards.size() - 1);
    int u = dist(mt_for_action);
    int v = dist(mt_for_action);
    while (u == v) {
      v = dist(mt_for_action);
    }
    return {u, v};
  }

  // ゲームが終了しているかどうか
  bool isDone() const
  {
    return _turn == MAX_TURN;
  }

  // 結果の出力を行う
  string outputCommand() {
    stringstream ss;
    ss << _commandHistory.size() << endl;
    for (auto command : _commandHistory) {
      ss << command.first + 1 << " " << command.second + 1 << endl;
    }
    return ss.str();
  }

  // 現在の実スコアを出力する
  double getScore() const {
    ull v = max( max((ull)_cards[0].omote, (ull)TARGET) - min((ull)_cards[0].omote, (ull)TARGET), max((ull)_cards[0].ura, (ull)TARGET) - min((ull)_cards[0].ura, (ull)TARGET) );
    long long score = 2000000 - 100000 * log10(static_cast<long double>(v));
    return score;
  }
};

bool operator<(const GameState &state1, const GameState &state2)
{
  return state1.nowScore < state2.nowScore;
}

pair<int,int> beamSearchAction(const GameState &initState, const int beamWidth, const int beamDepth)
{
  priority_queue<GameState> nowBeam;
  GameState bestState;

  nowBeam.push(initState);
  for (int t = 0; t < beamDepth; t++)
  {
      std::priority_queue<GameState> nextBeam;
      for (int i = 0; i < beamWidth; i++)
      {
          if (nowBeam.empty())
              break;
          GameState nowState = nowBeam.top();
          nowBeam.pop();
          auto legal_actions = nowState.legalActions();
          for (const auto &action : legal_actions)
          {
              GameState nextState = nowState;
              nextState.advance(action.first, action.second);
              if (t == 0)
                  nextState.first_action_ = action;
              nextBeam.push(nextState);
          }
      }
      nowBeam = nextBeam;
      bestState = nowBeam.top();
      if (bestState.isDone())
      {
        break;
      }
  }
  if (bestState.getScore() < initState.getScore()) {
    auto random_action = initState.randomAction();
    return random_action;
  }
  return bestState.first_action_;
}

pair<vector<unsigned long long>, vector<unsigned long long>> getFromInput() {
  int n;
  cin >> n;
  vector<unsigned long long> omotes(n), uras(n);
  for (int i = 0; i < n; i++) {
    cin >> omotes[i] >> uras[i];
  }
  return {omotes, uras};
}

int main() {
  // 本番用
  auto input = getFromInput();
  GameState state(input.first, input.second);
  while (!state.isDone()) {
    auto action = beamSearchAction(state, 20, 25);
    state.advance(action.first, action.second);
    // cerr<<state.nowScore<<endl;
    // cerr<<state.getScore()<<endl;
  }
  cout << state.outputCommand();
}
0