【ACWing】532. 货币系统（配数学证明）

题目地址：

https://www.acwing.com/problem/content/534/

在网友的国度中共有$n$种不同面额的货币，第$i$种货币的面额为$a[i]$，你可以假设每一种货币都有无穷多张。为了方便，我们把货币种数为$n$、面额数组为$a[1:n]$的货币系统记作$(n,a)$。在一个完善的货币系统中，每一个非负整数的金额$x$都应该可以被表示出，即对每一个非负整数$x$，都存在$n$个非负整数$t[i]$满足$a[i]\times t[i]$的和为$x$。然而，在网友的国度中，货币系统可能是不完善的，即可能存在金额$x$不能被该货币系统表示出。例如在货币系统$n=3,a=[2,5,9]$中，金额$1,3$就无法被表示出来。两个货币系统$(n,a)$和$(m,b)$是等价的，当且仅当对于任意非负整数$x$，它要么均可以被两个货币系统表出，要么不能被其中任何一个表出。现在网友们打算简化一下货币系统。他们希望找到一个货币系统$(m,b)$，满足$(m,b)$与原来的货币系统$(n,a)$等价，且$m$尽可能的小。他们希望你来协助完成这个艰巨的任务：找到最小的$m$。

输入格式：
输入文件的第一行包含一个整数$T$，表示数据的组数。接下来按照如下格式分别给出$T$组数据。每组数据的第一行包含一个正整数$n$。接下来一行包含$n$个由空格隔开的正整数$a[i]$。

输出格式：
输出文件共有$T$行，对于每组数据，输出一行一个正整数，表示所有与$(n,a)$等价的货币系统$(m,b)$中，最小的$m$。

数据范围：
$1\le n\le 100$
$1\le a[i]\le 25000$
$1\le T\le 20$

首先证明与$(n,a)$等价的且$m\le n$的最优解货币系统$(m,b)$，一定满足$b\subset a$（当然最优解也一定满足$m\le n$，因为$(n,a)$本身就是一组解），如若不然，则有$b_i\notin a$，由于等价，则有$b_i={\sum }_j{a}_j={\sum }_j{\sum }_k{c}_{jk}{b}_k,j>1,c_{jk}\ge 0$，由于$a$和$b$都是正整数数组，如果等式右边的$b_i$的系数是$0$，那么$b_i$就可以被别的$b$所表示，显然不是最优解，因为$b_i$可以被去掉；如果等式右边的$b_i$的系数是$1$（最多也只能是$1$了，否则就会造成左边小于右边），那会导致别的$b$的系数必须是$0$（否则也会造成左边小于右边），所以就有$\exists !j,a_j=b_i$（这里$\exists !$表示存在唯一的），矛盾。所以$b\subset a$，即我们只需要把$a$中那些“不取不行”的留下来即可。由于$a[i]$是正的，所以其只能被$a$中比它小的数表示出来。将$a$里的数从小到大排序，依次判断$a[i]$是否能被$a[j],j<i$表示出来，如果能，显然$a[i]$不是必需的；否则$a[j]$一定是必需的（因为如果不取$a[j]$，$a[j]$这个数就无法被表示，这就不等价了）。这个过程可以看成是做完全背包问题。代码如下：

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;

const int N = 110, M = 25010;
int n;
int a[N];
// f[i]表示i这个数字能否被货币系统所表示
bool f[M];

int main() {
  int T;
  cin >> T;
  while (T--) {
    cin >> n;
    for (int i = 1; i <= n; i++) cin >> a[i];
    sort(a + 1, a + n + 1);

	// 每次算完一个case，需要把f清零
    memset(f, 0, sizeof f);
    // 0是可以被表示的，可表示为每个货币都不取
    f[0] = true;

    int res = 0;
    // 从小到大枚举货币
    for (int i = 1; i <= n; i++) {
    	// 如果当前货币无法被比它小的货币表示，那当前货币是必需的，计数加1；
    	// 写法有点像筛素数
        if (!f[a[i]]) res++;
        // 完全背包问题体积要从小到大遍历
        for (int j = a[i]; j <= a[n]; j++)
            f[j] |= f[j - a[i]];
    }

    cout << res << endl;
  }
}

每个case时间复杂度$O(n\log n+n\max a_i)$，空间$O(n)$。