LightOJ - 1138 Trailing Zeroes (III)(n的阶乘结尾0的个数+二分)

传送门


题目大意

给出 q q q,我们需要求出最小的 n n n使得 n ! n! n!结尾有 q q q 0 0 0

题目大意

考虑 n n n的阶乘结尾 0 0 0的个数,实际上就是看最后乘了多少个 10 10 10,也就是将 n ! n! n!质因数分解得出的 2 x 5 y 2^x5^y 2x5y,结尾 0 0 0的个数就是 m i n { x , y } min\{x,y\} min{x,y},手推前几个阶乘,不难发现 x > y x >y x>y一定成立,因此我们只需考虑 5 5 5的个数。对于 n ! n! n!内含有多少质数的乘积,这是一个比较巧妙的数论知识:

  • 首先看第一层,也就是只含一个以上 5 5 5的数的个数,显然是 n / 5 n/5 n/5
  • 然后看第二层,因为阶乘是前 n n n个数相乘而不是考虑前面的每个数,那么第一层考虑的只含一个 5 5 5的数中也包含了含有两个 5 5 5以上的,因此含有两个以上 5 5 5的个数为 n / ( n / 5 ) n/(n/5) n/(n/5)
  • 以此类推,每次累加 n / 5 n/5 n/5,然后 n / = 5 n/=5 n/=5即可

随着 n n n的增大,显然其阶乘中含有 5 5 5的个数是非递减的,那么考虑二分,手推一下不难发现左边界为 1 1 1,右边界为 5 ∗ q 5*q 5q。还有就是手玩一下不难发现每 5 5 5个数其含有 5 5 5的个数是相等的,答案需要求最小的,实际上如果存在那么一定是 5 5 5的倍数,我们采用如下的二分形式:

while(l + 10 < r){
	....
}
代码

//
// Created by Happig on 2020/10/22
//
#include <bits/stdc++.h>
#include <unordered_map>
#include <unordered_set>

using namespace std;
#define fi first
#define se second
#define pb push_back
#define ins insert
#define Vector Point
#define ENDL "\n"
#define lowbit(x) (x&(-x))
#define mkp(x, y) make_pair(x,y)
#define mem(a, x) memset(a,x,sizeof a);
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
typedef pair<double, double> pdd;
const double eps = 1e-8;
const double pi = acos(-1.0);
const int inf = 0x3f3f3f3f;
const double dinf = 1e300;
const ll INF = 1e18;
const int Mod = 1e9 + 7;
const int maxn = 5e5 + 10;

int cal(ll n) {
    int ans = 0;
    while (n) {
        ans += n / 5;
        n /= 5;
    }
    return ans;
}

int main() {
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    //ios_base::sync_with_stdio(0), cin.tie(0), cout.tie(0);
    int T, kase = 0;
    ll q;
    cin >> T;
    while (T--) {
        cin >> q;
        ll l = 1, r = q * 5;
        while (l + 10 < r) {
            ll mid = (l + r) / 2;
            if (cal(mid) < q) {
                l = mid + 1;
            } else r = mid;
        }
        int ans = -1;
        for (int i = l; i <= r; i++) {
            if (cal(i) == q) {
                if (ans < 0) ans = i;
                else ans = min(ans, i);
            }
        }
        cout << "Case " << ++kase << ": ";
        if (ans < 0) cout << "impossible" << ENDL;
        else cout << ans << ENDL;
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值