二分与三分(备赛中)
创始人
2024-05-29 14:24:02
0

A.愤怒的牛儿

思路:找出最长距离,与最小距离,用二分法判断,如果当前距离满足放牛要求,就把距离区间l=mid+1;如果距离不合适就说明当前距离太大了,把区间变小r=mid-1;最后直到l
#include
using namespace std;
const int maxn=1e9+10;
int num[100010];
int n,m;
int judge(int dis)
{int dum=1;int xi=num[1];for(int i=2;i<=n;i++){if(num[i]-xi>=dis){dum++;xi=num[i];}if(dum==m)break;}return dum==m;}
int main()
{ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);cin>>n>>m;for(int i=1;i<=n;i++)cin>>num[i];sort(num+1,num+1+n);int l=num[1];int r=num[n]-num[1]+1;int k;while(l>1;if(judge(mid)){k=mid;l=mid+1;}else r=mid;}cout<

B.Best Cow Fences

思路:看到了平均值就知道应该要用前缀和了,我们要输出平均数最大的,当然就要考虑用二分了。我们依旧用二分法,每一次如果该平均值满足条件,我们就把区间l=mid;(这里为什么不l=mid+1;因为这样会跳过答案)。如何判断呢就是我们用前缀和,前缀和在开始的时候要把传过来的平均值直接减去,看看是否小于0即可。取最大平均值过程中是这样的,把前面最小的保存下来。然后用后面的减去那个最小值看看是否大于0即可。最好输出即可l或r。

#include
using namespace std;
const int maxn=1e5+10;
int n,m;
double num1[maxn];
double num2[maxn];
bool judge(double ave)
{for(int i=1;i<=n;i++)num2[i]=num2[i-1]+num1[i]-ave;//每次都减去平均值double mins=maxn;//筛出最小前缀值for(int i=m;i<=n;i++)//从最短距离开始{mins=min(num2[i-m],mins);if(num2[i]-mins>0)return true;}return false;
}
int main()
{cin>>n>>m;double l=0,r=0;for(int i=1;i<=n;i++){cin>>num1[i];r=max(num1[i],r);}while(r-l>1e-6){double mid=(l+r)/2;if(judge(mid))l=mid;else r=mid;}printf("%d\n",(int)(1000*r));
}

C.曲线

 思路:单峰函数用三分法。具体方法如下。

