結果

問題 No.1965 Heavier
ユーザー fumofumofunifumofumofuni
提出日時 2022-06-03 22:48:02
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 466 ms / 2,000 ms
コード長 3,891 bytes
コンパイル時間 2,397 ms
コンパイル使用メモリ 212,996 KB
最終ジャッジ日時 2025-01-29 17:54:45
ジャッジサーバーID
(参考情報)
judge5 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include<bits/stdc++.h>
using namespace std;
#define rep(i,n) for(ll i=0;i<n;i++)
#define repl(i,l,r) for(ll i=(l);i<(r);i++)
#define per(i,n) for(ll i=(n)-1;i>=0;i--)
#define perl(i,r,l) for(ll i=r-1;i>=l;i--)
#define fi first
#define se second
#define pb push_back
#define ins insert
#define pqueue(x) priority_queue<x,vector<x>,greater<x>>
#define all(x) (x).begin(),(x).end()
#define CST(x) cout<<fixed<<setprecision(x)
#define rev(x) reverse(x);
using ll=long long;
using vl=vector<ll>;
using vvl=vector<vector<ll>>;
using pl=pair<ll,ll>;
using vpl=vector<pl>;
using vvpl=vector<vpl>;
const ll MOD=1000000007;
const ll MOD9=998244353;
const int inf=1e9+10;
const ll INF=4e18;
const ll dy[8]={1,0,-1,0,1,1,-1,-1};
const ll dx[8]={0,1,0,-1,1,-1,1,-1};
template <typename T> inline bool chmax(T &a, T b) {
return ((a < b) ? (a = b, true) : (false));
}
template <typename T> inline bool chmin(T &a, T b) {
return ((a > b) ? (a = b, true) : (false));
}
//0-indexed,2
template <class T>
struct SegTree {
private:
int n;//
vector<T> data;// vector
T def; //
function<T(T, T)> operation; // 使
function<T(T, T)> change;// 使
T find(int a, int b) {
T val_left = def, val_right = def;
for (a += (n - 1), b += (n - 1); a < b; a >>= 1, b >>= 1)
{
if ((a & 1) == 0){
val_left = operation(val_left, data[a]);
}
if ((b & 1) == 0){
val_right = operation(data[--b],val_right);
}
}
return operation(val_left, val_right);
}
public:
// _n:, _def:, _operation:,
// _change:
SegTree(size_t _n, T _def, function<T(T, T)> _operation,
function<T(T, T)> _change=[](T a,T b){return b;})
: def(_def), operation(_operation), change(_change) {
n = 1;
while (n < _n) {
n *= 2;
}
data = vector<T>(2 * n - 1, def);
}
void set(int i, T x) { data[i + n - 1] = x; }
void build() {
for (int k=n-2;k>=0;k--) data[k] = operation(data[2*k+1],data[2*k+2]);
}
// i(0-indexed)x
void update(int i, T x) {
i += n - 1;
data[i] = change(data[i], x);
while (i > 0) {
i = (i - 1) / 2;
data[i] = operation(data[i * 2 + 1], data[i * 2 + 2]);
}
}
T all_prod(){
return data[0];
}
// [a, b)
T query(int a, int b) {
//return _query(a, b, 0, 0, n);
return find(a,b);
}
//
T operator[](int i) {
return data[i + n - 1];
}
};
int main(){
ll n;cin >> n;
vl a(n),b(n);rep(i,n)cin >> a[i];rep(i,n)cin >> b[i];
auto pls=b;
sort(all(pls));pls.erase(unique(all(pls)),pls.end());
auto m=pls.size();
vector<multiset<ll>> al(m),ar(m);
rep(i,m)al[i].insert(-INF),ar[i].insert(-INF);
SegTree<ll> stl(m,-INF,[](ll a,ll b){return max(a,b);}),str=stl;
ll ans=-n;
ll right=0;
for(ll l=0;l<n;l++){
while(right<n){
ll p=lower_bound(all(pls),b[right])-pls.begin();
auto lm=stl.query(0,p),rm=str.query(p,m);
if(lm<a[right]-b[right]&&rm<a[right]+b[right]){
al[p].insert(a[right]-b[right]);
ar[p].insert(a[right]+b[right]);
stl.update(p,*al[p].rbegin());
str.update(p,*ar[p].rbegin());
right++;
}
else break;
}
ans+=right-l;
//cout << l <<" " << right << endl;
ll p=lower_bound(all(pls),b[l])-pls.begin();
al[p].erase(al[p].find(a[l]-b[l]));
ar[p].erase(ar[p].find(a[l]+b[l]));
stl.update(p,*al[p].rbegin());
str.update(p,*ar[p].rbegin());
}
cout << ans << endl;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0