結果

問題 No.847 Divisors of Power
コンテスト
ユーザー claw88
提出日時 2019-07-05 22:44:58
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 3 ms / 2,000 ms
コード長 1,684 bytes
コンパイル時間 1,756 ms
コンパイル使用メモリ 177,296 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-06 22:41:43
合計ジャッジ時間 2,450 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
using namespace std;
#define mod 1000000007
#define dom 998244353
#define all(c) begin(c),end(c)
#define rep(i,n) repi(i,0,n)
#define repi(i,a,b) for(int i=(a);i<(b);i++)
#define rrep(i,a,b) for(int i=(b)-1;i>=(a);i--)
template <typename T> bool chmax(T &a, const T &b) { if (a < b) { a = b; return true; } return false; }
template <typename T> bool chmin(T &a, const T &b) { if (b < a) { a = b; return true; } return false; }
#define long long long
#define List vector
#define var auto
#define Count size()
#define Length size()
int dd[] = { 0, 1, 0, -1, 0 }; //→↓←↑
using pi = pair<long, long>;
long N, K, M;
vector<pi> P;
int ans = 0;

void dfs(int i, long now)
{
	if (i == P.size() - 1)
	{
		for (int j = 0; j <= P[i].second; j++)
		{
			ans++;
			now *= P[i].first;
			if (now > M)break;
		}

		//cout << now << " " << ans << endl;
		return;
	}

	for (int j = 0; j <= P[i].second; j++)
	{
		dfs(i + 1, now);
		now *= P[i].first;
		if (now > M)break;
	}
}
void solve()
{
	
	cin >> N >> K >> M;
	long W = N;
	if (N == 1)
	{
		cout << 1 << endl;
		return;
	}

	map<long, long> G;
	for (long i = 2; i * i <= W; i++)
	{
		while (W % i == 0)
		{
			G[i]++;
			W /= i;
		}
	}
	if (W != 1)
	{
		G[W]++;
	}

	for (auto v : G)
	{
		int cnt = 0;
		long mul = v.first;
		long now = 1;
		while (true)
		{
			now *= mul;
			if (now > M)
			{
				break;
			}
			cnt++;
		}
		G[v.first] = min((long)cnt, v.second * K);
		P.emplace_back(pi(v.first, G[v.first]));
	}
	//for (auto v : P)
	//{
	//	cout << v.first << " " << v.second << endl;
	//}

	dfs(0, 1);
	cout << ans << endl;
	
}

int main()
{
	cin.tie(0);
	ios::sync_with_stdio(false);
	solve();
	return 0;
}
0