結果

問題 No.826 連絡網
ユーザー Bwambocos
提出日時 2019-05-03 21:54:37
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 877 bytes
コンパイル時間 2,089 ms
コンパイル使用メモリ 180,732 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-12-31 17:53:17
合計ジャッジ時間 19,122 ms
ジャッジサーバーID
(参考情報)
judge1 / judge5
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 10 WA * 18 TLE * 2
権限があれば一括ダウンロードができます

ソースコード

diff #

#include "bits/stdc++.h"
#define in std::cin
#define out std::cout
#define rep(i,N) for(LL i=0;i<N;++i)
typedef long long int LL;

std::map<LL, LL> decomposit_prime(LL n)
{
	LL d = 2;
	std::map<LL, LL>res;
	while (d * d <= n)
	{
		if (n % d == 0)
		{
			++res[d];
			n /= d;
		}
		else ++d;
	}
	++res[n];
	return res;
}

LL gcd(LL a, LL b)
{
	if (a < b) std::swap(a, b);
	LL r = a % b;
	while (r != 0)
	{
		a = b;
		b = r;
		r = a % b;
	}
	return b;
}

LL lcm(LL a, LL b)
{
	return a / gcd(a, b) * b;
}

int main()
{
	LL N, P;
	in >> N >> P;

	auto primes = decomposit_prime(P);
	std::vector<LL>nums;
	for (auto i : primes) nums.push_back(i.first);
	std::sort(nums.begin(), nums.end());

	LL ans = 0;
	for (LL i = 2; i <= N; ++i)
	{
		auto ps = decomposit_prime(i);
		for (auto j : ps)
		{
			if (lcm(nums[0], j.first) <= N) ++ans;
			break;
		}
	}
	out << ans << std::endl;
}
0