oj.uz

Submission #942077

# Submission time Handle Problem Language Result Execution time Memory
942077 2024-03-10T08:20:53 Z MilosMilutinovic Harvest (JOI20_harvest) C++14
5 / 100
 5000 ms 32964 KB 
#include<bits/stdc++.h>
 
#define pb push_back
#define fi first
#define se second
#define mp make_pair
 
using namespace std;
 
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
typedef long double ld;
 
template <typename T> bool chkmin(T &x,T y){return x>y?x=y,1:0;}
template <typename T> bool chkmax(T &x,T y){return x<y?x=y,1:0;}
 
ll readint(){
	ll x=0,f=1; char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}

int n,m,l,c,q,T,root,tsz;
int a[200005],b[200005],prv[200005],deg[200005],rt[200005],fenw[200005],dfn[200005],dfo[200005],cnt[10000005],ls[10000005],rs[10000005];
ll d[200005],sum[10000005];
vector<pair<ll,int>> qs[200005];
vector<array<ll,4>> ev;
vector<ll> my[200005];
ll ans[200005];
vector<int> topo,ch[200005];
bool oncyc[200005];

ll calc(ll x){return x>=c?x:((c-x+l-1)/l)*1ll*l+x;}
ll divb(ll x,ll y){return x>=0?x-x%y:x-((y+(x%y))%y);}

int getprv(int x){
	int d=c%l;
	if(a[1]<=x-d){
		int l=1,r=n,p=1;
		while(l<=r){
			int mid=(l+r)/2;
			if(a[mid]<=x-d) l=mid+1,p=mid; 
			else r=mid-1;
		}
		return p;
	}else{
		int l=1,r=n,p=1;
		while(l<=r){
			int mid=(l+r)/2;
			if(x+::l-a[mid]>=d) l=mid+1,p=mid; 
			else r=mid-1;
		}
		return p;
	}
}

void change(int p,int v){
	for(;p<200005;p|=p+1) fenw[p]+=v;  
}

int ask(int p){
	int v=0;
	for(;p>=0;p=(p&(p+1))-1) v+=fenw[p];
	return v;
}

void dfs1(int x,ll d){
	dfn[x]=++T;
	for(auto&p:my[x]) ev.pb({d+p,0,dfn[x],0});
	for(auto&p:qs[x]) ev.pb({p.fi+d,1,x,p.se});
	for(auto&p:my[x]) my[rt[x]].pb(p+d);
	for(int v:ch[x]){
		rt[v]=rt[x];
		ll w=calc(a[x]<=a[v]?a[v]-a[x]:a[v]+l-a[x]);
		dfs1(v,d+w);
	}
	dfo[x]=T;
}

void modifyy(int&id,int l,int r,int p,ll v){
	if(!id) id=++tsz;
	sum[id]+=v;
	cnt[id]++;  
	if(l==r) return;
	int mid=(l+r)/2;
	if(p<=mid) modifyy(ls[id],l,mid,p,v);
	else modifyy(rs[id],mid+1,r,p,v);
}

void modifyx(int&id,int xl,int xr,int yl,int yr,int x,int y,ll v){
	if(!id) id=++tsz;
	modifyy(cnt[id],yl,yr,y,v);
	if(xl==xr) return;
	int mid=(xl+xr)/2;
	if(x<=mid) modifyx(ls[id],xl,mid,yl,yr,x,y,v);
	else modifyx(rs[id],mid+1,xr,yl,yr,x,y,v);
}

pair<ll,int> queryy(int id,int l,int r,int ql,int qr){
	if(ql<=l&&r<=qr) return mp(sum[id],cnt[id]);
	int mid=(l+r)/2;
	if(qr<=mid) return queryy(ls[id],l,mid,ql,qr);
	else if(ql>mid) return queryy(rs[id],mid+1,r,ql,qr);
	pair<ll,int> lv=queryy(ls[id],l,mid,ql,qr);
	pair<ll,int> rv=queryy(rs[id],mid+1,r,ql,qr);
	return mp(lv.fi+rv.fi,lv.se+rv.se);
}

pair<ll,int> queryx(int id,int xl,int xr,int yl,int yr,int x1,int x2,int y1,int y2){
	if(x1<=xl&&xr<=x2) return queryy(cnt[id],yl,yr,y1,y2);
	int mid=(xl+xr)/2;
	if(x2<=mid) return queryx(ls[id],xl,mid,yl,yr,x1,x2,y1,y2);
	else if(x1>mid) return queryx(rs[id],mid+1,xr,yl,yr,x1,x2,y1,y2);
	pair<ll,int> lv=queryx(ls[id],xl,mid,yl,yr,x1,x2,y1,y2);
	pair<ll,int> rv=queryx(rs[id],mid+1,xr,yl,yr,x1,x2,y1,y2);
	return mp(lv.fi+rv.fi,lv.se+rv.se);
}

