結果

問題 No.1361 [Zelkova 4th Tune *] QUADRUPLE-SEQUENCEの詩
ユーザー SSRS
提出日時 2021-01-22 22:42:54
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 3,326 bytes
コンパイル時間 2,419 ms
コンパイル使用メモリ 179,288 KB
実行使用メモリ 13,276 KB
最終ジャッジ日時 2024-12-28 03:57:12
合計ジャッジ時間 18,787 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 15 RE * 59
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
const long long INF = 1000000000000000000;
int main(){
  int K, L, M, N;
  long long S;
  cin >> K >> L >> M >> N >> S;
  vector<int> A(K);
  for (int i = 0; i < K; i++){
    cin >> A[i];
  }
  vector<int> B(L);
  for (int i = 0; i < L; i++){
    cin >> B[i];
  }
  vector<int> C(M);
  for (int i = 0; i < M; i++){
    cin >> C[i];
  }
  vector<int> D(N);
  for (int i = 0; i < N; i++){
    cin >> D[i];
  }
  vector<long long> X(K * L);
  for (int i = 0; i < K; i++){
    for (int j = 0; j < L; j++){
      X[i * L + j] = A[i] * B[j];
    }
  }
  vector<long long> Y(M * N);
  for (int i = 0; i < M; i++){
    for (int j = 0; j < N; j++){
      Y[i * N + j] = C[i] * D[j];
    }
  }
  sort(Y.begin(), Y.end());
  vector<long long> pos, neg;
  long long zero = 0;
  for (int i = 0; i < K * L; i++){
    if (X[i] > 0){
      pos.push_back(X[i]);
    }
    if (X[i] < 0){
      neg.push_back(X[i]);
    }
    if (X[i] == 0){
      zero++;
    }
  }
  sort(pos.begin(), pos.end());
  sort(neg.begin(), neg.end(), greater<long long>());
  int pcnt = pos.size();
  int ncnt = neg.size();
  long long tv = -INF;
  long long fv = INF;
  while (fv - tv > 1){
    long long mid = (tv + fv) / 2;
    long long cnt = 0;
    if (mid > 0){
      cnt += zero * (N * M);
    }
    int pc = 0;
    for (int i = 0; i < pcnt; i++){
      if (pc > 0){
        assert(pos[i] * Y[pc - 1] < mid);
      }
      while (pc < N * M){
        if (pos[i] * Y[pc] < mid){
          pc++;
        } else {
          break;
        }
      }
      cnt += pc;
    }
    int nc = 0;
    for (int i = 0; i < ncnt; i++){
      if (nc > 0){
        assert(neg[i] * Y[nc - 1] >= mid);
      }
      while (nc < N * M){
        if (neg[i] * Y[nc] >= mid){
          nc++;
        } else {
          break;
        }
      }
      cnt += N * M - nc;
    }
    if (cnt < S){
      tv = mid;
    } else {
      fv = mid;
    }
  }
  cout << tv << endl;
  for (int i = 0; i < K * L; i++){
    int id = -1;
    if (X[i] >= 0){
      int tv2 = -1;
      int fv2 = N * M;
      while (fv2 - tv2 > 1){
        int mid2 = (tv2 + fv2) / 2;
        if (X[i] * Y[mid2] <= tv){
          tv2 = mid2;
        } else {
          fv2 = mid2;
        }
      }
      if (tv2 != -1){
        if (X[i] * Y[tv2] == tv){
          id = tv2;
        }
      }
    } else {
      int tv2 = N * M;
      int fv2 = -1;
      while (tv2 - fv2 > 1){
        int mid2 = (tv2 + fv2) / 2;
        if (X[i] * Y[mid2] <= tv){
          tv2 = mid2;
        } else {
          fv2 = mid2;
        }
      }
      if (tv2 != N * M){
        if (X[i] * Y[tv2] == tv){
          id = tv2;
        }
      }
    }
    if (id != -1){
      bool ok1 = false;
      for (int j = 0; j < K; j++){
        for (int k = 0; k < L; k++){
          if (A[j] * B[k] == X[i]){
            cout << A[j] << ' ' << B[k] << ' ';
            ok1 = true;
            break;
          }
        }
        if (ok1){
          break;
        }
      }
      bool ok2 = false;
      for (int j = 0; j < M; j++){
        for (int k = 0; k < N; k++){
          if (C[j] * D[k] == Y[id]){
            cout << C[j] << ' ' << D[k] << endl;
            ok2 = true;
            break;
          }
        }
        if (ok2){
          break;
        }
      }
      break;
    }
  }
}
0