HDU 2062 Subset sequence

Description
Consider the aggregate An= { 1, 2, …, n }. For example, A1={1}, A3={1,2,3}. A subset sequence is defined as a array of a non-empty subset. Sort all the subset sequece of An in lexicography order. Your task is to find the m-th one. 
 
Input
The input contains several test cases. Each test case consists of two numbers n and m ( 0< n<= 20, 0< m<= the total number of the subset sequence of An ).
 
Output
For each test case, you should output the m-th subset sequence of An in one line.
 
Sample Input
1 1
2 1
2 2
2 3
2 4
3 10 
 
Sample Output
1
1
1 2
2
2 1
2 3 1 
 
   Problem Analyse
  考虑一个集合 An = { 1, 2, ..., n}。比如，A1={1}，A3={1,2,3}。我们称一个非空子集元素的排列为一个子集序列。对所有的子序列按字典顺序排序。你的任务就是给出第m个子序列。
Algorithm Analyse
  首先我们来看看An一共有多少个子集。
n=1时,只有{1}一个子集合
n=2时，就有:
{1}, {2}, 
{1, 2}, {2, 1}
4个子集合。


n=3时，有
{1}, {2}, {3}, 
{1, 2}, {1, 3}, {2, 1}, {2, 3}, {3, 1}, {3, 2}, 
{1, 2, 3}, {1, 3, 2}, {2, 1, 3}, {2, 3, 1}, {3, 1, 2}, {3, 2, 1}


也许你发现规律了。An子集合的个数为:
C1n·A11 + C2n·A22 + ... + Cnn·Ann
这个公式是对的。但我们换个角度看。
n=3时，有
{1}
{1, 2}
{1, 2, 3}
{1, 3}
{1, 3, 2}
{2}
{2, 1}
{2, 1, 3}
{2, 3}
{2, 3, 1}


{3}
{3, 1}
{3, 1, 2}
{3, 2}
{3, 2, 1}


不难发现，An可以按首数字分成n组，而每组里除了第一项，剩下的就是An-1的子集合了。
∴f(n) = n[f(n-1) + 1]
f(1) = 1
我们拿测试数据3 10来做个示范，解释一下怎么求解。
因为n=3，所以开始数组里1、2、3三个数。
我们知道，n=2时，有4种排列，所以上面n=3可以分成三组，每组5个(加上空集)。
因此第10个在第二组里。所以第一个是2，把2输出。原来的数组里删除2，变成1、3两个数。然后10 - (2 - 1) * 5 = 5，即它在第2组的第5个。
减去首个空集合，5 - 1 = 4 ≠ 0，表示2后面还有数字。
因为A1 = 1是，所以再第2组里又可以分成两组，每组2个(加上空集)。
所以，4在第2组，剩下的数组中，第二个元素是3，所以输出3。再把数组里的3删除，剩下1一个数。
然后4 - (2 - 1) * 2 = 2，既它是第2组的第2个。
减去首个空集，2 - 1 = 1 ≠ 0，表示2后面还有数字。
按上面的方法继续下去，直到n = 0 或 后面为空集为止。
最后输出数组里的第1个元素，就得到2 3 1，就是解了。


从上面的计算可以看出来，本题目的关键是先求的An中每一组的个数g(n)
不难得出:g(n) = f(n) / n
∵f(n) = n[f(n-1) + 1]
∴g(n) = n[f(n-1) + 1] / n = f(n-1) + 1
∵f(n-1) = (n-1) * g(n-1)
∴g(n) = (n-1) * g(n-1) + 1


代码如下：


[cpp] view plaincopy
#include <stdio.h>  
  
  
int main()  
{  
    int i,n,t;//n：一共多少元素<=20。t：所求子集位于分组后的第几组  
    __int64 m;//位于第几个子集  
    __int64 c[21]={0};//后面将子集分组后平均每组个数，如:c[2]表示n=2时的分组每组中子集数  
    int s[21];//后面将子集按字典序分组后每组的初始首元素,组数<=20  
  
  
    for (i=1;i<21;i++)  
        c[i]=c[i-1]*(i-1)+1;//推导出来的c[n]=(n-1)*c[n-1]+1  
    while (scanf("%d%I64d",&n,&m)!=EOF)  
    {  
        for(i=0;i<21;i++)  
            s[i]=i;//每循环一次就重新归位每组首元素  
        while (n>0&&m>0)  
        {  
            t=m/c[n]+(m%c[n]>0?1:0);  
            if(t>0)//得到第m个子集在分组后的第t组,若t>0  
            {  
                printf("%d",s[t]);  
                for(i=t;i<=n;i++)  
                    s[i]=s[i+1];//或s[i]+=1,我们发现:第n组中，第2个元素在第n个时变为它的下一个数  
                m-=((t-1)*c[n]+1);//减去(t-1组总子集数+1）,m变为表示在剩余子集中位于第几个  
                putchar(m==0?'\n':' ');  
            }  
            n--;//依次递减,直到执行上面的if代码或退出  
        }  
  
  
    }  
    return 0;  
}  
具体操作步骤如下：


