二分法51nod分割绳子

现在有N(1 <= N <= 1000)条绳子,
他们的长度分别为L1,L2,……,Ln(1 <= Li <= 10000),如果从他们中切割出K(1 <= K <= 1000)条长度相同的绳子,
这K条绳子每条最长能多长?

 收起

输入

共有两行,第一行包含两个正整数N和K,用一个空格分割;第二行包含N个数,一次表示N条绳子的长度,两数间用一个空格分隔,每条绳子的长度的小数不超过两位。

输出

仅包含一个数,表示所得K条绳子的最大长度。答案四舍五入保留小数点后两位

原题的样例似乎出了个问题,样例的精确的答案是2.005,四舍五入的话应该是2.01才对。

输入样例

4 11
8.02 7.43 4.57 5.39

输出样例

2.01

题目大概的意思并不是把绳子截成11段,有可能比11段要多,而必须要有11段相同的绳子,并且尽可能让绳子更长,本题用到了二分法,用来求绳子的长度。代码及解析如下:

#include <cmath>
#include <cstdio>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
double a,d,e; //我把几个变量设在这,因为设在主函数里,a在void里不能调用,所以就都设在这了
int b,k;
const double f=1e-4;//这个是让绳子长短的距离保持在0.0001以内,当然可以更小,但题目要求精度是0.01所以,这些就足够了
double c[100000];
void hs(double d,double e)
{ double mid;
    while(d+f<e)
    {
        mid=(d+e)/2.0;   //设中间值
        int g=0;
        for(int h=1;h<=a;++h)
        {
            g+=c[h]/mid;    //看中间值的绳子可以得到多少根
        }
        if(g>=k) d=mid; //大了这说明中间值小了,还可以更大,所以让小的那部分等于中间值,这样可以提高中间值的大小
        else e=mid;   //同上,减小中间值
    }
     printf("%.2f",e);   //ps!!!!这个地方不能写在主函数里,否则会输出e,这个函数里的e与主函数不一样!!

}
int main() {
cin>>a>>k;
for(b=1;b<=a;++b)
{
    cin>>c[b];
    d=min(d,c[b]);  //去最小,这时因为b放在了主函数外,顾b一直是0
    e=max(e,c[b]);  //取最大
}
hs(d,e);

    return 0;
}
 

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
二分法是一种用于求解函数零点的常见算法,它的基本思想是将一个区间逐步分割为两个子区间,直到最后找到一个足够小的区间,使得函数在该区间内的函数值与零点的误差小于给定的容忍度。在Matlab中,我们可以用递归的方法实现二分法分割区间,具体步骤如下: 1. 定义函数。首先要定义一个用于求解函数零点的函数,可以用Matlab自带的函数fzero(),也可以自己编写一个函数。 2. 定义区间。在主函数中,定义需要求解零点的区间[a,b]。 3. 分割区间。用while循环依次对区间进行分割,直到满足给定的精度要求。 4. 判断根的位置。在每个子区间中,判断根的位置,如果根在左边区间,则进入左边区间继续分割,否则进入右边区间分割。 5. 输出结果。当找到根之后,输出求解结果。 二分法分割区间的Matlab代码如下: function [x,fval,exitflag,output] = bisection(funfcn,a,b,options,varargin) if nargin < 4 || isempty(options), options = optimset; end if any(strcmp(fieldnames(options),'TolFun')) && strcmp(options.TolFun,'eps') options = optimset(options,'TolFun',1000*eps); end if (b-a)<options.TolX, x=a; fval=feval(funfcn,x,varargin{:}); exitflag=1; output=[]; return; end fa = feval(funfcn,a,varargin{:}); fb = feval(funfcn,b,varargin{:}); while (b-a)/2 > options.TolX c = (a+b)/2; fc = feval(funfcn,c,varargin{:}); if fc==0, x=c; fval=fc; exitflag=1; output=[]; return; end if sign(fc)==sign(fa), a = c; fa = fc; else b = c; fb = fc; end end x = (a+b)/2; fval = feval(funfcn,x,varargin{:}); exitflag = 1; output.iterations = NaN; output.funcCount = NaN; output.algorithm = 'bisection';

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值