二分捏啊哈+wqs二分

已于 2023-09-13 13:25:32 修改
阅读量178 收藏 1
点赞数
文章标签： 算法 c++ 图论
于 2022-10-27 16:30:38 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/cn_tigerhan/article/details/127540291
版权

嗯…这是一个难以开口的专题，不算难，通常是和其他的算法一起呃干一些难以描述的事？使用这个的前提是某一个东西具有单调性而且就是在单调闭区间中进行的。呃分为二分答案和正常二分两种？呃啊其实也差的不算太多，一个可能维护的是序列中的位置，另一个可能维护的是两个数中的某一个值？思路的话通常会将询问等操作，化成判断，枚举一个值看其是否成立，呃啊其实也就不是很难，做多了就懂了。
板子的话要背这个，不然太多太杂了：

while(l<=r)
{
	int mid=(l+r)/2;
	if(check(mid)) ans=mid,l=mid+1;
	else r=mid-1;
}

看一个！P2678 [NOIP2015 提高组] 跳石头来，给大家整个活：首先这个可以用贪心，不过嘛，欸，我选择更难的二分答案，那就要看判断函数啦。我们枚举一个可行的最短跳跃距离的最大值。那么小于这个的石头就要删除掉，并且令cnt++，若cnt>m 那么就说明不可行咯。好！好活当赏！呜呜呜我哭了难受啊。时间的话O（nlogn）吧。

#include<bits/stdc++.h>
using namespace std;
int n,m,L,ans;
int a[500001];
bool check(int lim)
{
	int cnt=0,s=0;
	for(int i=1;i<=n+1;i++)
	{
		if(lim>a[i]-s) cnt++;
		else s=a[i];
	}
	if(cnt<=m) return true;
	else return false;
}
int main()
{
	scanf("%d%d%d",&L,&n,&m);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	a[n+1]=L;int l=0,r=L;
	while(l<=r)	//好的枚举！
	{
		int mid=(l+r)/2;
		if(check(mid)) ans=mid,l=mid+1;
		else r=mid-1;
	}
	printf("%d",ans);
	return 0;
}

下一道也差不多吧P3853 [TJOI2007]路标设置这个和上面那道差不多啦，说的好啊嗯。

#include<bits/stdc++.h>
using namespace std;
int n,m,L,ans;
int a[1000001];
bool check(int lim)//好的更新方式
{
	int cnt=0,s=0;
	for(int i=1;i<=n;i++)
	{
		while(a[i]-s>lim) s+=lim,cnt++;
		s=max(s,a[i]); 
	}
	if(cnt>m) return true;
	else return false;
}
int main()
{
	scanf("%d%d%d",&L,&n,&m);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	int l=0,r=L;
	while(l<=r)//枚举的是区间最小最大长度 
	{
		int mid=(l+r)/2;
		if(check(mid)) l=mid+1;
		else ans=mid,r=mid-1;
	}
	printf("%d",ans);
	return 0;
}

下一题：P1281 书的复制呃这个的话一个小二分不说了，心痛破防累，没啥的，大概，这是，合理的：

#include<bits/stdc++.h>
using namespace std;
int n,m,l,r,ans;
int st[100001],ed[100001],a[100001];
bool check(int lim)
{
	int cnt=1,k=0;
	for(int i=1;i<=n;i++)
	{
		if(k+a[i]>lim) k=0,cnt++;
		k+=a[i];
	}
	if(cnt>m||k>lim) return false;
	else return true;
}
int main()
{
	memset(st,0,sizeof(st));memset(ed,0,sizeof(ed));
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]),r+=a[i];
	while(l<=r)
	{
		int mid=(l+r)/2;
		if(check(mid)) ans=mid,r=mid-1;
		else l=mid+1; 
	}	
	int now=n;
	for(int i=m;i>=1;i--)
	{
		int k=0;ed[i]=now;
		while(k+a[now]<=ans&&now>=1) st[i]=now,k+=a[now--];
	}
	for(int i=1;i<=m;i++) printf("%d %d\n",st[i],ed[i]);
	return 0;
}

有一道有小数数的二分，用我的板子不知道行不行，打打看看？P1570 KC 喝咖啡呃啊有点破防问题在于这道题没什么问题：

