7.2

真的，刷题日记是我写题的动力...

今天自习一整天，学个整体二分。

(P3834)[https://www.luogu.com.cn/problem/P3834] 板子

CODE

#include<bits/stdc++.h>
#define fst first
#define sec second
#define mkp(a,b) make_pair(a,b)
#define usetime() (double)clock () / CLOCKS_PER_SEC * 1000.0
using namespace std;
typedef long long LL;
typedef pair<int,int> pii;
const int maxn=2e5+5,inf=1e9+1;
void read(int& x){
	char c;
	bool f=0;
	while((c=getchar())<48) f|=(c==45);
	x=c-48;
	while((c=getchar())>47) x=(x<<3)+(x<<1)+c-48;
	x=(f ? -x : x);
}
struct BIT{
	int f[maxn];
	void add(int x,int p){
		for(int i=x;i<maxn;i+=i&(-i)) f[i]+=p;
	}
	int query(int x){
		int ans=0;
		for(int i=x;i>=1;i-=i&(-i)) ans+=f[i];
		return ans;
	}
}bt;
struct node{
	int x,y,k,id,tp;
};
int n,m;
node q[maxn<<1],q1[maxn<<1],q2[maxn<<1];
int ans[maxn];
void solve(int l,int r,int L,int R){
	if(l>r) return;
	if(L==R){
		for(int i=l;i<=r;i++){
			if(q[i].tp==2) ans[q[i].id]=L;
		}
		return;
	}
	int mid=(L+R)>>1;
	int p1=0,p2=0;
	for(int i=l;i<=r;i++){
		if(q[i].tp==1){
			if(q[i].x<=mid) bt.add(q[i].id,1),q1[++p1]=q[i];
			else q2[++p2]=q[i];
		}
		else{
			int res=bt.query(q[i].y)-bt.query(q[i].x-1);
			if(res>=q[i].k) q1[++p1]=q[i];
			else q[i].k-=res,q2[++p2]=q[i];
		}
	}
	for(int i=1;i<=p1;i++){
		if(q1[i].tp==1) bt.add(q1[i].id,-1);
	}
	for(int i=1;i<=p1;i++) q[l+i-1]=q1[i];
	for(int i=1;i<=p2;i++) q[l+p1+i-1]=q2[i];
	solve(l,l+p1-1,L,mid),solve(l+p1,r,mid+1,R);
}
int main(){
	read(n),read(m);
	for(int i=1;i<=n;i++){
		int x; read(x);
		q[i]=(node){x,-1,-1,i,1};
	}
	for(int i=1;i<=m;i++){
		int x,y,k; read(x),read(y),read(k);
		q[n+i]=(node){x,y,k,i,2};
	}
	solve(1,n+m,-inf,inf);
	for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
	return 0;
}
//^o^

https://www.luogu.com.cn/article/wff0kib6

这里学的，写得不错。

---

以下不是整体二分。

P3576

水紫，降绿我都没意见。

先倒过来，从根节点往叶子节点走，维护 \(l\) 和 \(r\) 表示当前节点的蚁群大小在该值域范围内，到达终点时大小会变为 \(k\)，显然是很好维护的，最后二分看有多少个点在值域范围内即可。

CODE

#include<bits/stdc++.h>
#define fst first
#define sec second
#define mkp(a,b) make_pair(a,b)
#define usetime() (double)clock () / CLOCKS_PER_SEC * 1000.0
using namespace std;
typedef long long LL;
typedef pair<int,int> pii;
const int maxn=1e6+5;
void read(int& x){
	char c;
	bool f=0;
	while((c=getchar())<48) f|=(c==45);
	x=c-48;
	while((c=getchar())>47) x=(x<<3)+(x<<1)+c-48;
	x=(f ? -x : x);
}
void read(LL& x){
	char c;
	bool f=0;
	while((c=getchar())<48) f|=(c==45);
	x=c-48;
	while((c=getchar())>47) x=(x<<3)+(x<<1)+c-48;
	x=(f ? -x : x);
}
int n,g,h1,h2;
LL k,a[maxn],ans=0;
vector<int> mp[maxn];
LL l[maxn],r[maxn];
void dfs(int u,int fa){
	//cout<<u<<endl;
	for(int v : mp[u]){
		if(v==fa) continue;
		l[v]=(mp[u].size()-1)*l[u],l[v]=min(l[v],(LL)1e9+3);
		r[v]=(mp[u].size()-1)*(r[u]+1)-1,r[v]=min(r[v],(LL)1e9+3);
		dfs(v,u);
	}
	if((int)mp[u].size()==1){
		//cout<<u<<' '<<l[u]<<' '<<r[u]<<endl;
		int li=lower_bound(a+1,a+g+1,l[u])-a,ri=upper_bound(a+1,a+g+1,r[u])-a;
		//cout<<"kkk::"<<li<<' '<<ri-1<<endl;
		ans+=(ri-li)*k;
	}
}
int main(){
	read(n),read(g),read(k);
	for(int i=1;i<=g;i++) read(a[i]);
	sort(a+1,a+g+1);
	read(h1),read(h2);
	for(int i=2;i<n;i++){
		int u,v; read(u),read(v);
		mp[u].push_back(v),mp[v].push_back(u);
	}
	l[0]=r[0]=k; mp[0].push_back(h1),mp[0].push_back(h2),mp[h1].push_back(0),mp[h2].push_back(0);
	dfs(0,0);
	printf("%lld",ans);
	return 0;
}
//^o^