Milk Measuring

题意：摘自NOCOW翻译（http://www.nocow.cn/index.php/Translate:USACO/milk4）

描述

农夫约翰要量取 Q（1 <= Q <= 20,000）夸脱（夸脱，quarts，容积单位——译者注） 他的最好的牛奶，并把它装入一个大瓶子中卖出。消费者要多少，他就给多少，从不有任何误差。

农夫约翰总是很节约。他现在在奶牛五金商店购买一些桶，用来从他的巨大的牛奶池中量出 Q 夸脱的牛奶。每个桶的价格一样。你的任务是计算出一个农夫约翰可以购买的最少的桶的集合，使得能够刚好用这些桶量出 Q 夸脱的牛奶。另外，由于农夫约翰必须把这些桶搬回家，对于给出的两个极小桶集合，他会选择“更小的”一个，即：把这两个集合按升序排序，比较第一个桶，选择第一个桶容积较小的一个。如果第一个桶相同，比较第二个桶，也按上面的方法选择。否则继续这样的工作，直到相比较的两个桶不一致为止。例如，集合 {3，5，7，100} 比集合 {3，6，7，8} 要好。

为了量出牛奶，农夫约翰可以从牛奶池把桶装满，然后倒进瓶子。他决不把瓶子里的牛奶倒出来或者把桶里的牛奶倒到别处。用一个容积为 1 夸脱的桶，农夫约翰可以只用这个桶量出所有可能的夸脱数。其它的桶的组合没有这么方便。

计算需要购买的最佳桶集，保证所有的测试数据都至少有一个解。

[编辑]格式

PROGRAM NAME: milk4

INPUT FORMAT

Line 1: 一个整数 Q

Line 2: 一个整数P（1 <= P <= 100），表示商店里桶的数量

Lines 3..P+2: 每行包括一个桶的容积（1 <= 桶的容积 <= 10000）

OUTPUT FORMAT

输出文件只有一行，由空格分开的整数组成：

为了量出想要的夸脱数，需要购买的最少的桶的数量，接着是：

一个排好序的列表（从小到大），表示需要购买的每个桶的容积

[编辑]SAMPLE INPUT (file milk4.in)

16
3
3
5
7

[编辑]SAMPLE OUTPUT (file milk4.out)

2 3 5

解题思路：

1. 用pails[0...P-1]存储桶的容量，然后升序排序
2. DFSID搜索桶的数量，从1到P
3. 在选定桶的数量后，按升序搜索所有的组合情况，搜索到正确情况直接exit(0)
4. 判断选定的组合是否满足条件：使用DP，设f[i] = j代表使用当前桶的组合能否得到i容量的奶（j = 0不行，j = 1行）。设桶的容量集合为seq[0...p]，边界f[0] = 1，那么如果存在j在[0, p]中，使f[i - seq[j]] = 1，那么f[i] = 1，否则为0
5. 除了使用DP，我想DFS应该也能解决。起点是容量Q，每个点的边都有P条，权值为seq[0...p]，每次下降将容量减去相应的边权，搜到0则当前组合满足条件。这样应该比DP要快

代码：

/*
ID: zc.rene1
LANG: C
PROG: milk4
 */

#include<stdio.h>
#include<stdlib.h>
#include<string.h>

#define MAX_P 100
#define MAX_Q 20000

FILE *fout;
int Q, P;
int pails[MAX_P];
int choice[MAX_P];
int choosen[MAX_P];
int f[MAX_Q + 1];

void GetInput(FILE *fin)
{
    int i;

    fscanf(fin, "%d %d", &Q, &P);
    for (i=0; i<P; i++)
    {
	fscanf(fin, "%d", &pails[i]);
    }
}

void SortPails(void)
{
    int i, j, temp;

    for (j=P-1; j>0; j--)
    {
	for (i=0; i<j; i++)
	{
	    if (pails[i] > pails[i + 1])
	    {
		temp = pails[i];
		pails[i] = pails[i + 1];
		pails[i + 1] = temp;
	    }
	}
    }
}

void print_seq(int dest)
{
    int i;

    fprintf(fout, "%d", dest);
    for (i=0; i<dest; i++)
    {
	fprintf(fout, " %d", pails[choice[i]]);
    }
    fprintf(fout, "\n");
    exit(0);
}


void TestChoice(int dest)
{
    int i, j, k;

    memset(f, 0, (MAX_Q + 1) * sizeof(int));
    f[0] = 1;

    for (i=1; i<=Q; i++)
    {
	for (j=0; j<dest; j++)
	{
	    k = i - pails[choice[j]];
	    if (k >= 0 && f[k] == 1)
	    {
		f[i] = 1;
		break;
	    }
	}
    }

    if (f[Q] == 1)
    {
	print_seq(dest);
    }
}

void DFS(int current, int dest)
{
    int i;

    if (current >= dest)
    {
	TestChoice(dest);
	return ;
    }

    if (current == 0)
    {
	i = 0;
    }
    else
    {
	i = choice[current - 1] + 1;
    }

    for (; i<P; i++)
    {
	if (choosen[i] == 0)
	{
	    choice[current] = i;
	    choosen[i] = 1;
	    DFS(current + 1, dest);
	    choosen[i] = 0;
	}
    }
}
void Calculate(void)
{
    int i;

    memset(choosen, 0, MAX_P * sizeof(int));
    for (i=1; i<=P; i++)
    {
	DFS(0, i);
    }
}

int main(void)
{
    FILE *fin;

    fin = fopen("milk4.in", "r");
    fout = fopen("milk4.out", "w");

    GetInput(fin);
    SortPails();
    Calculate();
    
    return 0;
}