結果
| 問題 |
No.847 Divisors of Power
|
| コンテスト | |
| ユーザー |
ate
|
| 提出日時 | 2019-11-21 01:11:21 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 4 ms / 2,000 ms |
| コード長 | 1,737 bytes |
| コンパイル時間 | 2,713 ms |
| コンパイル使用メモリ | 206,524 KB |
| 最終ジャッジ日時 | 2025-01-08 04:14:23 |
|
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 26 |
ソースコード
#include<bits/stdc++.h>
using namespace std;
#define rep(i,n) for(ll i=0;i<(n);++i)
using ll = long long;
using pll = pair<ll,ll>;
constexpr ll INF = (1LL<<60);
constexpr ll MOD = (1e9+7);
//constexpr ll MOD = (998244353);
template<class T> bool chmax(T &a,const T &b){if(a<b){a=b;return 1;}return 0;}
template<class T> bool chmin(T &a,const T &b){if(a>b){a=b;return 1;}return 0;}
void dump(){cout<<endl;}
template<class T,class... Ts> void dump(T&& h, Ts&&... t){cout<<h<<", ";dump(forward<Ts>(t)...);}
template<class T> istream &operator>>(istream&is,vector<T>&v){for(auto &elemnt:v)is>>elemnt;return is;}
template<class T,class U> istream &operator>>(istream&is,pair<T,U>&p){is>>p.first>>p.second;return is;}
template<class T>vector<T> make_vector(size_t a){return vector<T>(a);}
template<class T, class... Ts>auto make_vector(size_t a, Ts... ts){return vector<decltype(make_vector<T>(ts...))>(a, make_vector<T>(ts...));}
void solve1();void solve2();
int main(){
solve1();
return 0;
}
ll mpow(ll n,ll p){
ll ans=1LL;
while(p){
if(p&1)ans*=n;
n*=n;
p>>=1;
}
return ans;
}
void solve1(){
ll n,k,m;
cin>>n>>k>>m;
auto prime_factorization = [](ll n){
map<ll,ll> res;
for(ll p=2;p*p<=n;++p){
while(n%p==0){
n/=p;
res[p]++;
}
}
if(n!=1)res[n]++;
return res;
};
auto a = prime_factorization(n);
vector<pll> fac;
for(auto& [p,x]:a)fac.emplace_back(p,x*k);
auto dfs = [&](auto&& dfs,ll i,ll prod)->ll{
if(i==fac.size())return 1;
ll res=0;
rep(j,fac[i].second+1){
if(prod*mpow(fac[i].first,j)>m)break;
res += dfs(dfs,i+1,prod*mpow(fac[i].first,j));
}
return res;
};
ll ans = dfs(dfs,0,1);
cout<<ans<<endl;
}
ate