[POI 2007]Weights(拆位贪心)

题目链接

http://main.edu.pl/en/archive/oi/14/odw

题目大意

转自BZOJ
在byteotian公司搬家的时候，他们发现他们的大量的精密砝码的搬运是一件恼人的工作。公司有一些固定容量的容器可以装这些砝码。他们想装尽量多的砝码以便搬运，并且丢弃剩下的砝码。每个容器可以装的砝码数量有限制，但是他们能够装的总重量不能超过每个容器的限制。一个容器也可以不装任何东西。任何两个砝码都有一个特征，他们的中总有一个的重量是另外一个的整数倍，当然他们也可能相等。

思路

假设不同重量的砝码共有$top$种：$stack_1,stack_2...stack_{top}(stack_1<stack_2<...<stack_{top})$
然后我们可以把每个容器体积表示成这种奇怪进制的数字的形式

$$digit_{top}stack_{top}+...+digit_2stack_2+digit_1stack_1+t$$
这个t是对这个数字分解后得到的一个小数字，因为太小了，不够放任何一个砝码，用不上，后面直接无视掉它

我们把每个容器体积的奇怪进制的每一位加起来(t就不要算进来了)，得到一个新数字

比如说现在得到的这样一个奇怪进制的数字是2301
我们首先满足重量为$stack_1$的砝码，假设这种砝码只有一个，去掉以后这个奇怪进制数字就是2300
然后满足重量为$stack_2$的砝码，由于这一位为0，我们从左边与它最接近且不为0的那一位$i$借一个数位过来，这一位就增加数字$\frac {stack_i}{stack_{2}}$,可以发现这样做的话，结果一定不比不借位的答案差，因为借了一位以后，第i位要么可以满足的砝码个数不变，要么可以满足的砝码个数-1，而借了一位后却多放入了一个砝码，证毕。

直到出现当前位为0，而又找不到更高位来借位时，就说明不能放入新的砝码了。

我们再看为什么这种把每个容器的奇怪进制的每一位分别加起来，再统一做贪心是正确的。

如果当前位不需要借位，那么一定能从$a_1...a_n$中找到若干个$a_i$，使得这些$a_i$的奇怪进制的第$j$位非零，且这些$a_i$的第$j$位加起来是大于等于$stack_j$的个数的。

如果当前位需要从第x位借位，那么一定能从$a_1...a_n$中找到若干个$a_i$，使得这些$a_i$的奇怪进制的第$x$位非零，且这些$a_i$的第$x$位加起来再乘上$\frac{stack_x}{stack_j}$是大于等于$stack_j$的个数的，这样我们在实际上就是从这些$a_i$里分别借位来满足$stack_j$的需求的。

代码

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <algorithm>

#define MAXN 110000

using namespace std;

int n,m;
int a[MAXN],b[MAXN]; //a[i]=每个容器的容量，b[i]=每个物品的体积
int digit[MAXN]; //digit[i]=转换成奇怪进制的数字的第i位
int stack[MAXN],top=0; //保存每种砝码的重量

bool check(int pos) //检查从高位借一位到pos位上是否可行
{
    for(int i=pos+1;i<=top;i++)
        if(digit[i])
        {
            digit[i]--;
            digit[pos]+=stack[i]/stack[pos];
            return true;
        }
    return false;
}

int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);
    for(int i=1;i<=m;i++) scanf("%d",&b[i]);
    sort(b+1,b+m+1);
    for(int i=1;i<=m;i++)
        if(b[i]!=b[i-1]) stack[++top]=b[i];
    for(int i=1;i<=n;i++)
        for(int j=top;j>=1;j--)
        {
            digit[j]+=a[i]/stack[j];
            a[i]%=stack[j];
        }
    int now=1; //当前在从低到高第now位
    for(int i=1;i<=m;i++)
    {
        while(b[i]>stack[now]) now++;
        if(digit[now]) digit[now]--;
        else if(check(now)) digit[now]--;
        else
        {
            printf("%d\n",i-1);
            return 0;
        }
    }
    printf("%d\n",m);
    return 0;
}