连续邮资问题

/*   
连续邮资问题   
算法设计:   
该问题是设计最佳的邮票面值，用来表示最大的区间   
对于连续邮姿问题，用n元组x[1:n]表示n种不同的邮票面值并约定它们从小到大排列。   
整数r表示当前使用不超过m张邮票能贴出的最大连续邮资区间。 
x[1] = 1是唯一的选择。此时最大连续邮资区间是[1:m]   
接下来x[2]的可能取值范围是[2:m+1]，这个取值范围就决定了x[2]应该怎么选，   
只所以确定这样的范围是因为，至少要保证增加一个面值后，可以将区间增   
大1一般情况下，已选定x[1:i-1]，此时最大连续邮资区间是[1:r]，则接下来   
x[i]的可取值范围是[x[i-1]+1:r+1]。由此可以看出，在用回溯法时，可以用一棵树表示其   
解空间。该解空间树中各结点的度随x的不同取值而变化。   
   
设x[1～n]是从小到大排列的n种邮票面值。    
首先x[1]=1 是毫无疑问的，否则最小邮资 1 无法支付。    
面值为1的邮票构成的最大连续区间只能是[1,m] (只贴1张就是1……m张邮票全贴满是m)。   
如果这时候再给我们第二种面额的邮票。因为它需要比x[1]大，所以最少只能是2；   
又因为单纯由x[1]支付的邮资区间最大只到m，再往上就出现了断档，所以x[2]最大也只能取m+1   
（如果再大的话，m+1这个值就无法正常表示）    
同样的道理，假设前面k种邮票面值都已经有了，并且能构成[1, r]的连续邮资区间，   
那么第k+1种邮票的面值必须满足：    
x[k]+1 <= x[k+1] <= r+1    
算法就是要找到这么多种可能情况下的最优方案   
   
对于连续邮资问题，用n元组x[1：n]表示n种不同的邮票面值，并约定它们从小到大排列。   
x[1] = 1是惟一的选择。在未确定剩余其它n－1种邮票面值的情况下，只用x[1]一种邮票面值，   
所能得到的最大连续邮资区间是[1：m]。接下来，确定x[2]的取值范围，很明显，   
x[2]的值最小可以取到2，因为x[1] 这时已经为1了。那么最大可以取的值呢？应当是m+1，   
为什么？因为如果x[2]取m+2，则在这个方案中，邮资m+1不可能取得，   
（这个时候，x[2：n]的任何一个面值是不能取的，因为就算只取一张，总和也会至少为m+2，   
超过m+1了，反过来，如果只取x[1]面值，则就算把m张全取了，也只能凑到m）。所以，   
在搜索第2张邮票的时候，搜索范围是2 ~ m+2（m+2不剪取）。在一般的情况下，已选定x[1：i-1]，   
最大连续邮资区间是[1：r]时，接下来x[i]的可取值范围是[x[i-1]+1 ：r+1]。由此可以看出，   
[x[i-1]+1 ：r+1]（为限界函数）在用回溯法解连续邮资问题时，可用树表示其解空间，   
该解空间树中各结点的度随x的不同取值而变化。   
   
现在可以大致画出该解空间树的结构。   
现在的关键问题是，如何确定x3的取值范围？转换一下就是如何求得在x[1]，x[2]面值确定的情况下，   
用不超过m张邮票，所得到的最大连续邮资区间是多少？如果每一个这样的问题能够确定，   
则解空间树已经构造好了（每一个的儿子结点为r-x[i-1]），我们就可以通过深度优先遍历的方式得到最佳面值设计   
(每到达叶结点见比较更新)。   
   
    现在来看在x[1]，x[2]面值确定的情况下，用不超过m张邮票，所得到的最大连续邮资区间是多少？   
最大连续邮资区间总是以1作为起点的，所以我们用max来表示这个最大值，显然，max至少可以取到m，   
因为即使不用x[2]面值，只用x[1]面值的情况下，所能得到的最大值就已经是m了。现在来看在x[1]，   
x[2]面值确定的情况下，用不超过m张邮票，max+1能不能取到？    
假设在拼凑的过程中，x[2]面值的邮票取t张，则t>=0，t<=（max+1） /x[2] && t<=m   
现在x[2]面值的邮票张数已经确定，原问题转化为另一个子问题，即，在x[1]面值确定的情况下，   
用不超过m－t张邮票，max+1 - t*x[2] 能不能取到？ 。。。。。。   
到这一步就很容易理解了，因为x[1] 面值为1，所以如果要拼凑的值max+1- t*x[2] <= m-t的话，   
则只用取max+1- t*x[2]张就可以拉。相反，如果max+1- t*x[2] > m-t的话，则是不可以满足的。   
   
   
  在下面的回溯法描述中，递归函数Backtrack实现对整个解空间的回溯搜索。   
maxvalue记录当前已经找到的最大连续邮资区间，bestx是相应的当前最优解。   
数组y用来记录当前已经选定的邮票面值x[1:i]能贴出各种邮资所需的最少邮票数。   
   
  也就是说,y[k]是用不超过m张面值为x[1:i]的邮票，贴出邮资k所需的最少邮票张数。   
在函数Backtrack中，   
当i>n时，表示算法已经搜索到一个叶结点，得到一个新的邮票面值设计方案x[1:n]。如果该方案能贴出的   
最大连续邮资区间大于当前已经找到的最大连续邮资区间maxvalue，则更新当前最优值maxvalue和相应的最优解。   
当i <= n时，当前扩展结点z是解空间中的一个内部结点，在该结点处x[1:i-1]   
能贴出的最大最大邮资区间为r-1.因此在结点z处x[i]的可取范围是x[i-1]+1:r，   
从而，结点z有r-x[i-1]个儿子结点。算法对当前扩展结点z的每一个儿子结点，   
以深度优先的方式递归地对相应子树进行搜索   
   
  解空间是多叉树，孩子接点个数是每层都在变化的   

*/   

