如何将非常大的十进制数转换为阶乘数系统?

时间:2016-03-20 14:36:46

标签: c++ algorithm permutation combinatorics number-systems

我想将十进制数转换为阶乘数字系统。

我想这样做是为了找到最多100个元素的nth Lexicographic排列,例如。 A [87] = {1,2,3 ...,87}

我得到索引'n',我需要找到那个位置的词典排列。例如{1,2,3}的第二个排列是{1,3,2}

为此我试图使用阶乘数字系统。

以下链接提供有关转换方法的信息。

in the docs

如十进制所述(463)给出341010!在阶乘。

463÷1 = 463,余数为0

463÷2 = 231,其余为1

231÷3 = 77,余数为0

77÷4 = 19,余数1

19÷5 = 3,其余为4

3÷6 = 0,余数3

只有十进制数落在允许范围内时才能应用此方法,如unsigned long long int。

如果数字不能适合整数范围,该怎么办?

我的测试用例涉及的数字太大,以至于需要以字符串格式存储。(例如查找123456789012345678901234555623344数组排列[100] = {1,2,3,4,...... 100})

我正在尝试用c ++解决这个问题。

(在c ++中使用next_permutation()来获得给定的索引是一种代价高昂的方法并且需要花费很多时间。)

2 个答案:

答案 0 :(得分:0)

这是代码。你可以看到并问我任何困惑。

另外,我只提供了一个测试用例,并且我没有对代码进行详尽的测试。如果您发现任何错误,我将很乐意解决。

#include<bits/stdc++.h>
using namespace std;

#define MAX 1000000
#define MOD 1000000007
#define F first
#define S second
#define PB push_back
#define MP make_pair
#define V vector
#define I int
#define D double
#define B bool
#define pii pair<int,int>
#define LL long long

#define in(x) scanf("%d",&x)
#define in2(x,y) scanf("%d%d",&x,&y)
#define lin(x) scanf("%lld",&x)
#define lin2(x,y) scanf("%lld%lld",&x,&y)
#define FOR(i,a,b) for(i=a;i<b;i++)
#define all(v) v.begin(),v.end()

string q;    //this is the input
V<I> sol;    //this is final solution vector. (will be printed in reverse)

void divide(I n){    //function to divide a string by `n`
    string r = "";
    I i,k=0;
    FOR(i,0,q.length()){
        k *= 10;
        k += (q[i] - '0');
        I g = k / n;
        k = k % n;
        if((r.length()==0 && g!=0) || (r.length()>0)){
            r += (char)(g + '0');
        }
    }
    q = r;
    sol.PB(k);
}

I main(){
    cin>>q;
    I i;
    FOR(i,1,101){   //assuming 100 is the limit
        if(q.length()==0)
            break;
        divide(i);
    }
    //print the result
    for(i=sol.size()-1;i>=0;i--)
        //printf("%d ",sol[i]);
        cout<<sol[i]<<" ";
    printf("\n");
    return 0;
}

答案 1 :(得分:0)

虽然标题声明您正在寻找将十进制数转换为整数的方法,但我将回答您尝试解决的实际问题:如何获得N数组的Kth排列元件。

简单地说,你需要逐个数字地预测给定数组的Kth排列。事情的理论方面非常简单。假设您拥有数组A中的元素,并且存储有关是否在第二个数组S中使用每个元素的信息。当您为每个数字选择适当的值时,将更新S.结果将存储在数组R中。

有N!给定数组A中元素的排列。对于具有N个数字的数组,如果A中的最小元素被选为结果中最左边的数字R [0],则考虑有多少个排列。它(N-1)!,对吧?从#1到#(N-1)的排列!属于结果最左边的元素是A中最小元素的情况。排列#((N-1)!+ 1)到#(2 *(N-1)!)具有第二个最小值A作为R [0]。所以排列#((i-1)*(N-1)!+ 1)到#(i *(N-1)!)使用A中未使用且按字典顺序排列最小的数字作为R [0]。

在更广义的意义上,该值在Rth中使用,在第K个词典中,最小的排列是A [i],使得A [i]是第i个词典上最小的元素,因此不使用远和(i *(N-1-d)!+ 1)&lt; = k k <=((i + 1)*(N-1-) d)!)。

如果你遍历整个S,将需要 O(N)时间来找到合适的i值。我不确定你如何准确地实现它,但你也可以在S上进行二进制搜索,并在O(logN)时间内找到合适的i。

如果你有大K值,我认为你需要实现大整数乘法才能进行比较,但是如果我想到一个聪明的方法来解决这个问题,我会更新这部分答案此

一旦你选择了正确的i,你可以将A [i]指定为R [d]并继续寻找下一个数字。