#include<bits/stdc++.h>
using namespace std;
int n,m;
double l,r,ans;
struct node
{
	double v,c,dis;
};node e[100001];
bool cmp(const node &x,const node &y)
{
	return x.dis<y.dis;
}
bool check(double lim)
{
	double tot=0;
	for(int i=1;i<=n;i++) e[i].dis=(e[i].c*lim)-e[i].v;
	sort(e+1,e+n+1,cmp);
	for(int i=1;i<=m;i++) tot+=e[i].dis;
	if(tot<=0) return true;
	else return false ;
}
int main() 
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%lf",&e[i].v);
	for(int i=1;i<=n;i++) scanf("%lf",&e[i].c),r=max(r,e[i].v/e[i].c);
	while(l<=r)
	{
		double mid=(l+r)/2;
		if(check(mid)) l=mid+0.00001;
		else r=mid-0.00001,ans=mid;
	}
	printf("%.3lf",ans);
	return 0;
}

下一个吧P1577 切绳子嗯差不多的，不过有一个小小技巧，因为这个只要求答案保留到小数点后 2位(直接舍掉 2 位后的小数)。所以我们可以将这个转化为整数，后面再除二不过没什么用：

#include<bits/stdc++.h>
using namespace std;
int n,m;
double L[1000001];
int a[1000001];
int l=0,r=10000001,ans;
bool check(int lim)
{
	int cnt=0;
	for(int i=1;i<=n;i++) cnt+=a[i]/lim;
	if(cnt>=m) return true ;
	else return false ;
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%lf",&L[i]),a[i]=(int)(L[i]*100.0);
	while(l<=r)
	{
		int mid=(l+r)/2;
		if(mid==0) 
		{
			ans=mid;
			break;	
		} 
		if(check(mid)) l=mid+1,ans=mid;
		else r=mid-1;
	}
	printf("%.2lf",ans/100.0);
	return 0;
}

其实也没什么的，也就这样了，再做一两道题也就行了吧P1843 奶牛晒衣服：

#include<bits/stdc++.h>
using namespace std;
int n,m,a,b,ans;
int w[600001];
bool check(int lim)
{
	int cnt=0;
	for(int i=1;i<=n;i++) 
	{
		int k=w[i]-lim*a;
		if(k>0) 
		{
			if(k%b==0) cnt+=(k/b);
			else cnt+=(k/b)+1; 
		}
	}
	if(cnt<=lim) return true ;
	else return false ;
}
int main()
{
	scanf("%d%d%d",&n,&a,&b);
	int l=0,r=0;
	for(int i=1;i<=n;i++) scanf("%d",&w[i]),r=max(r,w[i]/a+1);
	while(l<=r)
	{
		int mid=(l+r)/2;
		if(check(mid)) r=mid-1,ans=mid;
		else l=mid+1;
	}
	printf("%d",ans);
	return 0;
}

下一题P2370 yyy2015c01 的 U 盘呃本来是dp+二分的不过后面发现可以直接dp，因为数据好像挺水，导致dp可以直接过哦。

#include<bits/stdc++.h>
using namespace std;
int n,m,p,s;
int w[100001],v[100001];
struct node
{
	int w,v;
};node e[100001];
bool cmp(const node &x,const node &y)
{
	return x.w<y.w;
}
int f[10001];
int main()
{
	scanf("%d%d%d",&n,&p,&s);
	for(int i=1;i<=n;i++) scanf("%d%d",&e[i].w,&e[i].v);
	sort(e+1,e+n+1,cmp);
	for(int i=1;i<=n;i++)
	{
		for(int j=s;j>=e[i].w;j--)
		{
			f[j]=max(f[j],f[j-e[i].w]+e[i].v);
			if(f[j]>=p) 
			{
				printf("%d",e[i].w);
				return 0;
			}
		}
	}
	printf("No Solution!");
	return 0;
}

那二分就写到这吧，不算难的东西就放下了再打一会儿倍增。好吧做倍增做到一道二分：P4403 [BJWC2008]秦腾与教学评估，考虑从奇偶数性入手，那么发现性质：奇数加偶数是奇数，偶数加偶数是偶数，那么既然题目保证了唯一性，就直接考虑用前缀和维护若一个点是奇数，则它后面的点都是，那么再套上一个二分答案就可以解决了。打一下？

