結果

問題 No.847 Divisors of Power
コンテスト
ユーザー kenta255
提出日時 2019-07-25 23:47:10
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 16 ms / 2,000 ms
コード長 1,395 bytes
コンパイル時間 1,964 ms
コンパイル使用メモリ 196,856 KB
最終ジャッジ日時 2025-01-07 07:31:16
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define P pair<ll,ll>
#define FOR(I,A,B) for(ll I = ll(A); I < ll(B); ++I)
#define FORR(I,A,B) for(ll I = ll((B)-1); I >= ll(A); --I)
#define TO(x,t,f) ((x)?(t):(f))
#define SORT(x) (sort(x.begin(),x.end())) // 0 2 2 3 4 5 8 9
#define POSL(x,v) (lower_bound(x.begin(),x.end(),v)-x.begin()) //xi>=v  x is sorted
#define POSU(x,v) (upper_bound(x.begin(),x.end(),v)-x.begin()) //xi>v  x is sorted
#define NUM(x,v) (POSU(x,v)-POSL(x,v))  //x is sorted
#define REV(x) (reverse(x.begin(),x.end())) //reverse
ll gcd(ll a,ll b){if(a<b)swap(a,b);if(a%b==0)return b;return gcd(b,a%b);}
ll lcm(ll a,ll b){ll c=gcd(a,b);return ((a/c)*(b/c)*c);}//saisyo kobaisu
#define NEXTP(x) next_permutation(x.begin(),x.end())
const ll INF=1e18+7;
const ll MOD=1e9+7;
#define out(a) cout<<fixed<<setprecision((a))

ll K,N,M,S,ans=0;
vector<P> p;

void dfs(ll m=1,ll ind=0,ll used=0){
	ans++;
	if(used<p[ind].second&&m*p[ind].first<=M){
		dfs(m*p[ind].first,ind,used+1);
	}
	FOR(i,ind+1,S){
		if(m*p[i].first<=M){
			dfs(m*p[i].first,i,1);
		}
	}
}

int main(){
	cin >> N >> K >> M;
	ll j = N;
	for(ll i=2;i*i<=j;++i){
		P x = {i,0};
		while(N%i==0){
			x.second+=K;
			N /= i;
		}
		if(x.second){
			p.push_back(x);
		}
	}
	if(N!=1){
		p.push_back({N,K});
	}
	S = p.size();
	if(S==0){
		cout << 1 << endl;
		return 0;
	}
	dfs();
	cout << ans << endl;
}
0