UVa 165 - Stamps

/*
    连续邮资问题  参考：王晓东《计算机算法设计与分析》
    方法一:直接枚举(h+k<=9)
*/

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
const int MAX_H = 9+5; // 最大邮票数
const int MAX_K = 9+5; // 最大面额数
const int MAX_V = 200; // 可组成面额最大值
int A[MAX_H]; //当前解
int h, k; //h是邮票数,k是面额种类, h+k<=9
int max_n; //最大连续面额数
int best[MAX_H]; //最优结果

int vis[MAX_V];
void get_max(int max_idx, int cur, int cur_n)
{
    if(cur == h) return;
    for(int j=0; j<=max_idx; j++) {
        vis[cur_n + A[j]] = 1;
        get_max(max_idx, cur+1, cur_n+A[j]);
    }
}

void dfs(int cur,int cur_n)
{
    if(cur == k) {
        if(cur_n > max_n) {
            max_n = cur_n;
            memcpy(best, A, sizeof(A));
        }
        return;
    }
    for(int i=A[cur-1]+1; i<=cur_n+1; i++) {
        A[cur] = i;
        memset(vis, 0, sizeof(vis));
        get_max(cur, 0, 0);
        int j;
        for(j=cur_n; j<MAX_V; j++)  if(!vis[j]) {
            break;
        }
        dfs(cur+1, j-1);
    }
}

int main(void) {
    #ifndef ONLINE_JUDGE
    freopen("in.txt", "r", stdin);
    freopen("out.txt", "w", stdout);
    #endif

    while(scanf("%d%d", &h, &k) == 2) {
        if(h==0 && k==0) break;
        A[0] = 1;
        max_n = 0;
        dfs(1, h);
        for(int i=0; i<k; i++) {
            printf("%3d", best[i]);
        }
        printf(" ->");
        printf("%3d\n", max_n);
    }
    return 0;
}

/*
    方法二：动态规划过程，用数组记录组成每种面额最少所需邮票数
    假设输入h=3，k=2, 初始化k=1，使用￥1 ：y[0]=0 y[1]=1 y[2]=2 y[3]=3 ...（y[3]=3表示生成面额3需要3张￥1）
    深度遍历下一层k=2,加入￥2，更新y数组，y[2]=1，y[4]=2, y[6]=3 .. （y[2]=1因为面额2的只要1张￥2）
*/
#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
const int MAX_H = 9+5; // 最大邮票数
const int MAX_K = 9+5; // 最大面额数
const int MAX_V = 200; // 可组成面额最大值
int A[MAX_H]; //当前解
int h, k; //h是邮票数,k是面额种类, h+k<=9
int max_n; //最大连续面额数
int best[MAX_H]; //最优结果
int y[MAX_V]; // y[i]表示目前组成面额i所需要最小邮票数

void dfs(int cur,int cur_n)
{
    if(cur == k) {
        if(cur_n > max_n) {
            max_n = cur_n;
            memcpy(best, A, sizeof(A));
        }
        return;
    }
    int buff_y[MAX_V];
    for(int i=A[cur-1]+1; i<=cur_n+1; i++) {
        A[cur] = i;
        memcpy(buff_y, y, sizeof(y));
        for(int j=0; j<=A[cur-1]*h; j++) if(y[j] < h) { //y[j]小于h，说明邮票数没有到达上届，还看以添加
            for(int n=1; n<=h-y[j]; n++) { // n表示添加i的个数
                if(y[j+n*i] > y[j]+n) {
                    y[j+n*i] = y[j]+n; // 目前组成面额j+n*i最少只要y[j]+n张邮票
                }
            }
        }
        int j;
        for(j=1; j<MAX_V; j++) if(y[j] == MAX_H) break;
        dfs(cur+1, j-1);
        memcpy(y, buff_y, sizeof(y));
    }
}

int main(void) {
    #ifndef ONLINE_JUDGE
    freopen("in.txt", "r", stdin);
    freopen("out.txt", "w", stdout);
    #endif

    while(scanf("%d%d", &h, &k) == 2) {
        if(h==0 && k==0) break;
        A[0] = 1;
        max_n = 0;
        int i;
        for(i=0; i<=h; i++) {
            y[i] = i;
        }
        while(i<MAX_V) y[i++] = MAX_H;
        dfs(1, h);
        for(int i=0; i<k; i++) {
            printf("%3d", best[i]);
        }
        printf(" ->");
        printf("%3d\n", max_n);
    }
    return 0;
}