//我们这题直接维护的是前缀和，所以我们可以直接拿这
//(min(lim,e[i])-s[i])/d[i]+1; 来维护一手1到lim的前缀和 注意要把一开始的点加上 
#include<bits/stdc++.h>
#define int long long 
using namespace std;
int n,m,t,ans;
int s[400000],e[400000],d[400000];
bool check(int lim)
{
	int sum=0;
	for(int i=1;i<=n;i++)
	{
		if(s[i]<=lim) sum+=(min(lim,e[i])-s[i])/d[i]+1;//s起始点与e边界咯，然后求出区间长度内的站着的人 
	} 
	if(sum%2==1) return true ;
	else return false;
}
signed main()
{
	scanf("%lld",&t);
	while(t--)
	{
		int l=1e9,r=0,pd=0,sum=0;
		scanf("%lld",&n);
		for(int i=1;i<=n;i++) 
		{
			scanf("%lld%lld%lld",&s[i],&e[i],&d[i]);
			l=min(l,s[i]);r=max(r,e[i]);//最大的左右端点 
			pd+=(e[i]-s[i])/d[i]+1;
		}
		if(pd%2==0) 
		{
			printf("Poor QIN Teng:(\n");
			continue ;
		}
		while(l<=r) 
		{
			int mid=(l+r)/2;
			if(check(mid)) ans=mid,r=mid-1;
			else l=mid+1;
		}
		for(int i=1;i<=n;i++)
		{
			if(ans<s[i]||ans>e[i]) continue ;
			if((ans-s[i])%d[i]==0) sum++;//若这个点刚好在s[i]上或者d[i]+s[i]上 
		}
		printf("%lld %lld\n",ans,sum);
	}
	return 0;
}

wqs二分学习笔记：
https://www.luogu.com.cn/blog/Flying2018/wqs-er-fen-min-ke-fu-si-ji-hu-xue-xi-bi-ji
https://www.luogu.com.cn/blog/daniu/wqs-er-fen
对于一个形如二次函数求最大值的东西，即其拥有单调性，我们可以使用wqs二分 二分斜率。
具体而言，对于一些给定的数值，在一定限制条件下，求解最大值。P5633 最小度限制生成树，这题的问题在于给定了一些边，并要求其中有一点有且仅有k条边，明显可以看出其有单调性。不妨尝试转化成wqs，最后再处理一些不合法情况即可。

#include<bits/stdc++.h>
#define int long long 
using namespace std;
int n,m,fa[1000001],len=0,s,k,check1=0;
struct pp
{
	int x,y,w,s;
};pp p[1000001];
int findfa(int x)
{
	if(fa[x]==x) return x;
	return fa[x]=findfa(fa[x]);
}
void ins(int x,int y,int w)
{
	int ss=0,now=++len;if(x==s||y==s) ss=1;
	p[now]={x,y,w,ss};check1+=ss;
	return ;
}
bool merge(int x,int y)
{
	int xx=findfa(x),yy=findfa(y);if(xx==yy) return false;
	fa[xx]=yy;return true;
}
bool cmp(const pp &x,const pp &y)
{
	if(x.w==y.w) return x.s>y.s;
	return x.w<y.w;
}
void remake(int val)
{
	for(int i=1;i<=len;i++) if(p[i].s) p[i].w+=val;
	return ;
}
bool check2()
{
	int now=n-1,now2=0;for(int i=1;i<=n;i++) fa[i]=i;remake(1e9);sort(p+1,p+m+1,cmp);
	for(int i=1;i<=m;i++) 
	{
		if(!merge(p[i].x,p[i].y)) continue ;
		now--,now2+=p[i].s;//printf("%d %d %d %d\n",p[i].x,p[i].y,now2,p[i].s);
	}remake(-1e9);//printf("%d %d ",now2,k);
	if((!now)&&now2<=k) return true;
	return false;
}
int check(int val)
{
	remake(val);for(int i=1;i<=n;i++) fa[i]=i;int now=n-1,sumval=0,now2=0;sort(p+1,p+m+1,cmp);
	for(int i=1;i<=m;i++)
	{
		int x=p[i].x,y=p[i].y;if(!merge(x,y)) continue ;
		now--,sumval+=p[i].w;if(x==s||y==s) now2++;
		if(!now) break ;
	}
	remake(-val);
	if(now||now2<k) return -1;return sumval;
}
signed main()
{
	scanf("%lld%lld%lld%lld",&n,&m,&s,&k);
	for(int i=1;i<=m;i++)
	{
		int x,y,w;scanf("%lld%lld%lld",&x,&y,&w);ins(x,y,w);
	}
	if(check1<k||!check2()) 
	{
//		printf("%d %d %d\n",check1,k,check2);
		printf("Impossible\n");return 0;
	}
	int l=-10000000ll,r=10000000ll,ans=-10000000ll;
	if(check(l)==-1) 
	{
		printf("Impossible\n");
		return 0;
	}
	while(l<=r)
	{
		int mid=(l+r)/2;
		if(check(mid)!=-1) ans=max(ans,mid),l=mid+1;
		else r=mid-1;
	}
	int sum=check(ans);
	printf("%lld",sum-ans*k);
	return 0;
}

