连续邮资问题(回溯+动态规划）

这个程序debug了好长时间....一个晚上都耗上面了

这个程序实际上每一部分并不复杂，但是动态规划那边边界错了好长时间。

题目：假设国家发行了n种不同面值的邮票，并且规定每张信封上最多只允许贴m张邮票。连续邮资问题要求对于给定的n和m的值，给出邮票面值的最佳设计，在1张信封上可贴出从邮资1开始，增量为1的最大连续邮资区间。例如，当n=5和m=4时，面值为(1,3,11,15,32)的5种邮票可以贴出邮资的最大连续邮资区间是1到70。

首先比较困难的这个题目看上去循环是没有上限的，也就是我们看似找不到循环边界，其实我们只需考虑可能生成的节点即可。

dfs(x,cur,max)其中x数组存储当前的解有哪些（如果你就是求具体数值，而不要求解的组成甚至可以省略），cur代表目前到第几种面值，当到n为止，max表示到目前为止的最大可到达邮资，而之后下一个x[cur+1]一定是取x[cur]+1~max+1，这个很好理解，max+1目前是访问不到的，那么我加入x[cur+1]后，要看看max更新到多少，而max是从max+1到

m*x[cur+1]（最多就是m个最大的那个邮票，最少肯定是要更新，不然要你何用...)，接下来的问题就是能否用n种邮票，面值都储存在x数组中,最多贴m张，表示出某个数,

方法是动态规划，状态转移方程dp[i,j]=min(dp[i-1][j-k*a[i]+k),其实有点像背包;（注意边界！！！！！！！)dp[i][0]=0(没有钱时0张邮票)dp[1][i]=i;(因为第一张邮票就是1,贴几块钱就是几张邮票)

话说我的程序跑的真慢，应该还能优化...

附一些数据

n= 5 m=4 {1,3,11,15,32} 最大邮资为70

n=5 m=5  {1,4,9,31,51} 最大邮资为126

n=4 m=2 {1,3,5,6}最大邮资为12

n=3 m=4 {1,5,8}最大邮资为26

n=6 m=4{1,4,9,16,38,49}最大邮资为108

n=5 m=6{1,7,12,43,52}最大邮资为216

详细的内容可以参考代码(代码打了很多注释),我前面转载的一片博客也不错，不过思路应该是差不多吧：

#include <stdio.h>#include <string.h>int n,m;//n为邮票种类,m为一封信上最多贴的邮票个数int Max;int ans[10000];//最终答案数组int min(int a,int b){    return a<b?a:b;}int panduan(int x[10000],int n,int sum)//能否用n种邮票，面值在x数组中,最多贴m张，表示出sum(是个动态规划问题，方法是求出dp[n][sum]看它是否小于sum,状态转移方程dp[i][j]=min(dp[i-1][j-k*x[i]]+k)(其中dp[i][j]表示用到第i种邮票，表示邮资为j的最少邮票 {     int i,j,k;     int dp[15][1005];     for (i=0;i<=n;i++)        dp[i][0]=0;       for (i=0;i<=sum;i++)        dp[1][i]=i;     for (i=2;i<=n;i++)        for (j=1;j<=sum;j++)     {         dp[i][j]=9999;         for (k=0;k<=j/x[i];k++)            dp[i][j]=min(dp[i][j],dp[i-1][j-x[i]*k]+k);     }     if (dp[n][sum]>m)        return 0;        return 1; }void DFS(int x[10000],int cur,int max){    int i,j,next;    if (cur==n)//如果已经得出了n种邮票    {        if (max>Max)//并且它的最大值已经大于当前最大邮资数        {            Max=max;            for (i=1;i<=cur;i++)                ans[i]=x[i];//更新答案数组        }         return;    }        for (next=x[cur]+1;next<=max+1;next++)//如果还没得到n中邮票，那么从x[cur]+1~max+1选一个作为下一个邮资，因为max+1没法表示，所以必定到max+1为止        {          x[cur+1]=next;//接下来是重点，用种类为cur+1,数目分别为x[1..cur+1]的邮票,最多使用m张，能否表示出大于max的某个数             for (i=max+1;i<=m*x[cur+1];i++)//这个数最少要为max+1(不然没有意义了),最多是x[cur+1]*m                if (panduan(x,cur+1,i)==0)//如果成立                break;                if (i>max+1)//如果至少让最大值更新了一次                DFS(x,cur+1,i-1);       }}  int main()  {      int i,j,max,cur;      int x[1000];//中间传递的数组，存储当前的邮票值的解      scanf("%d%d",&n,&m);      Max=0;      max=m;      cur=1;      x[cur]=1;      DFS(x,cur,max);//x存储当前的解,cur表示当前传递到第几种邮票,max表示目前能表示到的最大值      printf("%d\n",Max);      for (i=1;i<=n;i++)      printf("%d ",ans[i]);      return 0;  }