結果

問題 No.847 Divisors of Power
ユーザー ngtkanangtkana
提出日時 2020-03-12 14:30:30
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 99 ms / 2,000 ms
コード長 1,117 bytes
コンパイル時間 2,152 ms
コンパイル使用メモリ 204,256 KB
最終ジャッジ日時 2025-01-09 06:32:20
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#define endl enjoy_codeforces
using lint=long long;
int main(){
    std::cin.tie(nullptr);std::ios_base::sync_with_stdio(false);
    int n,m,k;std::cin>>n>>m>>k;
    std::vector<std::pair<int,lint>>fac;
    for(int p=2;p*p<=n;p++){
        if(n%p!=0)continue;
        for(int i=0;;i++){
            if(n%p!=0){
                fac.emplace_back(p,i);
                break;
            }
            n/=p;
        }
    }
    if(n!=1)fac.emplace_back(n,1);
    for(auto&&p:fac)p.second*=m;

    std::set<int,std::greater<>>set={1};
    for(auto[p,m]:fac){
        for(int i=0;i<std::min(m,lint(std::log(k)/std::log(p)+0.01));i++){
            for(auto it=set.lower_bound(k/p);it!=set.end();it++){
                set.insert(*it*p);
            }
        }
    }
    std::cout<<set.size()<<'\n';
}
/*
 * まずは N の K 乗の素因数分解を求めます。
 * 素因数ごとに見て、約数の集合を DP です。
 * p 倍、p*p 倍、p*p*p 倍、というふうにためします。
 * set.insert( 今見ている要素 * 掛け算したいもの )
 * を降順に行います。
 */
0