結果
問題 |
No.1496 不思議な数え上げ
|
ユーザー |
|
提出日時 | 2024-04-01 01:33:49 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 868 ms / 3,000 ms |
コード長 | 1,405 bytes |
コンパイル時間 | 2,066 ms |
コンパイル使用メモリ | 199,360 KB |
最終ジャッジ日時 | 2025-02-20 18:35:58 |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 39 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/segtree> using namespace atcoder; using namespace std; using ll = long long; ll INF = (1LL << 60); ll op(ll a, ll b) { return min(a, b); } ll e() { return INF; } int main() { int N;cin >> N; vector<ll> P(N + 1, 0); vector<ll>A(N + 1, 0); for (int i = 1;i <= N;i++)cin >> P[i]; for (int i = 1;i <= N;i++) cin >> A[i]; vector<ll> S(N + 1, 0); for (int i = 1;i <= N;i++) S[i] = S[i - 1] + P[i]; segtree<ll, op, e> seg(P); vector<ll> ans(N + 1, 0); for (int i = 1;i <= N;i++) { ll x = P[i]; int ok = i; int ng = 0; while (ok - ng > 1) { int m = (ok + ng) / 2; if (seg.prod(m, i + 1) != x) { ng = m; } else { ok = m; } } int L = ok; ok = i; ng = N + 1; while (ng - ok > 1) { int m = (ok + ng) / 2; if (seg.prod(i, m + 1) != x) { ng = m; } else { ok = m; } } int R = ok; if (i - L >= R - i) { for (int j = i;j <= R;j++) { auto itr = lower_bound(S.begin(), S.end(), S[j] - A[x]); int t = itr - S.begin(); if (t + 1 > i) continue; int st = max(t + 1, L); ans[x] += i - st + 1; } } else { for (int j = i;j >= L;j--) { auto itr = upper_bound(S.begin(), S.end(), S[j - 1] + A[x]); itr--; int t = itr - S.begin(); if (t < i) continue; int ed = min(t, R); ans[x] += ed - i + 1; } } } for (int x = 1;x <= N;x++) cout << ans[x] << endl; }