程序必需因素：


    1、每组子集的个数c[n];


    2、每组子集的首元素;


    3、所求子集位于当前分组后的第几组中t


    4、所求子集位于该组的第几个


主要递归步骤：


1、求出所在组t


2、输出所在组t的首元素s[t](同一组首元素相同)


3、将该子集的下一个元素到最后一个的值+1,注意这个规律：在第i组，首元素为i，删除首元素后，在第i个子集后首元素均变大+1.




程序步骤实例解说：


n=3，m=10时，有
{1}
{1, 2}
{1, 2, 3}
{1, 3}
{1, 3, 2}
{2}
{2, 1}
{2, 1, 3}
{2, 3}
{2, 3, 1}
{3}
{3, 1}
{3, 1, 2}
{3, 2}
{3, 2, 1}




1。求得t=2
先输出第2组首元素2，再去掉前面不需要的分组，和首元素,剩下唯一一组子集：
因此此时m-=((t-1)*c[n]+1)=4
//{}
{1}
{1, 3}
{3}
{3, 1}
此时的s[t~~n]均变大+1
2。然后再分成两组， t=m/c[n]+(m%c[n]>0?1:0)求得当前在第t=2组
输出第2组首元素3，再去掉前面不需要的分组，和首元素,剩下唯一一组子集
因此此时m-=((t-1)*c[n]+1)=1
//{}
{1}
3。然后剩最后一组， t=m/c[n]+(m%c[n]>0?1:0)求得当前在第t=1组
输出第1组首元素1，和首元素,剩下唯一一组子集
{}//空集
因此此时m-=((t-1)*c[n]+1)=0
最后退出。


题解二：
对于一个集合An= { 1, 2, …, n }，在n比较小的情况下，在纸上按字典顺序把所有子集排列一下。


以n=3，m=10举例：


bubuko.com,布布扣
1
1 2
1 2 3
1 3
1 3 2
2
2 1
2 1 3
2 3
2 3 1
3
3 1
3 1 2
3 2
3 2 1
n=3的情况
容易看出前5个打头的是1，紧接着5个子集打头的是2，最后5个开头的是3。


拿前五个来说，除了第一个，后面四个不看开头的1，后面的排列形式和n=2的子集的排列很相似。


f(n)代表集合An所有子集的个数，那么有递推关系：


f(n) = n * (f(n - 1) + 1)， f(1) = 1


 


这里数组taken的作用就是标记某个数是否被占用。


在这个例子里面，要求第一个数，计算(10 - 1) / 5 + 1 = 2。


表示这个数是所有未被占用的数里面从小到大第2个数，也就是2。


再计算一下余数r = (10 - 1) % 5等于4


如果r == 0说明后面的数没有了，跳出循环。


否则m = r;


继续下一轮循环


这里m == 4，计算第二个数 (4 - 1) / 2 + 1 == 2。


现在2已经被第一个数占用了，所以未被占用的第二个数就是3。


后面依次类推。


 


bubuko.com,布布扣
  //#define LOCAL
  #include <iostream>
  #include <cstdio>
  #include <cstring>
  using namespace std;
  
  int main(void)
  {
      #ifdef LOCAL
         freopen("2062in.txt", "r", stdin);
     #endif
 
     int n;
     bool taken[25];
     int b[25];
     long long m, a[25];
    a[1] = 1;
     for(int i = 2; i <= 20; ++i)
        a[i] = i * (a[i - 1] + 1);
 
     while(scanf("%d%I64d", &n, &m) == 2)
     {
         memset(taken, false, sizeof(taken));
         int i;
         long long r = 1;
         for(i = 1; i <= n; ++i)
        {
             b[i] = ((m - 1) / (a[n - i] + 1)) + 1;
             int j, k = 0;
             for(j = 1; j <= n; ++j)
             {
                 if(!taken[j])
                    ++k;
                 if(k == b[i])
                     break;
             }
             b[i] = j;
            taken[j] = true;
             r = (m - 1) % (a[n - i] + 1);
             if(r == 0)
                 break;
             m = r;
         }
        for(int j = 1; j < i; ++j)
             printf("%d ", b[j]);
        printf("%d\n", b[i]);
     }
     return 0;
 }