POJ2976 Dropping tests(二分+精度问题）

---恢复内容开始---

　　这个题就是大白P144页的一个变形，二分枚举x，对a[i]-x*b[i]从大到小进行排序，选取前n-k个判断和是否大于等于0，若满足条件，增大下限，否则，减小下限。

　　这个题因为精度问题wa了n次，后来干脆把a[i]和b[i]改成double就过了，循环终止条件写成while (abs(ub - lb) > 1e-4)比写成for(int i=0;i<100;i++)要好，既能减少时间消耗，又能满足精度。

　　lb和ub初始化时，ub必须等于1.0,否则会莫名其妙的wa

　　同时注意输出，一开始写成printf("%d\n",int(lb*100))wa了n次，正确的应该写成

printf("%.0f\n",lb*);

　　或者写成

cout << fixed << setprecision() << lb *  << endl;

　　AC代码

#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[]={,,-,},dy[]={,,,-};

const double INF=0x3f3f3f3f;

const ll mod=1e9+;

ll powmod(ll a,ll b) {ll res=;a%=mod;for(;b;b>>=){if(b&)res=res*a%mod;a=a*a%mod;}return res;}

const bool AC=true;

int n,k;

double a[]; //改成int很容易wa

double b[]; //改成int很容易wa

double c[];

bool cmp(int x,int y){

    return x>y;

}

bool C(double x){ //传入的为double

    double sum=;

    rep(i,,n){

        c[i]=a[i]-x*b[i];

    }

    sort(c,c+n,cmp);

    rep(i,,n-k){  //选取n-k个

        sum+=c[i];

    }

    return sum>=;

}

int main(){

    double lb,mid,ub;

    while(in2(n,k)==){

        if(n==&&k==) break;

        rep(i,,n){

            inlf(a[i]);

        }

        rep(i,,n){

            inlf(b[i]);

        }

        lb=;ub=;//开的过大可能会因为循环精度不够会莫名奇妙的wa

        for(int i=;i<;i++){

            mid=(lb+ub)/;

            if(C(mid)) lb=mid; //满足条件,增大下限

            else ub=mid;

        }

        printf("%.0f\n",lb*);//此处不能写成printf("%d\n",int(lb*100));

    }

    return ;

}

---恢复内容结束---