結果

問題 No.1623 三角形の制作
コンテスト
ユーザー siman
提出日時 2022-08-18 18:08:27
言語 C++17(clang)
(17.0.6 + boost 1.87.0)
結果
AC  
実行時間 895 ms / 2,000 ms
コード長 1,473 bytes
コンパイル時間 2,345 ms
コンパイル使用メモリ 144,748 KB
実行使用メモリ 6,144 KB
最終ジャッジ日時 2024-10-06 11:52:43
合計ジャッジ時間 19,962 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cassert>
#include <cmath>
#include <algorithm>
#include <iostream>
#include <iomanip>
#include <climits>
#include <map>
#include <queue>
#include <set>
#include <cstring>
#include <vector>

using namespace std;
typedef long long ll;

int main() {
  ll N;
  cin >> N;

  vector<int> R(N);
  vector<int> G(N);
  vector<int> B(N);
  map<int, ll> r_counter;
  map<int, ll> g_counter;
  map<int, ll> b_counter;
  for (int i = 0; i < N; ++i) {
    cin >> R[i];
    r_counter[R[i]]++;
  }
  for (int i = 0; i < N; ++i) {
    cin >> G[i];
    g_counter[G[i]]++;
  }
  for (int i = 0; i < N; ++i) {
    cin >> B[i];
    b_counter[B[i]]++;
  }
  sort(R.begin(), R.end());
  sort(G.begin(), G.end());
  sort(B.begin(), B.end());

  ll ans = N * N * N;
  ll dp[6010];
  memset(dp, 0, sizeof(dp));

  for (int i = 0; i < N; ++i) {
    int r = R[i];
    ll c1 = N - (ll) (lower_bound(G.begin(), G.end(), r + 1) - G.begin());
    ll c2 = N - (ll) (lower_bound(B.begin(), B.end(), r + 1) - B.begin());

    // fprintf(stderr, "c1: %lld, c2: %lld\n", c1, c2);
    ans -= c1 * N;
    ans -= c2 * N;
    ans += c1 * c2;
  }

  for (int v = 1; v <= 3000; ++v) {
    for (int w = 1; w <= 3000; ++w) {
      dp[v + w] += g_counter[v] * b_counter[w];
    }
  }

  for (int v = 1; v <= 3000; ++v) {
    if (r_counter[v] == 0) continue;
    for (int w = 1; w <= v; ++w) {
      if (dp[w] == 0) continue;
      ans -= r_counter[v] * dp[w];
    }
  }

  cout << ans << endl;

  return 0;
}
0