#include
using namespace std;
int num;
int num1;
double a[20005],b[20005],c[20005];
double judge(double x)//返回值
{double mas=INT_MIN;for(int i=1;i<=num1;i++){mas=max(mas,a[i]*x*x+b[i]*x+c[i]);//选出最大的那个值}return mas;
}
void solves()
{cin>>num1;for(int i=1;i<=num1;i++){cin>>a[i]>>b[i]>>c[i];}double left=0,right=1000;while(right>left+1e-4)//精度{double midleft=left+(right-left)/3;//三分左区间double midright=right-(right-left)/3;//三分右区间if(judge(midleft)>=judge(midright))//如果左高的话说明左离得远{left=midleft;//移动左区间}else right=midright;//否则移动右区间}printf("%.4f\n",judge(left));}
int main()
{ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);cin>>num;while(num--){solves();}
}

D.数列分段

 思路:依旧是二分思路,取全和做R,取1做L,移动区间即可。对于满足条件就把区间r=mid,不满足条件就l=mid+1;最后输出r。判断函数写法具体是,看看是否能在段数不大于m的情况,满足条件。

#include
using namespace std;
typedef long long ll;
const int maxn=1e5+10;
int n,m;
ll num[maxn];
int check(int x)
{ll sum=0;int cnt=1;for(int i=1;i<=n;i++){if(num[i]>x)return 0;sum+=num[i];if(sum>x){sum=num[i];cnt++;}}return m>=cnt;
}
int main()
{cin>>n>>m;int l=1,r=0;for(int i=1;i<=n;i++){cin>>num[i];r+=num[i];}while(l>1;if(check(mid))r=mid;else l=mid+1;}cout<

E.扩散

 思路:依旧是二分,但需要并查集的协助,我们设置最大时间R=1e9+10,最小时间L=0.我们开始二分算法,如果当前时间满足条件,就把时间区间R=mid;否则L=mid+1.最后输出R。判断函数具体思路是我们设置二重循环,每一次访问两个点看看能不能在当前时间内是否能扩散相交。具体条件是x,y坐标差值与两倍时间的大小关系。如果小于或等于则可以扩散到,反之则不可以。如果扩散到的话就,放到并查集中,如果没有就不用放。之后看看每个点是否只有一个祖先即可。

#include
using namespace std;
const int maxn=1e6+10;
int x[maxn],y[maxn];
int father[maxn];
int n;
int find(int x)
{if(father[x]==x)return x;else return find(father[x]);
}
bool check(int mid)
{int sum=0;for(int i=1;i<=n;i++)father[i]=i;for(int i=1;i<=n;i++){for(int i1=i+1;i1<=n;i1++){if(abs(x[i]-x[i1])+abs(y[i]-y[i1])<=2*mid){if(find(i)!=find(i1)){father[find(i)]=find(i1);}}}}for(int i=1;i<=n;i++){if(father[i]==i)sum++;if(sum==2)return false;}return true;}
int main()
{ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);cin>>n;for(int i=1;i<=n;i++){cin>>x[i]>>y[i];}int l=0,r=1e9+10;int ans;while(l>1;if(check(mid)){r=mid;}else l=mid+1;}cout<

H.灯泡

 

 思路:根据计算所得影子长度就两种情况

(1)当影子不会射到墙上时,影子长度为

(2)当影子射到了墙上时,影子长度为

看函数我们将x定义为等到人的距离。看函数是一个单峰函数。所以用三分法去做。具体如下

​
#include
using namespace std;
double H,h,D;
double check(double x)
{if(D-x>h*x/(H-h))//影子长度小于人到墙的距离{return h*x/(H-h);//返回第一种情况}else{return D-x+h-(H-h)*(D-x)/x;//反之影子在墙上,返回第二种情况。}
}
void solve()
{cin>>H>>h>>D;double l=0;double r=D;while(r-l>1e-6){double midl=l+(r-l)/3;double midr=r-(r-l)/3;if(check(midl)>=check(midr))//单峰函数移动法{r=midr;}else l=midl;}printf("%.3f\n",check(r));
}
int main()
{int num;cin>>num;while(num--){solve();}
}​

 G.传送带

 思路:这道题有点难度,想了好久,学习了好久才会的。共有三段路程,一段在AB上,一段在CD上,还有一段在路上。很简单看一下,那个最少时间满足单谷函数。所以我们用三分,但未知量有两个所以我们用三分套三分。具体如下

#include
using namespace std;
double Ax,Ay,Bx,By,Cx,Cy,Dx,Dy,p,q,r,ans;
double check(double x,double y)
{double P1x=Ax+(Bx-Ax)*x,P1y=Ay+(By-Ay)*x,P2x=Cx+(Dx-Cx)*y,P2y=Cy+(Dy-Cy)*y;return sqrt((P1x-Ax)*(P1x-Ax)+(P1y-Ay)*(P1y-Ay))/p+sqrt((Dx-P2x)*(Dx-P2x)+(Dy-P2y)*(Dy-P2y))/q+sqrt((P2x-P1x)*(P2x-P1x)+(P2y-P1y)*(P2y-P1y))/r;
}
double find2(double x,double l,double r)//在CD上三分
{double res=0.0;//返回的时间数while(r-l>=1e-6)//精度{double midl=l+(r-l)/3;double midr=r-(r-l)/3;double res1=check(x,midl);//按在AB中的比例,和CD中的比例返回时间double res2=check(x,midr);if(res1=1e-6){double midl=l+(r-l)/3;//在AB上三分double midr=r-(r-l)/3;double res1=find2(midl,0,1);//用三分出来的比例在CD上看返回来的值分布情况double res2=find2(midr,0,1);if(res1>Ax>>Ay>>Bx>>By>>Cx>>Cy>>Dx>>Dy>>p>>q>>r;//输入find1(0,1);//这里的0,1是在AB上的比例。printf("%.2f",ans);
}

相关内容

热门资讯

今天,宜爬“山”! 转自:珠海发布大家都说,父爱如山,母爱又何尝不是!它以山的姿态站在我们的人生坐标里,它用山的沉默垒出...
汽车车门把手将迎来国家标准,释... 转自:新华网  新华社北京5月10日电 题:汽车车门把手将迎来国家标准,释放什么信号?  新华社记者...
10日付费短剧热度大盘5650... 来源:新腕儿今天发布的是昨天(5月10日)短剧付费大盘日榜,整体大盘热度值为5650万。特别声明:以...
漂泊 53 年,苏联金星探测器... IT之家 5 月 11 日消息,俄罗斯国家航天集团宣布,近地空间危险情况自动预警系统监测显示,宇宙-...
加州州长再就关税批评特朗普政府 转自:海客新闻 【#加州州长再就关税批评特朗普政府#】#...
都有数|玩儿着就能挣钱?国家认... 转自:北京日报客户端人社部日前公示拟新增无人机群飞行规划员、装修管家等17个新职业,以及黄金鉴定估价...
我爱你,不止今天! 咿呀学语时当我第一次模糊地喊出“妈妈”你尚且青涩的眼眸中亮起柔软的光长大离家后一声“妈妈”成了电话里...
焦点访谈 | 突围拓展多元市场... 来源:央视网 央视网消息(焦点访谈):当前我国的外贸行业正面临前所未有的挑战。美国肆意挥舞关税大棒,...
雨水按下暂停键 阳光重返天气舞...     昨日19时许,南宁城区迎来降雨,体感稍凉。    本报记者宋延康 摄  本报讯(记者赵金玲)...
科普大篷车进乡村 科技体验零距...     科技辅导员为学生讲授科技知识。    本报通讯员肖重虎 摄  本报讯(记者凌剑伊 通讯员蒙蓓...
镇江人才吸引力升至全国37位 金山网讯 5月9日,智联招聘旗下专业研究机构智联研究院联合“泽平宏观”正式发布《中国城市人才吸引力排...
【加快落实“五大任务” 推动高... 来源:内蒙古新闻网-《内蒙古日报  5月7日,在位于乌兰察布市四子王旗忽鸡图乡小东营村灌木林退化修复...
时政微纪录丨习主席的俄罗斯时间 来源:央视新闻客户端 应俄罗斯联邦总统普京邀请,国家主席...
常见故障100%识别 我国货运... 记者今天(11日)从国家能源集团获悉,我国货运铁路首套智能巡检机器人在河北沧州投用,这对促进我国铁路...
中共太原市委组织部公示 现对市委研究拟任职干部予以公示。郭强,男,汉族,1981年10月生,在职大学,中共党员,现任太原物产...
苏丹一州首府遭无人机袭击,20... 苏丹政府当地时间10日说,苏丹北科尔多凡州首府欧拜伊德市当天遭无人机袭击,造成20名囚犯死亡。苏丹文...
陈钧书画金石题跋艺术作品展在铜...   本报讯 (记者 吴洋凯)5月7日,由铜川书画院主办,铜川书画公园文化艺术有限公司承办的“砚耕四十...
花季少女被限制高消费,原因竟然... 转自:上观新闻2024年的夏天,大一学生小赵攥着身份证站在高铁站售票厅,屏幕上“无法购买高铁票”的提...
新职业新工种,开启无限可能! 转自:新华网 早上叫醒你的不是闹钟,可能是睡眠健康管理师...
政策一经出台,便在金融市场激起... 2025 年 5 月 7 日,中国人民银行今日向市场释放重磅政策组合拳。潘功胜行长在国务院新闻办发布...