EOJ Monthly 2018.2 (Good bye 2017) A2. 坑爹的售票机【完全背包】

题目：http://acm.ecnu.edu.cn/contest/51/problem/A2/

因为所给纸币面值可以用来贪心，所以先用贪心算法得出买n张票的所需最少纸币数。

接着，因为没有规定总共可以买几次票，就是一个完全背包问题。背包容量为买票数目，每种买票方式为一个物品，方式所花纸币数目为物品价值，所买票数为物品重量。

但由于题中所给数据过大，无法创造如此庞大的dp[数组]，故我们采用类似于滚动数组的方法。

题解如下：改题为完全背包的变形题，我们可以求出价值为x需要的张数，从而我们可以知道单次买[1,k]张票需要的张数，到这里应该就是很裸的完全背包了吧，但是一看范围直接懵逼，1e9，直接弃疗，但是这里的k的取值最大为10，所以我们想到，当n很大时，其实我们一次取{1,2,3,···,k}的最小公倍数个数即可，这样我们把剩下的再次求完全背包即可，这里注意，我们取的时候要留一组，来和剩下的做一次完全背包，避免只剩1个或很少的这种极端情况。

OJ上的题解：

k 较小，因此 2520 张票的最优解一定为通过性价比最高的方式买。因此最终答案为 DP 算到 min(n,nmod2520+2520)（每种买票方式需要通过 DP 算的最多不超过 2520 张，若超过 2520 张则这 2520 张可以按性价比最优方式买，因此最多 DP 到 k!k

#include <bits/stdc++.h>


using namespace std;
typedef long long ll;

const ll INF = 99999999999999999LL;

struct {
    ll num;     ///个数
    ll s;       ///总价钱
    ll cnt;     ///张数
}a[11];

const int maxn = 1e6;
ll dp[maxn];

ll f(ll x) {
    ll res = 0;
    if(x >= 100) {
        res += x/100;
        x %= 100;
    }
    if(x >= 50) {
        res++;
        x -= 50;
    }
    while(x >= 20) {
        res++;
        x -= 20;
    }
    if(x >= 10) {
        res++;
        x -= 10;
    }
    if(x >= 5) {
        res++; x -= 5;
    }
    while (x >= 1) {
        res++;
        x -= 1;
    }
    return res;
}

ll gcd(ll a,ll b) {
    return b == 0? a : gcd(b,a%b);
}
int main() {
    ll n,k,p;cin>>n>>k>>p;
    for(int i = 1;i <= k;i++) {
        a[i].num = i;
        a[i].s = i*p;
        a[i].cnt = f(a[i].s);
    }
    ll lcm = 1;//maxLcm = 2520
    for(int i = 2;i <= k;i++) {
        lcm = lcm*i/gcd(lcm,i);
    }
    for(int i = 0;i < maxn;i++) dp[i] = INF;
    dp[0] = 0;
    for(int i = 1; i <= k; i++) {
        for(int j = a[i].num;j <= lcm;j++) {
            dp[j] = min(dp[j],dp[j - a[i].num] + a[i].cnt);
        }
    }
    ll mi = dp[lcm];
    ll res = 0;
    if(n >= 2*lcm) {
        res += (n/lcm - 1)*mi;
        n = n%lcm + lcm;
    }
    for(int i = 0;i < maxn;i++) dp[i] = INF;
    dp[0] = 0;
    for(int i = 1; i <= k; i++) {
        for(int j = a[i].num;j <= n;j++) {
            dp[j] = min(dp[j],dp[j - a[i].num] + a[i].cnt);
        }
    }
    cout<<dp[n]+res<<endl;
    return 0;
}

原帖地址：http://blog.csdn.net/nobleman__/article/details/79323495