Gym-101490H

/*************************************************************************
    > Author: MentalOmega
    > Mail: 965194745@qq.com
    > 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;
}