洛谷P3092 [USACO13NOV]No Change G （状压dp+二分）

传送门

题意

约翰到商场购物，他的钱包里有K(1 <= K <= 16)个硬币，面值的范围是1…100,000,000。

约翰想按顺序买 N个物品(1 <= N <= 100,000)，第i个物品需要花费c(i)块钱，(1 <= c(i) <= 10,000)。

在依次进行的购买N个物品的过程中，约翰可以随时停下来付款，每次付款只用一个硬币，支付购买的内容是从上一次支付后开始到现在的这些所有物品（前提是该硬币足以支付这些物品的费用）。不幸的是，商场的收银机坏了，如果约翰支付的硬币面值大于所需的费用，他不会得到任何找零。

请计算出在购买完N个物品后，约翰最多剩下多少钱。如果无法完成购买，输出-1。

分析

注意到k的范围，考虑状压。
f[i]表示当前状态能到购买到“最远” 的商品编号。
于是可以考虑用二分优化查找。
总体还是比较直观的。

代码

#include <bits/stdc++.h>

using namespace std;
//-----pre_def----
const double PI = acos(-1.0);
const int INF = 0x3f3f3f3f;
typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int, int> PII;
typedef pair<double, double> PDD;
#define fir(i, a, b) for (int i = (a); i <= (b); i++)
#define rif(i, a, b) for (int i = (a); i >= (b); i--)
#define endl '\n'
#define init_h memset(h, -1, sizeof h), idx = 0;
#define lowbit(x) x &(-x)

//---------------
const int N = 1e5 + 10;
int n, k;
LL c[20], d[N];
LL s[N];
int f[1 << 17]; //当前状态到了第几个物品
int func(int l, int r, int v)//找从第l+1个商品开始，面值为v的硬币最多能购买几个商品
{
    int temp = l;
    int tmp = s[l], res;
    while (l <= r)
    {
        int mid = l + r >> 1;
        if (s[mid] - tmp <= v)
            res = mid, l = mid + 1;
        else
            r = mid - 1;
    }
    return res - temp;
}

int main()
{
#ifndef ONLINE_JUDGE
    freopen("in.txt", "r", stdin);
    freopen("out.txt", "w", stdout);
    int StartTime = clock();
#endif
    scanf("%d%d", &k, &n);
    fir(i, 0, k - 1)
    {
        scanf("%lld", &c[i]);
    }
    fir(i, 1, n)
    {
        scanf("%lld", &d[i]);
        s[i] = s[i - 1] + d[i];
    }
    LL ans = 0;
    bool flag = 0;
    for (int i = 1; i < (1 << k); i++)
    {
        for (int j = 0; j < k; j++)
        {
            if ((i >> j) & 1)
            {
                f[i] = max(f[i], f[i - (1 << j)] + func(f[i - (1 << j)], n, c[j]));//这里把f[i - (1 << j)]敲成了f[i - (i << j)]一直re，调了好久才找到。。以至于代码写的很乱
            }
        }
        if (f[i] == n)
        {
            flag = 1;
            LL res = 0;
            for (int j = 0; j < k; j++)
            {
                if ((i >> j) & 1)
                {
                    continue;
                }
                res += c[j];
            }
            ans = max(ans, res);
        }
    }
    if (!flag)
        puts("-1");
    else
        cout << ans << endl;
#ifndef ONLINE_JUDGE
    printf("Run_Time = %d ms\n", clock() - StartTime);
#endif
    return 0;
}