康托逆展开式

关于康托展开式的问题我已经在第一篇文章提到了，需要的朋友可以点进去看看。
ACM题目之排列序数
康托逆展开式就是给你一个数组,然后让你求出其全排列第n大的序列.
这里借用百度百科给的大家解释下
例1 {1,2,3,4,5}的全排列，并且已经从小到大排序完毕
(1)找出第96个数
首先用96-1得到95
用95去除4! 得到3余23
有3个数比它小的数是4
所以第一位是4
用23去除3! 得到3余5
有3个数比它小的数是4但4已经在之前出现过了所以第二位是5（4在之前出现过，所以实际比5小的数是3个）
用5去除2!得到2余1
有2个数比它小的数是3，第三位是3
用1去除1!得到1余0
有1个数比它小的数是2，第二位是2
最后一个数只能是1
所以这个数是45321
(2)找出第16个数
首先用16-1得到15
用15去除4!得到0余15
用15去除3!得到2余3
用3去除2!得到1余1
用1去除1!得到1余0
有0个数比它小的数是1
有2个数比它小的数是3 但由于1已经在之前出现过了所以是4（因为1在之前出现过了所以实际比4小的数是2）
有1个数比它小的数是2 但由于1已经在之前出现过了所以是3（因为1在之前出现过了所以实际比3小的数是1）
有1个数比它小得数是2 但由于1，3，4已经在之前出现过了所以是5（因为1，3，4在之前出现过了所以实际比5小的数是1）
最后一个数只能是2
所以这个数是14352。
code:

//康托逆展开式
//求n个数的全排列中第m大的组合：
//n>20就需要用高精度来做
#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int fact[10]= {1};
int res[10];
int n,m;
int main()
{
    for(int i=1; i<=10; ++i)
        fact[i]=fact[i-1]*i;  //阶乘数组
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        bool *vis=new bool[n];
        memset(res,0,sizeof(res));
        memset(vis,false,sizeof(vis));
        //先判断下m是否比n!还要大
        if(m>fact[n]){cout<<"impossible"<<endl;continue;}
        int i,j,x,t;
        m--;
        for(i=1; i<=n; ++i)
        {
            x=m/fact[n-i];
            m%=fact[n-i];
            for(j=1,t=0; t<=x; j++)
                if(!vis[j])t++;
            j--;
            res[i-1]=j;
            vis[j]=true;
        }
        for(i=0; i<n; ++i)
            printf("%d",res[i]);
        printf("\n");
    }
    // cout << "Hello world!" << endl;
    return 0;
}