void clry(int id,int l,int r){
	if(!id) return;
	int mid=(l+r)/2;
	clry(ls[id],l,mid);
	clry(rs[id],mid+1,r);
	ls[id]=rs[id]=cnt[id]=sum[id]=0;
}

void clrx(int id,int xl,int xr,int yl,int yr){
	if(!id) return;
	clry(cnt[id],yl,yr);
	int mid=(xl+xr)/2;
	clrx(ls[id],xl,mid,yl,yr);
	clrx(rs[id],mid+1,xr,yl,yr);
	ls[id]=rs[id]=cnt[id]=sum[id]=0;
}

int main(){
	n=readint(); m=readint(); l=readint(); c=readint();
	for(int i=1;i<=n;i++) a[i]=readint();
	for(int i=1;i<=m;i++) b[i]=readint();
	q=readint();
	for(int i=1;i<=q;i++){
		int v=readint();
		ll t=readint();
		qs[v].pb(mp(t,i));
	}
	for(int i=1;i<=n;i++) prv[i]=getprv(a[i]);
	for(int i=1;i<=m;i++){
		int idx;
		if(a[1]<=b[i]){
			int l=1,r=n;
			while(l<=r){
				int mid=(l+r)/2;
				if(a[mid]<=b[i]) idx=mid,l=mid+1;
				else r=mid-1;
			}
		}else idx=n;
		my[idx].pb((b[i]-a[idx]+l)%l);
	}
	for(int i=1;i<=n;i++) deg[prv[i]]++;
	for(int i=1;i<=n;i++) if(deg[i]==0) topo.pb(i);
	for(int i=0;i<(int)topo.size();i++){
		int v=topo[i];
		deg[prv[v]]--;
		if(deg[prv[v]]==0) topo.pb(prv[v]);
	}
	for(int i=1;i<=n;i++) oncyc[i]=true;
	for(int i:topo) oncyc[i]=false;
	for(int i=1;i<=n;i++) ch[prv[i]].pb(i);
	for(int i=1;i<=n;i++){
		if(oncyc[i]||!oncyc[prv[i]]) continue;
		rt[i]=prv[i];
		ll w=calc(a[prv[i]]<=a[i]?a[i]-a[prv[i]]:a[i]+l-a[prv[i]]);
		T=0;
		ev.clear();
		dfs1(i,w);
		sort(ev.begin(),ev.end());
		for(auto&p:ev){
			if(p[1]==0){
				change(p[2],1);
			}else{
				ans[p[3]]=ask(dfo[p[2]])-ask(dfn[p[2]]-1);
			}
		}
		for(auto&p:ev) if(p[1]==0) change(p[2],-1);
	}
	for(int i=1;i<=n;i++){
		if(!oncyc[i]) continue;
		vector<int> f;
		int ver=i;
		while(oncyc[ver]){
			f.pb(ver);
			oncyc[ver]=false;
			ver=prv[ver];
		}
		int sz=(int)f.size();
		d[0]=0;
		for(int i=1;i<sz;i++){
			ll w=calc(a[f[i]]<=a[f[i-1]]?a[f[i-1]]-a[f[i]]:a[f[i-1]]+l-a[f[i]]);
			d[i]=d[i-1]+w;
		}
		ll len=d[sz-1]+calc(a[f[0]]<=a[f.back()]?a[f.back()]-a[f[0]]:a[f.back()]+l-a[f[0]]);
		if(f.size()==1) len=calc(l);
		vector<ll> xs1,xs2;
		for(int i=0;i<sz;i++){
			for(auto&p:my[f[i]]){
				xs1.pb(d[i]-p);
				xs2.pb(d[i]-p-divb(d[i]-p,len));
			}
			for(auto&p:qs[f[i]]){
				xs1.pb(p.fi-d[i]);
				xs1.pb(-p.fi+d[i]);
				xs2.pb(len-(p.fi-d[i]-divb(p.fi-d[i],len)));
			}
		}
		sort(xs1.begin(),xs1.end());
		xs1.erase(unique(xs1.begin(),xs1.end()),xs1.end());
		sort(xs2.begin(),xs2.end());
		xs2.erase(unique(xs2.begin(),xs2.end()),xs2.end());
		int sz1=(int)xs1.size(),sz2=(int)xs2.size();
		/*for(int i=0;i<sz;i++){
			for(auto&p:my[f[i]]){
				int idx1=(int)(lower_bound(xs1.begin(),xs1.end(),d[i]-p)-xs1.begin());
				int idx2=(int)(lower_bound(xs2.begin(),xs2.end(),d[i]-p-divb(d[i]-p,len))-xs2.begin());
				modifyx(root,0,sz1-1,0,sz2-1,idx1,idx2,divb(d[i]-p,len));
			}
			for(auto&p:qs[f[i]]){
				int idx1=(int)(lower_bound(xs1.begin(),xs1.end(),-(p.fi-d[i]))-xs1.begin());
				int idx2=(int)(lower_bound(xs2.begin(),xs2.end(),len-(p.fi-d[i]-divb(p.fi-d[i],len)))-xs2.begin());
				assert(idx1>=0&&idx1<sz1);
				pair<ll,int> tmp=queryx(root,0,sz1-1,0,sz2-1,idx1,sz1-1,0,sz2-1);
				ll s=(tmp.fi+tmp.se*divb(p.fi-d[i],len))/len+tmp.se;
				ans[p.se]+=s+queryx(root,0,sz1-1,0,sz2-1,idx1,sz1-1,idx2,sz2-1).se;
			}
		}*/
		vector<pair<ll,ll>> vec;
		for(int i=0;i<sz;i++){
			for(auto&p:my[f[i]]){
				vec.pb(mp(d[i]-p,divb(d[i]-p,len)));
			}
			for(auto&p:qs[f[i]]){
				ll cv=p.fi-d[i],cm=divb(cv,len);
				for(auto&pp:vec){
					if(cv+pp.fi<0) continue;
					ans[p.se]+=(pp.se+cm)/len+1+(pp.fi-pp.se+cv-cm>=len?1:0);
				}
			}
		}
		clrx(root,0,sz-1,0,sz2-1);
		root=tsz=0;
		for(int i=0;i<sz;i++){
			for(int j=0;j<i;j++){
				for(auto&p:my[f[i]]){
					for(auto&pp:qs[f[j]]){
						ans[pp.se]+=max(0ll,pp.fi-(i<=j?d[j]-d[i]:d[j]+len-d[i])-p+len)/len;
					}
				}
			}
		}
	}
	for(int i=1;i<=q;i++) printf("%lld\n",ans[i]);
}