#include<iostream>  
#include<Windows.h>  
using namespace std;  
   
class Stamp  
{  
friend int MaxStamp(int  ,int  ,int []);  
  
private:  
     // int Bound(int i);  
      void Backtrack(int i,int r);  
      int n;//邮票面值数  
      int m;//每张信封最多允许贴的邮票数  
      int maxvalue;//当前最优值  
      int maxint;//大整数  
      int maxl;//邮资上界  
      int *x;//当前解  
      int *y;//贴出各种邮资所需最少邮票数  
      int *bestx;//当前最优解  
   
};  
  
void Stamp::Backtrack(int i,int r)  
{  
    /*计算X[1:i]的最大连续邮资区间,考虑到直接递归的求解复杂度太高， 
    我们不妨尝试计算用不超过m张面值为x[1:i]的邮票贴出邮资k所需的最少邮票数y[k]。 
    通过y[k]可以很快推出r的值。事实上，y[k]可以通过递推在O(n)时间内解决*/  
	int j,k;
    for(j=0;j<=x[i-2]*(m-1);j++) //x[i-2]*(m-1)是第i-2层循环的一个上限，目的是找到r-1的值   
            if(y[j]<m)  
            {  
                for(k=1;k<=m-y[j];k++) //k是对表示j剩余的票数进行检查     
                {  
                    if(y[j]+k<y[j+x[i-1]*k])  
                //x[i-1]*k是k张邮票能表示的最大邮资     
                //+j表示增加了i邮资后能     
                //判断新增加的能表示的邮资需要多少   
                    {  
                        y[j+x[i-1]*k]=y[j]+k;  
                        //对第i-2层扩展一个x[i-1]后的邮资分布  
                    }  
                    
                }  
            }  
                //查看邮资范围扩大多少，然后查询y数组从而找到r   
                while(y[r]<maxint) //计算X[1:i]的最大连续邮资区间  
                {  
                    r++;  
                }  
                //搜索求出r-1的值，对应x[1:i-1]的在m张限制内的最大区间   
                if(i>n) // 如果到达发行邮票的张数，则更新最终结果值，并返回结果  
                {  
                    if(r-1>maxvalue)  // 用r计算可贴出的连续邮资最大值，而maxStamp存放最终结果  
                    {  
                        maxvalue=r-1;  
                        for(int j=1;j<=n;j++)  
                            bestx[j]=x[j];  // 用x[i]表示当前以确定的第i+1张邮票的面值，bestx保存最终结果  
                    }  
                        return;  
                }  
                int *z=new int[maxl+1];  
                for(k=1;k<=maxl;k++)  
                    z[k]=y[k];  
                //保留数据副本，以便返回上层时候能够恢复数据     
                //以上都是处理第i-1层及其之上的问题     
                for(j=x[i-1]+1;j<=r;j++) //在第i层有这么多的孩子结点供选择     
                {  
                    x[i]=j;  
                    Backtrack(i+1,r);//返回上层恢复信息  
                    for(int k=1;k<=maxl;k++)  
                        y[k]=z[k];  
                }  
                delete[] z;  
}  
  
  
  
int MaxStamp(int n,int m,int bestx[]){  
    Stamp X;  
    int maxint=32767;  
    int i,maxl=1500;  
      
    X.n=n;  
    X.m=m;  
    X.maxvalue=0;  
    X.maxint=maxint;  
    X.maxl=maxl;  
    X.bestx=bestx;  
    X.x=new int [n+1];  
    X.y=new int [maxl+1];  
      
    for(i=0;i<=n;i++)  
        X.x[i]=0;  
    for(i=1;i<=maxl;i++)  
        X.y[i]=maxint;  
    X.x[1]=1;  
    X.y[0]=0;  
    X.Backtrack(2,1);  
    cout<<"当前最优解:";  
    for(i=1;i<=n;i++)  
        cout<<bestx[i]<<"  ";  
    cout<<endl;  
    delete[] X.x;  
    delete [] X.y;  
    return X.maxvalue;  
}  
  
  
int main(){  
  
      
    int *bestx;  
    int n;  
    int m;  
    cout<<"请输入邮票面值数:";  
    cin>>n;  
    cout<<"请输入每张信封最多允许贴的邮票数:";  
    cin>>m;  
  
    bestx=new int[n+1];  
    for(int i=1;i<=n;i++)  
    bestx[i]=0;  
      
    cout<<"最大邮资:"<<MaxStamp(n,m,bestx)<<endl;  
    Sleep(5000);  
      
      return 0;
  
}