P2619 [国家集训队] Tree I,这题是上面的简化版，同样也是使用wqs的算法。

#include<bits/stdc++.h>
using namespace std;
int n,m,need,fa[100001];
struct pp
{
	int x,y,c,color;
};pp p[1000001];
int findfa(int x)
{
	if(fa[x]==x) return x;
	return fa[x]=findfa(fa[x]);
}
int check()
{
	int tot=0,num=0,rt=0;
	for(int i=1;i<=m;i++)
	{
		int x=p[i].x,y=p[i].y,fx=findfa(x),fy=findfa(y);
		if(fx!=fy) rt+=p[i].c,tot+=(p[i].color^1),fa[fx]=fy,num++;
		if(num==n-1)break;
	}
	if(tot>=need) return rt;
	return 0;
}
bool cmp(const pp &x,const pp &y)
{
	if(x.c==y.c) return x.color<y.color;
	return x.c<y.c;	
}
void remake(int val)
{
	for(int i=0;i<=n;i++) fa[i]=i;
	for(int i=1;i<=m;i++) if(!p[i].color) p[i].c+=val;
	sort(p+1,p+m+1,cmp);
	return ;
}
int main()
{
	scanf("%d%d%d",&n,&m,&need);
	int l=-111,r=111,ans;
	for(int i=1;i<=m;i++) scanf("%d%d%d%d",&p[i].x,&p[i].y,&p[i].c,&p[i].color);//l-=min(p[i].c,0),r+=max(0,p[i].c);
	while(l<=r)
	{
		int mid=(l+r)/2;
		remake(mid);
		if(check()) ans=mid,l=mid+1;
		else r=mid-1;
		remake(-mid);
	}
	remake(ans);
	printf("%d",check()-ans*need);
	return 0;
}
cn_tigerhan
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
WQS汇编语言试验指导.doc
04-22
WQS汇编语言试验指导。包含9个单独的实验指导，方便大家更好的掌握汇编语言
用ansys软件分析换能器代码（二）
10-29
- `voltvalue`、`tcbj`、`ouco`、`sdpd`、`dbj`、`bhw`、`wqs`、`wqx`、`hudu`、`dihw`、`sdps`、`sdpx`、`tchw`、`ys`、`js` 用于控制求解过程中的各种参数，如电压值、接触边界条件、输出控制等。 8. **求解准备...
二分问题题解
SSL_lwz的博客
08-23 123
思路：既然要求两个最大相减总和，那么对于一个答案x而言，每个i对它的贡献为，所以贪心取前i个最大的数，若贡献不为负数，即答案合法，注意精度1e-3，实际需要取到-0.02。
洛谷 P1570 KC喝咖啡【二分答案/最大化平均值模板】
weixin_33726313的博客
07-24 234
题目描述 话说KC和SH在福州的时候常常跑去85°C喝咖啡或者其他的一些什么东西。 这天，KC想要喝一杯咖啡，服务员告诉他，现在有n种调料，这杯咖啡只可以加入其中的m种(当然KC一定会加入m种，不会加入少于m种的调料)，根据加入的调料不同，制成这杯咖啡要用的时间也不同，得到的咖啡的美味度也不同。 KC在得知所有的n种调料后，作为曾经的化竞之神的他，马上就知道了所有调料消耗的时间ci以及调料的美味度...
洛谷 P1570 KC喝咖啡【二分】
韧锋zzZ的Blog
08-01 494
**题目描述 话说KC和SH在福州的时候常常跑去85°C喝咖啡或者其他的一些什么东西。 这天，KC想要喝一杯咖啡，服务员告诉他，现在有n种调料，这杯咖啡只可以加入其中的m种(当然KC一定会加入m种，不会加入少于m种的调料)，根据加入的调料不同，制成这杯咖啡要用的时间也不同，得到的咖啡的美味度也不同。 KC在得知所有的n种调料后，作为曾经的化竞之神的他，马上就知道了所有调料消耗的时间ci以及调...
【斜率优化dp+wqs二分】最小划分
andyc_03的博客
02-13 209
题意 给定一个长度为n的序列，把它划分成m段，每段的和为$w_i$ 并给定一个参数p，
WQS二分
Scar_Halo
12-31 550
WQS二分/带权二分/DP凸优化     模型：有 nnn 个物品，选择每一个都会有相应的权值，需要求出强制选 kkk 个物品时的最大/最小权值和     大致思路：     限制个数的最大/最小权值和可以用 n2n^2n2 的 DPDPDP 求出   &nb...
wqs二分（带权二分）总结
andyc_03的博客
04-10 1283
wqs二分一般用于解决这样一种问题： dp是凸函数，也就是多选一个物品，收获肯定更多 这种问题一般是有限制的，比如限制最多选择m个物品 且如果没有限制，将很容易计算 设fuf_ufu​表示选 u 个的收获，图像如下： 我们要找到这样一个(m,fm)(m,f_m)(m,fm​) 可以二分斜率 k ，得到斜率 k 的直线与图像相切的点，这个点相当于是截距最大的点，截距gu=fu−k∗ug_u=f_u-k*ugu​=fu​−k∗u，相当于每个物品都减少 k 的价值，得到的最大的贡献点对应的x坐标，如果x大了意味着
【学习笔记】WQS二分详解及常见理解误区解释
04-08 8164
绝对是良心的博客！新手入门绝对懂！！
tree【WQS二分+MST】
既然弱小，就只顾变强就是了
09-22 352
题目链接——洛谷（精确涉及到了WQS二分） BZOJ-2654（不推荐） 个人不推荐做BZOJ2654的这道题，因为那道题可以水过去，不用WQS二分也是可以的，可以直接二分答案，显然是没有这个好的。 先在这里讲一下什么是WQS二分吧，也是从网上看来的，一开始在做这道题的时候，想到的也是存在这种可能性，但是依然WA了几次： 先说题意：给你一个N个点M条边无向带权连通图，每条边是黑...
wqs.zip_保存文件
09-20
标题中的“wqs.zip_保存文件”表明这是一个与文件保存相关的压缩包，可能包含了某个程序或脚本的执行结果，该结果是计算和显示完全数的过程，并最终被保存到了一个文本文件中。描述中提到的任务涉及到计算数学概念...
通达信指标公式源码神奇分时主图可以精确.doc
07-30
在后面部分，我们可以看到更多的指标公式，如OPEN1、CLOSE1、DFZD、WQS1等，这些公式都是基于股票价格和交易量的计算，用于分析股票的买卖信号和趋势。 这个源码提供了一个完整的技术分析系统，包括各种指标公式和...
includeReplacer:wqs.jd.com域名下的页面片替换神器
05-08
这是一个简单的chrome插件，针对的是wqs.jd.com域名下的前端页面开发。因为经常会有页面片内嵌，本地预览经常会出问题，所以我制作了一个chrome小插件来解决本地预览的问题。如果有需要的同学可以拿过去针对自己的...
P4383 [八省联考2018]林克卡特树（树形dp+wqs二分）
Fighting_Peter的博客
02-28 303
[八省联考2018]林克卡特树 题目大意：给定一棵有负权边的树，现在必须恰好删去kkk条边，并加上恰好kkk条权值为000的边，要求最大化它的直径长度。 首先考虑删去KKK条边的效果：把整棵树变成k+1k+1k+1个连通块 然后用0权变把这几个连通块连起来，当k+1k+1k+1个连通块已经确定后有个显然的结论：用kkk条边把他们连接后的最大直径长度就是k+1k+1k+1个连通块的直径代数求和。 方法是将k+1k+1k+1条直径首位依次相连即可，显然不可能有比它更大的直径。这个方案引导我们转化问题：寻找k
排序算法(4)之快速排序(1)
wer24_25的博客
07-17 1093
排序算法(4)之快速排序(1),快速排序的概念和快速排序的三种实现方式.
(C++二叉树02) 翻转二叉树 对称二叉树 二叉树的深度
qq_44962949的博客
07-17 200
【代码】(C++二叉树02) 翻转二叉树 对称二叉树 二叉树的深度。
Leetcode—207. 课程表【中等】
最新发布
Atticus的博客
07-19 244
运行结果 之后我会持续更新，如果喜欢我的文章，请记得一键三连哦，点赞关注收藏，你的每一个赞每一份关注每一次收藏都将是我前进路上的无限动力 ！！！↖(▔▽▔)↗感谢支持！
WQS 二分是什么意思
07-28
WQS 二分是指在算法和数据结构中使用二分查找的方法。二分查找是一种高效的查找算法，它通过将查找区域分成两半来缩小查找范围，直到找到目标元素或确定目标元素不存在为止。它适用于有序数组或有序列表的查找，每次...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值