Gym-101490H

/*************************************************************************    > Author: MentalOmega    > Mail: [e-mail]    > Created Time: 2017年9月14日    > function:题目说的是给一个数n，求在B进制下的每一位的乘积等于n的条件下的最小的数 ************************************************************************/#include<bits/stdc++.h>const unsigned long long MAXN = 1ull<<63;using namespace std;map<long long,long long> dp;vector<int> bp;int main(){    if (fopen("in.txt", "r") != NULL)    {        freopen("in.txt", "r", stdin);        // freopen("out.txt", "w", stdout);    }    int b;    long long n;    cin>>b>>n;    if(n==1)    {        puts("1");        return 0;    }    long long tmp=n;    for(int i=2;i<b;i++)    {        if(n%i) continue;        bp.push_back(i);        while(tmp%i==0) tmp/=i;    }    if(tmp!=1)    {        puts("impossible");        return 0;    }    dp[n]=0;//定义的是数剩余n，分解成了满足情况的最小数0    for(auto it = dp.rbegin();it!=dp.rend();it++)//这种枚举方法利用到新状态肯定比现在的状态小，这样就能枚举到新的状态    {        long long u = it->first,v=it->second;        int last=v%b;        for(auto P = bp.rbegin();P!=bp.rend()&&*P>=last;P++)        {            int p=*P;            if(u%p!=0) continue;//剩余数u中都没有p这个因数可以分解的            if(((long double)v*b+p)>MAXN) continue;//如果选择这个作为分解数会导致超过最大值2^63那么也是非法的            long long nu=u/p,nv=v*b+p;            if(!dp.count(nu)) dp[nu]=nv;//判断大小就交给了这里            else if(dp[nu]>nv) dp[nu]=nv;        }    }    //其实这题就是个状压dp，用map来实现了，还有个long double伪存大数方法    cout<<dp[1];    return 0;}