Dropping tests POJ - 2976 ——01分数规划+二分 ——很奇妙的一个题 各个击破

博客目录

一、引言

    好久之前做过这道题,而且是提交了好多次AC代码,前几天又遇到类似的依然懵逼,不光我懵逼,一块做题的acmer也跟着懵逼了,而且大家都做过.......好了被钉在耻辱柱子上了。所以今天又回来补一下博客。

二、题意

vjudge原题传送门

大意:给n个数对(ai,bi),现在删掉不超过k个数对,使得最大。

三、思路

     看上去这道题有点像溶液浓度问题——贪心。但是你用公式化简溶液加入两种新的溶液的密度条件:(a+c)/(b+d) <(a+e)/(b+f)

化简结果为:af+c(b+f)<ad+e(b+d),跟a和b有关(也就是跟以前溶液的情况有关),现在的决策受到以前决策的影响,无法贪心。

     正解为:

设二分变量为ans,每次check  原式能否达到ans, 原式>=ans,也就是

ansa<=100*sigma(ai)/sigma(bi),去分母,然后移项:

ans*sigma(bi)-sigma(ai)<=0;

观察公式,形式为一次函数形式:f(ans)=a*ans+b <= 0的形式,所以可以用对ans二分求解。

而check函数:

将sigma合并

sigma(ans*bi-ai)<=0;

然后对ans*bi-ai排序,有序之后上述公式求和过程便是单调的了,我们check验证此时原式能否满足原式>=ans,所以我们很容易从大到小检查前n-k个求和是否大于等于0,大于等于0就是可以满足条件原式>=ans,所以ans小了应该往上偏移,反之亦然。

时间复杂度:二分ans:log100<7  乘以 [   check时排序了一下nlogn  加上  并列的从大到小遍历n]    结果为O(7nlogn)

四、AC代码

//#include<bits/stdc++.h>
#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<string.h>
using namespace std;
/*
二分
https://vjudge.net/problem/POJ-2976
取前n-k大的 
*/
typedef long long ll;
typedef struct
{
    ll x,y;
}Tone;
Tone co[1010];
ll n,k;
ll s;
double c[1010];
int check(double ans)
{
    int i;
    int p=n-1;
    double sum=0;

    for(i=0;i<n;i++)
        c[i]=co[i].x*100-ans*co[i].y;
    sort(c,c+n);
    for(i=0;i<s;i++)
    {
        sum+=c[p];
        p--;
    }
    return sum>=0;
}
int main()
{
    double m;
    ll i,k;
    double a,b;
    while(cin>>n>>k,n||k)
    {
        s=n-k;
        for(i=0;i<n;i++)
        {
            scanf("%I64d",&co[i].x);
        }
        for(i=0;i<n;i++)
        {
            scanf("%I64d",&co[i].y);
        }
        a=0;
        b=100;
        int x=200;//二分200次精度就差不多够了 
        while(x--)
        {
            m=(a+b)/2;
            if(check(m))
            {
                a=m;
            }
            else b=m;
        }
        printf("%.0f\n",a);
    }
}

猜你喜欢

转载自blog.csdn.net/GreyBtfly/article/details/81385011