結果

問題 No.1623 三角形の制作
コンテスト
ユーザー sak
提出日時 2021-06-20 02:08:23
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 190 ms / 2,000 ms
コード長 1,371 bytes
コンパイル時間 2,153 ms
コンパイル使用メモリ 196,596 KB
最終ジャッジ日時 2025-01-22 10:25:22
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
typedef pair<ll, ll> p_ll;

template<class T>
void debug(T itr1, T itr2) { auto now = itr1; while(now<itr2) { cout << *now << " "; now++; } cout << endl; }
#define repr(i,from,to) for (ll i=(ll)from; i<(ll)to; i++)
#define all(vec) vec.begin(), vec.end()
#define rep(i,N) repr(i,0,N)
#define per(i,N) for (ll i=(ll)N-1; i>=0; i--)
#define popcount __builtin_popcount

const ll LLINF = pow(2,61)-1;
const ll INF = pow(2,30)-1;

ll gcd(ll a, ll b) { if (a<b) swap(a,b); return b==0 ? a : gcd(b, a%b); }
ll lcm(ll a, ll b) { return a/gcd(a,b)*b; }

int main() {
  ll n; cin >> n;
  ll r[n]; rep(i,n) cin >> r[i];
  ll g[n]; rep(i,n) cin >> g[i];
  ll b[n]; rep(i,n) cin >> b[i];

  assert(1<=n && n<=2*pow(10,5));
  rep(i,n) assert(1<=r[i] && r[i]<=3*pow(10,3));
  rep(i,n) assert(1<=g[i] && g[i]<=3*pow(10,3));
  rep(i,n) assert(1<=b[i] && b[i]<=3*pow(10,3));

  ll MA = 3 * pow(10,3) + 10;
  vector<ll> cr(MA), cg(MA), cb(MA);
  rep(i,n) { cr[r[i]]++; cb[b[i]]++; cg[g[i]]++; }
  vector<ll> rui_cr(MA); rep(i,MA) rui_cr[i] = i==0 ? 0 : rui_cr[i-1] + cr[i];
  // debug(all(rui_cr));

  ll result = 0;
  rep(ig,MA) rep(ib,MA) {
    ll from = max(max(ig, ib), abs(ig-ib)+1), to = min(MA-1, ig+ib-1);
    result += (rui_cr[to] - rui_cr[from-1]) * cg[ig] * cb[ib];
  }
  cout << result << endl;
  return 0;
}
0