QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#297657#7872. 崩坏天际线zhouhuanyi0 4936ms21084kbC++145.8kb2024-01-04 21:54:262024-01-04 21:54:27

Judging History

你现在查看的是最新测评结果

  • [2024-01-04 21:54:27]
  • 评测
  • 测评结果:0
  • 用时:4936ms
  • 内存:21084kb
  • [2024-01-04 21:54:26]
  • 提交

answer

#include<iostream>
#include<cstdio>
#include<vector>
#include<algorithm>
#include<cmath>
#define N 100000
#define mod 998244353
using namespace std;
const int inv2=(mod+1)>>1;
int read()
{
	char c=0;
	int sum=0;
	while (c<'0'||c>'9') c=getchar();
	while ('0'<=c&&c<='9') sum=sum*10+c-'0',c=getchar();
	return sum;
}
void Adder(int &x,int d)
{
	x+=d;
	if (x>=mod) x-=mod;
	return;
}
void Adder2(int &x,int d)
{
	x+=d;
	if (x<0) x+=mod;
	return;
}
int MD(int x)
{
	return x>=mod?x-mod:x;
}
int MD2(int x)
{
	return x<0?x+mod:x;
}
struct reads
{
	int l,r,t;
};
reads tong[N+1];
struct rds
{
	int x,t;
};
rds st[N+1];
int n,q,szt,ans,pst[N+1],op[N+1],l[N+1],r[N+1],x[N+1],fa[N+1],dep[N+1],dep2[N+1],pw[N+1],invpw[N+1],block[N+1],length,leng,delta[N+1],sdelta[N+1],ls[N+1],rs[N+1],pv[N+1],nt[N+1],sps[N+1],num[N+1],num2[N+1],snum[N+1],snum2[N+1],dque[N+1],top,tongs[N+1],lengths,sz[N+1],hs[N+1],belong[N+1],sbelong[N+1],sq[N+1],cater,cl[N+1],tops[N+1],s[N+1],depth[N+1],sdepth[N+1],s2[N+1],dp[N+1],DSP[N+1],DSP2[N+1],DP[N+1],DP2[N+1],len;
vector<reads>v[N+1];
void dfs(int x,int l,int r)
{
	pv[x]=l,nt[x]=r,sz[x]=1;
	if (ls[x]) fa[ls[x]]=x,depth[ls[x]]=depth[x]+1,dfs(ls[x],l,x);
	if (rs[x]) fa[rs[x]]=x,depth[rs[x]]=depth[x]+1,dfs(rs[x],x,r);
	sz[x]=sz[ls[x]]+sz[rs[x]]+1,hs[x]=(sz[ls[x]]>sz[rs[x]])?ls[x]:rs[x];
	return;
}
void dfs2(int x)
{
	tops[hs[x]]=tops[x],dp[x]=pst[nt[x]]-pst[pv[x]];
	if (ls[x]^rs[x]^hs[x]) tops[ls[x]^rs[x]^hs[x]]=ls[x]^rs[x]^hs[x];
	if (ls[x]) dfs2(ls[x]),Adder(dp[x],1ll*dp[ls[x]]*inv2%mod);
	if (rs[x]) dfs2(rs[x]),Adder(dp[x],1ll*dp[rs[x]]*inv2%mod);
	return;
}
void dfs3(int x)
{
	if (ls[x]) sdepth[ls[x]]=sdepth[x]+1,dfs3(ls[x]);
	if (rs[x]) sdepth[rs[x]]=sdepth[x]+1,dfs3(rs[x]);
	return;
}
int lca(int x,int y)
{
	while (tops[x]!=tops[y])
	{
		if (depth[tops[x]]<depth[tops[y]]) swap(x,y);
		x=fa[tops[x]];
	}
	if (depth[x]>depth[y]) swap(x,y);
	return x;
}
int calc(int l,int r,int sl,int sr,int d)
{
	int res=1ll*(r-l)*invpw[d]%mod,rst=0,ds;
	if (sl<=sr)
	{
		ds=lca(sl,sr);
		Adder(rst,1ll*MD2(DP[sl]-DP[ds])*pw[dep[ds]-1]%mod);
		Adder(rst,1ll*MD2(DP2[sr]-DP2[ds])*pw[dep2[ds]-1]%mod);
		Adder(rst,1ll*MD2(1-invpw[dep[sl]-dep[ds]])*inv2%mod*MD2(-l)%mod);
		Adder(rst,1ll*MD2(DSP[sl]-DSP[ds])*pw[dep[ds]-1]%mod);
		Adder(rst,1ll*MD2(1-invpw[dep2[sr]-dep2[ds]])*inv2%mod*r%mod);
		Adder(rst,1ll*MD2(DSP2[sr]-DSP2[ds])*pw[dep2[ds]-1]%mod);
		Adder(rst,1ll*(r-l)*inv2%mod);
		Adder2(res,-1ll*rst*invpw[d]%mod);
	}
	return res;
}
int main()
{
	int ps=0,sl,sr;
	pw[0]=invpw[0]=1;
	for (int i=1;i<=N;++i) pw[i]=2ll*pw[i-1]%mod,invpw[i]=1ll*invpw[i-1]*inv2%mod;
	n=read(),q=read();
	for (int i=1;i<=q;++i)
	{
		op[i]=read();
		if (op[i]==1) l[i]=read(),r[i]=read(),tong[++length]=(reads){l[i],r[i],i};
		else x[i]=read(),st[++leng]=(rds){x[i],i};
	}
	reverse(tong+1,tong+length+1),reverse(st+1,st+leng+1),szt=leng;
	for (int i=1;i<=leng;++i) block[i]=(i-1)/szt+1;
	for (int i=1;i<=length;++i)
	{
		while (ps+1<=leng&&st[ps+1].t>=tong[i].t) ++ps;
		if (ps) v[block[ps]].push_back(tong[i]);
		else Adder(ans,tong[i].r-tong[i].l);
	}
	for (int i=1;i<=block[leng];++i)
		if (!v[i].empty())
		{
			if (i!=1)
			{
				lengths=top=cater=0;
				for (int j=1;j<=n;++j) delta[j]=0;
				for (int j=1;j<=(i-1)*szt;++j) delta[st[j].x]=st[j].t;
				for (int j=1;j<=n;++j)
					if (delta[j])
						tongs[++lengths]=j,ls[lengths]=rs[lengths]=hs[lengths]=fa[lengths]=0,sdelta[lengths]=delta[j],pst[lengths]=j;
				for (int j=1;j<=lengths;++j)
				{
					while (top&&sdelta[dque[top]]>sdelta[j]) ls[j]=dque[top],top--;
					if (top) rs[dque[top]]=j;
					dque[++top]=j;
				}
				depth[dque[1]]=1,dfs(dque[1],0,lengths+1),tops[dque[1]]=dque[1],pst[0]=1,pst[lengths+1]=n,dfs2(dque[1]);
				for (int j=lengths;j>=1;--j)
				{
					if (nt[j]!=lengths+1) dep[j]=dep[nt[j]]+1,DP[j]=MD(DP[nt[j]]+1ll*MD(1ll*dp[rs[j]]*inv2%mod+(pst[nt[j]]-pst[pv[j]]))*invpw[dep[j]]%mod),DSP[j]=MD(DSP[nt[j]]+1ll*invpw[dep[j]]*pst[pv[j]]%mod);
				    else dep[j]=1,DP[j]=1ll*MD(1ll*dp[rs[j]]*inv2%mod+(pst[nt[j]]-pst[pv[j]]))*inv2%mod,DSP[j]=1ll*invpw[dep[j]]*pst[pv[j]]%mod;
				}
				for (int j=1;j<=lengths;++j)
				{
					if (pv[j]!=0) dep2[j]=dep2[pv[j]]+1,DP2[j]=MD(DP2[pv[j]]+1ll*MD(1ll*dp[ls[j]]*inv2%mod+(pst[nt[j]]-pst[pv[j]]))*invpw[dep2[j]]%mod),DSP2[j]=MD(DSP2[pv[j]]+1ll*invpw[dep2[j]]*MD2(-pst[nt[j]])%mod);
				    else dep2[j]=1,DP2[j]=1ll*MD(1ll*dp[ls[j]]*inv2%mod+(pst[nt[j]]-pst[pv[j]]))*inv2%mod,DSP2[j]=1ll*invpw[dep2[j]]*MD2(-pst[nt[j]])%mod;
				}
			}
			for (int j=1;j<=n;++j) belong[j]=0;
			for (int j=(i-1)*szt+1;j<=min(i*szt,leng);++j) num[st[j].x]=lower_bound(tongs+1,tongs+lengths+1,st[j].x)-tongs-1,num2[st[j].x]=lower_bound(tongs+1,tongs+lengths+1,st[j].x+1)-tongs,belong[st[j].x]=st[j].t;
			for (int j=1;j<=n;++j)
				if (belong[j])
					sq[++cater]=j,snum[cater]=num[j],snum2[cater]=num2[j],sbelong[cater]=belong[j],sps[cater]=j;
			for (int j=0;j<v[i].size();++j)
			{
				sl=lower_bound(tongs+1,tongs+lengths+1,v[i][j].l+1)-tongs,sr=lower_bound(tongs+1,tongs+lengths+1,v[i][j].r)-tongs-1,top=len=0;
				for (int k=1;k<=cater;++k) ls[k]=rs[k]=0;
				for (int k=1;k<=cater;++k)
					if (v[i][j].l<sq[k]&&sq[k]<v[i][j].r&&sbelong[k]>=v[i][j].t)
					{
						while (top&&sbelong[dque[top]]>sbelong[k]) ls[k]=dque[top],top--;
						if (top) rs[dque[top]]=k;
						cl[++len]=dque[++top]=k;
					}
				if (top)
				{
					sdepth[dque[1]]=1,dfs3(dque[1]);
					Adder(ans,calc(v[i][j].l,sps[cl[1]],sl,snum[cl[1]],sdepth[cl[1]]));
					for (int k=1;k<=len-1;++k) Adder(ans,calc(sps[cl[k]],sps[cl[k+1]],snum2[cl[k]],snum[cl[k+1]],max(sdepth[cl[k]],sdepth[cl[k+1]])));
					Adder(ans,calc(sps[cl[len]],v[i][j].r,snum2[cl[len]],sr,sdepth[cl[len]]));
				}
				else Adder(ans,calc(v[i][j].l,v[i][j].r,sl,sr,0));
			}
		}
    printf("%d\n",ans);
	return 0;
}

