結果

問題 No.1383 Numbers of Product
ユーザー carrot46carrot46
提出日時 2021-02-07 21:19:54
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 1,237 ms / 2,000 ms
コード長 1,916 bytes
コンパイル時間 1,848 ms
コンパイル使用メモリ 179,732 KB
実行使用メモリ 68,464 KB
最終ジャッジ日時 2023-08-14 22:11:51
合計ジャッジ時間 34,053 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge14
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,380 KB
testcase_01 AC 2 ms
4,376 KB
testcase_02 AC 1 ms
4,380 KB
testcase_03 AC 1 ms
4,380 KB
testcase_04 AC 1 ms
4,380 KB
testcase_05 AC 2 ms
4,380 KB
testcase_06 AC 1 ms
4,376 KB
testcase_07 AC 709 ms
43,080 KB
testcase_08 AC 693 ms
42,568 KB
testcase_09 AC 859 ms
51,944 KB
testcase_10 AC 1,146 ms
67,556 KB
testcase_11 AC 1,177 ms
66,672 KB
testcase_12 AC 1,210 ms
67,452 KB
testcase_13 AC 1,140 ms
64,764 KB
testcase_14 AC 947 ms
55,472 KB
testcase_15 AC 740 ms
45,400 KB
testcase_16 AC 1,134 ms
67,008 KB
testcase_17 AC 1 ms
4,380 KB
testcase_18 AC 1 ms
4,376 KB
testcase_19 AC 1 ms
4,376 KB
testcase_20 AC 2 ms
4,380 KB
testcase_21 AC 1 ms
4,384 KB
testcase_22 AC 1 ms
4,376 KB
testcase_23 AC 1 ms
4,380 KB
testcase_24 AC 1 ms
4,380 KB
testcase_25 AC 1 ms
4,380 KB
testcase_26 AC 1 ms
4,380 KB
testcase_27 AC 1 ms
4,376 KB
testcase_28 AC 15 ms
4,384 KB
testcase_29 AC 903 ms
53,200 KB
testcase_30 AC 1,160 ms
68,372 KB
testcase_31 AC 1 ms
4,376 KB
testcase_32 AC 228 ms
16,704 KB
testcase_33 AC 2 ms
4,384 KB
testcase_34 AC 2 ms
4,376 KB
testcase_35 AC 1,105 ms
62,576 KB
testcase_36 AC 1,040 ms
59,880 KB
testcase_37 AC 1,230 ms
68,276 KB
testcase_38 AC 1,225 ms
68,200 KB
testcase_39 AC 1,218 ms
68,196 KB
testcase_40 AC 1,233 ms
68,272 KB
testcase_41 AC 1,231 ms
68,248 KB
testcase_42 AC 1,227 ms
68,192 KB
testcase_43 AC 1,221 ms
68,184 KB
testcase_44 AC 1,213 ms
68,464 KB
testcase_45 AC 1,222 ms
68,180 KB
testcase_46 AC 1 ms
4,380 KB
testcase_47 AC 1,157 ms
66,932 KB
testcase_48 AC 1,188 ms
67,200 KB
testcase_49 AC 1,155 ms
68,028 KB
testcase_50 AC 1,237 ms
67,828 KB
testcase_51 AC 1 ms
4,376 KB
testcase_52 AC 2 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
//#include <chrono>
//#pragma GCC optimize("Ofast")
using namespace std;
#define reps(i,s,n) for(int i = s; i < n; i++)
#define rep(i,n) reps(i,0,n)
#define Rreps(i,n,e) for(int i = n - 1; i >= e; --i)
#define Rrep(i,n) Rreps(i,n,0)
#define ALL(a) a.begin(), a.end()

using ll = long long;
using vec = vector<ll>;
using mat = vector<vec>;

ll N,M,H,W,Q,K,A,B;
string S;
using P = pair<ll, ll>;
const ll INF = (1LL<<60);

template<class T> bool chmin(T &a, const T b){
    if(a > b) {a = b; return true;}
    else return false;
}
template<class T> bool chmax(T &a, const T b){
    if(a < b) {a = b; return true;}
    else return false;
}
template<class T> void my_printv(std::vector<T> v,bool endline = true){
    if(!v.empty()){
        for(std::size_t i{}; i<v.size()-1; ++i) std::cout<<v[i]<<" ";
        std::cout<<v.back();
    }
    if(endline) std::cout<<std::endl;
}

ll mfact(ll c, ll p){
    ll res = c, temp = c;
    rep(i, p) {
        temp += K;
        if(INF / temp < res) return INF;
        else res *= temp;
    }
    return res;
}

ll calc(ll n){
    //number of f(X) = 1 if X <= n
    ll lb = 0, ub = (ll)2e+9;
    while(ub - lb > 1){
        ll cen = (ub + lb) / 2;
        (mfact(cen, 1) <= n ? lb : ub) = cen;
    }
    return lb;
}

bool if_b1(ll n){
    return calc(n) - calc(n - 1) > 0;
}

int main() {
    cin.tie(nullptr);
    ios::sync_with_stdio(false);
    cin>>N>>K>>M;
    map<ll, int> f;
    for(ll b = 2; b <= 20; ++b){
        for(ll a = 1, key = mfact(a, b); key <= N; ++a){
            key = mfact(a, b);
            if(key <= N) {
                ++f[key];
            }
        }
    }
    map<int, int> num;
    num[1] = calc(N);
    for(auto &p : f){
        if(p.first > N) break;
        if(if_b1(p.first)){
            --num[1];
            ++num[p.second + 1];
        }else{
            ++num[p.second];
        }
    }
    cout<<num[M]<<endl;
}
0