結果

問題 No.1361 [Zelkova 4th Tune *] QUADRUPLE-SEQUENCEの詩
ユーザー SSRSSSRS
提出日時 2021-01-22 22:18:13
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 2,966 bytes
コンパイル時間 1,798 ms
コンパイル使用メモリ 176,664 KB
実行使用メモリ 15,652 KB
最終ジャッジ日時 2024-12-28 02:11:00
合計ジャッジ時間 77,852 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 57 TLE * 17
権限があれば一括ダウンロードができます

ソースコード

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());
  long long tv = -INF;
  long long fv = INF;
  while (fv - tv > 1){
    long long mid = (tv + fv) / 2;
    long long cnt = 0;
    for (int i = 0; i < K * L; i++){
      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] < mid){
            tv2 = mid2;
          } else {
            fv2 = mid2;
          }
        }
        cnt += fv2;
      } else {
        int tv2 = N * M;
        int fv2 = -1;
        while (tv2 - fv2 > 1){
          int mid2 = (tv2 + fv2) / 2;
          if (X[i] * Y[mid2] < mid){
            tv2 = mid2;
          } else {
            fv2 = mid2;
          }
        }
        cnt += N * M - tv2;
      }
    }
    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