結果
| 問題 | No.318 学学学学学 |
| コンテスト | |
| ユーザー |
mugen_1337
|
| 提出日時 | 2020-06-07 13:35:16 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 422 ms / 2,000 ms |
| コード長 | 6,569 bytes |
| コンパイル時間 | 2,760 ms |
| コンパイル使用メモリ | 217,588 KB |
| 最終ジャッジ日時 | 2025-01-10 23:48:52 |
|
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 26 |
ソースコード
#include<bits/stdc++.h>
using namespace std;
#define ALL(x) x.begin(),x.end()
#define rep(i,n) for(int i=0;i<(n);i++)
#define debug(v) cout<<#v<<":";for(auto x:v){cout<<x<<' ';}cout<<endl;
#define INF 1000000000
#define mod 1000000007
using ll=long long;
const ll LINF=1001002003004005006ll;
int dx[]={1,0,-1,0};
int dy[]={0,1,0,-1};
// ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return true;}return false;}
template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return true;}return false;}
class Beats{
const ll inf=LLONG_MAX;
vector<ll> mx,smx,mxc,mi,smi,mic,sum,lval,ladd;
int sz;
// smx[k]<x<mx[k]の部分の更新
void update_node_max(int k,ll x){
sum[k]+=(x-mx[k])*mxc[k];
if(mx[k]==mi[k]) mx[k]=x,mi[k]=x;
else if(mx[k]==smi[k]) mx[k]=x,smi[k]=x;
else mx[k]=x;
if(lval[k]!=inf and x<lval[k]) lval[k]=x;
}
void update_node_min(int k,ll x){
sum[k]+=(x-mi[k])*mic[k];
if(mx[k]==mi[k]) mx[k]=x,mi[k]=x;
else if(smx[k]==mi[k]) smx[k]=x,mi[k]=x;
else mi[k]=x;
if(lval[k]!=inf and lval[k]<x) lval[k]=x;
}
void update_node_add(int k,ll len,ll x){
mx[k]+=x;
if(smx[k]!=-inf)smx[k]+=x;
mi[k]+=x;
if(smi[k]!=inf) smi[k]+=x;
sum[k]+=x*len;
if(lval[k]!=inf)lval[k]+=x;
else ladd[k]+=x;
}
// 遅延分を処理,node-kは今の値,子は古い値かもしれない
void push(int k,ll len){
if(k>=sz) return ;
if(lval[k]!=inf){
update_all(2*k,len/2,lval[k]);
update_all(2*k+1,len/2,lval[k]);
lval[k]=inf;
return ;
}
if(ladd[k]!=0){
update_node_add(2*k,len/2,ladd[k]);
update_node_add(2*k+1,len/2,ladd[k]);
ladd[k]=0;
}
// 子が古い値なら更新.ここがポイント
if(mx[2*k]>mx[k]) update_node_max(2*k,mx[k]);
if(mx[2*k+1]>mx[k])update_node_max(2*k+1,mx[k]);
if(mi[2*k]<mi[k]) update_node_min(2*k,mi[k]);
if(mi[2*k+1]<mi[k])update_node_min(2*k+1,mi[k]);
}
// 子からnode-kの更新
void update_from_children(int k){
sum[k]=sum[2*k]+sum[2*k+1];
if(mx[2*k]<mx[2*k+1]){
mx[k]=mx[2*k+1];mxc[k]=mxc[2*k+1];
smx[k]=max(mx[2*k],smx[2*k+1]);
}else if(mx[2*k]>mx[2*k+1]){
mx[k]=mx[2*k];mxc[k]=mxc[2*k];
smx[k]=max(smx[2*k],mx[2*k+1]);
}else{
mx[k]=mx[2*k];mxc[k]=mxc[2*k]+mxc[2*k+1];
smx[k]=max(smx[2*k],smx[2*k+1]);
}
if(mi[2*k]<mi[2*k+1]){
mi[k]=mi[2*k];mic[k]=mic[2*k];
smi[k]=min(smi[2*k],mi[2*k+1]);
}else if(mi[2*k]>mi[2*k+1]){
mi[k]=mi[2*k+1];mic[k]=mic[2*k+1];
smi[k]=min(mi[2*k],smi[2*k+1]);
}else{
mi[k]=mi[2*k];mic[k]=mic[2*k]+mic[2*k+1];
smi[k]=min(smi[2*k],smi[2*k+1]);
}
}
void update_all(int k,ll len,ll x){
mx[k]=x,smx[k]=-inf,mi[k]=x,smi[k]=inf;
mxc[k]=len,mic[k]=len;
sum[k]=x*len;
lval[k]=x,ladd[k]=0;
}
public:
Beats(int n){
sz=1;
while(sz<n)sz<<=1;
mx.resize(2*sz,-inf);smx.resize(2*sz,-inf);mxc.resize(2*sz,0);
mi.resize(2*sz,inf);smi.resize(2*sz,inf);mic.resize(2*sz,0);
sum.resize(2*sz,0);ladd.resize(2*sz,0);lval.resize(2*sz,inf);
}
void set(int k,ll x){
mx[k+sz]=x;mxc[k+sz]=1;
mi[k+sz]=x;mic[k+sz]=1;
sum[k+sz]=x;
}
void build(){
for(int i=sz-1;i;i--) update_from_children(i);
}
// [a,b)を更新,[l,r)探索区間,k-node
void chmin(int a,int b,ll x,int k=1,int l=0,int r=-1){
if(r==-1) r=sz;
if(r<=a or b<=l or mx[k]<=x) return ;
if(a<=l and r<=b and smx[k]<x){
update_node_max(k,x);
return ;
}
push(k,r-l);
chmin(a,b,x,2*k,l,(l+r)/2);
chmin(a,b,x,2*k+1,(l+r)/2,r);
update_from_children(k);
}
void chmax(int a,int b,ll x,int k=1,int l=0,int r=-1){
if(r==-1) r=sz;
if(r<=a or b<=l or mi[k]>=x) return ;
if(a<=l and r<=b and smi[k]>x){
update_node_min(k,x);
return ;
}
push(k,r-l);
chmax(a,b,x,2*k,l,(l+r)/2);
chmax(a,b,x,2*k+1,(l+r)/2,r);
update_from_children(k);
}
void add(int a,int b,ll x,int k=1,int l=0,int r=-1){
if(r==-1) r=sz;
if(r<=a or b<=l) return ;
if(a<=l and r<=b){
update_node_add(k,r-l,x);
return ;
}
push(k,r-l);
add(a,b,x,2*k,l,(l+r)/2);
add(a,b,x,2*k+1,(l+r)/2,r);
update_from_children(k);
}
void update(int a,int b,ll x,int k=1,int l=0,int r=-1){
if(r==-1)r=sz;
if(r<=a or b<=l) return ;
if(a<=l and r<=b){
update_all(k,r-l,x);
return ;
}
push(k,r-l);
update(a,b,x,2*k,l,(l+r)/2);
update(a,b,x,2*k+1,(l+r)/2,r);
update_from_children(k);
}
ll query_sum(int a,int b,int k=1,int l=0,int r=-1){
if(r==-1)r=sz;
if(r<=a or b<=l)return 0;
if(a<=l and r<=b) return sum[k];
push(k,r-l);
ll lsum=query_sum(a,b,2*k,l,(l+r)/2);
ll rsum=query_sum(a,b,2*k+1,(l+r)/2,r);
return lsum+rsum;
}
ll query_min(int a,int b,int k=1,int l=0,int r=-1){
if(r==-1)r=sz;
if(b<=l or r<=a) return inf;
if(a<=l and r<=b) return mi[k];
push(k,r-l);
ll lmin=query_min(a,b,2*k,l,(l+r)/2);
ll rmin=query_min(a,b,2*k+1,(l+r)/2,r);
return min(lmin,rmin);
}
ll query_max(int a,int b,int k=1,int l=0,int r=-1){
if(r==-1)r=sz;
if(b<=l or r<=a) return -inf;
if(a<=l and r<=b) return mx[k];
push(k,r-l);
ll lmax=query_max(a,b,2*k,l,(l+r)/2);
ll rmax=query_max(a,b,2*k+1,(l+r)/2,r);
return max(lmax,rmax);
}
};
signed main(){
cin.tie(0);
ios::sync_with_stdio(0);
int n;cin>>n;
map<int,int> l,r;
set<int> st;
rep(i,n){
int x;cin>>x;
st.insert(x);
if(!l.count(x)) l[x]=i;
else l[x]=min(l[x],i);
if(!r.count(x)) r[x]=i;
else r[x]=max(r[x],i);
}
Beats seg(n);
seg.build();
for(auto x:st) seg.update(l[x],r[x]+1,x);
rep(i,n) cout<<seg.query_max(i,i+1)<<(i+1==n?"\n":" ");
return 0;
}
mugen_1337