詳細信息

Subtask #1:

score: 0
Wrong Answer

Test #1:

score: 0
Wrong Answer
time: 0ms
memory: 21084kb

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:

169902410

result:

wrong answer 1st lines differ - expected: '855279801', found: '169902410'

Subtask #2:

score: 0
Skipped

Dependency #1:

0%

Subtask #3:

score: 0
Time Limit Exceeded

Test #13:

score: 40
Accepted
time: 3850ms
memory: 19064kb

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:

733099543

result:

ok single line: '733099543'

Test #14:

score: 0
Accepted
time: 2896ms
memory: 20516kb

input:

49951 43686
1 21796 23464
1 29304 46959
1 5034 41719
1 7779 35334
1 27566 36486
1 20347 26165
1 12508 30387
1 18363 20335
1 8540 21417
1 5728 49086
1 46038 47603
1 10371 15910
1 27293 43572
1 18915 45279
1 7388 48342
1 6802 43746
1 4361 40049
1 41177 43375
1 23287 48354
1 37097 41733
1 2406 11638
1 ...

output:

792296531

result:

ok single line: '792296531'

Test #15:

score: 0
Accepted
time: 4936ms
memory: 18564kb

input:

49914 43874
1 8935 40963
1 4425 44317
1 1769 45855
1 2436 40257
1 1778 47216
1 383 42149
1 5398 40732
1 1079 43346
1 6578 41660
1 9689 45985
1 6131 42681
1 8862 47431
1 3979 46189
1 6456 43485
1 2028 46574
1 3802 47787
1 6990 41659
1 9221 41204
1 2271 43554
1 8018 45280
1 9344 43917
1 6623 41152
1 7...

output:

831211412

result:

ok single line: '831211412'

Test #16:

score: -40
Time Limit Exceeded

input:

50000 50000
1 1310 49344
1 5755 44255
1 3582 41465
1 6800 42160
1 1651 44584
1 7967 44410
1 3116 48795
1 1855 41120
1 27 42294
1 2455 49629
1 4196 42487
1 7070 44542
1 136 42053
1 5715 44222
1 8794 43115
1 4048 45579
1 635 46703
1 9246 41055
1 3678 41276
1 4871 41715
1 1659 44679
1 1639 46392
1 2479...

output:


result:


Subtask #4:

score: 0
Skipped

Dependency #1:

0%