QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#282208 | #7872. 崩坏天际线 | sjc061031 | 0 | 0ms | 0kb | C++14 | 7.0kb | 2023-12-11 16:25:40 | 2023-12-11 16:25:40 |
answer
#include <bits/stdc++.h>
#include<ext/pb_ds/assoc_container.hpp>
#include<ext/pb_ds/hash_policy.hpp>
#define MOD 998244353
using namespace __gnu_pbds;
using namespace std;
struct custom_hash {
static uint64_t splitmix64(uint64_t x) {
x += 0x9e3779b97f4a7c15;
x = (x ^ (x >> 30)) * 0xbf58476d1ce4e5b9;
x = (x ^ (x >> 27)) * 0x94d049bb133111eb;
return x ^ (x >> 31);
}
size_t operator()(uint64_t x) const {
static const uint64_t FIXED_RANDOM = chrono::steady_clock::now().time_since_epoch().count();
return splitmix64(x + FIXED_RANDOM);
}
};
int n,q,ans,block,pw[50010],flag[50010],l[50010],r[50010],x[50010],f[50010],g[50010],h[50010],b[50010],c[50010];
int dat[200010],tot,rt[50010][2],L[5000010],R[5000010],tag[5000010],sum[5000010];
vector<pair<int,int> > v0[50010],v1[50010];
const int rev2=499122177;
inline void add1(int &x,int y)
{
x+=y;
if(x>=MOD) x-=MOD;
}
inline void add2(int &x,int y)
{
x+=y;
if(x<0) x+=MOD;
}
inline void build1(int o,int l,int r)
{
dat[o]=q+1;
if(l<r){
int mid=(l+r)/2;
build1(o*2,l,mid);
build1(o*2+1,mid+1,r);
}
}
inline void update1(int o,int l,int r,int x,int v)
{
if(l==r){
dat[o]=v;
return;
}
int mid=(l+r)/2;
if(x<=mid) update1(o*2,l,mid,x,v);
else update1(o*2+1,mid+1,r,x,v);
dat[o]=min(dat[o*2],dat[o*2+1]);
}
inline int query1(int o,int l,int r,int x,int y)
{
if(r<x||l>y) return q+1;
if(x<=l&&r<=y) return dat[o];
int mid=(l+r)/2,res=q+1;
res=min(res,query1(o*2,l,mid,x,y));
res=min(res,query1(o*2+1,mid+1,r,x,y));
return res;
}
inline void pushdown(int o)
{
if(tag[o]!=1){
if(L[o]){
tag[L[o]]=1ll*tag[L[o]]*tag[o]%MOD;sum[L[o]]=1ll*sum[L[o]]*tag[o]%MOD;
}
if(R[o]){
tag[R[o]]=1ll*tag[R[o]]*tag[o]%MOD;sum[R[o]]=1ll*sum[R[o]]*tag[o]%MOD;
}
tag[o]=1;
}
}
inline int split0(int o,int l,int r,int x)
{
if(!o) return 0;
int root=++tot;tag[root]=1;
if(l==r){
sum[root]=sum[o];sum[o]=0;
return root;
}
int mid=(l+r)/2;
pushdown(o);
if(x<=mid){
L[root]=split0(L[o],l,mid,x);R[root]=0;
}
else{
L[root]=L[o];L[o]=0;R[root]=split0(R[o],mid+1,r,x);
}
sum[o]=0;
if(L[o]) add1(sum[o],sum[L[o]]);
if(R[o]) add1(sum[o],sum[R[o]]);
sum[root]=0;
if(L[root]) add1(sum[root],sum[L[root]]);
if(R[root]) add1(sum[root],sum[R[root]]);
return root;
}
inline int split1(int o,int l,int r,int x)
{
if(!o) return 0;
int root=++tot;tag[root]=1;
if(l==r){
sum[root]=sum[o];sum[o]=0;
return root;
}
int mid=(l+r)/2;
pushdown(o);
if(x<=mid){
R[root]=R[o];R[o]=0;L[root]=split1(L[o],l,mid,x);
}
else{
R[root]=split1(R[o],mid+1,r,x);L[root]=0;
}
sum[o]=0;
if(L[o]) add1(sum[o],sum[L[o]]);
if(R[o]) add1(sum[o],sum[R[o]]);
sum[root]=0;
if(L[root]) add1(sum[root],sum[L[root]]);
if(R[root]) add1(sum[root],sum[R[root]]);
return root;
}
inline void update(int o,int l,int r,int x,int v)
{
if(l==r){
add1(sum[o],v);
return;
}
int mid=(l+r)/2;
pushdown(o);
if(x<=mid){
if(!L[o]) L[o]=++tot,tag[L[o]]=1;
update(L[o],l,mid,x,v);
}
else{
if(!R[o]) R[o]=++tot,tag[R[o]]=1;
update(R[o],mid+1,r,x,v);
}
sum[o]=0;
if(L[o]) add1(sum[o],sum[L[o]]);
if(R[o]) add1(sum[o],sum[R[o]]);
}
inline void calc0(int o,int l,int r,int x)
{
if(!o) return;
if(l==r){
add1(ans,1ll*sum[o]*(x-l)%MOD);
return;
}
int mid=(l+r)/2;
pushdown(o);
calc0(L[o],l,mid,x);
calc0(R[o],mid+1,r,x);
}
inline void calc1(int o,int l,int r,int x)
{
if(!o) return;
if(l==r){
add1(ans,1ll*sum[o]*(l-x)%MOD);
return;
}
int mid=(l+r)/2;
pushdown(o);
calc1(L[o],l,mid,x);
calc1(R[o],mid+1,r,x);
}
signed main()
{
freopen("test.in","r",stdin);
freopen("test.out","w",stdout);
ios::sync_with_stdio(false);cin.tie(0);
cin>>n>>q;
pw[0]=1;
for(int i=1;i<=n;i++) pw[i]=1ll*pw[i-1]*rev2%MOD;
for(int i=1;i<=q;i++){
cin>>flag[i];
if(flag[i]==1){
cin>>l[i]>>r[i];
}
else{
cin>>x[i];
}
}
for(int i=1;i<=q;i++){
int ll=i,rr=q,Sum=0,cnt=0;
for(int j=ll;j<=q;j++){
if(flag[j]==1) Sum++,cnt++;
else Sum+=cnt;
if(Sum>n){
rr=j-1;break;
}
}
i=rr;
build1(1,1,n);
for(int j=q;j>=rr+1;j--){
if(flag[j]==2){
update1(1,1,n,x[j],j);
v0[j].clear();v1[j].clear();
}
}
vector<pair<pair<int,int>,int> > nw;
vector<int> u;
for(int j=rr;j>=ll;j--){
if(flag[j]==2){
g[x[j]]=j;
bool flag=false;
for(int k=0;k<(int)u.size();k++){
if(x[j]==u[k]){
flag=true;break;
}
if(x[j]<u[k]){
u.insert(u.begin()+k,x[j]);
flag=true;break;
}
}
if(!flag) u.push_back(x[j]);
}
else{
int id=0;
b[0]=0;
int cnt=0;
for(int k=0;k<(int)u.size();k++) if(l[j]<u[k]&&u[k]<r[j]){
while(id>0&&h[id]>g[u[k]]) id--;
h[++id]=g[u[k]];b[++cnt]=id;
}
id=0;
c[0]=0;cnt=0;
for(int k=(int)u.size()-1;k>=0;k--) if(l[j]<u[k]&&u[k]<r[j]){
while(id>0&&h[id]>g[u[k]]) id--;
h[++id]=g[u[k]];c[++cnt]=id;
}
int pre=l[j],lvl=0;
for(int k=0;k<(int)u.size();k++) if(l[j]<u[k]&&u[k]<r[j]){
nw.push_back(make_pair(make_pair(pre,u[k]),pw[b[lvl]+c[cnt-lvl]]));
pre=u[k];lvl++;
}
nw.push_back(make_pair(make_pair(pre,r[j]),pw[b[lvl]+c[cnt-lvl]]));
}
}
for(int j=0;j<(int)nw.size();j++){
int loc=q+1,l=nw[j].first.first,r=nw[j].first.second,val=nw[j].second;
if(l<r-1) loc=query1(1,1,n,l+1,r-1);
if(loc==q+1) add1(ans,1ll*(r-l)*val%MOD);
else{
v0[loc].push_back(make_pair(l,1ll*val*rev2%MOD));
v1[loc].push_back(make_pair(r,1ll*val*rev2%MOD));
}
}
set<int> s;
set<int>::iterator it;
s.insert(1);s.insert(n);
for(int j=rr+1;j<=q;j++){
if(flag[j]==2){
if(s.find(x[j])!=s.end()) continue;
s.insert(x[j]);
it=s.find(x[j]);
int pre,nxt,u=x[j];
it--;pre=(*it);it++;
it++;nxt=(*it);it--;
f[u]=1ll*f[pre]*rev2%MOD;f[pre]=1ll*f[pre]*rev2%MOD;
rt[u][0]=split0(rt[nxt][0],1,n,u-1);
tag[rt[u][0]]=1ll*tag[rt[u][0]]*rev2%MOD;sum[rt[u][0]]=1ll*sum[rt[u][0]]*rev2%MOD;
add1(f[u],sum[rt[u][0]]);
rt[u][1]=split1(rt[pre][1],1,n,u+1);
tag[rt[u][1]]=1ll*tag[rt[u][1]]*rev2%MOD;sum[rt[u][1]]=1ll*sum[rt[u][1]]*rev2%MOD;
add1(f[pre],sum[rt[u][1]]);
for(int k=0;k<(int)v0[j].size();k++){
if(!rt[u][0]) rt[u][0]=++tot;
update(rt[u][0],1,n,v0[j][k].first,v0[j][k].second);
}
for(int k=0;k<(int)v1[j].size();k++){
if(!rt[u][1]) rt[u][1]=++tot;
update(rt[u][1],1,n,v1[j][k].first,v1[j][k].second);
}
}
}
vector<int> vec;
for(it=s.begin();it!=s.end();it++) vec.push_back(*it);
for(int j=0;j<(int)vec.size()-1;j++){
add1(ans,1ll*f[vec[j]]*(vec[j+1]-vec[j])%MOD);
}
for(int j=0;j<(int)vec.size();j++){
calc0(rt[vec[j]][0],1,n,vec[j]);
calc1(rt[vec[j]][1],1,n,vec[j]);
}
for(int j=0;j<(int)vec.size();j++) f[vec[j]]=0,rt[vec[j]][0]=0,rt[vec[j]][1]=0;
for(int j=1;j<=tot;j++) L[j]=0,R[j]=0,tag[j]=0,sum[j]=0;
tot=0;
}
cout<<ans<<'\n';
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Subtask #1:
score: 0
Dangerous Syscalls
Test #1:
score: 0
Dangerous Syscalls
input:
500 500 1 119 258 2 134 2 417 2 176 2 61 2 60 2 110 1 335 336 1 96 111 2 202 1 138 344 2 358 2 134 1 29 54 1 73 381 1 179 495 2 490 2 418 2 186 2 183 1 168 340 2 78 1 15 27 2 373 1 245 498 1 372 495 2 244 2 63 1 174 490 2 282 2 417 1 272 408 1 109 134 2 303 2 345 1 238 401 1 36 480 1 21 483 2 10 1 3...
output:
result:
Subtask #2:
score: 0
Skipped
Dependency #1:
0%
Subtask #3:
score: 0
Dangerous Syscalls
Test #13:
score: 0
Dangerous Syscalls
input:
50000 50000 1 24367 33007 1 14396 42256 1 6375 22327 1 7892 42501 1 10100 37998 1 6284 48524 1 7357 18164 1 16200 46424 1 18972 34131 1 16849 32591 1 1917 3018 1 19897 30272 1 45044 45753 1 18999 25448 1 5167 31033 1 6182 35335 1 7270 37270 1 12651 39965 1 28896 38022 1 13853 35426 1 35516 48244 1 1...
output:
result:
Subtask #4:
score: 0
Skipped
Dependency #1:
0%