1. 程式人生 > >求一個數的階乘在 m 進制下末尾 0 的個數

求一個數的階乘在 m 進制下末尾 0 的個數

long out ret pri bsp init 做的 題意 true

技術分享圖片

技術分享圖片

題意 :

  求一個數 n 的階層在 m 進制下末尾 0 的個數

思路分析 :

  如果是 10 進制地話我們是很容易知道怎麽做的,數一下其對 5 約數地個數即可,但是換成 m 進制的話就需要先將 m 分解質因數,然後然後看 n! 下因數個數最少的是幾個,即是最終答案。

代碼示例 :

#define ll long long
const ll maxn = 1e6+5;
const ll mod = 1e9+7;
const double eps = 1e-9;
const double pi = acos(-1.0);
const ll inf = 0x3f3f3f3f;

ll n, b;
ll prime[maxn];
vector<ll>ve;

void init(){
    for(ll i = 2; i <= 1000000; i++){
        if (!prime[i]){
            ve.push_back(i);
            for(ll j = 2*i; j <= 1000000; j += i){
                prime[j] = 1;
            }
        }
    }
}
ll get(ll pp, ll x){
    ll res = 0;
    
    while(x){
        res += x/pp;
        x /= pp;
    }
    return res;
}

ll cnt[maxn], num[maxn];
void solve(){
    ll f = b;
    for(ll i = 0; i < ve.size(); i++){
        if (f == 1) break;
        while(f%ve[i] == 0){
            cnt[ve[i]]++;
            f /= ve[i];
        }        
    }
    for(ll i = 0; i < ve.size(); i++){
        if (cnt[ve[i]]){
            num[ve[i]] = get(ve[i], n);
        }
    }
    ll ans = 1e18+10;
    if (f != 1) ans = min(ans, get(f, n));
    for(ll i = 0; i < ve.size(); i++){
        if (cnt[ve[i]]){
            ans = min(ans, num[ve[i]]/cnt[ve[i]]);
        }
    }
    cout << ans << endl;
}

int main() {
    //freopen("in.txt", "r", stdin);
    //freopen("out.txt", "w", stdout);
    
    cin >> n >> b;
    init();
    solve();
    return 0;
}

求一個數的階乘在 m 進制下末尾 0 的個數