题意:摘自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
解题思路:
- 用pails[0...P-1]存储桶的容量,然后升序排序
- DFSID搜索桶的数量,从1到P
- 在选定桶的数量后,按升序搜索所有的组合情况,搜索到正确情况直接exit(0)
- 判断选定的组合是否满足条件:使用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
- 除了使用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;
}