Compilation message

harvest.cpp: In function 'int main()':
harvest.cpp:222:7: warning: unused variable 'sz1' [-Wunused-variable]
  222 |   int sz1=(int)xs1.size(),sz2=(int)xs2.size();
      |       ^~~

Subtask #1
5.0 / 5.0

# Verdict Execution time Memory Grader output
1 Correct 6 ms 21080 KB Output is correct
2 Correct 6 ms 21084 KB Output is correct
3 Correct 86 ms 21428 KB Output is correct
4 Correct 10 ms 25628 KB Output is correct
5 Correct 7 ms 26140 KB Output is correct
6 Correct 9 ms 26292 KB Output is correct
7 Correct 7 ms 25996 KB Output is correct
8 Correct 6 ms 25468 KB Output is correct
9 Correct 7 ms 25276 KB Output is correct
10 Correct 7 ms 25664 KB Output is correct
11 Correct 7 ms 25560 KB Output is correct
12 Correct 34 ms 21336 KB Output is correct
13 Correct 58 ms 21340 KB Output is correct
14 Correct 41 ms 21340 KB Output is correct
15 Correct 7 ms 25848 KB Output is correct
16 Correct 7 ms 25884 KB Output is correct
17 Correct 7 ms 25884 KB Output is correct
18 Correct 7 ms 25880 KB Output is correct
19 Correct 7 ms 25884 KB Output is correct
20 Correct 7 ms 25884 KB Output is correct

Subtask #2
0 / 20.0

# Verdict Execution time Memory Grader output
1 Execution timed out 5021 ms 32964 KB Time limit exceeded
2 Halted 0 ms 0 KB -

Subtask #3
0 / 75.0

# Verdict Execution time Memory Grader output
1 Correct 6 ms 21080 KB Output is correct
2 Correct 6 ms 21084 KB Output is correct
3 Correct 86 ms 21428 KB Output is correct
4 Correct 10 ms 25628 KB Output is correct
5 Correct 7 ms 26140 KB Output is correct
6 Correct 9 ms 26292 KB Output is correct
7 Correct 7 ms 25996 KB Output is correct
8 Correct 6 ms 25468 KB Output is correct
9 Correct 7 ms 25276 KB Output is correct
10 Correct 7 ms 25664 KB Output is correct
11 Correct 7 ms 25560 KB Output is correct
12 Correct 34 ms 21336 KB Output is correct
13 Correct 58 ms 21340 KB Output is correct
14 Correct 41 ms 21340 KB Output is correct
15 Correct 7 ms 25848 KB Output is correct
16 Correct 7 ms 25884 KB Output is correct
17 Correct 7 ms 25884 KB Output is correct
18 Correct 7 ms 25880 KB Output is correct
19 Correct 7 ms 25884 KB Output is correct
20 Correct 7 ms 25884 KB Output is correct
21 Execution timed out 5021 ms 32964 KB Time limit exceeded
22 Halted 0 ms 0 KB -