以下是实现此解决方案的代码段。它很冗长,但大多数只是大整数实现。算法的要点实际上少于30行。我只是想提供一个有效的代码,以便您可以自己测试一下。

#include <iostream>
#include <cstdio>
#include <string>
#include <algorithm>
#define NLIMIT 100
#define ASIZELIMIT 101
#define BIGINTBUCKETSLIMIT 100
#define BUCKETCAPACITY 1000000000
#define DIGITSPERBUCKET 9

using namespace std;

/* sufficient big integer implementation */
class BigInt
{
    /*
     * Note that BIGINTBUCKETSLIMIT should be high enough so that
     * the values given as input does not cause overflow
     * or access violation from the last bucket in operations
     * multiply and subtract.
     */
public:
    long long buckets[BIGINTBUCKETSLIMIT];

    BigInt() {
        for(int i=0; i < BIGINTBUCKETSLIMIT; ++i) {
            buckets[i] = 0LL;
        }
    }

    BigInt(int initialValue) {
        for(int i=0; i < BIGINTBUCKETSLIMIT; ++i)
        {
            buckets[i] = initialValue % BUCKETCAPACITY;
            initialValue /= BUCKETCAPACITY;
        }
    }

    void multiply(int val) {
        for(int i= BIGINTBUCKETSLIMIT - 1; i >= 0; --i)
            buckets[i] = buckets[i] * val;

        for(int i=0; i < BIGINTBUCKETSLIMIT - 1; ++i) {
            buckets[i+1] += buckets[i] / BUCKETCAPACITY;
            buckets[i] = buckets[i] % BUCKETCAPACITY;
        }
    }

    void subtract(BigInt B) {
        for(int i= 0; i < BIGINTBUCKETSLIMIT; ++i) {
            buckets[i] = buckets[i] - B.buckets[i];
            if(buckets[i] < 0LL) {
                buckets[i] += BUCKETCAPACITY;
                buckets[i+1]--;
            }
        }
    }

    const BigInt & operator=(const BigInt &B) {
        for(int i=0; i < BIGINTBUCKETSLIMIT; ++i)
            buckets[i] = B.buckets[i];
        return *this;
    }

    bool operator<(const BigInt &B) {
        for(int i=BIGINTBUCKETSLIMIT-1; i >= 0; --i)
            if(buckets[i] != B.buckets[i])
                return buckets[i] < B.buckets[i];
        return false;
    }

    void importFromStr(string &src)
    {
        long long buffer = 0, j = 0;
        for(int i=src.size() - 1; i >= 0; i -= DIGITSPERBUCKET) {
            buffer = 0;
            for(int k=max(0, i - DIGITSPERBUCKET + 1); k <= i; ++k) {
                buffer = buffer * 10 + (src[k] - '0');
            }
            buckets[j++] = buffer;
        }
    }
};

BigInt factorials[ASIZELIMIT];

void preprocessFactorials(int n)
{
    factorials[0] = BigInt(1);
    for(int i=1; i <= n; ++i) {
        factorials[i] = factorials[i-1];
        factorials[i].multiply(i);
    }
}

void findKthPermutation(int N, int A[], BigInt K, int result[]) {
    BigInt tmpBigInt;

    bool S[ASIZELIMIT];
    for(int i=0; i < N; ++i)
        S[i] = true;
    K.subtract(BigInt(1));
    preprocessFactorials(N);

    for(int d=0; d < N; ++d) {
        for(int i=0, j=0; i < N; ++i) {
            if(S[i]) {
                tmpBigInt = factorials[N-1-d];
                tmpBigInt.multiply(j+1);
                if(K < tmpBigInt) {
                    result[d] = A[i];
                    S[i] = 0;
                    tmpBigInt = factorials[N-1-d];
                    tmpBigInt.multiply(j);
                    K.subtract(tmpBigInt);
                    break;
                }
                ++j;
            }
        }
    }
}

int main() {
    string k;
    BigInt K;
    int N;
    int A[ASIZELIMIT], R[ASIZELIMIT];

    cin >> N >> k;
    for(int i=0; i < N; ++i)
        cin >> A[i];
    K.importFromStr(k);

    sort(A, A+N);
    findKthPermutation(N, A, K, R);

    cout << R[0];
    for(int i=1; i < N; ++i)
        cout << " " << R[i];
    cout << endl;

    return 0;
}

你可以很容易地观察函数findKthPermutation和我的BigInt类中的2个循环,无论K如何,实现都在O(N 3 )中工作。虽然我不知道你的确切性能需求,如N <= 100,它可能足够有效。如果它没有你想要的那么高效,我的第一个建议是优化使用一些其他数据结构在S中存储信息,这些数据结构可以在O(logN)时间内为每个数字d寻找适当的i值。

最后,请注意,此解决方案假定A不包含重复元素,因为插入了可能的排列的词典编辑