目录
1、一元三次方程求解
#include<iostream>
using namespace std;
double a,b,c,d;
double f(double x){
return a*x*x*x+b*x*x+c*x+d;
}
int main(){
double l,r,m,x1,x2;
int s=0;//统计点的个数
scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
for(int i=-100;i<100;i++){
l=i;
r=i+1;//两点之间找答案
x1=f(l);
x2=f(r);
if(!x1){
printf("%.2lf ",l);
s++;
}
//如果左端点为0了,左端点找到答案了,不能判段右端点回合下面的二分重复
if(x1*x2<0){//中间有解
while(r-l>0.001)//确保精度
{
m=(r+l)/2;
if(f(m)*f(r)<=0){
l=m;
}else{
r=m;
}
}
printf("%.2lf ",l);
s++;
}
if(s==3) break;
}
return 0;
}
2、数的分段
#include<iostream>
#include<algorithm>
using namespace std;
const int N=1e5+10;//最多1e5个数
int n,m,x,f[N],p,sum,maxn;//p是分段的点位
int l,r,mid;
int main(){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++){
scanf("%d",&x);
f[i]=f[i-1]+x;//数组全部组成一段,总和为最大的
maxn=max(maxn,x);//最小为把数组全单独一段,最大的数就是最大的
}
int l=maxn,r=f[n];
while(l<r){
mid=(l+r)>>1;
p=1,sum=1;//当前总共分了1段
for(int i=2;i<=n;i++){
if(f[i]-f[p-1]>mid) p=i,sum++;
//p到i-1的和大于mid说明当前区间和大于要求了,p取到当前位置,但是p是属于下一个区间的
}
if(sum>m) l=mid+1;//搜后半段
else r=mid;
}
cout<<r<<endl;
return 0;
}
3、营救
代码不解释了,用的Kruskal算法,排序一下,只要s到t联通了输出当前的w就是最大的那个路径了
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1e4+10,M=1e5+10;
int n,m,s,t;
int p[N];
struct Edge{
int a,b,w;
}edges[M];
bool cmp(Edge A,Edge B){
return A.w<B.w;
}
int find(int x){
if(p[x]!=x) p[x]=find(p[x]);
return p[x];
}
int main(){
scanf("%d%d%d%d",&n,&m,&s,&t);
for(int i=0;i<m;i++){
int a,b,c;
scanf("%d%d%d",&a,&b,&c);
edges[i]={a,b,c};
}
for(int i=1;i<=n;i++) p[i]=i;
sort(edges,edges+m,cmp);
for(int i=0;i<m;i++){
int a=edges[i].a,b=edges[i].b,w=edges[i].w;
a=find(a),b=find(b);
if(a!=b) p[a]=b;
if(find(s)==find(t)){
cout<<w<<endl;
break;
}
}
return 0;
}
4、KC喝咖啡
#include<iostream>
#include<algorithm>
using namespace std;
const int N=210;
struct coffee{
int c,v;//花费和美味度
double div;
}co[N];
int n,m;
double l,r,mid;
bool cmp(coffee A,coffee B){//按照性价比排序
return A.div>B.div;
}
bool check(double x){
for(int i=1;i<=n;i++){
co[i].div=co[i].v-x*co[i].c;
}//计算当前x下的性价比
sort(co+1,co+1+n,cmp);
double tot=0;
for(int i=0;i<=m;i++){//取前m个
tot+=co[i].div;
}
return tot>=0;
}
int main(){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++) scanf("%d",&co[i].v);
for(int i=1;i<=n;i++) {
scanf("%d",&co[i].c);
r=max(r,co[i].v*1.0/co[i].c);//有边界为最不划算的情况
}
while(r-l>=1e-5){
mid=(l+r)/2;
if(check(mid)) l=mid;//>=说明ans取小了
else r=mid;
}
printf("%.3lf",l);
return 0;
}
5、切绳子
二分边界真的很重要!!!边界磨半天了
#include<iostream>
#include<algorithm>
using namespace std;
const int N=1e4+10;
int n,k;
double c;
int arr[N];
int l,r,mid;
bool check(double x){
int res=0;
for(int i=1;i<=n;i++) res+=arr[i]/x;
return res>=k;
}
int main(){
scanf("%d%d",&n,&k);
for(int i=1;i<=n;i++){
scanf("%lf",&c);
arr[i]=c*100;
r=max(r,arr[i]);
}
while(l<r){
mid=l+r+1>>1;
if(check(mid)) l=mid;
else r=mid-1;
}
printf("%.2lf",r*0.01);
return 0;
}
6、逛画展
这题有必要说下数据的问题,我通常会设到l<r,但是做这题有一个数据卡住了,如果起点到终点才相同的话,判段到最后会直接退出了,并且没有一次进入到ansl和ansr的更新。
所以这里看了下题解,设到l<=r,但是r要取到-1。l要取到+1,防止两个相同进入死循环,还有就是mid这里l+r+1>>1,也是防止死循环的,不过这里既然两个边界不取就没有问题了,加不加都一样。
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
const int N=1e6+10,M=2e3+10;
int n,m;
int arr[N],flag[M];
int ansl,ansr;
bool check(int x){
int cnt=0;
memset(flag,0,sizeof flag);//将记录数组初始化为0
for(int i=1;i<=x;i++){
if(!flag[arr[i]]) cnt++;
flag[arr[i]]++;
}//找到1到x中有多少个不相同数
if(cnt==m){//如果直接相同了更新区间
ansl=1,ansr=x;
return true;
}
int p=1;//滑动窗口,因为x固定,所以左边-1右边+1
while(x<n){
flag[arr[p]]--;
if(!flag[arr[p++]]) cnt--;
if(!flag[arr[++x]]) cnt++;
flag[arr[x]]++;
if(cnt==m){
ansl=p,ansr=x;
return true;
}
}
return false;
}
int main(){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++) scanf("%d",arr+i);
int l=1,r=n;
while(l<=r){
int mid=l+r+1>>1;
if(check(mid)) r=mid-1;
else l=mid+1;
}
cout<<ansl<<" "<<ansr<<endl;
return 0;
}