POJ 1742 Coins 【多重背包】【可行性问题】

POJ 1742 Coins 

大致意思： 给出硬币面额及每种硬币的个数，求从1到m能凑出面额的个数。  Input 多组数据，每组数据前两个数字为n，m。n表示硬币种类数，m为最大面额，之后前n个数为每种硬币的面额，后n个数为相应每种硬币的个数。 （n<=100,m<=100000，面额<=100000，每种个数<=1000） Output 面额数 Sample Input 3 10 1 2 4 2 1 1 2 5 1 4 2 1 0 0 Sample Output 8 4 --------------------------------------------------- 第一组数据，可以组成面额为1.2.3.4.5.6.7.8的·，都在m=10之内，共8种 第二组数据，可以组成面额为1.2.4.5.6的·，但是6不在m=5之内，满足的共4种

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#define max_n 100010
using namespace std; 
int a[max_n], b[max_n], dp[max_n], vis[max_n];

int main() {
    int n, m, k, ans;
    while(scanf("%d %d", &n, &m) && (n + m)) {
        for(int i = 0; i < n; i++)
            scanf("%d", &a[i]);
        for(int i = 0; i < n; i++)
            scanf("%d", &b[i]);
        memset(vis, 0, sizeof(vis));
        ans = 0,vis[0] = 1;
        for(int i = 0; i < n; i++) {
            memset(dp, 0, sizeof(dp));
            for(int j = a[i]; j <= m; j++) {
                if(!vis[j] && vis[j - a[i]] && dp[j - a[i]] + 1 <= b[i]) {
                    vis[j] = 1;
                    dp[j] = dp[j - a[i]] + 1; //表示a[i]的使用次数
                    ans++;
                }
            }
        }
        printf("%d\n", ans);
    }
    return 0;
}