結果

問題 No.847 Divisors of Power
コンテスト
ユーザー ldsyb
提出日時 2019-07-05 23:35:21
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 955 bytes
コンパイル時間 1,921 ms
コンパイル使用メモリ 177,828 KB
実行使用メモリ 13,640 KB
最終ジャッジ日時 2024-10-06 23:18:32
合計ジャッジ時間 5,566 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample -- * 4
other AC * 5 WA * 10 TLE * 1 -- * 10
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

int main()
{
	int64_t n, k, m;
	cin >> n >> k >> m;

	if (n == 1)
	{
		cout << 1 << endl;
		return 0;
	}

	vector<int64_t> pf;

	for (int64_t i = 2; i * i <= n; i++)
	{
		while (n % i == 0)
		{
			pf.emplace_back(i);
			n /= i;
		}
	}
	if (n != 1)
	{
		pf.emplace_back(n);
	}

	sort(pf.begin(), pf.end());

	vector<pair<int64_t, int64_t>> v;
	v.push_back({pf[0], 1});
	for (int64_t i = 1; i < pf.size(); i++)
	{
		if (v.rbegin()->first == pf[i])
		{
			v.rbegin()->second++;
		}
		else
		{
			v.push_back({pf[i], 1});
		}
	}
	for (auto &p : v)
	{
		p.second *= k;
	}

	auto f = [&](auto &g, int64_t now = 1, int64_t i = 0) {
		if (m < now)
		{
			return int64_t(0);
		}
		if (i == v.size())
		{
			return int64_t(now <= m);
		}

		int64_t r = 0;

		for (int64_t j = 0, l = 1; j <= v[i].second; j++, l *= v[i].first)
		{
			r += g(g, now * l, i + 1);
		}

		return r;
	};

	cout << f(f) << endl;

	return 0;
}
0