POJ 1064 Cable master（二分查找+精度）（神坑题）

Posted 2020-07-02 xgboost

POJ 1064 Cable master 

　　一开始把 int C（double x) 里面写成了  int C（int x) ，莫名奇妙竟然过了样例，交了以后直接就wa。

    后来发现又把二分查找的判断条件写错了，wa了n次，当 c（mid）＜＝ｋ时，令ｕｂ＝ｍｉｄ，这个判断是错的，因为要找到最大切割长度，当满足这个条件时，可能已经不是最大长度了，此时还继续缩小区间，自然就ｗａ了，（从大到小递减，第一次满足这个条件的值，就是最大的值），正确的判断是当 c（mid）＜ｋ时，令ｕｂ＝ｍｉｄ，这样循环１００次时，自然能达到精度，１００次循环能达到的精度是１０^（－３０）；

　　然后就是精度问题了

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <ctime>
#include <iostream>
#include <algorithm>
#include <string>
#include <vector>
#include <deque>
#include <list>
#include <set>
#include <map>
#include <stack>
#include <queue>
#include <numeric>
#include <iomanip>
#include <bitset>
#include <sstream>
#include <fstream>
using namespace std;
#define rep(i,a,n) for (int i=a;i<n;i++)
#define per(i,a,n) for (int i=n-1;i>=a;i--)
#define in(n) scanf("%d",&(n))
#define in2(x1,x2) scanf("%d%d",&(x1),&(x2))
#define inll(n) scanf("%I64d",&(n))
#define inll2(x1,x2) scanf("%I64d%I64d",&(x1),&(x2))
#define inlld(n) scanf("%lld",&(n))
#define inlld2(x1,x2) scanf("%lld%lld",&(x1),&(x2))
#define inf(n) scanf("%f",&(n))
#define inf2(x1,x2) scanf("%f%f",&(x1),&(x2))
#define inlf(n) scanf("%lf",&(n))
#define inlf2(x1,x2) scanf("%lf%lf",&(x1),&(x2))
#define inc(str) scanf("%c",&(str))
#define ins(str) scanf("%s",(str))
#define out(x) printf("%d\n",(x))
#define out2(x1,x2) printf("%d %d\n",(x1),(x2))
#define outf(x) printf("%f\n",(x))
#define outlf(x) printf("%lf\n",(x))
#define outlf2(x1,x2) printf("%lf %lf\n",(x1),(x2));
#define outll(x) printf("%I64d\n",(x))
#define outlld(x) printf("%lld\n",(x))
#define outc(str) printf("%c\n",(str))
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define SZ(x) ((int)(x).size())
#define mem(X,Y) memset(X,Y,sizeof(X));
typedef vector<int> vec;
typedef long long ll;
typedef pair<int,int> P;
const int dx[4]={1,0,-1,0},dy[4]={0,1,0,-1};
const double INF=0x3f3f3f3f;
const ll mod=1e9+7;
ll powmod(ll a,ll b) {ll res=1;a%=mod;for(;b;b>>=1){if(b&1)res=res*a%mod;a=a*a%mod;}return res;}
const bool AC=true;

double a[10005];
int n,k; //k为应该切割的绳子段数
int C(double x){ //计算切割成长度为x的绳子时能得到的绳子的段数(注意x为double)
    int ans=0;
    rep(i,0,n){
        ans+=floor(a[i]/x);
    }
    return ans;
}
int main(){
    in2(n,k);
    rep(i,0,n){
        inlf(a[i]);
    }
    double lb,mid,ub;
    lb=0.0,ub=INF;
    for(int i=0;i<100;i++){ //100次循环可以达到10^(-30)精度
        mid=(lb+ub)/2.0; //二分查找
        if(C(mid)<k) ub=mid; //切割的段数不够，要减小切割长度，此处不是C(mid)<=k
        //(因为要找到第一次等于k的值，此时若等于k,mid可能不是最大的值)
        else lb=mid; //切割段数满足，长度过小，要增大切割长度
    }
    printf("%.2f\n",floor(ub*100)/100);//一定要加floor函数,不加floor函数是四舍五入
    return 0;
}

 

　　由于ｕｂ这个值是稍大于ｘ的，最后用floor(ub*100)/100进行截取，直接截取不加floor函数是四舍五入。