結果

問題 No.847 Divisors of Power
ユーザー sugarrrsugarrr
提出日時 2020-04-07 19:15:17
言語 C++17(clang)
(17.0.6 + boost 1.83.0)
結果
AC  
実行時間 342 ms / 2,000 ms
コード長 2,684 bytes
コンパイル時間 1,601 ms
コンパイル使用メモリ 135,176 KB
実行使用メモリ 40,600 KB
最終ジャッジ日時 2023-08-20 14:38:11
合計ジャッジ時間 11,954 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 262 ms
39,496 KB
testcase_01 AC 276 ms
40,372 KB
testcase_02 AC 260 ms
39,084 KB
testcase_03 AC 277 ms
39,632 KB
testcase_04 AC 251 ms
39,348 KB
testcase_05 AC 241 ms
39,236 KB
testcase_06 AC 245 ms
39,776 KB
testcase_07 AC 263 ms
40,600 KB
testcase_08 AC 264 ms
39,556 KB
testcase_09 AC 256 ms
40,436 KB
testcase_10 AC 257 ms
39,808 KB
testcase_11 AC 276 ms
39,324 KB
testcase_12 AC 270 ms
40,080 KB
testcase_13 AC 290 ms
39,132 KB
testcase_14 AC 291 ms
39,032 KB
testcase_15 AC 342 ms
39,124 KB
testcase_16 AC 326 ms
39,364 KB
testcase_17 AC 308 ms
39,768 KB
testcase_18 AC 310 ms
39,792 KB
testcase_19 AC 308 ms
39,604 KB
testcase_20 AC 305 ms
39,316 KB
testcase_21 AC 316 ms
39,284 KB
testcase_22 AC 316 ms
39,504 KB
testcase_23 AC 308 ms
39,532 KB
testcase_24 AC 300 ms
39,048 KB
testcase_25 AC 304 ms
39,920 KB
testcase_26 AC 290 ms
39,504 KB
testcase_27 AC 288 ms
39,596 KB
testcase_28 AC 267 ms
39,800 KB
testcase_29 AC 275 ms
39,856 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

//#include <bits/stdc++.h>
#include "bits/stdc++.h"
using namespace std;
typedef long long ll;
//#include "boost/multiprecision/cpp_int.hpp"
//typedef boost::multiprecision::cpp_int LL;
typedef long double dd;
#define i_7 (ll)(1E9+7)
//#define i_7 998244353
#define i_5 i_7-2

ll mod(ll a){
    ll c=a%i_7;
    if(c>=0)return c;
    return c+i_7;
}
typedef pair<ll,ll> l_l;
typedef pair<dd,dd> d_d;
ll inf=(ll)1E16;
#define rep(i,l,r) for(ll i=l;i<=r;i++)
#define pb push_back
ll max(ll a,ll b){if(a<b)return b;else return a;}
ll min(ll a,ll b){if(a>b)return b;else return a;}
void Max(ll &pos,ll val){pos=max(pos,val);}//Max(dp[n],dp[n-1]);
void Min(ll &pos,ll val){pos=min(pos,val);}
void Add(ll &pos,ll val){pos=mod(pos+val);}
dd EPS=1E-9;
#define fastio ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);
#define fi first
#define se second
#define endl "\n"
#define SORT(v) sort(v.begin(),v.end())
#define ERASE(v) v.erase(unique(v.begin(),v.end()),v.end())
#define POSL(v,x) (lower_bound(v.begin(),v.end(),x)-v.begin())
#define POSU(v,x) (upper_bound(v.begin(),v.end(),x)-v.begin())
//template<class T>void max(T a,T b){if(a<b)return b;else return a;}
//template<class T>void min(T a,T b){if(a>b)return b;else return a;}
//template<class T>bool Max(T&a, T b){if(a < b){a = b;return 1;}return 0;}
//template<class T>bool Min(T&a, T b){if(a > b){a = b;return 1;}return 0;}

//////////////////////////



#define MAX_N (ll)2E7///変更！！！！！！！
vector<ll>prime;
bool is_prime[MAX_N+1];
void precal(){
    rep(i,0,MAX_N)is_prime[i]=true;
    is_prime[0]=is_prime[1]=false;
    rep(i,2,MAX_N){
        if(is_prime[i]){
            prime.pb(i);
            for(ll j=2*i;j<=MAX_N;j+=i)is_prime[j]=false;
        }
    }
}
vector<l_l> fact(ll n){
    vector<l_l>v;
    for(ll p:prime){
        if(p*p>n)break;
        ll c=0;
        while(n%p==0){
            n/=p;
            c++;
        }
        if(c)v.pb(l_l(p,c));//pという素因数がc個入っている
    }
    if(n!=1){
        ll vsize=v.size();
        if(vsize>=1&&v[vsize-1].first==n)v[vsize-1].second++;
        else v.pb(l_l(n,1));
    }
    if(v.size()==0){
        v.pb(l_l(1,1));
    }
    return v;
}
ll n,k,m;
ll vs;
vector<l_l>v;
ll ans=0;
void dfs(ll pos,ll res){
    
    if(pos==vs){
        if(res<=m)ans++;
        return;
    }
    
    rep(i,0,v[pos].se){
        ll a=res*powl(v[pos].fi,i);
        if(a<=m)dfs(pos+1,a);
        else return;
    }
}

int main(){fastio
    precal();
    cin>>n>>k>>m;
    if(n==1){cout<<1;return 0;}
    v=fact(n);
    
    vs=v.size();
    rep(i,0,vs-1){
        v[i].se=min(31,v[i].se*k);
    }
    dfs(0,1);
    cout<<ans<<endl;
    
    return 0;
}
0