货币系统(AcWing, NOIP 2018提高组)

题目描述:  

题目链接:

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

思路:相当经典的动态规划问题，有一定的难度

一开始我往同余方程ax+by=d,去考虑该同余方程有非负整数解的时候的d的取值范围，然后分析了半天感觉没啥很好的结果。

实际上，本题应该是先考虑以下3个性质:

性质1: a1,a2,...,an 一定都可以被数组b表示出来

性质2: 最优数组b中的任意一个数bi，一定都不能被数组b中的其他数给表示出来

性质3: 最优数组b中的所有元素一定都是属于数组a的。(i.e. 数组b是由数组a中取了若干个元素出来而得到的)

性质1的证明: 由于ai可以被数组a表示出来，而数组b和数组a是等价的，因此性质1成立

性质2的证明: 反证法，若数组b中的某一个数bi可以被数组b中其他的数给表示出来，那么bi就没有存在的必要了，因为任何用到bi的地方，我可以用其他数给他替代掉。那这样，把bi去掉，数组b仍然满足要求，与数组b的最优性矛盾！ 

性质3的证明: 这是本题的关键，很有意思的结论。反证法，假设某个bi不属于数组a，首先由于数组b和数组a的等价性，且bi可以被数组b表示出来，可知: 数字bi可以被数组a表示出来。i.e.

bi=ai1*t1+ai2*t2+...aik*tk， 其中ti是正整数。 由于bi不属于数组a，那么必然有: bi> ai1, ai2,...,aik

注意到，数字ai1, ai2,...,aik可由数组b表示，且aij的表式中一定不含bi(因为bi比aij大)，故，aij可以用数组b中除了bi以外的其他数表示。 进一步可知: bi可以用数组b中除了bi以外的其他数表示。这与性质2矛盾！

故bi一定属于数组a，性质3得证。

有了以上3个性质，本题就比较容易了。

完全背包问题，用f[j]表示数字j是否能被表示，然后从小到大遍历数组a。对于当前的ai，如果f[a[i]]=1,那就说明数字a[i]可以用之前的数表示出来，那a[i]就没必要选进数组b了；

如果f[a[i]]=0, 说明数组a[i] 不能用前面的数表示出来，a[i]必须要选进数组b(否则，数组b无法表示a[i])

总的时间复杂度O(nT*25000)

 代码:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int N=25010;
int f[N];
int a[N];
int n;
int t;


int main(void){
	
	scanf("%d",&t);
	while(t--){
		scanf("%d",&n);
		int maxnum=0;
		for(int i=1;i<=n;i++){
			scanf("%d",&a[i]);
			maxnum=max(maxnum,a[i]);
		} 
		sort(a+1,a+1+n);
		memset(f,0,sizeof(f));
		int ans=0;
		f[0]=1;
		for(int i=1;i<=n;i++){
			if(f[a[i]]) continue;
			ans++;
			for(int j=a[i];j<=maxnum;j++){
				f[j]=f[j]|f[j-a[i]];
			}
		}
		printf("%d\n",ans);
	}
	
	return 0;
}